OFFSET
1,1
COMMENTS
See comments in A006549.
FORMULA
a(n) = A006549(n) + 1. - Michel Marcus, Aug 25 2013
MATHEMATICA
Join[{2}, Select[Range[0, 10^6], PrimePowerQ[#]&&PrimePowerQ[#-1]&]] (* James C. McMahon, Oct 05 2024 *)
PROG
(PARI) isokp(n) = (n==1) || isprime(n) || (ispower(n, , &k) && isprime(k));
isok(n) = isokp(n) && isokp(n-1); \\ Michel Marcus, Aug 25 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Greg Huber, Jul 13 2006
STATUS
approved