A340793 -id:A340793 - cp's OEIS frontend

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).

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

Offset: 1

Clark Kimberling

Row sums of triangle A128489. E.g., a(5) = 15 = (10 + 3 + 1 + 1), sum of row 4 terms of triangle A128489. - Gary W. Adamson, Jun 03 2007
Row sums of triangle A134867. - Gary W. Adamson, Nov 14 2007
a(10^4) = 82256014, a(10^5) = 8224740835, a(10^6) = 822468118437, a(10^7) = 82246711794796; see A072692. - M. F. Hasler, Nov 22 2007
Equals row sums of triangle A158905. - Gary W. Adamson, Mar 29 2009
n is prime if and only if a(n) - a(n-1) - 1 = n. - Omar E. Pol, Dec 31 2012
Also the alternating row sums of A236104. - Omar E. Pol, Jul 21 2014
a(n) is also the total number of parts in all partitions of the positive integers <= n into equal parts. - Omar E. Pol, Apr 30 2017
a(n) is also the total area of the terraces of the stepped pyramid with n levels described in A245092. - Omar E. Pol, Nov 04 2017
a(n) is also the area under the Dyck path described in the n-th row of A237593 (see example). - Omar E. Pol, Sep 17 2018
From Omar E. Pol, Feb 17 2020: (Start)
Convolution of A340793 and A000027.
Convolved with A340793 gives A000385. (End)
a(n) is also the number of cubic cells (or cubes) in the n-th level starting from the top of the stepped pyramid described in A245092. - Omar E. Pol, Jan 12 2022

			From _Omar E. Pol_, Aug 20 2021: (Start)
For n = 6 the sum of all divisors of the first six positive integers is [1] + [1 + 2] + [1 + 3] + [1 + 2 + 4] + [1 + 5] + [1 + 2 + 3 + 6] = 1 + 3 + 4 + 7 + 6 + 12 = 33, so a(6) = 33.
On the other hand the area under the Dyck path of the 6th diagram as shown below is equal to 33, so a(6) = 33.
Illustration of initial terms:                        _ _ _ _
                                        _ _ _        |       |_
                            _ _ _      |     |       |         |_
                  _ _      |     |_    |     |_ _    |           |
          _ _    |   |_    |       |   |         |   |           |
    _    |   |   |     |   |       |   |         |   |           |
   |_|   |_ _|   |_ _ _|   |_ _ _ _|   |_ _ _ _ _|   |_ _ _ _ _ _|
.
    1      4        8          15           21             33         (End)

Hardy and Wright, "An introduction to the theory of numbers", Oxford University Press, fifth edition, p. 266.

Daniel Mondot, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Vaclav Kotesovec, Plot of (a(n) - Pi^2*n^2/12) / (n*log(n)^(2/3)) for n = 2..100000.
P. L. Patodia (pannalal(AT)usa.net), PARI program for A072692 and A024916.
Peter Polm, C# program for A024916.
A. Walfisz, Weylsche Exponentialsummen in der neueren Zahlentheorie, ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik, Volume 44, Issue 12, page 607, 1964.

Partial sums of A000203.
Cf. A056550, A104471(2*n-1, n), A123229, A128489, A000217, A134867, A072692, A158905, A237593, A245092, A006218, A222548, A092406, A160664.
Cf. A000385, A010766, A340793.

a024916 n = sum $ map (\k -> k * div n k) [1..n]
-- Reinhard Zumkeller, Apr 20 2015

[(&+[DivisorSigma(1, k): k in [1..n]]): n in [1..60]]; // G. C. Greubel, Mar 15 2019

A024916 := proc(n)
    add(numtheory[sigma](k),k=0..n) ;
end proc: # Zerinvary Lajos, Jan 11 2009
# second Maple program:
a:= proc(n) option remember; `if`(n=0, 0,
      numtheory[sigma](n)+a(n-1))
    end:
seq(a(n), n=1..100);  # Alois P. Heinz, Sep 12 2019

Table[Plus @@ Flatten[Divisors[Range[n]]], {n, 50}] (* Alonso del Arte, Mar 06 2006 *)
Table[Sum[n - Mod[n, m], {m, n}], {n, 50}] (* Roger L. Bagula and Gary W. Adamson, Oct 06 2006 *)
a[n_] := Sum[DivisorSigma[1, k], {k, n}]; Table[a[n], {n, 51}] (* Jean-François Alcover, Dec 16 2011 *)
Accumulate[DivisorSigma[1,Range[60]]] (* Harvey P. Dale, Mar 13 2014 *)

A024916(n)=sum(k=1,n,n\k*k) \\ M. F. Hasler, Nov 22 2007

A024916(z) = { my(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); } \\ See the link for a nicely formatted version. - P. L. Patodia (pannalal(AT)usa.net), Jan 11 2008

