OFFSET
0,3
MATHEMATICA
k = 1; lst = {0}; While[k < 30000, If[PrimeQ[2^k + k^2 + 2], AppendTo[lst, k]]; k += 2]; lst (* Robert G. Wilson v, Jun 08 2010 *)
Join[{0}, Select[Range[1, 21000, 2], PrimeQ[2^#+#^2+2]&]] (* Harvey P. Dale, Dec 18 2011 *)
PROG
(Magma) [n: n in [0..2200]|IsPrime(2^n+n^2+2)]
(PARI) is(n)=ispseudoprime(2^n+n^2+2) \\ Charles R Greathouse IV, Jun 13 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Vincenzo Librandi, May 31 2010
EXTENSIONS
a(17) & a(18) from Robert G. Wilson v, Jun 08 2010
STATUS
approved