cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 31-40 of 228 results. Next

A221529 Triangle read by rows: T(n,k) = A000203(k)*A000041(n-k), 1 <= k <= n.

Original entry on oeis.org

1, 1, 3, 2, 3, 4, 3, 6, 4, 7, 5, 9, 8, 7, 6, 7, 15, 12, 14, 6, 12, 11, 21, 20, 21, 12, 12, 8, 15, 33, 28, 35, 18, 24, 8, 15, 22, 45, 44, 49, 30, 36, 16, 15, 13, 30, 66, 60, 77, 42, 60, 24, 30, 13, 18, 42, 90, 88, 105, 66, 84, 40, 45, 26, 18, 12, 56, 126, 120, 154, 90, 132, 56, 75, 39, 36, 12, 28
Offset: 1

Views

Author

Omar E. Pol, Jan 20 2013

Keywords

Comments

Since A000203(k) has a symmetric representation, both T(n,k) and the partial sums of row n can be represented by symmetric polycubes. For more information see A237593 and A237270. For another version see A245099. - Omar E. Pol, Jul 15 2014
From Omar E. Pol, Jul 10 2021: (Start)
The above comment refers to a symmetric tower whose terraces are the symmetric representation of sigma(i), for i = 1..n, starting from the top. The levels of these terraces are the partition numbers A000041(h-1), for h = 1 to n, starting from the base of the tower, where n is the length of the largest side of the base.
The base of the tower is the symmetric representation of A024916(n).
The height of the tower is equal to A000041(n-1).
The surface area of the tower is equal to A345023(n).
The volume (or the number of cubes) of the tower equals A066186(n).
The volume represents the n-th term of the convolution of A000203 and A000041, that is A066186(n).
Note that the terraces that are the symmetric representation of sigma(n) and the terraces that are the symmetric representation of sigma(n-1) both are unified in level 1 of the structure. That is because the first two partition numbers A000041 are [1, 1].
The tower is an object of the family of the stepped pyramid described in A245092.
T(n,k) can be represented with a set of A237271(k) right prisms of height A000041(n-k) since T(n,k) is the total number of cubes that are exactly below the parts of the symmetric representation of sigma(k) in the tower.
T(n,k) is also the sum of all divisors of all k's that are in the first n rows of triangle A336811, or in other words, in the first A000070(n-1) terms of the sequence A336811. Hence T(n,k) is also the sum of all divisors of all k's in the n-th row of triangle A176206.
The mentioned property is due to the correspondence between divisors and parts explained in A338156: all divisors of the first A000070(n-1) terms of A336811 are also all parts of all partitions of n.
Therefore the set of all partitions of n >= 1 has an associated tower.
The partial column sums of A340583 give this triangle showing the growth of the structure of the tower.
Note that the convolution of A000203 with any integer sequence S can be represented with a symmetric tower or structure of the same family where its terraces are the symmetric representation of sigma starting from the top and the heights of the terraces starting from the base are the terms of the sequence S. (End)

Examples

			Triangle begins:
------------------------------------------------------
    n| k    1   2   3   4   5   6   7   8   9  10
------------------------------------------------------
    1|      1;
    2|      1,  3;
    3|      2,  3,  4;
    4|      3,  6,  4,  7;
    5|      5,  9,  8,  7,  6;
    6|      7, 15, 12, 14,  6, 12;
    7|     11, 21, 20, 21, 12, 12,  8;
    8|     15, 33, 28, 35, 18, 24,  8, 15;
    9|     22, 45, 44, 49, 30, 36, 16, 15, 13;
   10|     30, 66, 60, 77, 42, 60, 24, 30, 13, 18;
...
The sum of row 10 is [30 + 66 + 60 + 77 + 42 + 60 + 24 + 30 + 13 + 18] = A066186(10) = 420.
.
For n = 10 the calculation of the row 10 is as follows:
    k    A000203         T(10,k)
    1       1   *  30   =   30
    2       3   *  22   =   66
    3       4   *  15   =   60
    4       7   *  11   =   77
    5       6   *   7   =   42
    6      12   *   5   =   60
    7       8   *   3   =   24
    8      15   *   2   =   30
    9      13   *   1   =   13
   10      18   *   1   =   18
                 A000041
.
From _Omar E. Pol_, Jul 13 2021: (Start)
For n = 10 we can see below three views of two associated polycubes called here "prism of partitions" and "tower". Both objects contain the same number of cubes (that property is valid for n >= 1).
        _ _ _ _ _ _ _ _ _ _
  42   |_ _ _ _ _          |
       |_ _ _ _ _|_        |
       |_ _ _ _ _ _|_      |
       |_ _ _ _      |     |
       |_ _ _ _|_ _ _|_    |
       |_ _ _ _        |   |
       |_ _ _ _|_      |   |
       |_ _ _ _ _|_    |   |
       |_ _ _      |   |   |
       |_ _ _|_    |   |   |
       |_ _    |   |   |   |
       |_ _|_ _|_ _|_ _|_  |                             _
  30   |_ _ _ _ _        | |                            | | 30
       |_ _ _ _ _|_      | |                            | |
       |_ _ _      |     | |                            | |
       |_ _ _|_ _ _|_    | |                            | |
       |_ _ _ _      |   | |                            | |
       |_ _ _ _|_    |   | |                            | |
       |_ _ _    |   |   | |                            | |
       |_ _ _|_ _|_ _|_  | |                           _|_|
  22   |_ _ _ _        | | |                          |   |  22
       |_ _ _ _|_      | | |                          |   |
       |_ _ _ _ _|_    | | |                          |   |
       |_ _ _      |   | | |                          |   |
       |_ _ _|_    |   | | |                          |   |
       |_ _    |   |   | | |                          |   |
       |_ _|_ _|_ _|_  | | |                         _|_ _|
  15   |_ _ _ _      | | | |                        | |   |  15
       |_ _ _ _|_    | | | |                        | |   |
       |_ _ _    |   | | | |                        | |   |
       |_ _ _|_ _|_  | | | |                       _|_|_ _|
  11   |_ _ _      | | | | |                      | |     |  11
       |_ _ _|_    | | | | |                      | |     |
       |_ _    |   | | | | |                      | |     |
       |_ _|_ _|_  | | | | |                     _| |_ _ _|
   7   |_ _ _    | | | | | |                    |   |     |   7
       |_ _ _|_  | | | | | |                   _|_ _|_ _ _|
   5   |_ _    | | | | | | |                  | | |       |   5
       |_ _|_  | | | | | | |                 _| | |_ _ _ _|
   3   |_ _  | | | | | | | |               _|_ _|_|_ _ _ _|   3
   2   |_  | | | | | | | | |           _ _|_ _|_|_ _ _ _ _|   2
   1   |_|_|_|_|_|_|_|_|_|_|          |_ _|_|_|_ _ _ _ _ _|   1
.
             Figure 1.                       Figure 2.
         Front view of the                 Lateral view
        prism of partitions.               of the tower.
