OFFSET
1,3
MATHEMATICA
Flatten[{0, Select[Range[1000], PrimeQ[16*#^# + 1] &]}] (* Vaclav Kotesovec, Mar 25 2018 *)
PROG
(PARI) for(n=0, 1000, if(isprime(16*n^n+1), print1(n", ")))
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Seiichi Manyama, Mar 23 2018
EXTENSIONS
a(10) from Michael S. Branicky, Apr 16 2023
STATUS
approved