RosettaCodeData/Task/Perfect-totient-numbers/Mathematica/perfect-totient-numbers.math
2023-07-01 13:44:08 -04:00

9 lines
222 B
Text

ClearAll[PerfectTotientNumberQ]
PerfectTotientNumberQ[n_Integer] := Total[Rest[Most[FixedPointList[EulerPhi, n]]]] == n
res = {};
i = 0;
While[Length[res] < 20,
i++;
If[PerfectTotientNumberQ[i], AppendTo[res, i]]
]
res