OFFSET
1,2
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 *)
PROG
(PARI) isok(k) = my(f = factor(k)); !(sum(i = 0, 2, sigma(f, i)) % k); \\ Amiram Eldar, Feb 16 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Nov 10 2000
EXTENSIONS
a(26)-a(34) from Donovan Johnson, Jun 08 2011
a(35)-a(37) from Amiram Eldar, Feb 16 2026
STATUS
approved
