OFFSET
1,2
COMMENTS
sigma_2(n) is the sum of the squares of the divisors of n (A001157).
MATHEMATICA
Do[If[Mod[DivisorSigma[2, n]+EulerPhi[n]^2, n]==0, Print[n]], {n, 1, 2*10^6}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Jun 09 2000
EXTENSIONS
a(21)-a(30) from Donovan Johnson, Feb 13 2010
STATUS
approved