OFFSET
1,4
LINKS
MAPLE
a:= n-> (m-> m*add(i[2]*numtheory[pi](i[1])/i[1], i=ifactors(m)[2]))
(mul(ithprime(i[2]), i=ifactors(n)[2])):
seq(a(n), n=1..120); # Alois P. Heinz, Apr 28 2022
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Apr 28 2022
STATUS
approved