OFFSET
1,1
COMMENTS
All terms are odd since 16*(2*m)^(2*m) - 1 = (2^(m+2)*m^m - 1)*(2^(m+2)*m^m + 1). - Altug Alkan, Mar 23 2018
a(5), if it exists, is greater than 5000. - Vaclav Kotesovec, Mar 25 2018
a(5), if it exists, is greater than 20000. - Michael S. Branicky, Sep 03 2024
MATHEMATICA
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
STATUS
approved