OFFSET
1,2
COMMENTS
a(22) > 200000. - Giovanni Resta, Mar 23 2014
a(22) > 400000. - Robert Price, Sep 01 2014
MAPLE
A100358:=n->`if`(isprime(2^n+n^3+1), n, NULL): seq(A100358(n), n=0..10^3); # Wesley Ivan Hurt, Sep 01 2014
MATHEMATICA
{ta={{0}}, tb={{0}}}; Do[g=n; s=2^n+n^3+1; If[PrimeQ[s], Print[n]; ta=Append[ta, n]; tb=Append[tb, s]], {n, 1, 10000}]; {ta, tb, g}
PROG
(PARI)
for(n=1, 10^5, if(ispseudoprime(2^n+n^3+1), print1(n, ", "))) \\ Derek Orr, Sep 01 2014
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Labos Elemer, Nov 19 2004
EXTENSIONS
a(17)-a(21) from Giovanni Resta, Mar 23 2014
STATUS
approved