%I #7 Mar 30 2012 17:34:48
%S 1,2,6,182,330,850,2150,5964,23028,85800,168150,172250,179850,308700,
%T 391300,484692,507725,1098750,1248148,1377075,3105218,4723500,9547650,
%U 17353350,22058792,37465350,47094425,54658450,110019000,144820200
%N Numbers n such that n | Sigma_2(n) + Phi(n)^2.
%C sigma_2(n) is the sum of the squares of the divisors of n (A001157).
%t Do[If[Mod[DivisorSigma[2, n]+EulerPhi[n]^2, n]==0, Print[n]], {n, 1, 2*10^6}]
%K nonn
%O 1,2
%A _Robert G. Wilson v_, Jun 09 2000
%E a(21)-a(30) from _Donovan Johnson_, Feb 13 2010
|