OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
PROG
(PARI) list(lim)=my(v=List(), t); lim+=.5; for(e=2, log(lim)\log(3), forprime(p=3, lim^(1/e), ispower(t=p^e+2, , &t); if(isprime(t), listput(v, p^e)))); vecsort(Vec(v))
\\ Charles R Greathouse IV, Apr 30 2012
(PARI) list(lim)=my(v=List()); if(lim>=25, listput(v, 25)); lim+=.5; for(e=2, log(lim)\log(3), forprime(p=3, lim^(1/e), if(isprime(p^e+2), listput(v, p^e)))); vecsort(Vec(v))
/* This second program assumes A076427(2) = 1 but is about a hundred times faster. I proved that it is correct up to 10^20 without this assumption. */
\\ Charles R Greathouse IV, Apr 30 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Sep 10 2002
EXTENSIONS
More terms from Sascha Kurz, Jan 30 2003
STATUS
approved