OFFSET
1,3
PROG
(PARI) lista(nn) = {my(m, q, r); print1(0); for(k=1, 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(", ", k); r=m))); } \\ Jinyuan Wang, Dec 05 2020
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Klaus Brockhaus, Apr 14 2007
EXTENSIONS
a(16) inserted by and more terms from Jinyuan Wang, Dec 04 2020
STATUS
approved