|
|
OFFSET
| 1,2
|
|
|
COMMENTS
| No more terms to 10^13. [Charles R Greathouse IV, Mar 21, 2011]
|
|
|
MATHEMATICA
| a(1) = 1; s = 2^4; Do[s = s + Prime[2n]^4+Prime[2n+1]^4; If[ Mod[s, 2n+1] == 0, Print[2n+1]], {n, 1, 20000000}]
|
|
|
PROG
| (PARI) s=0; n=0; forprime(p=2, 4e9, s+=p^4; if(s%n++==0, print1(n", "))) \\ Charles R Greathouse IV, Mar 21, 2011
|
|
|
CROSSREFS
| Cf. A085450, A122142, A122103, A098999, A007504, A045345, A024450, A111441, A122102, A122140.
Sequence in context: A096628 A051024 A177087 * A025514 A204393 A178273
Adjacent sequences: A125904 A125905 A125906 * A125908 A125909 A125910
|
|
|
KEYWORD
| bref,hard,more,nonn
|
|
|
AUTHOR
| Alexander Adamchuk (alex(AT)kolmogorov.com), Feb 04 2007
|
|