OFFSET
1,2
COMMENTS
Fixed point is always a prime or a power of prime: fixed points are terms of A000961.
LINKS
Nathaniel Johnston, Table of n, a(n) for n = 1..10000
EXAMPLE
n=10!=3628800:list to fixed point={3628800,369,50,27}.
MATHEMATICA
ffi[x_] := Flatten[FactorInteger[x]] lf[x_] := Length[FactorInteger[x]] ba[x_] := Table[Part[ffi[x], 2*w-1], {w, 1, lf[x]}] ep[x_] := Table[Part[ffi[x], 2*w], {w, 1, lf[x]}] sex[x_] := Apply[Plus, ba[x]^ep[x]] Table[FixedPoint[sex, w], {w, 1, 128}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Apr 08 2003
STATUS
approved