OFFSET
0,5
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
FORMULA
a(n) = A003415^9(n).
MAPLE
d:= n-> n*add(i[2]/i[1], i=ifactors(n)[2]):
A:= proc(n, k) option remember; `if`(k=0, n, d(A(n, k-1))) end:
a:= n-> A(n, 9):
seq(a(n), n=0..70);
PROG
(Python)
from sympy import factorint
def A258649(n):
for _ in range(9):
if n <= 1: return 0
n = sum((n*e//p for p, e in factorint(n).items()))
return n # Chai Wah Wu, Nov 03 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jun 06 2015
STATUS
approved