OFFSET
1,2
MATHEMATICA
Do[If[Mod[DivisorSigma[2, n]+EulerPhi[n]^2, n]==0, Print[n]], {n, 1, 2*10^6}]
PROG
(PARI) isok(k) = {my(f = factor(k)); !((sigma(f, 2) + eulerphi(f)^2) % k); } \\ Amiram Eldar, Feb 16 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Jun 09 2000
EXTENSIONS
a(21)-a(30) from Donovan Johnson, Feb 13 2010
a(31)-a(34) from Amiram Eldar, Feb 16 2026
STATUS
approved
