OFFSET
0,2
FORMULA
A333609(a(n)) = n.
MATHEMATICA
f[p_, e_] := p^(2^(-1 + Position[Reverse @ IntegerDigits[e, 2], 1])); iphi[1] = 1; iphi[n_] := Times @@ (Flatten@(f @@@ FactorInteger[n]) - 1); numiter[n_] := Length @ NestWhileList[iphi, n, # != 1 &] - 1; n = 0; seq = {}; Do[If[numiter[k] == n, AppendTo[seq, k]; n++], {k, 1, 1000}]; seq
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Mar 28 2020
STATUS
approved