.
.                                      _ _ _ _ _ _ _ _ _ _
                                      |   | | | | | | | |_|   1
                                      |   | | | | | |_|_ _|   2
                                      |   | | | |_|_  |_ _|   3
                                      |   | |_|_    |_ _ _|   4
                                      |   |_ _  |_  |_ _ _|   5
                                      |_ _    |_  |_ _ _ _|   6
                                          |_    | |_ _ _ _|   7
                                            |_  |_ _ _ _ _|   8
                                              |           |   9
                                              |_ _ _ _ _ _|  10
.
                                             Figure 3.
                                             Top view
                                           of the tower.
.
Figure 1 is a two-dimensional diagram of the partitions of 10 in colexicographic order (cf. A026792, A211992). The area of the diagram is 10*42 = A066186(10) = 420. Note that the diagram can be interpreted also as the front view of a right prism whose volume is 1*10*42 = 420 equaling the volume and the number of cubes of the tower that appears in the figures 2 and 3.
Note that the shape and the area of the lateral view of the tower are the same as the shape and the area where the 1's are located in the diagram of partitions. In this case the mentioned area equals A000070(10-1) = 97.
The connection between these two associated objects is a representation of the correspondence divisor/part described in A338156. See also A336812.
The sum of the volumes of both objects equals A220909.
For the connection with the table of A338156 see also A340035. (End)
		

Crossrefs

Programs

  • Mathematica
    nrows=12; Table[Table[DivisorSigma[1,k]PartitionsP[n-k],{k,n}],{n,nrows}] // Flatten (* Paolo Xausa, Jun 17 2022 *)
  • PARI
    T(n,k)=sigma(k)*numbpart(n-k) \\ Charles R Greathouse IV, Feb 19 2013

Formula

T(n,k) = sigma(k)*p(n-k) = A000203(k)*A027293(n,k).
T(n,k) = A245093(n,k)*A027293(n,k).

A296508 Irregular triangle read by rows: T(n,k) is the size of the subpart that is adjacent to the k-th peak of the largest Dyck path of the symmetric representation of sigma(n), or T(n,k) = 0 if the mentioned subpart is already associated to a previous peak or if there is no subpart adjacent to the k-th peak, with n >= 1, k >= 1.

Original entry on oeis.org

1, 3, 2, 2, 7, 0, 3, 3, 11, 1, 0, 4, 0, 4, 15, 0, 0, 5, 3, 5, 9, 0, 9, 0, 6, 0, 0, 6, 23, 5, 0, 0, 7, 0, 0, 7, 12, 0, 12, 0, 8, 7, 1, 0, 8, 31, 0, 0, 0, 0, 9, 0, 0, 0, 9, 35, 2, 0, 2, 0, 10, 0, 0, 0, 10, 39, 0, 3, 0, 0, 11, 5, 0, 5, 0, 11, 18, 0, 0, 0, 18, 0, 12, 0, 0, 0, 0, 12, 47, 13, 0, 0, 0, 0, 13, 0, 5, 0, 0, 13
Offset: 1

Views

Author

Omar E. Pol, Feb 10 2018

Keywords

Comments

Conjecture: row n is formed by the odd-indexed terms of the n-th row of triangle A280850 together with the even-indexed terms of the same row but listed in reverse order. Examples: the 15th row of A280850 is [8, 8, 7, 0, 1] so the 15th row of this triangle is [8, 7, 1, 0, 8]. The 75th row of A280850 is [38, 38, 21, 0, 3, 3, 0, 0, 0, 21, 0] so the 75h row of this triangle is [38, 21, 3, 0, 0, 0, 21, 0, 3, 0, 38].
For the definition of "subparts" see A279387.
For more information about the mentioned Dyck paths see A237593.
T(n,k) could be called the "charge" of the k-th peak of the largest Dyck path of the symmetric representation of sigma(n).
The number of zeros in row n is A238005(n). - Omar E. Pol, Sep 11 2021

Examples

			Triangle begins (rows 1..28):
   1;
   3;
   2,  2;
   7,  0;
   3,  3;
  11,  1,  0;
   4,  0,  4;
  15,  0,  0;
   5,  3,  5;
   9,  0,  9,  0;
   6,  0,  0,  6;
  23,  5,  0,  0;
   7,  0,  0,  7;
  12,  0, 12,  0;
   8,  7,  1,  0,  8;
  31,  0,  0,  0,  0;
   9,  0,  0,  0,  9;
  35,  2,  0,  2,  0;
  10,  0,  0,  0, 10;
  39,  0,  3,  0,  0;
  11,  5,  0,  5,  0, 11;
  18,  0,  0,  0, 18,  0;
  12,  0,  0,  0,  0, 12;
  47, 13,  0,  0,  0,  0;
  13,  0,  5,  0,  0, 13;
  21,  0,  0,  0  21,  0;
  14,  6,  0,  6,  0, 14;
  55,  0,  0,  1,  0,  0,  0;
  ...
