A319383 Numbers k such that phi(k)^phi(k) == 1 (mod k^2).

1, 2, 19043, 289627, 6674419, 49865347, 185014655

Offset

1,2

Comments

All terms are cyclic numbers (A003277).

The next term, if it exists, is > 10^10. - Vaclav Kotesovec, Oct 23 2018

a(8) > 10^12, if it exists. - Giovanni Resta, Oct 25 2018

Links

Table of n, a(n) for n=1..7.

Mathematica

Select[Range[20000], Divisible[EulerPhi[#]^EulerPhi[#] - 1, #^2] &] (* Vaclav Kotesovec, Oct 21 2018 *)
Join[{1}, Select[Range[1851*10^5], With[{c=EulerPhi[#]}, PowerMod[c, c, #^2] == 1&]]] (* Harvey P. Dale, Oct 09 2020 *)

Prog

(PARI) isok(n) = Mod(eulerphi(n), n^2)^eulerphi(n)==1;
for(n=1, 10000000, if(isok(n), print1(n, ", ")))

Crossrefs

Cf. A000010, A003277, A063439, A077816.

Sequence in context: A260252 A324437 A158344 * A124364 A173156 A214598

Adjacent sequences: A319380 A319381 A319382 * A319384 A319385 A319386

Keyword

nonn,more

Author

Altug Alkan, Sep 18 2018

Status

approved