use Prime::Factor;
my \π = lazy 0, |(1..*).hyper.map: -> \t { t * [*] t.&prime-factors.squish.map: 1 - 1/* }
my \ππ = Nil, |(3, *+2 β¦ *).grep: -> \p { p == sum π[p], { π[$_] } β¦ 1 };
put "The first twenty Perfect totient numbers:\n", ππ[1..20];
The first twenty Perfect totient numbers:
3 9 15 27 39 81 111 183 243 255 327 363 471 729 2187 2199 3063 4359 4375 5571