A301522 - OEIS

%I #16 Apr 16 2023 20:26:11

%S 0,1,3,6,12,60,85,87,612,9073

%N Numbers k such that 16*k^k + 1 is prime.

%C Next term, if it exists, is greater than 5000. - _Vaclav Kotesovec_, Mar 25 2018

%t Flatten[{0, Select[Range[1000], PrimeQ[16*#^# + 1] &]}] (* _Vaclav Kotesovec_, Mar 25 2018 *)

%o (PARI) for(n=0, 1000, if(isprime(16*n^n+1), print1(n", ")))

%Y Cf. A110932, A160600, A301519.

%K nonn,more

%O 1,3

%A _Seiichi Manyama_, Mar 23 2018

%E a(10) from _Michael S. Branicky_, Apr 16 2023