|
|
A059919
|
|
Generalized Fermat numbers: 3^(2^n)+1, n >= 0.
|
|
16
|
|
|
|
OFFSET
|
0,1
|
|
COMMENTS
|
Generalized Fermat numbers (Ribenboim (1996))
F_n(a) := F_n(a,1) = a^(2^n) + 1, a >= 2, n >= 0, can't be prime if a is odd (as is the case for this sequence). - Daniel Forgues, Jun 19-20 2011
|
|
LINKS
|
|
|
FORMULA
|
a(0) = 4; a(n) = (a(n-1)-1)^2 + 1, n >= 1.
a(n) = 2*a(n-1)*a(n-2)*...*a(1)*a(0) + 2, n >= 0, where for n = 0, we get 2*(empty product, i.e., 1) + 2 = 4 = a(0).
The above formula implies the GCD of any pair of terms is 2, which means that the terms of (3^(2^n)+1)/2 (A059917) are pairwise coprime. - Daniel Forgues, Jun 20 & 22 2011
|
|
EXAMPLE
|
a(0) = 3^(2^0)+1 = 3^1+1 = 4 = 2*(1)+2 = 2*(empty product)+2;
a(1) = 3^(2^1)+1 = 3^2+1 = 10 = 2*(4)+2;
a(2) = 3^(2^2)+1 = 3^4+1 = 82 = 2*(4*10)+2;
a(3) = 3^(2^3)+1 = 3^8+1 = 6562 = 2*(4*10*82)+2;
a(4) = 3^(2^4)+1 = 3^16+1 = 43046722 = 2*(4*10*82*6562)+2;
a(5) = 3^(2^5)+1 = 3^32+1 = 1853020188851842 = 2*(4*10*82*6562*43046722)+2;
|
|
MAPLE
|
|
|
MATHEMATICA
|
|
|
PROG
|
(PARI) { for (n=0, 11, write("b059919.txt", n, " ", 3^(2^n) + 1); ) } \\ Harry J. Smith, Jun 30 2009
|
|
CROSSREFS
|
Cf. A000215 (Fermat numbers: 2^(2^n) + 1, n >= 0).
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|