OFFSET
1,2
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..258
MAPLE
a:= proc(n) option remember; local k;
for k from 1+ a(n-1) while
not isprime(1+(k^37-k)/(k^2-1)) do od; k
end: a(1):=1:
seq(a(n), n=1..30); # Alois P. Heinz, Jun 26 2014
MATHEMATICA
Do[If[PrimeQ[1+n+n^3+n^5+n^7+n^9+n^11+n^13+n^15+n^17+n^19+n^21+n^23 +n^25 +n^27 +n^29+n^31+n^33+n^35], Print[n]], {n, 1, 2400}]
Select[Range[5000], PrimeQ[Total[#^Range[1, 35, 2]] + 1] &] (* Vincenzo Librandi, Jun 27 2014 *)
PROG
(PARI) for(n=1, 10^4, if(ispseudoprime(sum(i=0, 17, n^(2*i+1))+1), print1(n, ", "))) \\ Derek Orr, Jun 24 2014
(Magma) [n: n in [0..4000] | IsPrime(1+n*(1+n^2)*(1+n^4+n^8)*(1+n^12+n^24))]; // Vincenzo Librandi, Jun 27 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Dec 13 2006
EXTENSIONS
a(45) and beyond from Derek Orr, Jun 24 2014
STATUS
approved