OFFSET
1,2
COMMENTS
LINKS
MATHEMATICA
With[{m = 40000}, Position[Accumulate[Table[(-1)^n * DivisorSigma[2, n], {n, 1, m}]]/Range[m], _?IntegerQ] // Flatten]
PROG
(PARI) list(lim) = my(s = 0); for(k = 1, lim, s += (-1)^k * sigma(k, 2); if(!(s % k), print1(k, ", ")));
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Jan 06 2025
STATUS
approved