OFFSET
1,1
COMMENTS
All a(n) == 0 or 2 (mod 6). - Robert Israel, Jul 13 2014
LINKS
Vincenzo Librandi and Robert Israel, Table of n, a(n) for n = 1..2500 (first 210 terms from Vincenzo Librandi)
MAPLE
f:= unapply(1 + sum(n^(2*j+1), j=0..28), n):
select(isprime @ f, [seq(2*i, i=1..1000)]); # Robert Israel, Jul 13 2014
MATHEMATICA
Select[Range[13000], PrimeQ[Total[#^Range[1, 57, 2]] + 1] &]
PROG
(Magma) [n: n in [0..13000] | IsPrime(s) where s is 1+&+[n^i: i in [1..57 by 2]]];
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Jun 27 2014
STATUS
approved
