|
|
A245466
|
|
a(n) = sigma_1(1) + sigma_2(2) + sigma_3(3) + ... + sigma_n-1(n-1) + sigma_n(n).
|
|
7
|
|
|
1, 6, 34, 307, 3433, 50883, 874427, 17717436, 405157609, 10414924259, 295726594871, 9214021138217, 312089127730471, 11424774176377721, 449318695089164129, 18896344248070459234, 846136606134407223412, 40192694877626586149007, 2018612350537940175272987
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Let sigma_k(n) represent the sum of the k-th powers of the divisors of n.
Then a(n) = Sum_{k=1..n} sigma_k(k), the partial sums of sigma_k(k) for k from 1 to n.
|
|
LINKS
|
|
|
FORMULA
|
a(n) = Sum_{k=1..n} sigma_k(k).
a(n) = n + Sum_{d=2..n} (d^(d*(floor(n/d)+1))-d^d)/(d^d-1). - Chayim Lowen, Aug 04 2015
|
|
EXAMPLE
|
a(1) = 1 because sigma_1(1) = sigma(1) = 1.
a(2) = 6: sigma_1(1) + sigma_2(2) = 1 + (1^2 + 2^2) = 6.
a(3) = 34: sigma_1(1) + sigma_2(2) + sigma_3(3) = 6 + (1^3 + 3^3) = 34.
a(4) = 307: sigma_1(1) + ... + sigma_4(4) = 34 + (1^4 + 2^4 + 4^4) = 307.
|
|
MAPLE
|
B:= [seq(numtheory:-sigma[n](n), n=1..100)]:
|
|
MATHEMATICA
|
Table[Sum[DivisorSigma[k, k], {k, n}], {n, 20}]
Accumulate[Table[DivisorSigma[n, n], {n, 20}]] (* Harvey P. Dale, Apr 10 2018 *)
|
|
PROG
|
(PARI)
a(n) = sum(i=1, n, sigma(i, i))
vector(50, n, a(n)) \\ Derek Orr, Jul 27 2014
(Magma) [&+[DivisorSigma(i, i): i in [1..n]]: n in [1..20]]; // Bruno Berselli, Jul 29 2014
(Magma) [n eq 1 select 1 else Self(n-1)+ DivisorSigma(n, n): n in [1..20]]; // Vincenzo Librandi, Aug 05 2015
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|