OFFSET
1,2
COMMENTS
The refactorable prime powers A000961.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..102
MATHEMATICA
f[p_, max_] := Module[{k = 1, s = {}, r}, While[(r = p^(p^k - 1)) < max, AppendTo[s, r]; k++]; s]; seq[max_] := Module[{s = {1}, p = 2, s1}, While[(s1 = f[p, max]) != {}, s = Join[s, s1]; p = NextPrime[p]]; Union[s]]; seq[10^30] (* Amiram Eldar, Aug 01 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Franz Vrabec, Oct 17 2005
EXTENSIONS
More terms from Robert G. Wilson v, Oct 21 2005
a(18)-a(19) from Amiram Eldar, Aug 01 2024
STATUS
approved