|
|
A152585
|
|
Generalized Fermat numbers: 12^(2^n) + 1, n >= 0.
|
|
13
|
|
|
|
OFFSET
|
0,1
|
|
COMMENTS
|
There appears to be no divisibility rule for this sequence.
13 is the only prime up to 12^(2^15)+1.
|
|
LINKS
|
|
|
FORMULA
|
a(0) = 13; a(n)=(a(n-1)-1)^2 + 1, n >= 1.
a(n) = 11*a(n-1)*a(n-2)*...*a(1)*a(0) + 2, n >= 0, where for n = 0, we get 11*(empty product, i.e., 1)+ 2 = 13 = a(0). This implies that the terms, all odd, are pairwise coprime. - Daniel Forgues, Jun 20 2011
|
|
EXAMPLE
|
a(0) = 12^1+1 = 13 = 11(1)+2 = 11(empty product)+2.
a(1) = 12^2+1 = 145 = 11(13)+2.
a(2) = 12^4+1 = 20737 = 11(13*145)+2.
a(3) = 12^8+1 = 429981697 = 11(13*145*20737)+2.
a(4) = 12^16+1 = 184884258895036417 = 11(13*145*20737*429981697)+2.
a(5) = 12^32+1 = 34182189187166852111368841966125057 = 11(13*145*20737*429981697*184884258895036417)+2.
|
|
MATHEMATICA
|
|
|
PROG
|
(PARI) g(a, n) = if(a%2, b=2, b=1); for(x=0, n, y=a^(2^x)+b; print1(y", "))
(Python)
|
|
CROSSREFS
|
Cf. A000215 (Fermat numbers: 2^(2^n)+1, n >= 0).
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|