|
Search: id:A024916
|
|
|
| A024916 |
|
a(n) = sum_{k=1..n} k*floor(n/k); also sum_{k=1..n} sigma(k) where sigma(n) = sum of divisors of n (A000203). |
|
+0
26
|
|
| 1, 4, 8, 15, 21, 33, 41, 56, 69, 87, 99, 127, 141, 165, 189, 220, 238, 277, 297, 339, 371, 407, 431, 491, 522, 564, 604, 660, 690, 762, 794, 857, 905, 959, 1007, 1098, 1136, 1196, 1252, 1342, 1384, 1480, 1524, 1608, 1686, 1758, 1806, 1930, 1987, 2080, 2152
(list; graph; listen)
|
|
|
OFFSET
|
1,2
|
|
|
COMMENT
|
a(10^4) = 82256014, a(10^5) = 8224740835, a(10^6) = 822468118437, a(10^7) = 82246711794796; see A072692. - M. F. Hasler (Maximilian.Hasler(AT)gmail.com), Nov 22 2007
|
|
REFERENCES
|
For the asymptotic formula see Theorem 324 in Hardy and Wright "An introduction to the theory of numbers", Oxford university press, fifth edition, p. 266
|
|
LINKS
|
T. D. Noe, Table of n, a(n) for n = 1..1000
P. L. Patodia (pannalal(AT)usa.net), PARI program for A072692 and A024916
|
|
FORMULA
|
a(n)=n^2-A004125(n); asymptotically a(n)=n^2*Pi^2/12+O(n*Log(n)) - Benoit Cloitre (benoit7848c(AT)orange.fr), Apr 28 2002
G.f.: 1/(1-x)*sum(k>=1, x^k/(1-x^k)^2). - Benoit Cloitre (benoit7848c(AT)orange.fr), Apr 23 2003
a(n) = Sum[n - Mod[n, m], {m, 1, n}] - Roger Bagula and Gary Adamson (rlbagulatftn(AT)yahoo.com), Oct 06 2006
Row sums of triangle A130541. E.g. a(5) = 15 = (10 + 3 + 1 + 1), sum of row 4 terms of triangle A130541. - Gary W. Adamson (qntmpkt(AT)yahoo.com), Jun 03 2007
Row sums of triangle A134867 - Gary W. Adamson (qntmpkt(AT)yahoo.com), Nov 14 2007
|
|
MATHEMATICA
|
Table[Plus @@ Flatten[Divisors[Range[n]]], {n, 50}] - Alonso Delarte (alonso.delarte(AT)gmail.com), Mar 06 2006
Table[Sum[n - Mod[n, m], {m, 1, n}], {n, 1, 50}] - Roger Bagula and Gary Adamson (rlbagulatftn(AT)yahoo.com), Oct 06 2006
|
|
PROGRAM
|
(PARI) A024916(n)=sum(k=1, n, n\k*k) - M. F. Hasler (Maximilian.Hasler(AT)gmail.com), Nov 22 2007
(PARI) A024916(z) = { local(s, u, d, n, a, p); s = z*z; u = sqrtint(z); p = 2; for(d=1, u, n = z\d - z\(d+1); if(n<=1, p=d; break(), a = z%d; s -= (2*a+(n-1)*d)*n/2); ); u = z\p; for(d=2, u, s -= z%d); return(s); } - P L Patodia (pannalal(AT)usa.net), Jan 11 2008
|
|
CROSSREFS
|
Cf. A056550, A104471(2*n-1, n).
Cf. A123229, A130541, A000217, A134867, A072692.
Adjacent sequences: A024913 A024914 A024915 this_sequence A024917 A024918 A024919
Sequence in context: A136403 A071422 A113902 this_sequence A102216 A001182 A122247
|
|
KEYWORD
|
nonn,easy,nice
|
|
AUTHOR
|
Clark Kimberling (ck6(AT)evansville.edu)
|
|
|
Search completed in 0.002 seconds
|
