OFFSET
1,1
FORMULA
Solutions of k^4 mod sigma(k) = phi(k).
EXAMPLE
sigma(76) = 140 and 76^4 mod 140 = 36 = phi(76).
MAPLE
with(numtheory): op(select(n->n^4 mod sigma(n)=phi(n), [$1..2670974]));
MATHEMATICA
Select[Range[2, 41*10^5], PowerMod[#, 4, DivisorSigma[1, #]]==EulerPhi[#]&] (* The program generates the first 20 terms of the sequence. *) (* Harvey P. Dale, Jul 09 2022 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paolo P. Lava, Feb 18 2019
EXTENSIONS
a(23)-a(32) from Giovanni Resta, Feb 19 2019
STATUS
approved