OFFSET
1,1
EXAMPLE
PROG
(PARI) lista(nn) = {my(m, q, r); for(k=2, nn, if(k%4 && !(k%2 && isprime((k+1)/2)), q=4+k%2; while(!ispower(q, , &m) || !isprime((q+k)/2), q+=2); if(m>r, print1(m, ", "); r=m))); } \\ Jinyuan Wang, Dec 04 2020
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Klaus Brockhaus, Apr 14 2007
EXTENSIONS
a(16) inserted by and a(19)-a(26) from Jinyuan Wang, Dec 04 2020
STATUS
approved