|
| |
|
|
A297402
|
|
a(n) = gcd_{k=1..n} (prime(k+1)^n-1)/2.
|
|
2
|
|
|
|
1, 4, 1, 8, 1, 4, 1, 16, 1, 4, 1, 8, 1, 4, 1, 32, 1, 4, 1, 8, 1, 4, 1, 16, 1, 4, 1, 8, 1, 4, 1, 64, 1, 4, 1, 8, 1, 4, 1, 16, 1, 4, 1, 8, 1, 4, 1, 32, 1, 4, 1, 8, 1, 4, 1, 16, 1, 4, 1, 8, 1, 4, 1, 128, 1, 4, 1, 8, 1, 4, 1, 16, 1, 4, 1, 8, 1, 4, 1, 32, 1, 4, 1, 8, 1, 4, 1, 16, 1, 4, 1, 8, 1, 4, 1, 64, 1, 4, 1, 8
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
1,2
|
|
|
COMMENTS
|
If p is an odd prime and p^n is the length of the odd leg of a primitive Pythagorean triangle it constrains the other leg and hypotenuse to be (p^(2n)-1)/2 and (p^(2n)+1)/2. The resulting triangle has a semiperimeter of p^n(p^n+1)/2, an area of (p^n-1)p^n(p^n+1)/4 and an inradius of (p^n-1)/2. a(n) equals the GCD of the inradius terms (p^n-1)/2 for at least the first n odd primes.
Conjecture: a(n) equals the GCD of the inradius terms (p^n-1)/2 for all odd primes, i.e. a(n) = GCD_{k=1..oo} (prime(k+1)^n-1)/2.
All terms are powers of 2. Proof: suppose p | a(n) for some odd prime p. Then p | (p^n - 1) / 2 and so p | (p^n - 1) which isn't the case.
If n is odd then a(n) = 1. Proof: 2 | (p^k - 1) for all k and odd primes p. 3^n - 1 = 3 * 9^k - 1 = 3 - 1 = 2 (mod 4), so 3^n - 1 is of the form 2*m for some odd m, hence the GCD of all (p^n - 1) / 2 is 1 for odd n. (End)
|
|
|
LINKS
|
|
|
|
FORMULA
|
It appears that for k > 0, a(2^k) = 2^(k+1).
Multiplicative with a(2^e) = 2^(e+1), a(p^e) = 1 for odd prime p. - Andrew Howroyd, Jul 25 2018
It appears that for m > 0, a(2m-1) = 1 (proved in comments) and a(2m) = 2^(k+1) where k is the exponent of the even prime in the prime factorization of 2m. - Frank M Jackson, Jul 28 2018
Dirichlet g.f.: zeta(s) * (1 + 1/2^s + 1/(2^(s-1) - 1)).
Sum_{k=1..n} a(k) ~ (n/log(2)) * (log(n) + gamma + log(2) - 1), where gamma is Euler's constant (A001620). (End)
|
|
|
EXAMPLE
|
a(4)=8 because for n=4 and for the first 4 odd primes {3, 5, 7, 11}, the term (p^n-1)/2 gives {40, 312, 1200, 7320} with a GCD of 8.
|
|
|
MATHEMATICA
|
a[n_] := GCD @@ Array[(Prime[# +1]^n -1)/2 &, n]; Array[a, 90] (* slightly modified by Robert G. Wilson v, Jan 01 2018 *)
a[n_] := If[EvenQ[n], 2^(FactorInteger[n][[1]][[2]] + 1), 1]; Array[a, 90] (* Frank M Jackson, Jul 28 2018 *)
|
|
|
PROG
|
(PARI) a(n) = gcd(vector(n, i, (prime(i+1)^n-1)/2)) \\ Iain Fox, Dec 29 2017
|
|
|
CROSSREFS
|
|
|
|
KEYWORD
|
nonn,easy,mult
|
|
|
AUTHOR
|
|
|
|
STATUS
|
approved
|
| |
|
|
