OFFSET
1,2
COMMENTS
The powers of primes correspond to A000961.
Apparently, a(n) > 0 for all n > 0.
LINKS
Rémy Sigrist, Table of n, a(n) for n = 1..10000
EXAMPLE
PROG
(PARI) a(n) = for (k=0, oo, if (omega(n) <= 1, return (n), my (f=factor(n)); n += f[#f~, 1]))
CROSSREFS
KEYWORD
nonn
AUTHOR
Rémy Sigrist, Jan 05 2019
STATUS
approved