OFFSET
1,2
PROG
(PARI) a(n) = { my (f=factor(n)); #Set(apply (k -> prod (i=1, #f~, f[i, 1]^(f[i, 2]\k)), [1..1+if (n==1, 0, vecmax(f[, 2]~))])) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Rémy Sigrist, Apr 19 2020
STATUS
approved