OFFSET
1,2
MATHEMATICA
Do[ If[ Mod[ DivisorSigma[3, n] + DivisorSigma[2, n] + DivisorSigma[1, n] + DivisorSigma[0, n], n] == 0, Print[n]], {n, 1, 6.6 10^6}]
PROG
(PARI) isok(k) = {my(f = factor(k)); !(sum(i = 0, 3, sigma(f, i)) % k); } \\ Amiram Eldar, Feb 16 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Nov 11 2000
EXTENSIONS
a(22)-a(33) from Donovan Johnson, Mar 06 2011
a(34)-a(36) from Amiram Eldar, Feb 16 2026
STATUS
approved
