OFFSET
1,2
COMMENTS
MATHEMATICA
Do[ If[ Mod[ DivisorSigma[ 2, n] + DivisorSigma[ 1, n] + DivisorSigma[ 0, n], n] == 0, Print[n]], {n, 1, 10^7}]
With[{c=Total[Table[DivisorSigma[x, #], {x, 0, 2}]]}, Select[Range[800000], Mod[c, #]==0&]] (* Harvey P. Dale, May 23 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Nov 10 2000
EXTENSIONS
a(26)-a(34) from Donovan Johnson, Jun 08 2011
STATUS
approved