OFFSET
0,5
COMMENTS
a(n) is zero for all prime n.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..80
MATHEMATICA
dn[0]=0; dn[1]=0; dn[n_] := Module[{f=Transpose[FactorInteger[n]]}, If[PrimeQ[n], 1, Plus@@(n*f[[2]]/f[[1]])]]; Table[Nest[dn, n, n], {n, 50}]
PROG
(Python)
from sympy import factorint
def A185232(n):
for _ in range(n):
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
José María Grau Ribas, Jan 24 2012
STATUS
approved