OFFSET
1,3
COMMENTS
Next term, if it exists, is greater than 10000. - Vaclav Kotesovec, Mar 25 2018
MAPLE
a:=n->`if`(isprime(4^(2*k)+2*k+1), k, NULL): seq(a(k), k=0..5000); # Muniru A Asiru, Mar 25 2018
MATHEMATICA
Flatten[{0, Select[Range[1000], PrimeQ[16^# + 2*# + 1] &]}] (* Vaclav Kotesovec, Mar 25 2018 *)
PROG
(PARI) for(n=0, 500, if(isprime(4^(2*n)+2*n+1), print1(n", ")))
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Seiichi Manyama, Mar 25 2018
EXTENSIONS
a(7) from Muniru A Asiru, Mar 25 2018
a(8) from Vaclav Kotesovec, Mar 25 2018
STATUS
approved