OFFSET
1,2
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..500
Amiram Eldar, Table of n, a(n), A333611(a(n)) for n=1..500
MATHEMATICA
f[p_, e_] := p^(2^(-1 + Position[Reverse @ IntegerDigits[e, 2], 1])); iphi[1] = 1; iphi[n_] := Times @@ (Flatten@(f @@@ FactorInteger[n]) - 1); s[n_] := Plus @@ NestWhileList[iphi, n, # != 1 &] - n; seq = {}; smax = -1; Do[s1 = s[n]; If[s1 > smax, smax = s1; AppendTo[seq, n]], {n, 1, 10^5}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 28 2020
STATUS
approved