1,2
Do[If[Mod[DivisorSigma[1, n]-DivisorSigma[0, n]+EulerPhi[n], n]==0, Print[n]], {n, 1, 10^7}]
Sequence in context: A015217 A113094 A076113 * A007257 A045487 A022376
Adjacent sequences: A056010 A056011 A056012 * A056014 A056015 A056016
nonn
Robert G. Wilson v (rgwv(AT)rgwv.com), Jul 24 2000