OFFSET
1,2
COMMENTS
MATHEMATICA
ephi[n_] := Times @@ EulerPhi[FactorInteger[n][[;; , 2]]]; s[n_] := Plus @@ FixedPointList[ephi, n] - n - 1; seq = {}; smax = -1; Do[s1 = s[n]; If[s1 > smax, smax = s1; AppendTo[seq, n]], {n, 1, 5000}]; seq
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Feb 25 2020
STATUS
approved