OFFSET
1,3
COMMENTS
a(13) > 2*10^5. - Robert Price, Mar 06 2014
MATHEMATICA
Select[Range[0, 5000], PrimeQ[3^# + #^3 + 1] & ]
PROG
(Magma) [n: n in [0..800] | IsPrime(3^n+n^3+1)];
(PARI) is(n)=isprime(3^n+n^3+1) \\ Charles R Greathouse IV, Mar 06 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Vincenzo Librandi, Sep 03 2012
EXTENSIONS
a(12) from Robert Price, Mar 06 2014
STATUS
approved