OFFSET
1,4
COMMENTS
LINKS
Rémy Sigrist, Table of n, a(n) for n = 1..10000
Rémy Sigrist, Illustration of initial terms
FORMULA
EXAMPLE
See Links section.
PROG
(PARI) a(n) = { if (n==1, 1, my (f=factor(n), w=#f~, p=f[w, 1], e=f[w, 2]); if (e==1, n/p, n*p^(a(e)-e))) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Rémy Sigrist, Apr 23 2021
STATUS
approved