A024916(n)={my(s=0,d=1,q=n);while(dPeter Polm, Aug 18 2014

A024916(n)={ my(s=n^2, r=sqrtint(n), nd=n, D); for(d=1, r, (1>=D=nd-nd=n\(d+1)) && (r=d-1) && break; s -= n%d*D+(D-1)*D\2*d); s - sum(d=2, n\(r+1), n%d)} \\ Slightly optimized version of Patodia's code. - M. F. Hasler, Apr 18 2015
(C#) See Polm link.

def A024916(n): return sum(k*(n//k) for k in range(1,n+1)) # Chai Wah Wu, Dec 17 2021

from math import isqrt
def A024916(n): return (-(s:=isqrt(n))**2*(s+1) + sum((q:=n//k)*((k<<1)+q+1) for k in range(1,s+1)))>>1 # Chai Wah Wu, Oct 21 2023

[sum(sigma(k) for k in (1..n)) for n in (1..60)] # G. C. Greubel, Mar 15 2019

From Benoit Cloitre, Apr 28 2002: (Start)
a(n) = n^2 - A004125(n).
Asymptotically a(n) = n^2*Pi^2/12 + O(n*log(n)). (End)
G.f.: (1/(1-x))*Sum_{k>=1} x^k/(1-x^k)^2. - Benoit Cloitre, Apr 23 2003
a(n) = Sum_{m=1..n} (n - (n mod m)). - Roger L. Bagula and Gary W. Adamson, Oct 06 2006
a(n) = n^2*Pi^2/12 + O(n*log(n)^(2/3)) [Walfisz]. - Charles R Greathouse IV, Jun 19 2012
a(n) = A000217(n) + A153485(n). - Omar E. Pol, Jan 28 2014
a(n) = A000292(n) - A076664(n), n > 0. - Omar E. Pol, Feb 11 2014
a(n) = A078471(n) + A271342(n). - Omar E. Pol, Apr 08 2016
a(n) = (1/2)*(A222548(n) + A006218(n)). - Ridouane Oudra, Aug 03 2019
From Greg Dresden, Feb 23 2020: (Start)
a(n) = A092406(n) + 8, n>3.
a(n) = A160664(n) - 1, n>0. (End)
a(2*n) = A326123(n) + A326124(n). - Vaclav Kotesovec, Aug 18 2021
a(n) = Sum_{k=1..n} k * A010766(n,k). - Georg Fischer, Mar 04 2022

A066186 Sum of all parts of all partitions of n.

Original entry on oeis.org

0, 1, 4, 9, 20, 35, 66, 105, 176, 270, 420, 616, 924, 1313, 1890, 2640, 3696, 5049, 6930, 9310, 12540, 16632, 22044, 28865, 37800, 48950, 63336, 81270, 104104, 132385, 168120, 212102, 267168, 334719, 418540, 520905, 647172, 800569, 988570, 1216215, 1493520

Offset: 0

Wouter Meeussen, Dec 15 2001

Sum of the zeroth moments of all partitions of n.
Also the number of one-element transitions from the integer partitions of n to the partitions of n-1 for labeled parts with the assumption that any part z is composed of labeled elements of amount 1, i.e., z = 1_1 + 1_2 + ... + 1_z. Then one can take from z a single element in z different ways. E.g., for n=3 to n=2 we have A066186(3) = 9 and [111] --> [11], [111] --> [11], [111] --> [11], [12] --> [111], [12] --> [111], [12] --> [2], [3] --> 2, [3] --> 2, [3] --> 2. For the unlabeled case, one can take a single element from z in only one way. Then the number of one-element transitions from the integer partitions of n to the partitions of n-1 is given by A000070. E.g., A000070(3) = 4 and for the transition from n=3 to n=2 one has [111] --> [11], [12] --> [11], [12] --> [2], [3] --> [2]. - Thomas Wieder, May 20 2004
Also sum of all parts of all regions of n (Cf. A206437). - Omar E. Pol, Jan 13 2013
From Omar E. Pol, Jan 19 2021: (Start)
Apart from initial zero this is also as follows:
Convolution of A000203 and A000041.
Convolution of A024916 and A002865.
For n >= 1, a(n) is also the number of cells in a symmetric polycube in which the terraces are the symmetric representation of sigma(k), for k = n..1, (cf. A237593) starting from the base and located at the levels A000041(0)..A000041(n-1) respectively. The polycube looks like a symmetric tower (cf. A221529). A dissection is a three-dimensional spiral whose top view is described in A239660. The growth of the volume of the polycube represents each convolution mentioned above. (End)
From Omar E. Pol, Feb 04 2021: (Start)
a(n) is also the sum of all divisors of all positive integers in a sequence with n blocks where the m-th block consists of A000041(n-m) copies of m, with 1 <= m <= n. The mentioned divisors are also all parts of all partitions of n.
Apart from initial zero this is also the convolution of A340793 and A000070. (End)

			a(3)=9 because the partitions of 3 are: 3, 2+1 and 1+1+1; and (3) + (2+1) + (1+1+1) = 9.
a(4)=20 because A000041(4)=5 and 4*5=20.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
F. G. Garvan, Higher-order spt functions, Adv. Math. 228 (2011), no. 1, 241-265, alternate copy. - From _N. J. A. Sloane_, Jan 02 2013
F. G. Garvan, Higher-order spt functions, arXiv:1008.1207 [math.NT], 2010.
T. J. Osler, A. Hassen and T. R. Chandrupatia, Surprising connections between partitions and divisors, The College Mathematics Journal, Vol. 38. No. 4, Sep. 2007, 278-287 (see p. 287).
Omar E. Pol, Illustration of a(10), prism and tower, each polycube contains 420 cubes.
Omar E. Pol, Illustration of initial terms of A066186 and of A139582 (n>=1)

Cf. A000041, A093694, A000070, A132825, A001787 (same for ordered partitions), A277029, A000203, A221529, A237593, A239660.
First differences give A138879. - Omar E. Pol, Aug 16 2013
Row sums of triangles A138785, A181187, A245099, A337209, A339106, A340423, A340424, A221529, A302246, A338156, A340035, A340056, A340057, A346741. - Omar E. Pol, Aug 02 2021

a066186 = sum . concat . ps 1 where
   ps _ 0 = [[]]
   ps i j = [t:ts | t <- [i..j], ts <- ps t (j - t)]
-- Reinhard Zumkeller, Jul 13 2013

with(combinat): a:= n-> n*numbpart(n): seq(a(n), n=0..50); # Zerinvary Lajos, Apr 25 2007

PartitionsP[ Range[0, 60] ] * Range[0, 60]

a(n)=numbpart(n)*n \\ Charles R Greathouse IV, Mar 10 2012

from sympy import npartitions
def A066186(n): return n*npartitions(n) # Chai Wah Wu, Oct 22 2023

[n*Partitions(n).cardinality() for n in range(41)] # Peter Luschny, Jul 29 2014

a(n) = n * A000041(n). - Omar E. Pol, Oct 10 2011
G.f.: x * (d/dx) Product_{k>=1} 1/(1-x^k), i.e., derivative of g.f. for A000041. - Jon Perry, Mar 17 2004 (adjusted to match the offset by Geoffrey Critzer, Nov 29 2014)
Equals A132825 * [1, 2, 3, ...]. - Gary W. Adamson, Sep 02 2007
a(n) = A066967(n) + A066966(n). - Omar E. Pol, Mar 10 2012
a(n) = A207381(n) + A207382(n). - Omar E. Pol, Mar 13 2012
a(n) = A006128(n) + A196087(n). - Omar E. Pol, Apr 22 2012
a(n) = A220909(n)/2. - Omar E. Pol, Jan 13 2013
a(n) = Sum_{k=1..n} A000203(k)*A000041(n-k), n >= 1. - Omar E. Pol, Jan 20 2013
a(n) = Sum_{k=1..n} k*A036043(n,n-k+1). - L. Edson Jeffery, Aug 03 2013
a(n) = Sum_{k=1..n} A024916(k)*A002865(n-k), n >= 1. - Omar E. Pol, Jul 13 2014
a(n) ~ exp(Pi*sqrt(2*n/3))/(4*sqrt(3)) * (1 - (sqrt(3/2)/Pi + Pi/(24*sqrt(6))) / sqrt(n)). - Vaclav Kotesovec, Oct 24 2016
a(n) = Sum_{k=1..n} A340793(k)*A000070(n-k), n >= 1. - Omar E. Pol, Feb 04 2021

a(0) added by Franklin T. Adams-Watters, Jul 28 2014

A244050 Partial sums of A243980.

Original entry on oeis.org

4, 20, 52, 112, 196, 328, 492, 716, 992, 1340, 1736, 2244, 2808, 3468, 4224, 5104, 6056, 7164, 8352, 9708, 11192, 12820, 14544, 16508, 18596, 20852, 23268, 25908, 28668, 31716, 34892, 38320, 41940, 45776, 49804, 54196, 58740, 63524, 68532, 73900

Offset: 1

Omar E. Pol, Jun 18 2014

nonn

a(n) is also the volume of a special stepped pyramid with n levels related to the symmetric representation of sigma. Note that starting at the top of the pyramid, the total area of the horizontal regions at the n-th level is equal to A239050(n), and the total area of the vertical regions at the n-th level is equal to 8*n.
From Omar E. Pol, Sep 19 2015: (Start)
Also, consider that the area of the central square in the top of the pyramid is equal to 1, so the total area of the horizontal regions at the n-th level starting from the top is equal to sigma(n) = A000203(n), and the total area of the vertical regions at the n-th level is equal to 2*n.
Also note that this stepped pyramid can be constructed with four copies of the stepped pyramid described in A245092 back-to-back (one copy in every quadrant). (End)
From Omar E. Pol, Jan 20 2021: (Start)
Convolution of A000203 and the nonzero terms of A008586.
Convolution of A074400 and the nonzero terms of A005843.
Convolution of A340793 and the nonzero terms of A046092.
Convolution of A239050 and A000027.
(End)

			From _Omar E. Pol_, Aug 29 2015: (Start)
Illustration of the top view of the stepped pyramid with 16 levels. The pyramid is formed of 5104 unit cubes:
.                 _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
.                |  _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _  |
.                | |_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _| |
.             _ _| |  _ _ _ _ _ _ _ _ _ _ _ _ _ _  | |_ _
.           _|  _ _| |_ _ _ _ _ _ _ _ _ _ _ _ _ _| |_ _  |_
.         _|  _|  _| |  _ _ _ _ _ _ _ _ _ _ _ _  | |_  |_  |_
.        |  _|   |_ _| |_ _ _ _ _ _ _ _ _ _ _ _| |_ _|   |_  |
.   _ _ _| |  _ _|     |  _ _ _ _ _ _ _ _ _ _  |     |_ _  | |_ _ _
.  |  _ _ _|_| |      _| |_ _ _ _ _ _ _ _ _ _| |_      | |_|_ _ _  |
.  | | |  _ _ _|    _|_ _|  _ _ _ _ _ _ _ _  |_ _|_    |_ _ _  | | |
.  | | | | |  _ _ _| |  _| |_ _ _ _ _ _ _ _| |_  | |_ _ _  | | | | |
.  | | | | | | |  _ _|_|  _|  _ _ _ _ _ _  |_  |_|_ _  | | | | | | |
.  | | | | | | | | |  _ _|   |_ _ _ _ _ _|   |_ _  | | | | | | | | |
.  | | | | | | | | | | |  _ _|  _ _ _ _  |_ _  | | | | | | | | | | |
.  | | | | | | | | | | | | |  _|_ _ _ _|_  | | | | | | | | | | | | |
.  | | | | | | | | | | | | | | |  _ _  | | | | | | | | | | | | | | |
.  | | | | | | | | | | | | | | | |   | | | | | | | | | | | | | | | |
.  | | | | | | | | | | | | | | | |_ _| | | | | | | | | | | | | | | |
.  | | | | | | | | | | | | | |_|_ _ _ _|_| | | | | | | | | | | | | |
.  | | | | | | | | | | | |_|_  |_ _ _ _|  _|_| | | | | | | | | | | |
.  | | | | | | | | | |_|_    |_ _ _ _ _ _|    _|_| | | | | | | | | |
.  | | | | | | | |_|_ _  |_  |_ _ _ _ _ _|  _|  _ _|_| | | | | | | |
.  | | | | | |_|_ _  | |_  |_ _ _ _ _ _ _ _|  _| |  _ _|_| | | | | |
.  | | | |_|_ _    |_|_ _| |_ _ _ _ _ _ _ _| |_ _|_|    _ _|_| | | |
.  | |_|_ _ _  |     |_  |_ _ _ _ _ _ _ _ _ _|  _|     |  _ _ _|_| |
.  |_ _ _  | |_|_      | |_ _ _ _ _ _ _ _ _ _| |      _|_| |  _ _ _|
.        | |_    |_ _  |_ _ _ _ _ _ _ _ _ _ _ _|  _ _|    _| |
.        |_  |_  |_  | |_ _ _ _ _ _ _ _ _ _ _ _| |  _|  _|  _|
.          |_  |_ _| |_ _ _ _ _ _ _ _ _ _ _ _ _ _| |_ _|  _|
.            |_ _  | |_ _ _ _ _ _ _ _ _ _ _ _ _ _| |  _ _|
.                | |_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _| |
.                | |_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _| |
.                |_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _|
.
Note that the above diagram contains a hidden pattern, simpler, which emerges from the front view of every corner of the stepped pyramid.
For more information about the hidden pattern see A237593 and A245092.
(End)

Robert G. Wilson v, Table of n, a(n) for n = 1..10000 (first 7342 terms from Robert Price)
Omar E. Pol, Illustration of a(11) = 1736, Perspective view of the stepped pyramid with 11 levels which contains 1736 unit cubes.

Cf. A000203, A024916, A046092, A008586, A175254, A196020, A235791, A236104, A237048, A237270, A237271, A237591, A237593, A239050, A239660, A239931, A239932, A239933, A239934, A243980, A245092, A262626, A340793.

[4*(&+[(n-k+1)*DivisorSigma(1,k): k in [1..n]]): n in [1..40]]; // G. C. Greubel, Apr 07 2019

a[n_] := 4 Sum[(n - k + 1) DivisorSigma[1, k], {k, n}]; Array[a, 40] (* Robert G. Wilson v, Aug 06 2018 *)
Nest[Accumulate,4*DivisorSigma[1,Range[50]],2] (* Harvey P. Dale, Sep 07 2022 *)

a(n) = 4*sum(k=1, n, sigma(k)*(n-k+1)); \\ Michel Marcus, Aug 07 2018

from math import isqrt
def A244050(n): return (((s:=isqrt(n))**2*(s+1)*((s+1)*((s<<1)+1)-6*(n+1))>>1) + sum((q:=n//k)*(-k*(q+1)*(3*k+(q<<1)+1)+3*(n+1)*((k<<1)+q+1)) for k in range(1,s+1))<<1)//3 # Chai Wah Wu, Oct 22 2023

[4*sum(sigma(k)*(n-k+1) for k in (1..n)) for n in (1..40)] # G. C. Greubel, Apr 07 2019

a(n) = 4*A175254(n).

A138879 Sum of all parts of the last section of the set of partitions of n.

Original entry on oeis.org

1, 3, 5, 11, 15, 31, 39, 71, 94, 150, 196, 308, 389, 577, 750, 1056, 1353, 1881, 2380, 3230, 4092, 5412, 6821, 8935, 11150, 14386, 17934, 22834, 28281, 35735, 43982, 55066, 67551, 83821, 102365, 126267, 153397, 188001, 227645, 277305, 334383

Offset: 1

Omar E. Pol, Apr 30 2008

nonn

Row sums of the triangles A135010, A138121, A138151 and others related to the section model of partitions (see A135010 and A138121).
From Omar E. Pol, Jan 20 2021: (Start)
Convolution of A000203 and A002865.
Convolution of A340793 and A000041.
Row sums of triangles A339278, A340426, A340583. (End)
a(n) is also the sum of all divisors of all terms of n-th row of A336811. These divisors are also all parts in the last section of the set of partitions of n. - Omar E. Pol, Jul 27 2021
Row sums of A336812. - Omar E. Pol, Aug 03 2021

			a(6)=31 because the parts of the last section of the set of partitions of 6 are (6), (3,3), (4,2), (2,2,2), (1), (1), (1), (1), (1), (1), (1), so the sum is a(6) = 6 + 3 + 3 + 4 + 2 + 2 + 2 + 2 + 1 + 1 + 1 + 1 + 1 + 1 + 1 = 31.
From _Omar E. Pol_, Aug 13 2013: (Start)
Illustration of initial terms:
.                                           _ _ _ _ _ _
.                                          |_ _ _ _ _ _|
.                                          |_ _ _|_ _ _|
.                                          |_ _ _ _|_ _|
.                               _ _ _ _ _  |_ _|_ _|_ _|
.                              |_ _ _ _ _|           |_|
.                     _ _ _ _  |_ _ _|_ _|           |_|
.                    |_ _ _ _|         |_|           |_|
.             _ _ _  |_ _|_ _|         |_|           |_|
.       _ _  |_ _ _|       |_|         |_|           |_|
.   _  |_ _|     |_|       |_|         |_|           |_|
.  |_|   |_|     |_|       |_|         |_|           |_|
.
.   1    3      5        11         15           31
.
(End)
On the other hand for n = 6 the 6th row of triangle A336811 is [6, 4, 3, 2, 2, 1, 1] and the sum of all divisors of these terms is [1 + 2 + 3 + 6] + [1 + 2 + 4] + [1 + 3] + [1 + 2] + [1 + 2] + [1] + [1] = 31, so a(6) = 31. - _Omar E. Pol_, Jul 27 2021

Vaclav Kotesovec, Table of n, a(n) for n = 1..10000

Cf. A000041, A000203, A002865, A066186, A133041, A135010, A138121, A138135 - A138138, A138151, A138880, A139100, A237593, A336811, A336812, A338156, A339278, A340035, A340426, A340583, A340793.

A066186 := proc(n) n*combinat[numbpart](n) ; end proc:
A138879 := proc(n) A066186(n)-A066186(n-1) ; end proc:
seq(A138879(n),n=1..80) ; # R. J. Mathar, Jan 27 2011

Table[PartitionsP[n]*n - PartitionsP[n-1]*(n-1), {n, 1, 50}] (* Vaclav Kotesovec, Oct 21 2016 *)

for(n=1, 50, print1(numbpart(n)*n - numbpart(n - 1)*(n - 1),", ")) \\ Indranil Ghosh, Mar 19 2017

from sympy.ntheory import npartitions
print([npartitions(n)*n - npartitions(n - 1)*(n - 1) for n in range(1, 51)]) # Indranil Ghosh, Mar 19 2017

a(n) = A000041(n)*n - A000041(n-1)*(n-1) = A138880(n) + A000041(n-1).
a(n) = A066186(n) - A066186(n-1), for n>=1.
a(n) ~ exp(Pi*sqrt(2*n/3)) * Pi/(12*sqrt(2*n)) * (1 - (72 + 13*Pi^2) / (24*Pi*sqrt(6*n)) + (7/12 + 3/(2*Pi^2) + 217*Pi^2/6912)/n - (15*sqrt(3/2)/(16*Pi) + 115*Pi/(288*sqrt(6)) + 4069*Pi^3/(497664*sqrt(6)))/n^(3/2)). - Vaclav Kotesovec, Oct 21 2016, extended Jul 06 2019
G.f.: x*(1 - x)*f'(x), where f(x) = Product_{k>=1} 1/(1 - x^k). - Ilya Gutkovskiy, Apr 13 2017

a(34) corrected by R. J. Mathar, Jan 27 2011

A175254 a(n) = Sum_{k<=n} A000203(k)*(n-k+1), where A000203(m) is the sum of divisors of m.

Original entry on oeis.org

1, 5, 13, 28, 49, 82, 123, 179, 248, 335, 434, 561, 702, 867, 1056, 1276, 1514, 1791, 2088, 2427, 2798, 3205, 3636, 4127, 4649, 5213, 5817, 6477, 7167, 7929, 8723, 9580, 10485, 11444, 12451, 13549, 14685, 15881, 17133, 18475, 19859, 21339, 22863, 24471, 26157

Offset: 1

Jaroslav Krizek, Mar 14 2010

Partial sums of A024916. - Omar E. Pol, Jul 03 2014
a(n) is also the volume of the stepped pyramid with n levels described in A245092. - Omar E. Pol, Aug 12 2015
Also the alternating row sums of A262612. - Omar E. Pol, Nov 23 2015
From Omar E. Pol, Jan 20 2021: (Start)
Convolution of A000203 and A000027.
Convolution of A340793 and the nonzero terms of A000217.
Antidiagonal sums of A319073.
Row sums of A274824. (End)
Row sums of A345272. - Omar E. Pol, Jun 14 2021
Also the alternating row sums of A353690. - Omar E. Pol, Jun 05 2022

			For n = 4: a(4) = sigma(1)*4 + sigma(2)*3 + sigma(3)*2 + sigma(4)*1 = 1*4 + 3*3 + 4*2 + 7*1 = 28.

Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 2209 terms from Indranil Ghosh)

Cf. A000203, A000217, A006218, A024916, A072481, A237593, A244050, A245092, A262612, A274824, A319073, A340793, A345272, A353690.

Cf. A143128, A353908.

b:= proc(n) option remember; `if`(n<1, [0$2],
      (p-> p+[numtheory[sigma](n), p[1]])(b(n-1)))
    end:
a:= n-> b(n+1)[2]:
seq(a(n), n=1..45);  # Alois P. Heinz, Oct 07 2021

Table[Sum[DivisorSigma[1, k] (n - k + 1), {k, n}], {n, 45}] (* Michael De Vlieger, Nov 24 2015 *)

a(n) = sum(x=1, n, sigma(x)*(n-x+1)) \\ Michel Marcus, Mar 18 2013

from math import isqrt
def A175254(n): return (((s:=isqrt(n))**2*(s+1)*((s+1)*(2*s+1)-6*(n+1))>>1) + sum((q:=n//k)*(-k*(q+1)*(3*k+2*q+1)+3*(n+1)*(2*k+q+1)) for k in range(1,s+1)))//6 # Chai Wah Wu, Oct 21 2023

Conjecture: a(n) = Sum_{k=0..n} A006218(n-k). - R. J. Mathar, Oct 17 2012
a(n) = A000330(n) - A072481(n). - Omar E. Pol, Aug 12 2015
a(n) ~ Pi^2*n^3/36. - Vaclav Kotesovec, Sep 25 2016
G.f.: (1/(1 - x)^2)*Sum_{k>=1} k*x^k/(1 - x^k). - Ilya Gutkovskiy, Jan 03 2017
a(n) = Sum_{k=1..n} Sum_{i=1..k} k - (k mod i). - Wesley Ivan Hurt, Sep 13 2017
a(n) = A244050(n)/4. - Omar E. Pol, Jan 22 2021
a(n) = (n+1)*A024916(n) - A143128(n). - Vaclav Kotesovec, May 11 2022

Corrected by Jaroslav Krizek, Mar 17 2010

More terms from Michel Marcus, Mar 18 2013

A000385 Convolution of A000203 with itself.

Original entry on oeis.org

1, 6, 17, 38, 70, 116, 185, 258, 384, 490, 686, 826, 1124, 1292, 1705, 1896, 2491, 2670, 3416, 3680, 4602, 4796, 6110, 6178, 7700, 7980, 9684, 9730, 12156, 11920, 14601, 14752, 17514, 17224, 21395, 20406, 24590, 24556, 28920, 27860, 34112, 32186, 38674, 37994, 43980, 42136, 51646, 47772, 56749, 55500, 64316, 60606, 73420, 67956, 80500, 77760, 88860, 83810, 102284, 92690, 108752, 105236, 120777, 112672, 135120, 123046, 145194, 138656, 157512, 146580, 177515, 159396, 185744, 179122

Offset: 1

N. J. A. Sloane

Convolution of A340793 and A024916. - Omar E. Pol, Feb 17 2021

N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Jacques Touchard, On prime numbers and perfect numbers, Scripta Math., 129 (1953), 35-39.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Nicolae Ciprian Bonciocat, Congruences for the Convolution of Divisor sum function, Bull. Greek Math. Soc., Vol 47 (2003), pp. 19-29.
MathOverflow, Searching for a proof for a series identity.
Srinivasa Ramanujan, On certain arithmetical functions, Transactions of the Cambridge Philosophical Society, 22, No.9 (1916), 169- 184 (see Table IV, line 1).
Jacques Touchard, On prime numbers and perfect numbers, Scripta Math., 129 (1953), 35-39. [Annotated scanned copy]

Cf. A000203, A024916, A001158, A340793.

Column k=2 of A319083 (shifted).

a000385 n = sum $ zipWith (*) sigmas $ reverse sigmas where
   sigmas = take n a000203_list
-- Reinhard Zumkeller, Sep 20 2011

f:= n -> 5/12*numtheory:-sigma[3](n+1)-(5+6*n)/12*numtheory:-sigma(n+1):
map(f, [$1..100]); # Robert Israel, Sep 17 2018

a[n_] := Sum[DivisorSigma[1, k] DivisorSigma[1, n-k+1], {k, 1, n}];
Array[a, 100] (* Jean-François Alcover, Aug 01 2018 *)

a(n) = sum(k=1, n, sigma(k)*sigma(n-k+1)); \\ Michel Marcus, Nov 10 2016

a(n) = my(f = factor(n+1)); (5 * sigma(f, 3) - (6*n + 5) * sigma(f)) / 12; \\ Amiram Eldar, Jan 04 2025

from sympy import factorint
def A000385(n):
    f = factorint(n+1).items()
    return(5*prod((p**(3*(e+1))-1)//(p**3-1) for p,e in f)-(5+6*n)*prod((p**(e+1)-1)//(p-1) for p, e in f))//12 # Chai Wah Wu, Jul 25 2024

a(n) = Sum_{k=1..n} A000203(k)*A000203(n-k+1).
G.f.: (1/x)*(Sum_{k>=1} k*x^k/(1 - x^k))^2. - Ilya Gutkovskiy, Nov 10 2016
a(5*n+1)==0 (mod 5) and a(7*n+6)==0 (mod 7). See Bonciocat link. - Michel Marcus, Nov 10 2016
a(n) = (5/12)*A001158(n+1) - ((5+6*n)/12)*A000203(n+1). - Robert Israel, Sep 17 2018
Sum_{k=1..n} a(k) ~ Pi^4 * n^4 / 864. - Vaclav Kotesovec, Apr 02 2019

More terms from Sean A. Irvine, Nov 14 2010

A182738 Partial sums of A066186.

Original entry on oeis.org

1, 5, 14, 34, 69, 135, 240, 416, 686, 1106, 1722, 2646, 3959, 5849, 8489, 12185, 17234, 24164, 33474, 46014, 62646, 84690, 113555, 151355, 200305, 263641, 344911, 449015, 581400, 749520, 961622, 1228790, 1563509, 1982049

Offset: 1

Omar E. Pol, Jan 22 2011

a(n) is also the volume of a three-dimensional version of the section model of partitions: the 3D illustrations in A135010 show boxes with face areas of 1 X 1, 2 X 2, 3 X 3, 4 X 5, 5 X 7 units along the m and p(m) axis, which is sequence A066186. Assuming that the boxes are 1 unit deep, the total volume of all boxes up to layer n is a(n). See the first two links.
From Omar E. Pol, Jan 20 2021: (Start)
a(n) is the sum of all parts of all partitions of all positive integers <= n.
Convolution of A000203 and A000070.
Convolution of A024916 and A000041.
Convolution of A175254 and A002865.
Convolution of A340793 and A014153.
Row sums of triangles A340527, A340531, A340579.
Consider a symmetric tower (a polycube) in which the terraces are the symmetric representation of sigma (n..1) respectively starting from the base (cf. A237270, A237593). The total area of the terraces equals A024916(n), the same as the area of the base.
The levels of the terraces starting from the base are the first n terms of A000070, that is A000070(0)..A000070(n-1), hence the differences between two successive levels give the partition numbers A000041, that is A000041(0)..A000041(n-1).
a(n) is the volume (or the total number of unit cubes) of the polycube.
That is due to the correspondence between divisors and partitions (cf. A336811).
The symmetric tower is a member of the family of the pyramid described in A245092.
The growth of the volume of the polycube represents every convolution mentioned above. (End)

			a(6) = 135 because the volume V(6) = p(1) + 2*p(2) + 3*p(3) + 4*p(4) + 5*p(5) + 6*p(6) = 1 + 2*2 + 3*3 + 4*5 + 5*7 + 6*11 = 1 + 4 + 9 + 20 + 35 + 66 = 135 where p(n) = A000041(n).

Vaclav Kotesovec, Table of n, a(n) for n = 1..10000
Omar E. Pol, Illustration of a(6) = 135, the polycube contains 135 unit cubes.
Omar E. Pol, Illustration of a(9) = 686, the polycube contains 686 unit cubes.

Cf. A000041, A000070, A000203, A002865, A014153, A024916, A066186, A135010, A175254, A237270, A237593, A245092, A336811, A340527, A340531, A340579, A340793.

With[{no=35},Accumulate[PartitionsP[Range[no]]Range[no]]] (* Harvey P. Dale, Feb 02 2011 *)

a(n) = n*A000070(n) - A014153(n-1). - Vaclav Kotesovec, Jun 23 2015
a(n) ~ sqrt(n) * exp(Pi*sqrt(2*n/3)) / (Pi*2^(3/2)) * (1 + (11*Pi/(24*sqrt(6)) - sqrt(6)/Pi)/sqrt(n) + (73*Pi^2/6912 - 3/16)/n). - Vaclav Kotesovec, Jun 23 2015, extended Nov 04 2016
G.f.: x*f'(x)/(1 - x), where f(x) = Product_{k>=1} 1/(1 - x^k). - Ilya Gutkovskiy, Apr 10 2017

A000441 a(n) = Sum_{k=1..n-1} ksigma(k)sigma(n-k).

Original entry on oeis.org

0, 1, 9, 34, 95, 210, 406, 740, 1161, 1920, 2695, 4116, 5369, 7868, 9690, 13640, 16116, 22419, 25365, 34160, 38640, 50622, 55154, 73320, 77225, 100100, 107730, 135576, 141085, 182340, 184760, 233616, 243408, 297738, 301420, 385110, 377511, 467210, 478842

Offset: 1

N. J. A. Sloane

nonn

Apart from initial zero this is the convolution of A340793 and A143128. - Omar E. Pol, Feb 16 2021

N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Jacques Touchard, On prime numbers and perfect numbers, Scripta Math., 129 (1953), 35-39.

Vaclav Kotesovec, Table of n, a(n) for n = 1..1000
Jacques Touchard, On prime numbers and perfect numbers, Scripta Math., 129 (1953), 35-39. [Annotated scanned copy]

Cf. A000385, A000477, A000499, A259692, A259693, A259694, A259695, A259696.
Cf. A000203 (sigma_1), A001158 (sigma_3), A064987.
Cf. A143128, A340793.

S:=(n,e)->add(k^e*sigma(k)*sigma(n-k),k=1..n-1);
f:=e->[seq(S(n,e),n=1..30)];f(1); # N. J. A. Sloane, Jul 03 2015

a[n_] := Sum[k*DivisorSigma[1, k]*DivisorSigma[1, n-k], {k, 1, n-1}]; Array[a, 40] (* Jean-François Alcover, Feb 08 2016 *)

a(n) = sum(k=1, n-1, k*sigma(k)*sigma(n-k)); \\ Michel Marcus, Feb 02 2014

a(n) = my(f = factor(n)); ((n - 6*n^2) * sigma(f) + 5*n * sigma(f, 3)) / 24; \\ Amiram Eldar, Jan 04 2025

from sympy import divisor_sigma
def A000441(n): return (n*(1-6*n)*divisor_sigma(n)+5*n*divisor_sigma(n,3))//24 # Chai Wah Wu, Jul 25 2024

Convolution of A000203 with A064987. - Sean A. Irvine, Nov 14 2010
G.f.: x*f(x)*f'(x), where f(x) = Sum_{k>=1} k*x^k/(1 - x^k). - Ilya Gutkovskiy, Apr 28 2018
a(n) = (n/24 - n^2/4)*sigma_1(n) + (5*n/24)*sigma_3(n). - Ridouane Oudra, Sep 17 2020
Sum_{k=1..n} a(k) ~  Pi^4 * n^5 / 2160. - Vaclav Kotesovec, May 09 2022

More terms from Sean A. Irvine, Nov 14 2010

a(1)=0 prepended by Michel Marcus, Feb 02 2014

A341062 Sequence whose partial sums give A000005.

Original entry on oeis.org

1, 1, 0, 1, -1, 2, -2, 2, -1, 1, -2, 4, -4, 2, 0, 1, -3, 4, -4, 4, -2, 0, -2, 6, -5, 1, 0, 2, -4, 6, -6, 4, -2, 0, 0, 5, -7, 2, 0, 4, -6, 6, -6, 4, 0, -2, -2, 8, -7, 3, -2, 2, -4, 6, -4, 4, -4, 0, -2, 10, -10, 2, 2, 1, -3, 4, -6, 4, -2, 4, -6, 10, -10, 2, 2, 0, -2, 4, -6, 8, -5, -1, -2, 10, -8, 0, 0, 4, -6, 10

Offset: 1

Omar E. Pol, Feb 04 2021

Essentially a duplicate of A051950.
Convolved with A000041 gives A138137.
Convolved with A000027 gives the nonzero terms of A006218.
Convolved with A000070 gives the nonzero terms of A006128.
Convolved with A014153 gives the nonzero terms of A284870.
Convolved with A036469 gives the nonzero terms of A305082.
Convolved with the nonzero terms of A006218 gives A055507.
Convolved with the nonzero terms of A000217 gives the nonzero terms of A078567.

1 together with A051950.
Cf. A000005 (partial sums).
Cf. A000027, A000041, A000070, A000217, A006128, A006218, A014153, A036469, A055507, A078567, A138137, A284870, A305082, A340793.

Join[{1}, Differences[Table[DivisorSigma[0, n], {n, 1, 90}]]] (* Amiram Eldar, Feb 06 2021 *)

a(n) = A051950(n) for n > 1.

A191831 a(n) = Sum_{i+j=n, i,j >= 1} tau(i)*sigma(j), where tau() = A000005(), sigma() = A000203().

Original entry on oeis.org

0, 1, 5, 12, 24, 39, 60, 87, 113, 158, 189, 249, 286, 372, 402, 516, 545, 696, 709, 886, 912, 1125, 1110, 1401, 1348, 1674, 1654, 1992, 1906, 2390, 2226, 2735, 2648, 3141, 2926, 3705, 3346, 4069, 3898, 4604, 4223, 5282, 4707, 5757, 5426, 6326, 5754, 7269, 6324, 7669, 7230, 8468, 7556, 9456, 8240, 10018, 9320, 10748, 9621, 12246

Offset: 1

N. J. A. Sloane, Jun 17 2011

This is Andrews's D_{0,1}(n).
From Omar E. Pol, Dec 08 2021: (Start)
Zero together with the convolution of A000005 and A000203.
Zero together with the convolution of A341062 and A024916.
Zero together with the convolution of the nonzero terms of A006218 and A340793.
a(n) is also the volume of a symmetric polycube which belongs to the family of symmetric polycubes that represent the convolution of A000203 with any other integer sequence, n >= 1. (End)

George E. Andrews, Stacked lattice boxes, Ann. Comb. 3 (1999), 115-130.

Cf. A000005, A000203 A006218, A024916, A086718, A237593, A340793, A341062.

with(numtheory); D01:=n->add(tau(j)*sigma(n-j),j=1..n-1);
[seq(D01(n),n=1..60)];

Table[Sum[DivisorSigma[0, j] DivisorSigma[1, n - j], {j, n - 1}], {n, 60}] (* Michael De Vlieger, Jan 01 2017 *)

a(n)=sum(i=1,n-1,numdiv(i)*sigma(n-i)) \\ Charles R Greathouse IV, Feb 19 2013

G.f.: (Sum_{k>=1} x^k/(1 - x^k))*(Sum_{k>=1} k*x^k/(1 - x^k)). - Ilya Gutkovskiy, Jan 01 2017

cp's OEIS Frontend

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).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

Python

Python

Sage

Formula

A066186 Sum of all parts of all partitions of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Sage

Formula

Extensions

A244050 Partial sums of A243980.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Python

Sage

Formula

A138879 Sum of all parts of the last section of the set of partitions of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A175254 a(n) = Sum_{k<=n} A000203(k)*(n-k+1), where A000203(m) is the sum of divisors of m.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

A000441 a(n) = Sum_{k=1..n-1} ksigma(k)sigma(n-k).