For n = 15 we have that the 14th row of triangle A237593 is [8, 3, 1, 2, 2, 1, 3, 8] and the 15th row of the same triangle is [8, 3, 2, 1, 1, 1, 1, 2, 3, 8], so the diagram of the symmetric representation of sigma(15) is constructed in the third quadrant as shown below in Figure 1:
.    _                                  _
.   | |                                | |
.   | |                                | |
.   | |                                | |
. 8 | |                                | |
.   | |                                | |
.   | |                                | |
.   | |                                | |
.   |_|_ _ _                           |_|_ _ _
.         | |_ _                      8      | |_ _
.         |_    |                            |_ _  |
.           |_  |_                          7  |_| |_
.          8  |_ _|                           1  |_ _|
.                 |                             0    |
.                 |_ _ _ _ _ _ _ _                   |_ _ _ _ _ _ _ _
.                 |_ _ _ _ _ _ _ _|                  |_ _ _ _ _ _ _ _|
.                         8                         8
.
.   Figure 1. The symmetric            Figure 2. After the dissection
.   representation of sigma(15)        of the symmetric representation
.   has three parts of size 8          of sigma(15) into layers of
.   because every part contains        width 1 we can see four subparts,
.   8 cells, so the 15th row of        so the 15th row of this triangle is
.   triangle A237270 is [8, 8, 8].     [8, 7, 1, 0, 8]. See also below.
.
Illustration of first 50 terms (rows 1..16 of triangle) in an irregular spiral which can be find in the top view of the pyramid described in A244050:
.
.               12 _ _ _ _ _ _ _ _
.                 |  _ _ _ _ _ _ _|_ _ _ _ _ _ _ 7
.                 | |             |_ _ _ _ _ _ _|
.              0 _| |                           |
.               |_ _|9 _ _ _ _ _ _              |_ _ 0
.         12 _ _|     |  _ _ _ _ _|_ _ _ _ _ 5      |_ 0
.    0 _ _ _| |    0 _| |         |_ _ _ _ _|         |
.     |  _ _ _|  9 _|_ _|                   |_ _ 3    |_ _ _ 7
.     | |    0 _ _| |   11 _ _ _ _          |_  |         | |
.     | |     |  _ _|  1 _|  _ _ _|_ _ _ 3    |_|_ _ 5    | |
.     | |     | |    0 _|_| |     |_ _ _|         | |     | |
.     | |     | |     |  _ _|           |_ _ 3    | |     | |
.     | |     | |     | |    3 _ _        | |     | |     | |
.     | |     | |     | |     |  _|_ 1    | |     | |     | |
.    _|_|    _|_|    _|_|    _|_| |_|    _|_|    _|_|    _|_|    _
.   | |     | |     | |     | |         | |     | |     | |     | |
.   | |     | |     | |     |_|_ _     _| |     | |     | |     | |
.   | |     | |     | |    2  |_ _|_ _|  _|     | |     | |     | |
.   | |     | |     |_|_     2    |_ _ _|  0 _ _| |     | |     | |
.   | |     | |    4    |_               7 _|  _ _|0    | |     | |
.   | |     |_|_ _     0  |_ _ _ _        |  _|    _ _ _| |     | |
.   | |    6      |_      |_ _ _ _|_ _ _ _| |  0 _|  _ _ _|0    | |
.   |_|_ _ _     0  |_   4        |_ _ _ _ _|  _|  _| |    _ _ _| |
.  8      | |_ _   0  |                     15|  _|  _|   |  _ _ _|
.         |_ _  |     |_ _ _ _ _ _            | |_ _|  0 _| |      0
.        7  |_| |_    |_ _ _ _ _ _|_ _ _ _ _ _| |    5 _|  _|
.          1  |_ _|  6            |_ _ _ _ _ _ _|  _ _|  _|  0
.            0    |                             23|  _ _|  0
.                 |_ _ _ _ _ _ _ _                | |    0
.                 |_ _ _ _ _ _ _ _|_ _ _ _ _ _ _ _| |
.                8                |_ _ _ _ _ _ _ _ _|
.                                                    31
.
The diagram contains 30 subparts equaling A060831(16), the total number of partitions of all positive integers <= 16 into consecutive parts.
For the construction of the spiral see A239660.
From _Omar E. Pol_, Nov 26 2020: (Start)
Also consider the infinite double-staircases diagram defined in A335616 (see the theorem). For n = 15 the diagram with first 15 levels looks like this:
.
Level                         "Double-staircases" diagram
.                                          _
1                                        _|1|_
2                                      _|1 _ 1|_
3                                    _|1  |1|  1|_
4                                  _|1   _| |_   1|_
5                                _|1    |1 _ 1|    1|_
6                              _|1     _| |1| |_     1|_
7                            _|1      |1  | |  1|      1|_
8                          _|1       _|  _| |_  |_       1|_
9                        _|1        |1  |1 _ 1|  1|        1|_
10                     _|1         _|   | |1| |   |_         1|_
11                   _|1          |1   _| | | |_   1|          1|_
12                 _|1           _|   |1  | |  1|   |_           1|_
13               _|1            |1    |  _| |_  |    1|            1|_
14             _|1             _|    _| |1 _ 1| |_    |_             1|_
15            |1              |1    |1  | |1| |  1|    1|              1|
.
Starting from A196020 and after the algorithm described n A280850 and the conjecture applied to the above diagram we have a new diagram as shown below:
.
Level                             "Ziggurat" diagram
.                                          _
6                                         |1|
7                            _            | |            _
8                          _|1|          _| |_          |1|_
9                        _|1  |         |1   1|         |  1|_
10                     _|1    |         |     |         |    1|_
11                   _|1      |        _|     |_        |      1|_
12                 _|1        |       |1       1|       |        1|_
13               _|1          |       |         |       |          1|_
14             _|1            |      _|    _    |_      |            1|_
15            |1              |     |1    |1|    1|     |              1|
.
The 15th row
of A249351:   [1,1,1,1,1,1,1,1,0,0,0,1,1,1,2,1,1,1,0,0,0,1,1,1,1,1,1,1,1]
The 15th row
of A237270:   [              8,            8,            8              ]
The 15th row
of this seq:  [              8,      7,    1,    0,      8              ]
The 15th row
of A280851:   [              8,      7,    1,            8              ]
.
(End)
		

Crossrefs

Row sums give A000203.
Row n has length A003056(n).
Column k starts in row A000217(k).
Nonzero terms give A280851.
The number of nonzero terms in row n is A001227(n).
The triangle with n rows contain A060831(n) nonzero terms.

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

Views

Author

Jaroslav Krizek, Mar 14 2010

Keywords

Comments

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

Examples

			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.
		

Crossrefs

Programs

  • Maple
    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
  • Mathematica
    Table[Sum[DivisorSigma[1, k] (n - k + 1), {k, n}], {n, 45}] (* Michael De Vlieger, Nov 24 2015 *)
  • PARI
    a(n) = sum(x=1, n, sigma(x)*(n-x+1)) \\ Michel Marcus, Mar 18 2013
    
  • Python
    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

Formula

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

Extensions

Corrected by Jaroslav Krizek, Mar 17 2010
More terms from Michel Marcus, Mar 18 2013

A153485 Sum of all aliquot divisors of all positive integers <= n.

Original entry on oeis.org

0, 1, 2, 5, 6, 12, 13, 20, 24, 32, 33, 49, 50, 60, 69, 84, 85, 106, 107, 129, 140, 154, 155, 191, 197, 213, 226, 254, 255, 297, 298, 329, 344, 364, 377, 432, 433, 455, 472, 522, 523, 577, 578, 618, 651, 677, 678, 754, 762, 805, 826
Offset: 1

Views

Author

Omar E. Pol, Dec 27 2008

Keywords

Comments

a(n) is also the sum of first n terms of A000203, minus n-th triangular number.
n is prime if and only if a(n) - a(n-1) = 1. - Omar E. Pol, Dec 31 2012
Also the alternating row sums of A236540. - Omar E. Pol, Jun 23 2014
Sum of the areas of all x X z rectangles with x and y integers, x + y = n, x <= y and z = floor(y/x). - Wesley Ivan Hurt, Dec 21 2020
Apart from the symmetric representation of a(n) given in the Example section we have that a(n) can be represented with an arrowhead-shaped polygon formed by two zig-zag paths and the Dyck path described in the n-th row of A237593 as shown in the Links section. - Omar E. Pol, Jun 13 2022

Examples

			Assuming that a(1) = 0, for n = 6 the aliquot divisors of the first six positive integers are [0], [1], [1], [1, 2], [1], [1, 2, 3], so a(6) = 0 + 1 + 1 + 1 + 2 + 1 + 1 + 2 + 3 = 12.
