OFFSET
1,2
COMMENTS
FORMULA
(Sum_{j=1..k} sigma_4(j)) mod k = A064604(k) mod k = 0.
EXAMPLE
Adding 4th-power divisor-sums for j = 1..7 gives 1+17+82+273+626+1394+2402 = 4795 which is divisible by 7, so 7 is a term and the integer quotient is 655.
MATHEMATICA
k = 1; lst = {}; s = 0; While[k < 1000000001, s = s + DivisorSigma[4, k]; If[ Mod[s, k] == 0, AppendTo[lst, k]; Print@ k]; k++]; lst (* _Robert G.Wilson v_, Aug 25 2011 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Labos Elemer, Sep 24 2001
EXTENSIONS
a(13)-a(18) from Donovan Johnson, Jun 21 2010
a(19)-a(23) from Amiram Eldar, Jan 18 2024
STATUS
approved