OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MATHEMATICA
Select[Range[7000], PrimeQ[Total[#^Range[1, 17, 2]] + 1] &]
PROG
(Magma) [n: n in [0..2500] | IsPrime(s) where s is 1+&+[n^i: i in [1..17 by 2]]];
(PARI) isok(n) = isprime(1 + n + n^3 + n^5 + n^7 + n^9 + n^11 + n^13 + n^15 + n^17); \\ Michel Marcus, Jun 27 2014
(Sage)
i, n = var('i, n')
[n for n in (1..2000) if is_prime(1+(n^(2*i+1)).sum(i, 0, 8))] # Bruno Berselli, Jun 27 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Jun 27 2014
STATUS
approved