From _Omar E. Pol_, Mar 27 2021: (Start)
The following diagrams show a square dissection into regions that are the symmetric representation of A000203, A004125, A244048 and this sequence.
In order to construct every diagram we use the following rules:
At stage 1 in the first quadrant of the square grid we draw the symmetric representation of sigma(n) using the two Dyck paths described in the rows n and n-1 of A237593.
At stage 2 we draw a pair of orthogonal line segments (if it's necessary) such that in the drawing appears totally formed a square n X n. The area of the region that is above the symmetric representation of sigma(n) equals A004125(n).
At stage 3 we draw a zig-zag path with line segments of length 1 from (0,n-1) to (n-1,0) such that appears a staircase with n-1 steps. The area of the region (or regions) that is below the symmetric representation of sigma(n) and above the staircase equals A244048(n).
At stage 4 we draw a copy of the symmetric representation of A004125(n) rotated 180 degrees such that one of its vertices is the point (0,0). a(n) is the area of the region (or regions) that is above of this region and below the staircase.
Illustration for n = 1..6:
.                                                                    _ _ _ _ _ _
.                                                     _ _ _ _ _     |_ _ _  |_ R|
.                                        _ _ _ _ R   |_ _S_|  R|    | |_T | S |_|
.                             _ _ _ R   |_ _  |_|    | |_  |_ _|    |   |_|_ _  |
.                    _ _     |_S_|_|    | |_|_S |    |_U_|_T | |    |_  U |_T | |
.             _ S   |_ S|   U|_|_|S|    |_ U|_| |    |   | |_|S|    | |_    |_| |
.            |_|    |_|_|    |_|_|_|    |_|_ _|_|    |_V_|_U_|_|    |_V_|_ _ _|_|
.                  U        V   U       V
.
n:            1       2         3           4             5               6
R: A004125    0       0         1           1             4               3
S: A000203    1       3         4           7             6              12
T: A244048    0       0         1           2             5               6
U: a(n)       0       1         2           5             6              12
V: A004125    0       0         1           1             4               3
.
Illustration for n = 7..9:
.                                                      _ _ _ _ _ _ _ _ _
.                                _ _ _ _ _ _ _ _      |_ _ _S_ _|       |
.            _ _ _ _ _ _ _      |_ _ _ _  |     |     | |_      |_ _ R  |
.           |_ _S_ _|     |     | |_    | |_ R  |     |   |_    |_ S|   |
.           | |_    |_ R  |     |   |_  |_S |_ _|     |     |_  T |_|_ _|
.           |   |_  T |_ _|     |     |_T |_ _  |     |_ _    |_      | |
.           |_ _  |_    | |     |_ _  U |_    | |     |   |  U  |_    | |
.           |   |_U |_  |S|     |   |_    |_  | |     |   |_ _    |_  |S|
.           |  V  |   |_| |     |  V  |     |_| |     |  V    |     |_| |
.           |_ _ _|_ _ _|_|     |_ _ _|_ _ _ _|_|     |_ _ _ _|_ _ _ _|_|
.
n:                 7                    8                      9
R: A004125         8                    8                     12
S: A000203         8                   15                     12
T: A244048        12                   13                     20
U: a(n)           13                   20                     24
V: A004125         8                    8                     12
.
Illustration for n = 10..12:
.                                                         _ _ _ _ _ _ _ _ _ _ _ _
.                              _ _ _ _ _ _ _ _ _ _ _     |_ _ _ _ _ _  |         |
.     _ _ _ _ _ _ _ _ _ _     |_ _ _S_ _ _|         |    | |_        | |_ _   R  |
.    |_ _ _S_ _  |       |    | |_        |      R  |    |   |_      |     |_    |
.    | |_      | |_  R   |    |   |_      |_        |    |     |_    |_  S   |   |
.    |   |_    |_ _|_    |    |     |_      |_      |    |       |_    |_    |_ _|
.    |     |_      | |_ _|    |       |_   T  |_ _ _|    |         |_ T  |_ _ _  |
.    |       |_ T  |_ _  |    |_ _ _    |_        | |    |_ _        |_        | |
.    |_ _      |_      | |    |     |_ U  |_      | |    |   |    U    |_      | |
.    |   |_ U    |_    |S|    |       |_    |_    |S|    |   |_          |_    | |
.    |     |_      |_  | |    |         |     |_  | |    |     |_ _        |_  | |
.    |  V    |       |_| |    |  V      |       |_| |    |  V      |         |_| |
.    |_ _ _ _|_ _ _ _ _|_|    |_ _ _ _ _|_ _ _ _ _|_|    |_ _ _ _ _|_ _ _ _ _ _|_|
.
n:            10                         11                          12
R: A004125    13                         22                          17
S: A000203    18                         12                          28
T: A244048    24                         32                          33
U: a(n)       32                         33                          49
V: A004125    13                         22                          17
.
Note that in the diagrams the symmetric representation of A244048(n+1) is the same as the symmetric representation of a(n) rotated 180 degrees.
The diagrams for n = 11 and n = 12 both are copies from the diagrams that are in A244048 dated Jun 24 2014.
[Another way for the illustration of this sequence which is visible in the pyramid described in A245092 will be added soon.]
(End)
		

Crossrefs

Programs

  • Mathematica
    f[n_] := Sum[ DivisorSigma[1, m] - m, {m, n}]; Array[f, 60] (* Robert G. Wilson v, Jun 30 2014 *)
    Accumulate@ Table[DivisorSum[n, # &, # < n &], {n, 51}] (* or *)
    Table[Sum[k Floor[(n - k)/k], {k, n}], {n, 51}] (* Michael De Vlieger, Apr 02 2017 *)
  • PARI
    a(n) = sum(k=1, n, sigma(k)-k); \\ Michel Marcus, Jan 22 2017
    
  • Python
    from math import isqrt
    def A153485(n): return (-n*(n+1)-(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

Formula

a(n) = A024916(n) - A000217(n).
a(n) = A000217(n-1) - A004125(n). - Omar E. Pol, Jan 28 2014
a(n) = A000290(n) - A000203(n) - A024816(n) - A004125(n) = A024816(n+1) - A004125(n+1). - Omar E. Pol, Jun 23 2014
G.f.: (1/(1 - x))*Sum_{k>=1} k*x^(2*k)/(1 - x^k). - Ilya Gutkovskiy, Jan 22 2017
a(n) = Sum_{k=1..n} k * floor((n-k)/k). - Wesley Ivan Hurt, Apr 02 2017
a(n) ~ n^2 * (Pi^2/12 - 1/2). - Vaclav Kotesovec, Dec 21 2020
a(n) = A000290(n) - A000217(n) - A004125(n). - Omar E. Pol, Feb 26 2021
a(n) = A244048(n+1). - Omar E. Pol, Mar 28 2021

Extensions

Better name from Omar E. Pol, Jan 28 2014, Jun 23 2014

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

Views

Author

Keywords

Comments

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

References

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

Crossrefs

Column k=2 of A319083 (shifted).

Programs

  • Haskell
    a000385 n = sum $ zipWith (*) sigmas $ reverse sigmas where
       sigmas = take n a000203_list
    -- Reinhard Zumkeller, Sep 20 2011
    
  • Maple
    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
  • Mathematica
    a[n_] := Sum[DivisorSigma[1, k] DivisorSigma[1, n-k+1], {k, 1, n}];
    Array[a, 100] (* Jean-François Alcover, Aug 01 2018 *)
  • PARI
    a(n) = sum(k=1, n, sigma(k)*sigma(n-k+1)); \\ Michel Marcus, Nov 10 2016
    
  • PARI
    a(n) = my(f = factor(n+1)); (5 * sigma(f, 3) - (6*n + 5) * sigma(f)) / 12; \\ Amiram Eldar, Jan 04 2025
    
  • Python
    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

Formula

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

Extensions

More terms from Sean A. Irvine, Nov 14 2010

A336812 Irregular triangle read by rows T(n,k), n >= 1, k >= 1, in which row n is constructed replacing every term of row n of A336811 with its divisors.

Original entry on oeis.org

1, 1, 2, 1, 3, 1, 1, 2, 4, 1, 2, 1, 1, 5, 1, 3, 1, 2, 1, 1, 1, 2, 3, 6, 1, 2, 4, 1, 3, 1, 2, 1, 2, 1, 1, 1, 7, 1, 5, 1, 2, 4, 1, 3, 1, 3, 1, 2, 1, 2, 1, 1, 1, 1, 1, 2, 4, 8, 1, 2, 3, 6, 1, 5, 1, 2, 4, 1, 2, 4, 1, 3, 1, 3, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 3, 9, 1, 7, 1, 2, 3, 6
Offset: 1

Views

Author

Omar E. Pol, Nov 20 2020

Keywords

Comments

Here we introduce a new type of table which shows the correspondence between divisors and partitions. More precisely the table shows the corresponce between all parts of the last section of the set of partitions of n and all divisors of all terms of the n-th row of A336811, with n >= 1. The mentionded parts and the mentioned divisors are the same numbers (see Example section).
For an equivalent table showing the same kind of correspondence for all partitions of all positive integers see the supersequence A338156.

Examples

			Triangle begins:
  [1];
  [1, 2];
  [1, 3],       [1];
  [1, 2, 4],    [1, 2],    [1];
  [1, 5],       [1, 3],    [1, 2], [1],    [1];
  [1, 2, 3, 6], [1, 2, 4], [1, 3], [1, 2], [1, 2], [1], [1];
  ...
For n = 6 the 6th row of A336811 is [6, 4, 3, 2, 2, 1, 1] so replacing every term with its divisors we have {[1, 2, 3, 6], [1, 2, 4], [1, 3], [1, 2], [1, 2], [1], [1]} the same as the 6th row of this triangle.
Also, if the sequence is written as an irregular tetrahedron so the first six slices are:
  -------------
  [1],
  -------------
  [1, 2];
  -------------
  [1, 3],
  [1];
  -------------
  [1, 2, 4],
  [1, 2],
  [1];
  -------------
  [1, 5],
  [1, 3],
  [1, 2],
  [1],
  [1];
  -------------
  [1, 2, 3, 6],
  [1, 2, 4],
  [1, 3],
  [1, 2],
  [1, 2],
  [1],
  [1];
  -------------
The above slices appear in the lower zone of the following table which shows the correspondence between the mentioned divisors and the parts of the last section of the set of partitions of the positive integers.
The table is infinite. It is formed by three zones as follows:
The upper zone shows the last section of the set of partitions of every positive integer.
The lower zone shows the same numbers but arranged as divisors in accordance with the slices of the tetrahedron mentioned above.
Finally the middle zone shows the connection between the upper zone and the lower zone.
For every positive integer the numbers in the upper zone are the same numbers as in the lower zone.
|---|---------|-----|-------|---------|-----------|-------------|---------------|
| n |         |  1  |   2   |    3    |     4     |      5      |       6       |
|---|---------|-----|-------|---------|-----------|-------------|---------------|
|   |         |     |       |         |           |             |  6            |
| P |         |     |       |         |           |             |  3 3          |
| A |         |     |       |         |           |             |  4 2          |
| R |         |     |       |         |           |             |  2 2 2        |
| T |         |     |       |         |           |  5          |    1          |
| I |         |     |       |         |           |  3 2        |      1        |
| T |         |     |       |         |  4        |    1        |      1        |
| I |         |     |       |         |  2 2      |      1      |        1      |
| O |         |     |       |  3      |    1      |      1      |        1      |
| N |         |     |  2    |    1    |      1    |        1    |          1    |
| S |         |  1  |    1  |      1  |        1  |          1  |            1  |
|---|---------|-----|-------|---------|-----------|-------------|---------------|
.
|---|---------|-----|-------|---------|-----------|-------------|---------------|
|   | A207031 |  1  |  2 1  |  3 1 1  |  6 3 1 1  |  8 3 2 1 1  | 15 8 4 2 1 1  |
| L |         |  |  |  |/|  |  |/|/|  |  |/|/|/|  |  |/|/|/|/|  |  |/|/|/|/|/|  |
| I | A182703 |  1  |  1 1  |  2 0 1  |  3 2 0 1  |  5 1 1 0 1  |  7 4 2 1 0 1  |
| N |         |  *  |  * *  |  * * *  |  * * * *  |  * * * * *  |  * * * * * *  |
| K | A002260 |  1  |  1 2  |  1 2 3  |  1 2 3 4  |  1 2 3 4 5  |  1 2 3 4 5 6  |
|   |         |  =  |  = =  |  = = =  |  = = = =  |  = = = = =  |  = = = = = =  |
|   | A207383 |  1  |  1 2  |  2 0 3  |  3 4 0 4  |  5 2 3 0 5  |  7 8 6 4 0 6  |
|---|---------|-----|-------|---------|-----------|-------------|---------------|
.
|---|---------|-----|-------|---------|-----------|-------------|---------------|
|   | A027750 |  1  |  1 2  |  1   3  |  1 2   4  |  1       5  |  1 2 3     6  |
| D |---------|-----|-------|---------|-----------|-------------|---------------|
| I | A027750 |     |       |  1      |  1 2      |  1   3      |  1 2   4      |
| V |---------|-----|-------|---------|-----------|-------------|---------------|
| I | A027750 |     |       |         |  1        |  1 2        |  1   3        |
| S |---------|-----|-------|---------|-----------|-------------|---------------|
| O | A027750 |     |       |         |           |  1          |  1 2          |
| R | A027750 |     |       |         |           |  1          |  1 2          |
| S |---------|-----|-------|---------|-----------|-------------|---------------|
|   | A027750 |     |       |         |           |             |  1            |
|   | A027750 |     |       |         |           |             |  1            |
|---|---------|-----|-------|---------|-----------|-------------|---------------|
.
Note that every row in the lower zone lists A027750.
The "section" is the simpler substructure of the set of partitions of n that has this property in the three zones.
Also the lower zone for every positive integer can be constructed using the first n terms of A002865. For example: for n = 6 we consider the first 6 terms of A002865 (that is [1, 0, 1, 1, 2, 2]) and then the 6th slice is formed by a block with the divisors of 6, no block with the divisors of 5, one block with the divisors of 4, one block with the divisors of 3, two blocks with the divisors of 2 and two blocks with the divisors of 1.
Note that the lower zone is also in accordance with the tower (a polycube) described in A221529 in which its terraces are the symmetric representation of sigma starting from the top (cf. A237593) and the heights of the mentioned terraces are the partition numbers A000041 starting from the base.
The tower has the same volume (also the same number of cubes) equal to A066186(n) as a prism of partitions of size 1*n*A000041(n).
The above table shows the growth step by step of both the prism of partitions and its associated tower since the number of parts in the last section of the set of partitions of n is equal to A138137(n) equaling the number of divisors in the n-th slice of the lower table and equaling the same the number of terms in the n-th row of triangle. Also the sum of all parts in the last section of the set of partitions of n is equal to A138879(n) equaling the sum of all divisors in the n-th slice of the lower table and equaling the sum of the n-th row of triangle.
		

Crossrefs

Programs

  • Mathematica
    A336812[row_]:=Flatten[Table[ConstantArray[Divisors[row-m],PartitionsP[m]-PartitionsP[m-1]],{m,0,row-1}]];
    Array[A336812,10] (* Generates 10 rows *) (* Paolo Xausa, Feb 16 2023 *)

A098198 Decimal expansion of Pi^4/36 = zeta(2)^2.

Original entry on oeis.org

2, 7, 0, 5, 8, 0, 8, 0, 8, 4, 2, 7, 7, 8, 4, 5, 4, 7, 8, 7, 9, 0, 0, 0, 9, 2, 4, 1, 3, 5, 2, 9, 1, 9, 7, 5, 6, 9, 3, 6, 8, 7, 7, 3, 7, 9, 7, 9, 6, 8, 1, 7, 2, 6, 9, 2, 0, 7, 4, 4, 0, 5, 3, 8, 6, 1, 0, 3, 0, 1, 5, 4, 0, 4, 6, 7, 4, 2, 2, 1, 1, 6, 3, 9, 2, 2, 7, 4, 0, 8, 9, 8, 5, 4, 2, 4, 9, 7, 9, 3, 0, 8, 2, 4, 7
Offset: 1

Views

Author

Labos Elemer, Sep 21 2004

Keywords

Examples

			2.70580808427784547879000924135291975693687737979... = 2*A152649 = A013661^2.
		

Crossrefs

Programs

Formula

Decimal expansion of limit of q(n)= A024916(n)/A002088(n) = SummatorySigma / SummatoryTotient.
Equals Sum_{n>=1} A000005(n)/n^2. - R. J. Mathar, Dec 18 2010
Equals 10*Sum_{n>=2} (psi(n)+gamma)/n^3. - Jean-François Alcover, Feb 25 2013
Equals Zeta(4)*10/4 = A013662/0.4 = 1/A227929. - R. J. Mathar, Jul 20 2025
Equals 10 * zeta(3,1) = 10 * Sum_{n >= 1} 1/n Sum_{k >= n+1} 1/k^3 = 10 * Sum_{n >= 1} 1/n^3 * Sum_{k = 1..n-1} 1/k. - Peter Bala, Aug 07 2025

A237590 a(n) is the total number of regions (or parts) after n-th stage in the diagram of the symmetries of sigma described in A236104.

Original entry on oeis.org

1, 2, 4, 5, 7, 8, 10, 11, 14, 16, 18, 19, 21, 23, 26, 27, 29, 30, 32, 33, 37, 39, 41, 42, 45, 47, 51, 52, 54, 55, 57, 58, 62, 64, 67, 68, 70, 72, 76, 77, 79, 80, 82, 84, 87, 89, 91, 92, 95, 98, 102, 104, 106, 107, 111, 112, 116, 118, 120, 121, 123, 125, 130, 131, 135, 136, 138, 140, 144, 147, 149, 150, 152, 154
Offset: 1

Views

Author

Omar E. Pol, Mar 31 2014

Keywords

Comments

The total area (or total number of cells) of the diagram after n stages is equal to A024916(n), the sum of all divisors of all positive integers <= n.
Note that the region between the virtual circumscribed square and the diagram is a symmetric polygon whose area is equal to A004125(n), see example.
For more information see A237593 and A237270.
a(n) is also the total number of terraces of the stepped pyramid with n levels described in A245092. - Omar E. Pol, Apr 20 2016

Examples

			Illustration of initial terms:
.                                                         _ _ _ _
.                                           _ _ _        |_ _ _  |_
.                               _ _ _      |_ _ _|       |_ _ _|   |_
.                     _ _      |_ _  |_    |_ _  |_ _    |_ _  |_ _  |
.             _ _    |_ _|_    |_ _|_  |   |_ _|_  | |   |_ _|_  | | |
.       _    |_  |   |_  | |   |_  | | |   |_  | | | |   |_  | | | | |
.      |_|   |_|_|   |_|_|_|   |_|_|_|_|   |_|_|_|_|_|   |_|_|_|_|_|_|
.
.
.       1      2        4          5            7              8
.
For n = 6 the diagram contains 8 regions (or parts), so a(6) = 8.
The sum of all divisors of all positive integers <= 6 is [1] + [1+2] + [1+3] + [1+2+4] + [1+5] + [1+2+3+6] = 33. On the other hand after 6 stages the sum of all parts of the diagram is [1] + [3] + [2+2] + [7] + [3+3] + [12] = 33, equaling the sum of all divisors of all positive integers <= 6.
Note that the region between the virtual circumscribed square and the diagram is a symmetric polygon whose area is equal to A004125(6) = 3.
From _Omar E. Pol_, Dec 25 2020: (Start)
Illustration of the diagram after 29 stages (contain 215 vertices, 268 edges and 54 regions or parts):
._ _ _ _ _ _ _ _ _ _ _ _ _ _ _
|_ _ _ _ _ _ _ _ _ _ _ _ _ _ _|
|_ _ _ _ _ _ _ _ _ _ _ _ _ _  |
|_ _ _ _ _ _ _ _ _ _ _ _ _ _| |
|_ _ _ _ _ _ _ _ _ _ _ _ _  | |
|_ _ _ _ _ _ _ _ _ _ _ _ _| | |
|_ _ _ _ _ _ _ _ _ _ _ _  | | |_ _ _
|_ _ _ _ _ _ _ _ _ _ _ _| | |_ _ _  |
|_ _ _ _ _ _ _ _ _ _ _  | | |_ _  | |_
|_ _ _ _ _ _ _ _ _ _ _| | |_ _ _| |_  |_
|_ _ _ _ _ _ _ _ _ _  | |       |_ _|   |_
|_ _ _ _ _ _ _ _ _ _| | |_ _    |_  |_ _  |_ _
|_ _ _ _ _ _ _ _ _  | |_ _ _|     |_  | |_ _  |
|_ _ _ _ _ _ _ _ _| | |_ _  |_      |_|_ _  | |
|_ _ _ _ _ _ _ _  | |_ _  |_ _|_        | | | |_ _ _ _ _ _
|_ _ _ _ _ _ _ _| |     |     | |_ _    | |_|_ _ _ _ _  | |
|_ _ _ _ _ _ _  | |_ _  |_    |_  | |   |_ _ _ _ _  | | | |
|_ _ _ _ _ _ _| |_ _  |_  |_ _  | | |_ _ _ _ _  | | | | | |
|_ _ _ _ _ _  | |_  |_  |_    | |_|_ _ _ _  | | | | | | | |
|_ _ _ _ _ _| |_ _|   |_  |   |_ _ _ _  | | | | | | | | | |
|_ _ _ _ _  |     |_ _  | |_ _ _ _  | | | | | | | | | | | |
|_ _ _ _ _| |_      | |_|_ _ _  | | | | | | | | | | | | | |
|_ _ _ _  |_ _|_    |_ _ _  | | | | | | | | | | | | | | | |
|_ _ _ _| |_  | |_ _ _  | | | | | | | | | | | | | | | | | |
|_ _ _  |_  |_|_ _  | | | | | | | | | | | | | | | | | | | |
|_ _ _|   |_ _  | | | | | | | | | | | | | | | | | | | | | |
|_ _  |_ _  | | | | | | | | | | | | | | | | | | | | | | | |
|_ _|_  | | | | | | | | | | | | | | | | | | | | | | | | | |
|_  | | | | | | | | | | | | | | | | | | | | | | | | | | | |
|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|_|
.
(End)
		

Crossrefs

Partial sums of A237271.
Compare with A060831 (analog for the diagram that contains subparts).

Programs

  • Mathematica
    (* total number of parts in the first n symmetric representations *)
    (* Function a237270[] is defined in A237270 *)
    (* variable "previous" represents the sum from 1 through m-1 *)
    a237590[previous_,{m_,n_}]:=Rest[FoldList[Plus[#1,Length[a237270[#2]]]&,previous,Range[m,n]]]
    a237590[n_]:=a237590[0,{1,n}]
    a237590[78] (* data *)
    (* Hartmut F. W. Hoft, Jul 07 2014 *)

Formula

a(n) = A317109(n) - A294723(n) + 1 (Euler's formula). - Omar E. Pol, Jul 21 2018

Extensions

Definition clarified by Omar E. Pol, Jul 21 2018

A339278 Irregular triangle read by rows T(n,k), (n >= 1, k >= 1), in which the partition number A000041(n-1) is the length of row n and every column k is A000203, the sum of divisors function.

Original entry on oeis.org

1, 3, 4, 1, 7, 3, 1, 6, 4, 3, 1, 1, 12, 7, 4, 3, 3, 1, 1, 8, 6, 7, 4, 4, 3, 3, 1, 1, 1, 1, 15, 12, 6, 7, 7, 4, 4, 3, 3, 3, 3, 1, 1, 1, 1, 13, 8, 12, 6, 6, 7, 7, 4, 4, 4, 4, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 18, 15, 8, 12, 12, 6, 6, 7, 7, 7, 7, 4, 4, 4, 4, 3, 3, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Omar E. Pol, Nov 29 2020

Keywords

Comments

The sum of row n equals A138879(n), the sum of all parts in the last section of the set of partitions of n.
T(n,k) is also the number of cubic cells (or cubes) added at the n-th stage in the k-th level starting from the base in the tower described in A221529, assuming that the tower is an object under construction (see the example). - Omar E. Pol, Jan 20 2022

Examples

			Triangle begins:
   1;
   3;
   4,  1;
   7,  3,  1;
   6,  4,  3, 1, 1;
  12,  7,  4, 3, 3, 1, 1;
   8,  6,  7, 4, 4, 3, 3, 1, 1, 1, 1;
  15, 12,  6, 7, 7, 4, 4, 3, 3, 3, 3, 1, 1, 1, 1;
  13,  8, 12, 6, 6, 7, 7, 4, 4, 4, 4, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1;
...
From _Omar E. Pol_, Jan 13 2022: (Start)
Illustration of the first six rows of triangle showing the growth of the symmetric tower described in A221529:
    Level k: 1              2         3        4       5      6     7
Stage
  n   _ _ _ _ _ _ _ _
     |            _  |
  1  |           |_| |
     |_ _ _ _ _ _ _ _|
     |          _    |
     |         | |_  |
  2  |         |_ _| |
     |_ _ _ _ _ _ _ _|_ _ _ _ _ _
     |        _      |        _  |
     |       | |     |       |_| |
  3  |       |_|_ _  |           |
     |         |_ _| |           |
     |_ _ _ _ _ _ _ _|_ _ _ _ _ _|_ _ _ _ _
     |      _        |      _    |      _  |
     |     | |       |     | |_  |     |_| |
  4  |     | |_      |     |_ _| |         |
     |     |_  |_ _  |           |         |
     |       |_ _ _| |           |         |
     |_ _ _ _ _ _ _ _|_ _ _ _ _ _|_ _ _ _ _|_ _ _ _ _ _ _ _
     |    _          |    _      |    _    |    _  |    _  |
     |   | |         |   | |     |   | |_  |   |_| |   |_| |
     |   | |         |   |_|_ _  |   |_ _| |       |       |
  5  |   |_|_        |     |_ _| |         |       |       |
     |       |_ _ _  |           |         |       |       |
     |       |_ _ _| |           |         |       |       |
     |_ _ _ _ _ _ _ _|_ _ _ _ _ _|_ _ _ _ _|_ _ _ _|_ _ _ _|_ _ _ _ _ _
     |  _            |  _        |  _      |  _    |  _    |  _  |  _  |
     | | |           | | |       | | |     | | |_  | | |_  | |_| | |_| |
     | | |           | | |_      | |_|_ _  | |_ _| | |_ _| |     |     |
     | | |_ _        | |_  |_ _  |   |_ _| |       |       |     |     |
  6  | |_    |       |   |_ _ _| |         |       |       |     |     |
     |   |_  |_ _ _  |           |         |       |       |     |     |
     |     |_ _ _ _| |           |         |       |       |     |     |
     |_ _ _ _ _ _ _ _|_ _ _ _ _ _|_ _ _ _ _|_ _ _ _|_ _ _ _|_ _ _|_ _ _|
.
Every cell in the diagram of the symmetric representation of sigma represents a cubic cell or cube.
For n = 6 and k = 3 we add four cubes at 6th stage in the third level of the structure of the tower starting from the base so T(6,3) = 4.
For n = 9 another connection with the tower is as follows:
First we take the columns from the above triangle and build a new triangle in which all columns start at row 1 as shown below:
.
   1,  1,  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
   3,  3,  3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3;
   4,  4,  4, 4, 4, 4, 4, 4, 4, 4, 4;
   7,  7,  7, 7, 7, 7, 7;
   6,  6,  6, 6, 6;
  12, 12, 12;
   8,  8;
  15;
  13;
.
Then we rotate the triangle by 90 degrees as shown below:
                                       _
  1;                                  | |
  1;                                  | |
  1;                                  | |
  1;                                  | |
  1;                                  | |
  1;                                  | |
  1;                                  |_|_
  1, 3;                               |   |
  1, 3;                               |   |
  1, 3;                               |   |
  1, 3;                               |_ _|_
  1, 3, 4;                            |   | |
  1, 3, 4;                            |   | |
  1, 3, 4;                            |   | |
  1, 3, 4;                            |_ _|_|_
  1, 3, 4, 7;                         |     | |
  1, 3, 4, 7;                         |_ _ _| |_
  1, 3, 4, 7, 6;                      |     |   |
  1, 3, 4, 7, 6;                      |_ _ _|_ _|_
  1, 3, 4, 7, 6, 12;                  |_ _ _ _| | |_
  1, 3, 4, 7, 6, 12, 8;               |_ _ _ _|_|_ _|_ _
  1, 3, 4, 7, 6, 12, 8, 15; 13;       |_ _ _ _ _|_ _|_ _|
.
                                         Lateral view
                                         of the tower
.                                      _ _ _ _ _ _ _ _ _
                                      |_| | | | | | |   |
                                      |_ _|_| | | | |   |
                                      |_ _|  _|_| | |   |
                                      |_ _ _|    _|_|   |
                                      |_ _ _|  _|    _ _|
                                      |_ _ _ _|     |
                                      |_ _ _ _|  _ _|
                                      |         |
                                      |_ _ _ _ _|
.
                                           Top view
                                         of the tower
.
The sum of the m-th row of the new triangle equals A024916(j) where j is the length of the m-th row, equaling the number of cubic cells in the m-th level of the tower. For example: the last row of triangle has 9 terms and the sum of the last row is 1 + 3 + 4 + 7 + 6 + 12 + 8 + 15 + 13 = A024916(9) = 69, equaling the number of cubes in the base of the tower. (End)
		

Crossrefs

Sum of divisors of A336811.
Row n has length A000041(n-1).
Every column gives A000203.
The length of the m-th block in row n is A187219(m), m >= 1.
Row sums give A138879.
Cf. A337209 (another version).
Cf. A272172 (analog for the stepped pyramid described in A245092).

Programs

  • Mathematica
    A339278[rowmax_]:=Table[Flatten[Table[ConstantArray[DivisorSigma[1,n-m],PartitionsP[m]-PartitionsP[m-1]],{m,0,n-1}]],{n,rowmax}];
    A339278[15] (* Generates 15 rows *) (* Paolo Xausa, Feb 17 2023 *)
  • PARI
    f(n) = numbpart(n-1);
    T(n, k) = {if (k > f(n), error("invalid k")); if (k==1, return (sigma(n))); my(s=0); while (k <= f(n-1), s++; n--;); sigma(1+s);}
    tabf(nn) = {for (n=1, nn, for (k=1, f(n), print1(T(n,k), ", ");); print;);} \\ Michel Marcus, Jan 13 2021
    
  • PARI
    A339278(rowmax)=vector(rowmax,n,concat(vector(n,m,vector(numbpart(m-1)-numbpart(m-2),i,sigma(n-m+1)))));
    A339278(15) \\ Generates 15 rows \\ Paolo Xausa, Feb 17 2023

Formula

a(m) = A000203(A336811(m)).
T(n,k) = A000203(A336811(n,k)).

A222548 a(n) = Sum_{k=1..n} floor(n/k)^2.

Original entry on oeis.org

1, 5, 11, 22, 32, 52, 66, 92, 115, 147, 169, 219, 245, 289, 333, 390, 424, 496, 534, 612, 672, 740, 786, 898, 957, 1037, 1113, 1219, 1277, 1413, 1475, 1595, 1687, 1791, 1883, 2056, 2130, 2246, 2354, 2526, 2608, 2792, 2878, 3040, 3190, 3330, 3424, 3662, 3773
Offset: 1

Views

Author

Benoit Cloitre, Feb 24 2013

Keywords

Comments

a(n) is the number of common divisors of integers 1<=i,j<=n over all ordered pairs (i,j). - Geoffrey Critzer, Jan 15 2015

References

  • J. V. Uspensky and M. A. Heaslet, Elementary Number Theory, McGraw-Hill, NY, 1939, p. 98.

Crossrefs

Similar sequences for Sum_{k=1..n} floor(n/k)^m: A006218 (m=1), this sequence (m=2), A318742 (m=3), A318743 (m=4), A318744 (m=5).

Programs

  • Magma
    [&+[Floor(n/k)^2:k in [1..n] ]: n in [1..40]]; // Marius A. Burtea, Jul 16 2019
    
  • Mathematica
    Table[Sum[Floor[n/k]^2, {k, n}], {n, 50}] (* T. D. Noe, Feb 26 2013 *)
    Table[nn = n;Total[Level[Table[Table[DivisorSigma[0, GCD[i, j]], {i, 1, nn}], {j, 1, nn}], {2}]], {n, 1, 49}] (* Geoffrey Critzer, Jan 15 2015 *)
    Table[Sum[2*DivisorSigma[1, k] - DivisorSigma[0, k], {k, 1, n}], {n, 1, 50}] (* Vaclav Kotesovec, Sep 02 2018 *)
  • PARI
    a(n)=sum(k=1,n,(n\k)^2)
    
  • Python
    from math import isqrt
    def A222548(n): return -(s:=isqrt(n))**3 + sum((q:=n//k)*((k<<1)+q-1) for k in range(1,s+1)) # Chai Wah Wu, Oct 21 2023

Formula

a(n) = zeta(2)*n^2 + O(n log n).
a(n) = 2*A024916(n) - A006218(n). - Vaclav Kotesovec, Sep 02 2018
G.f.: (1/(1 - x)) * Sum_{k>=1} (2*k - 1) * x^k/(1 - x^k). - Ilya Gutkovskiy, Jul 16 2019
a(n) = Sum_{d=1..n} (2*d-1)*floor(n/d). [Uspensky and Heaslet] - Michael Somos, Feb 16 2020
a(n) = Sum_{k=1..n} Sum_{d|k} floor(n/d). - Ridouane Oudra, Jul 16 2020
a(n) = Sum_{i=1..n} Sum_{j=1..n} tau(gcd(i,j)). - Ridouane Oudra, Nov 23 2021
Previous Showing 31-40 of 228 results. Next