

A297402


a(n) = gcd_{k=1..n} (prime(k+1)^n1)/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^n1)p^n(p^n+1)/4 and an inradius of (p^n1)/2. a(n) equals the GCD of the inradius terms (p^n1)/2 for at least the first n odd primes.
Conjecture: a(n) equals the GCD of the inradius terms (p^n1)/2 for all odd primes, i.e. a(n) = GCD_{k=1..oo} (prime(k+1)^n1)/2.
From David A. Corneth, Dec 29 2017: (Start)
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)
This is the even bisection of A059159.  Rémy Sigrist, Dec 30 2017


LINKS

Frank M Jackson, Table of n, a(n) for n = 1..10000


FORMULA

It appears that for k > 0, a(2^k) = 2^(k+1).
a(n) = A006519(2n) for even n and a(n) = 1 for odd n.  David A. Corneth, Dec 29 2017
a(n) = A074723(n)/2.  Iain Fox, Dec 30 2017
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(2m1) = 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


EXAMPLE

a(4)=8 because for n=4 and for the first 4 odd primes {3, 5, 7, 11}, the term (p^n1)/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)^n1)/2)) \\ Iain Fox, Dec 29 2017
(PARI) a(n)=if(n%2, 1, 2)<<valuation(n, 2) \\ Charles R Greathouse IV, Jan 06 2018


CROSSREFS

Cf. A006519, A059159, A074723, A127922, A171977, A202318.
Sequence in context: A329998 A143320 A089146 * A239666 A101512 A085994
Adjacent sequences: A297399 A297400 A297401 * A297403 A297404 A297405


KEYWORD

nonn,easy,mult


AUTHOR

Frank M Jackson, Dec 29 2017


STATUS

approved



