A301522 Numbers k such that 16*k^k + 1 is prime.
0, 1, 3, 6, 12, 60, 85, 87, 612, 9073
Offset: 1
Programs
-
Mathematica
Flatten[{0, Select[Range[1000], PrimeQ[16*#^# + 1] &]}] (* Vaclav Kotesovec, Mar 25 2018 *)
-
PARI
for(n=0, 1000, if(isprime(16*n^n+1), print1(n", ")))
Extensions
a(10) from Michael S. Branicky, Apr 16 2023