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 11-20 of 3898 results. Next

A336811 Irregular triangle read by rows T(n,k) in which the length of row n equals the partition number A000041(n-1) and every column k gives the positive integers A000027, with n >= 1 and k >= 1.

Original entry on oeis.org

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

Views

Author

Omar E. Pol, Nov 20 2020

Keywords

Comments

In other words: row n lists A028310(n-1) blocks where the m-th block consists of A187219(m) copies of n - m + [m=1], with n >= 1 and m >= 1, where [] is the Iverson bracket. [Corrected by Paolo Xausa, Feb 10 2023]
All divisors of all terms in row n are also all parts in the last section of the set of partitions of n.
Thus all divisors of all terms of the first n rows of triangle are also all parts of all partitions of n. In other words: all divisors of the first A000070(n-1) terms of the sequence are also all parts of all partitions of n. - Omar E. Pol, Jun 19 2021
From Omar E. Pol, Jul 31 2021: (Start)
The number of k's in row n is equal to A002865(n-k), 1 <= k <= n.
The number of terms >= k in row n is equal to A000041(n-k), 1 <= k <= n.
The number of k's in the first n rows (or in the first A000070(n-1) terms of the sequence) is equal to A000041(n-k), 1 <= k <= n.
The number of terms >= k in the first n rows (or in the first A000070(n-1) terms of the sequence) is equal to A000070(n-k), 1 <= k <= n.
First n rows of triangle (or first A000070(n-1) terms of the sequence) in nonincreasing order give the n-th row of A176206. (End)

Examples

			Triangle begins:
1;
2;
3, 1;
4, 2, 1;
5, 3, 2, 1, 1;
6, 4, 3, 2, 2, 1, 1;
7, 5, 4, 3, 3, 2, 2, 1, 1, 1, 1;
8, 6, 5, 4, 4, 3, 3, 2, 2, 2, 2, 1, 1, 1, 1;
9, 7, 6, 5, 5, 4, 4, 3, 3, 3, 3, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1;
...
For n = 6, by definition the length of row 6 is A000041(6-1) = A000041(5) = 7, so the row 6 of triangle has seven terms. Since every column lists the positive integers A000027 so the row 6 is [6, 4, 3, 2, 2, 1, 1].
Then we have that the divisors of the numbers of the 6th row are:
.
6th row of the triangle ---------->   6 4 3 2 2 1 1
                                      3 2 1 1 1
                                      2 1
                                      1
.
There are seven 1's, four 2's, two 3's, one 4 and one 6.
In total there are 7 + 4 + 2 + 1 + 1 = 15 divisors.
On the other hand the last section of the set of the partitions of 6 can be represented in several ways, five of them as shown below:
._ _ _ _ _ _
|_ _ _      |       6    6                  6                       6
|_ _ _|_    |     3 3    3 3              3   3                     3   3
|_ _    |   |     4 2    4 2            4       2                     4     2
|_ _|_ _|_  |   2 2 2    2 2 2        2   2       2                 2 2   2
          | |       1      1                        1                           1
          | |       1        1                        1                       1
          | |       1        1                          1                   1
          | |       1          1                          1               1
          | |       1          1                            1           1
          | |       1            1                            1       1
          |_|       1              1                            1   1
.
   Figure 1.  Figure 2.  Figure 3.        Figure 4.                   Figure 5.
.
In every figure there are seven 1's, four 2's, two 3's, one 4 and one 6, as shown also the 6th row of A182703.
In total there are 7 + 4 + 2 + 1 + 1 = A138137(6) = 15 parts in every figure.
Figure 5 is an arrangement that shows the correspondence between divisors and parts since the columns give the divisors of the terms of 6th row of triangle.
Finally we can see that all divisors of all numbers in the 6th row of the triangle are the same positive integers as all parts in the last section of the set of the partitions of 6.
Example edited by _Omar E. Pol_, Aug 10 2021
		

Crossrefs

Row sums give A000070.
Row n has length A000041(n-1).
Every column k gives A000027.
Companion of A176206.

Programs

  • Mathematica
    A336811[row_]:=Flatten[Table[ConstantArray[row-m,PartitionsP[m]-PartitionsP[m-1]],{m,0,row-1}]];
    Array[A336811,10] (* Generates 10 rows *) (* Paolo Xausa, Feb 10 2023 *)
  • PARI
    f(n) = numbpart(n-1);
    T(n, k) = {if (k > f(n), error("invalid k")); if (k==1, return (n)); my(s=0); while (k <= f(n-1), s++; n--;); 1+s;}
    tabf(nn) = {for (n=1, nn, for (k=1, f(n), print1(T(n,k), ", ");); print;);} \\ Michel Marcus, Jan 13 2021

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

A053445 Second differences of partition numbers A000041.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 3, 1, 4, 2, 7, 3, 10, 7, 14, 11, 22, 17, 32, 28, 45, 43, 67, 63, 95, 96, 134, 139, 192, 199, 269, 287, 373, 406, 521, 566, 718, 792, 983, 1092, 1346, 1496, 1827, 2045, 2465, 2772, 3323, 3733, 4449, 5016, 5929, 6696, 7882, 8897, 10426
Offset: 0

Views

Author

Alford Arnold, Jan 12 2000

Keywords

Comments

First differences of 0 1 1 2 2 4 4 7 8 12 14 21 24 34 41 55 ... (A002865).
For n > 2, a(n-2) is the number of partitions of n with all parts > 1 and with the largest part occurring more than once. The list of partitions counted begins 22 (so a(2) = 1); 33, 222 (so a(4) = 2); 44, 332, 2222 (so a(6) = 3); 333; 55, 442, 3322, 22222; 443, 3332; 66, 552, 444, 4422, 3333, 33222, 222222; 553, 4432, 33322; ...
a(n) is the number of certain level-n quasi-primary states of a quotient space of certain Verma modules. See the Furlan et al. reference p. 67. - Wolfdieter Lang, Apr 25 2003

Examples

			a(8) = 7 - 4 = 3; the corresponding partitions are 44, 332 and 2222.
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 115, 3.

Crossrefs

Programs

  • Magma
    m:=58; S:=[ NumberOfPartitions(n): n in [0..m] ]; [ S[n+2]-2*S[n+1]+S[n]: n in [1..m-2] ]; // Klaus Brockhaus, Jun 09 2009
    
  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
           b(n, i-1)+`if`(i>n, 0, b(n-i, i))))
        end:
    a:= n-> b(n$2) -2*b(n+1$2) +b(n+2$2):
    seq(a(n), n=0..80);  # Alois P. Heinz, May 19 2014
    # alternative Maple program:
    P:= [seq(combinat:-numbpart(n),n=0..1002)]:
    seq(P[i]-2*P[i+1]+P[i+2],i=1..1001); # Robert Israel, Dec 15 2014
  • Mathematica
    Table[(PartitionsP[n+2]-PartitionsP[n+1])-(PartitionsP[n+1]-PartitionsP[n]),{n,0,42}] (* Vladimir Joseph Stephan Orlovsky, Apr 23 2008 *)
    Differences[Table[ PartitionsP[n], {n, 0, 56}], 2] (* Jean-François Alcover, Sep 07 2011 *)
    Differences[PartitionsP[Range[0,60]],2] (* Harvey P. Dale, Jan 29 2016 *)
  • PARI
    lista(nn) = {v = vector(nn, n, numbpart(n-1)); dv = vector(#v-1, n, v[n+1] - v[n]); vector(#dv-1, n, dv[n+1] - dv[n]);} \\ Michel Marcus, Dec 15 2014
    
  • Python
    from sympy import npartitions
    def A053445(n): return npartitions(n+2)-(npartitions(n+1)<<1)+npartitions(n) # Chai Wah Wu, Mar 30 2023

Formula

a(n) = A002865(n+2) - A002865(n+1).
a(n) ~ exp(sqrt(2*n/3)*Pi)*Pi^2/(24*sqrt(3)*n^2) * (1 + (23*Pi/(24*sqrt(6)) - 3*sqrt(6)/Pi)/sqrt(n) + (625*Pi^2/6912 + 45/(2*Pi^2) - 115/24)/n). - Vaclav Kotesovec, Nov 04 2016
G.f.: 1/x - 1/x^2 + ((1 - x)/x^2)*Product_{k>=2} 1/(1 - x^k). - Ilya Gutkovskiy, Feb 28 2017

Extensions

More terms from James Sellers, Feb 02 2000
Start of sequence changed, Apr 25 2003

A338156 Irregular triangle read by rows in which row n lists n blocks, where the m-th block consists of A000041(m-1) copies of the divisors of (n - m + 1), with 1 <= m <= n.

Original entry on oeis.org

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

Views

Author

Omar E. Pol, Oct 14 2020

Keywords

Comments

In other words: in row n replace every term of n-th row of A176206 with its divisors.
The terms in row n are also all parts of all partitions of n.
As in A336812 here we introduce a new type of table which shows the correspondence between divisors and partitions. More precisely the table shows the correspondence between all divisors of all terms of the n-th row of A176206 and all parts of all partitions of n, with n >= 1. Both the mentionded divisors and the mentioned parts are the same numbers (see Example section). That is because all divisors of the first A000070(n-1) terms of A336811 are also all parts of all partitions of n.
For an equivalent table for all parts of the last section of the set of partitions of n see the subsequence A336812. The section is the smallest substructure of the set of partitions in which appears the correspondence divisor/part.
From Omar E. Pol, Aug 01 2021: (Start)
The terms of row n appears in the triangle A346741 ordered in accordance with the successive sections of the set of partitions of n.
The terms of row n in nonincreasing order give the n-th row of A302246.
The terms of row n in nondecreasing order give the n-th row of A302247.
For the connection with the tower described in A221529 see also A340035. (End)

Examples

			Triangle begins:
  [1];
  [1,2],   [1];
  [1,3],   [1,2],   [1],   [1];
  [1,2,4], [1,3],   [1,2], [1,2], [1],   [1],   [1];
  [1,5],   [1,2,4], [1,3], [1,3], [1,2], [1,2], [1,2], [1], [1], [1], [1], [1];
  ...
For n = 5 the 5th row of A176206 is [5, 4, 3, 3, 2, 2, 2, 1, 1, 1, 1, 1] so replacing every term with its divisors we have the 5th row of this triangle.
Also, if the sequence is written as an irregular tetrahedron so the first six slices are:
  [1],
  -------
  [1, 2],
  [1],
  -------
  [1, 3],
  [1, 2],
  [1],
  [1];
  ----------
  [1, 2, 4],
  [1, 3],
  [1, 2],
  [1, 2],
  [1],
  [1],
  [1];
  ----------
  [1, 5],
  [1, 2, 4],
  [1, 3],
  [1, 3],
  [1, 2],
  [1, 2],
  [1, 2],
  [1],
  [1],
  [1],
  [1],
  [1];
.
The above slices appear in the lower zone of the following table which shows the correspondence between the mentioned divisors and all parts of all partitions of the positive integers.
The table is infinite. It is formed by three zones as follows:
The upper zone shows the partitions of every positive integer in colexicographic order (cf. A026792, A211992).
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       |
|---|---------|-----|-------|---------|------------|---------------|
| P |         |     |       |         |            |               |
| A |         |     |       |         |            |               |
| R |         |     |       |         |            |               |
| T |         |     |       |         |            |  5            |
| I |         |     |       |         |            |  3  2         |
| T |         |     |       |         |  4         |  4  1         |
| I |         |     |       |         |  2  2      |  2  2  1      |
| O |         |     |       |  3      |  3  1      |  3  1  1      |
| N |         |     |  2    |  2 1    |  2  1 1    |  2  1  1 1    |
| S |         |  1  |  1 1  |  1 1 1  |  1  1 1 1  |  1  1  1 1 1  |
----|---------|-----|-------|---------|------------|---------------|
.
|---|---------|-----|-------|---------|------------|---------------|
|   | A181187 |  1  |  3 1  |  6 2 1  | 12  5 2 1  | 20  8  4 2 1  |
|   |         |  |  |  |/|  |  |/|/|  |  |/ |/|/|  |  |/ | /|/|/|  |
| L | A066633 |  1  |  2 1  |  4 1 1  |  7  3 1 1  | 12  4  2 1 1  |
| I |         |  *  |  * *  |  * * *  |  *  * * *  |  *  *  * * *  |
| N | A002260 |  1  |  1 2  |  1 2 3  |  1  2 3 4  |  1  2  3 4 5  |
| K |         |  =  |  = =  |  = = =  |  =  = = =  |  =  =  = = =  |
|   | A138785 |  1  |  2 2  |  4 2 3  |  7  6 3 4  | 12  8  6 4 5  |
|   |         |  |  |  |\|  |  |\|\|  |  |\ |\|\|  |  |\ |\ |\|\|  |
|   | A206561 |  1  |  4 2  |  9 5 3  | 20 13 7 4  | 35 23 15 9 5  |
|---|---------|-----|-------|---------|------------|---------------|
.
|---|---------|-----|-------|---------|------------|---------------|
|   | A027750 |  1  |  1 2  |  1   3  |  1  2   4  |  1         5  |
|   |---------|-----|-------|---------|------------|---------------|
|   | A027750 |     |  1    |  1 2    |  1    3    |  1  2    4    |
|   |---------|-----|-------|---------|------------|---------------|
| D | A027750 |     |       |  1      |  1  2      |  1     3      |
| I | A027750 |     |       |  1      |  1  2      |  1     3      |
| V |---------|-----|-------|---------|------------|---------------|
| I | A027750 |     |       |         |  1         |  1  2         |
| S | A027750 |     |       |         |  1         |  1  2         |
| O | A027750 |     |       |         |  1         |  1  2         |
| R |---------|-----|-------|---------|------------|---------------|
| S | A027750 |     |       |         |            |  1            |
|   | A027750 |     |       |         |            |  1            |
|   | A027750 |     |       |         |            |  1            |
|   | A027750 |     |       |         |            |  1            |
|   | A027750 |     |       |         |            |  1            |
|---|---------|-----|-------|---------|------------|---------------|
.
Note that every row in the lower zone lists A027750.
Also the lower zone for every positive integer can be constructed using the first n terms of the partition numbers. For example: for n = 5 we consider the first 5 terms of A000041 (that is [1, 1, 2, 3, 5]) then the 5th slice is formed by a block with the divisors of 5, one block with the divisors of 4, two blocks with the divisors of 3, three blocks with the divisors of 2, and five 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 correspondence between the prism of partitions and its associated tower since the number of parts in all partitions of n is equal to A006128(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 of all partitions of n is equal to A066186(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

Nonzero terms of A340031.
Row n has length A006128(n).
The sum of row n is A066186(n).
The product of row n is A007870(n).
Row n lists the first n rows of A336812 (a subsequence).
The number of parts k in row n is A066633(n,k).
The sum of all parts k in row n is A138785(n,k).
The number of parts >= k in row n is A181187(n,k).
The sum of all parts >= k in row n is A206561(n,k).
The number of parts <= k in row n is A210947(n,k).
The sum of all parts <= k in row n is A210948(n,k).

Programs

  • Mathematica
    A338156[rowmax_]:=Table[Flatten[Table[ConstantArray[Divisors[n-m],PartitionsP[m]],{m,0,n-1}]],{n,rowmax}];
    A338156[10] (* Generates 10 rows *) (* Paolo Xausa, Jan 12 2023 *)
  • PARI
    A338156(rowmax)=vector(rowmax,n,concat(vector(n,m,concat(vector(numbpart(m-1),i,divisors(n-m+1))))));
    A338156(10) \\ Generates 10 rows - Paolo Xausa, Feb 17 2023

A134264 Coefficients T(j, k) of a partition transform for Lagrange compositional inversion of a function or generating series in terms of the coefficients of the power series for its reciprocal. Enumeration of noncrossing partitions and primitive parking functions. T(n,k) for n >= 1 and 1 <= k <= A000041(n-1), an irregular triangle read by rows.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 1, 4, 2, 6, 1, 1, 5, 5, 10, 10, 10, 1, 1, 6, 6, 3, 15, 30, 5, 20, 30, 15, 1, 1, 7, 7, 7, 21, 42, 21, 21, 35, 105, 35, 35, 70, 21, 1, 1, 8, 8, 8, 4, 28, 56, 56, 28, 28, 56, 168, 84, 168, 14, 70, 280, 140, 56, 140, 28, 1, 1, 9, 9, 9, 9, 36, 72
Offset: 1

Views

Author

Tom Copeland, Jan 14 2008

Keywords

Comments

Coefficients are listed in Abramowitz and Stegun order (A036036).
Given an invertible function f(t) analytic about t=0 (or a formal power series) with f(0)=0 and Df(0) not equal to 0, form h(t) = t / f(t) and let h_n denote the coefficient of t^n in h(t).
Lagrange inversion gives the compositional inverse about t=0 as g(t) = Sum_{j>=1} ( t^j * (1/j) * Sum_{permutations s with s(1) + s(2) + ... + s(j) = j - 1} h_s(1) * h_s(2) * ... * h_s(j) ) = t * T(1,1) * h_0 + Sum_{j>=2} ( t^j * Sum_{k=1..(# of partitions for j-1)} T(j,k) * H(j-1,k ; h_0,h_1,...) ), where H(j-1,k ; h_0,h_1,...) is the k-th partition for h_1 through h_(j-1) corresponding to n=j-1 on page 831 of Abramowitz and Stegun (ordered as in A&S) with (h_0)^(j-m)=(h_0)^(n+1-m) appended to each partition subsumed under n and m of A&S.
Denoting h_n by (n') for brevity, to 8th order in t,
g(t) = t * (0')
+ t^2 * [ (0') (1') ]
+ t^3 * [ (0')^2 (2') + (0') (1')^2 ]
+ t^4 * [ (0')^3 (3') + 3 (0')^2 (1') (2') + (0') (1')^3 ]
+ t^5 * [ (0')^4 (4') + 4 (0')^3 (1') (3') + 2 (0')^3 (2')^2 + 6 (0')^2 (1')^2 (2') + (0') (1')^4 ]
+ t^6 * [ (0')^5 (5') + 5 (0')^4 (1') (4') + 5 (0')^4 (2') (3') + 10 (0')^3 (1')^2 (3') + 10 (0')^3 (1') (2')^2 + 10 (0')^2 (1')^3 (2') + (0') (1')^5 ]
+ t^7 * [ (0')^6 (6') + 6 (0')^5 (1') (5') + 6 (0')^5 (2') (4') + 3 (0')^5 (3')^2 + 15 (0')^4 (1')^2 (4') + 30 (0')^4 (1') (2') (3') + 5 (0')^4 (2')^3 + 20 (0')^3 (1')^3 (3') + 30 (0')^3 (1')^2 (2')^2 + 15 (0')^2 (1')^4 (2') + (0') (1')^6]
+ t^8 * [ (0')^7 (7') + 7 (0')^6 (1') (6') + 7 (0')^6 (2') (5') + 7 (0')^6 (3') (4') + 21 (0')^5 (1')^2* (5') + 42 (0')^5 (1') (2') (4') + 21 (0')^5 (1') (3')^2 + 21 (0')^5 (2')^2 (3') + 35 (0')^4 (1')^3 (4') + 105 (0)^4 (1')^2 (2') (3') + 35 (0')^4 (1') (2')^3 + 35 (0')^3 (1')^4 (3') + 70 (0')^3 (1')^3 (2')^2 + 21 (0')^2 (1')^5 (2') + (0') (1')^7 ]
+ ..., where from the formula section, for example, T(8,1',2',...,7') = 7! / ((8 - (1'+ 2' + ... + 7'))! * 1'! * 2'! * ... * 7'!) are the coefficients of the integer partitions (1')^1' (2')^2' ... (7')^7' in the t^8 term.
A125181 is an extended, reordered version of the above sequence, omitting the leading 1, with alternate interpretations.
If the coefficients of partitions with the same exponent for h_0 are summed within rows, A001263 is obtained, omitting the leading 1.
From identification of the elements of the inversion with those on page 25 of the Ardila et al. link, the coefficients of the irregular table enumerate non-crossing partitions on [n]. - Tom Copeland, Oct 13 2014
From Tom Copeland, Oct 28-29 2014: (Start)
Operating with d/d(1') = d/d(h_1) on the n-th partition polynomial Prt(n;h_0,h_1,..,h_n) in square brackets above associated with t^(n+1) generates n * Prt(n-1;h_0,h_1,..,h_(n-1)); therefore, the polynomials are an Appell sequence of polynomials in the indeterminate h_1 when h_0=1 (a special type of Sheffer sequence).
Consequently, umbrally, [Prt(.;1,x,h_2,..) + y]^n = Prt(n;1,x+y,h_2,..); that is, Sum_{k=0..n} binomial(n,k) * Prt(k;1,x,h_2,..) * y^(n-k) = Prt(n;1,x+y,h_2,..).
Or, e^(x*z) * exp[Prt(.;1,0,h_2,..) * z] = exp[Prt(.;1,x,h_2,..) * z]. Then with x = h_1 = -(1/2) * d^2[f(t)]/dt^2 evaluated at t=0, the formal Laplace transform from z to 1/t of this expression generates g(t), the comp. inverse of f(t), when h_0 = 1 = df(t)/dt eval. at t=0.
I.e., t / (1 - t*(x + Prt(.;1,0,h_2,..))) = t / (1 - t*Prt(.;1,x,h_2,..)) = g(t), interpreted umbrally, when h_0 = 1.
(End)
Connections to and between arrays associated to the Catalan (A000108 and A007317), Riordan (A005043), Fibonacci (A000045), and Fine (A000957) numbers and to lattice paths, e.g., the Motzkin, Dyck, and Łukasiewicz, can be made explicit by considering the inverse in x of the o.g.f. of A104597(x,-t), i.e., f(x) = P(Cinv(x),t-1) = Cinv(x) / (1 + (t-1)*Cinv(x)) = x*(1-x) / (1 + (t-1)*x*(1-x)) = (x-x^2) / (1 + (t-1)*(x-x^2)), where Cinv(x) = x*(1-x) is the inverse of C(x) = (1 - sqrt(1-4*x)) / 2, a shifted o.g.f. for the Catalan numbers, and P(x,t) = x / (1+t*x) with inverse Pinv(x,t) = -P(-x,t) = x / (1-t*x). Then h(x,t) = x / f(x,t) = x * (1+(t-1)Cinv(x)) / Cinv(x) = 1 + t*x + x^2 + x^3 + ..., i.e., h_1=t and all other coefficients are 1, so the inverse of f(x,t) in x, which is explicitly in closed form finv(x,t) = C(Pinv(x,t-1)), is given by A091867, whose coefficients are sums of the refined Narayana numbers above obtained by setting h_1=(1')=t in the partition polynomials and all other coefficients to one. The group generators C(x) and P(x,t) and their inverses allow associations to be easily made between these classic number arrays. - Tom Copeland, Nov 03 2014
From Tom Copeland, Nov 10 2014: (Start)
Inverting in x with t a parameter, let F(x;t,n) = x - t*x^(n+1). Then h(x) = x / F(x;t,n) = 1 / (1-t*x^n) = 1 + t*x^n + t^2*x^(2n) + t^3*x^(3n) + ..., so h_k vanishes unless k = m*n with m an integer in which case h_k = t^m.
Finv(x;t,n) = Sum_{j>=0} {binomial((n+1)*j,j) / (n*j + 1)} * t^j * x^(n*j + 1), which gives the Catalan numbers for n=1, and the Fuss-Catalan sequences for n>1 (see A001764, n=2). [Added braces to disambiguate the formula. - N. J. A. Sloane, Oct 20 2015]
This relation reveals properties of the partitions and sums of the coefficients of the array. For n=1, h_k = t^k for all k, implying that the row sums are the Catalan numbers. For n = 2, h_k for k odd vanishes, implying that there are no blocks with only even-indexed h_k on the even-numbered rows and that only the blocks containing only even-sized bins contribute to the odd-row sums giving the Fuss-Catalan numbers for n=2. And so on, for n > 2.
These relations are reflected in any combinatorial structures enumerated by this array and the partitions, such as the noncrossing partitions depicted for a five-element set (a pentagon) in Wikipedia.
(End)
From Tom Copeland, Nov 12 2014: (Start)
An Appell sequence possesses an umbral inverse sequence (cf. A249548). The partition polynomials here, Prt(n;1,h_1,...), are an Appell sequence in the indeterminate h_1=u, so have an e.g.f. exp[Prt(.;1,u,h_2...)*t] = e^(u*t) * exp[Prt(.;1,0,h2,...)*t] with umbral inverses with an e.g.f e^(-u*t) / exp[Prt(.;1,0,h2,...)*t]. This makes contact with the formalism of A133314 (cf. also A049019 and A019538) and the signed, refined face partition polynomials of the permutahedra (or their duals), which determine the reciprocal of exp[Prt(.,0,u,h2...)*t] (cf. A249548) or exp[Prt(.;1,u,h2,...)*t], forming connections among the combinatorics of permutahedra and the noncrossing partitions, Dyck paths and trees (cf. A125181), and many other important structures isomorphic to the partitions of this entry, as well as to formal cumulants through A127671 and algebraic structures of Lie algebras. (Cf. relationship of permutahedra with the Eulerians A008292.)
(End)
From Tom Copeland, Nov 24 2014: (Start)
The n-th row multiplied by n gives the number of terms in the homogeneous symmetric monomials generated by [x(1) + x(2) + ... + x(n+1)]^n under the umbral mapping x(m)^j = h_j, for any m. E.g., [a + b + c]^2 = [a^2 + b^2 + c^2] + 2 * [a*b + a*c + b*c] is mapped to [3 * h_2] + 2 * [3 * h_1^2], and 3 * A134264(3) = 3 *(1,1)= (3,3) the number of summands in the two homogeneous polynomials in the square brackets. For n=3, [a + b + c + d]^3 = [a^3 + b^3 + ...] + 3 [a*b^2 + a*c^2 + ...] + 6 [a*b*c + a*c*d + ...] maps to [4 * h_3] + 3 [12 * h_1 * h_2] + 6 [4 * (h_1)^3], and the number of terms in the brackets is given by 4 * A134264(4) = 4 * (1,3,1) = (4,12,4).
The further reduced expression is 4 h_3 + 36 h_1 h_2 + 24 (h_1)^3 = A248120(4) with h_0 = 1. The general relation is n * A134264(n) = A248120(n) / A036038(n-1) where the arithmetic is performed on the coefficients of matching partitions in each row n.
Abramowitz and Stegun give combinatorial interpretations of A036038 and relations to other number arrays.
This can also be related to repeated umbral composition of Appell sequences and topology with the Bernoulli numbers playing a special role. See the Todd class link.
(End)
These partition polynomials are dubbed the Voiculescu polynomials on page 11 of the He and Jejjala link. - Tom Copeland, Jan 16 2015
See page 5 of the Josuat-Verges et al. reference for a refinement of these partition polynomials into a noncommutative version composed of nondecreasing parking functions. - Tom Copeland, Oct 05 2016
(Per Copeland's Oct 13 2014 comment.) The number of non-crossing set partitions whose block sizes are the parts of the n-th integer partition, where the ordering of integer partitions is first by total, then by length, then lexicographically by the reversed sequence of parts. - Gus Wiseman, Feb 15 2019
With h_0 = 1 and the other h_n replaced by suitably signed partition polynomials of A263633, the refined face partition polynomials for the associahedra of normalized A133437 with a shift in indices are obtained (cf. In the Realm of Shadows). - Tom Copeland, Sep 09 2019
Number of primitive parking functions associated to each partition of n. See Lemma 3.8 on p. 28 of Rattan. - Tom Copeland, Sep 10 2019
With h_n = n + 1, the d_k (A006013) of Table 2, p. 18, of Jong et al. are obtained, counting the n-point correlation functions in a quantum field theory. - Tom Copeland, Dec 25 2019
By inspection of the diagrams on Robert Dickau's website, one can see the relationship between the monomials of this entry and the connectivity of the line segments of the noncrossing partitions. - Tom Copeland, Dec 25 2019
Speicher has examples of the first four inversion partition polynomials on pp. 22 and 23 with his k_n equivalent to h_n = (n') here with h_0 = 1. Identifying z = t, C(z) = t/f(t) = h(t), and M(z) = f^(-1)(t)/t, then statement (3), on p. 43, of Theorem 3.26, C(z M(z)) = M(z), is equivalent to substituting f^(-1)(t) for t in t/f(t), and statement (4), M(z/C(z)) = C(z), to substituting f(t) for t in f^(-1)(t)/t. - Tom Copeland, Dec 08 2021
Given a Laurent series of the form f(z) = 1/z + h_1 + h_2 z + h_3 z^2 + ..., the compositional inverse is f^(-1)(z) = 1/z + Prt(1;1,h_1)/z^2 + Prt(2;1,h_1,h_2)/z^3 + ... = 1/z + h_1/z^2 + (h_1^2 + h_2)/z^3 + (h_1^3 + 3 h_1 h_2 + h_3)/z^4 + (h_1^4 + 6 h_1^2 h_2 + 4 h_1 h_3 + 2 h_2^2 + h_4)/z^5 + ... for which the polynomials in the numerators are the partition polynomials of this entry. For example, this formula applied to the q-expansion of Klein's j-invariant / function with coefficients A000521, related to monstrous moonshine, gives the compositional inverse with the coefficients A091406 (see He and Jejjala). - Tom Copeland, Dec 18 2021
The partition polynomials of A350499 'invert' the polynomials of this entry giving the indeterminates h_n. A multinomial formula for the coefficients of the partition polynomials of this entry, equivalent to the multinomial formula presented in the first four sentences of the formula section below, is presented in the MathOverflow question referenced in A350499. - Tom Copeland, Feb 19 2022

Examples

			1) With f(t) = t / (t-1), then h(t) = -(1-t), giving h_0 = -1, h_1 = 1 and h_n = 0 for n>1. Then g(t) = -t - t^2 - t^3 - ... = t / (t-1).
2) With f(t) = t*(1-t), then h(t) = 1 / (1-t), giving h_n = 1 for all n. The compositional inverse of this f(t) is g(t) = t*A(t) where A(t) is the o.g.f. for the Catalan numbers; therefore the sum over k of T(j,k), i.e., the row sum, is the Catalan number A000108(j-1).
3) With f(t) = (e^(-a*t)-1) / (-a), h(t) = Sum_{n>=0} Bernoulli(n) * (-a*t)^n / n! and g(t) = log(1-a*t) / (-a) = Sum_{n>=1} a^(n-1) * t^n / n. Therefore with h_n = Bernoulli(n) * (-a)^n / n!, Sum_{permutations s with s(1)+s(2)+...+s(j)=j-1} h_s(1) * h_s(2) * ... * h_s(j) = j * Sum_{k=1..(# of partitions for j-1)} T(j,k) * H(j-1,k ; h_0,h_1,...) = a^(j-1). Note, in turn, Sum_{a=1..m} a^(j-1) = (Bernoulli(j,m+1) - Bernoulli(j)) / j for the Bernoulli polynomials and numbers, for j>1.
4) With f(t,x) = t / (x-1+1/(1-t)), then h(t,x) = x-1+1/(1-t), giving (h_0)=x and (h_n)=1 for n>1. Then g(t,x) = (1-(1-x)*t-sqrt(1-2*(1+x)*t+((x-1)*t)^2)) / 2, a shifted o.g.f. in t for the Narayana polynomials in x of A001263.
5) With h(t)= o.g.f. of A075834, but with A075834(1)=2 rather than 1, which is the o.g.f. for the number of connected positroids on [n] (cf. Ardila et al., p. 25), g(t) is the o.g.f. for A000522, which is the o.g.f. for the number of positroids on [n]. (Added Oct 13 2014 by author.)
6) With f(t,x) = x / ((1-t*x)*(1-(1+t)*x)), an o.g.f. for A074909, the reverse face polynomials of the simplices, h(t,x) = (1-t*x) * (1-(1+t)*x) with h_0=1, h_1=-(1+2*t), and h_2=t*(1+t), giving as the inverse in x about 0 the o.g.f. (1+(1+2*t)*x-sqrt(1+(1+2*t)*2*x+x^2)) / (2*t*(1+t)*x) for signed A033282, the reverse face polynomials of the Stasheff polytopes, or associahedra. Cf. A248727. (Added Jan 21 2015 by author.)
7) With f(x,t) = x / ((1+x)*(1+t*x)), an o.g.f. for the polynomials (-1)^n * (1 + t + ... + t^n), h(t,x) = (1+x) * (1+t*x) with h_0=1, h_1=(1+t), and h_2=t, giving as the inverse in x about 0 the o.g.f. (1-(1+t)*x-sqrt(1-2*(1+t)*x+((t-1)*x)^2)) / (2*x*t) for the Narayana polynomials A001263. Cf. A046802. (Added Jan 24 2015 by author.)
From _Gus Wiseman_, Feb 15 2019: (Start)
Triangle begins:
   1
   1
   1   1
   1   3   1
   1   4   2   6   1
   1   5   5  10  10  10   1
   1   6   6   3  15  30   5  20  30  15   1
   1   7   7   7  21  42  21  21  35 105  35  35  70  21   1
Row 5 counts the following non-crossing set partitions:
  {{1234}}  {{1}{234}}  {{12}{34}}  {{1}{2}{34}}  {{1}{2}{3}{4}}
            {{123}{4}}  {{14}{23}}  {{1}{23}{4}}
            {{124}{3}}              {{12}{3}{4}}
            {{134}{2}}              {{1}{24}{3}}
                                    {{13}{2}{4}}
                                    {{14}{2}{3}}
(End)
		

References

  • A. Nica and R. Speicher (editors), Lectures on the Combinatorics of Free Probability, London Mathematical Society Lecture Note Series: 335, Cambridge University Press, 2006 (see in particular, Eqn. 9.14 on p. 141, enumerating noncrossing partitions).

Crossrefs

(A001263,A119900) = (reduced array, associated g(x)). See A145271 for meaning and other examples of reduced and associated.
Other orderings are A125181 and A306438.
Cf. A119900 (e.g.f. for reduced W(x) with (h_0)=t and (h_n)=1 for n>0).
Cf. A248927 and A248120, "scaled" versions of this Lagrange inversion.
Cf. A091867 and A125181, for relations to lattice paths and trees.
Cf. A249548 for use of Appell properties to generate the polynomials.
Cf. A133314, A049019, A019538, A127671, and A008292 for relations to permutahedra, Eulerians.
Cf. A006013.

Programs

  • Mathematica
    Table[Binomial[Total[y],Length[y]-1]*(Length[y]-1)!/Product[Count[y,i]!,{i,Max@@y}],{n,7},{y,Sort[Sort/@IntegerPartitions[n]]}] (* Gus Wiseman, Feb 15 2019 *)
  • PARI
    C(v)={my(n=vecsum(v), S=Set(v)); n!/((n-#v+1)!*prod(i=1, #S, my(x=S[i]); (#select(y->y==x, v))!))}
    row(n)=[C(Vec(p)) | p<-partitions(n-1)]
    { for(n=1, 7, print(row(n))) } \\ Andrew Howroyd, Feb 01 2022

Formula

For j>1, there are P(j,m;a...) = j! / [ (j-m)! (a_1)! (a_2)! ... (a_(j-1))! ] permutations of h_0 through h_(j-1) in which h_0 is repeated (j-m) times; h_1, repeated a_1 times; and so on with a_1 + a_2 + ... + a_(j-1) = m.
If, in addition, a_1 + 2 * a_2 + ... + (j-1) * a_(j-1) = j-1, then each distinct combination of these arrangements is correlated with a partition of j-1.
T(j,k) is [ P(j,m;a...) / j ] for the k-th partition of j-1 as described in the comments.
For example from g(t) above, T(5,4) = (5! / ((5-3)! * 2!)) / 5 = 6 for the 4th partition under n=5-1=4 with m=3 parts in A&S.
From Tom Copeland, Sep 30 2011: (Start)
Let W(x) = 1/(df(x)/dx)= 1/{d[x/h(x)]/dx}
= [(h_0)-1+:1/(1-h.*x):]^2 / {(h_0)-:[h.x/(1-h.x)]^2:}
= [(h_0)+(h_1)x+(h_2)x^2+...]^2 / [(h_0)-(h_2)x^2-2(h_3)x^3-3(h_4)x^4-...], where :" ": denotes umbral evaluation of the expression within the colons and h. is an umbral coefficient.
Then for the partition polynomials of A134264,
Poly[n;h_0,...,h_(n-1)]=(1/n!)(W(x)*d/dx)^n x, evaluated at x=0, and the compositional inverse of f(t) is g(t) = exp(t*W(x)*d/dx) x, evaluated at x=0. Also, dg(t)/dt = W(g(t)), and g(t) gives A001263 with (h_0)=u and (h_n)=1 for n>0 and A000108 with u=1.
(End)
From Tom Copeland, Oct 20 2011: (Start)
With exp(x* PS(.,t)) = exp(t*g(x)) = exp(x*W(y)d/dy) exp(t*y) eval. at y=0, the raising (creation) and lowering (annihilation) operators defined by R PS(n,t) = PS(n+1,t) and L PS(n,t) = n*PS(n-1,t) are
R = t*W(d/dt) = t*((h_0) + (h_1)d/dt + (h_2)(d/dt)^2 + ...)^2 / ((h_0) - (h_2)(d/dt)^2 - 2(h_3)(d/dt)^3 - 3(h_4)(d/dt)^4 + ...), and
L = (d/dt)/h(d/dt) = (d/dt) 1/((h_0) + (h_1)*d/dt + (h_2)*(d/dt)^2 + ...)
Then P(n,t) = (t^n/n!) dPS(n,z)/dz eval. at z=0 are the row polynomials of A134264. (Cf. A139605, A145271, and link therein to Mathemagical Forests for relation to planted trees on p. 13.)
(End)
Using the formalism of A263634, the raising operator for the partition polynomials of this array with h_0 = 1 begins as R = h_1 + h_2 D + h_3 D^2/2! + (h_4 - h_2^2) D^3/3! + (h_5 - 5 h_2 h_3) D^4/4! + (h_6 + 5 h_2^3 - 7 h_3^2 - 9 h_2 h_4) D^5/5! + (h_7 - 14 h_2 h_5 + 56 h_2^2 h_3) D^6/6! + ... with D = d/d(h_1). - Tom Copeland, Sep 09 2016
Let h(x) = x/f^{-1}(x) = 1/[1-(c_2*x+c_3*x^2+...)], with c_n all greater than zero. Then h_n are all greater than zero and h_0 = 1. Determine P_n(t) from exp[t*f^{-1}(x)] = exp[x*P.(t)] with f^{-1}(x) = x/h(x) expressed in terms of the h_n (cf. A133314 and A263633). Then P_n(b.) = 0 gives a recursion relation for the inversion polynomials of this entry a_n = b_n/n! in terms of the lower order inversion polynomials and P_j(b.)P_k(b.) = P_j(t)P_k(t)|{t^n = b_n} = d{j,k} >= 0 is the coefficient of x^j/j!*y^k/k! in the Taylor series expansion of the formal group law FGL(x,y) = f[f^{-1}(x)+f^{-1}(y)]. - Tom Copeland, Feb 09 2018
A raising operator for the partition polynomials with h_0 = 1 regarded as a Sheffer Appell sequence in h_1 is described in A249548. - Tom Copeland, Jul 03 2018

Extensions

Added explicit t^6, t^7, and t^8 polynomials and extended initial table to include the coefficients of t^8. - Tom Copeland, Sep 14 2016
Title modified by Tom Copeland, May 28 2018
More terms from Gus Wiseman, Feb 15 2019
Title modified by Tom Copeland, Sep 10 2019

A129129 An irregular triangular array of natural numbers read by rows, with shape sequence A000041(n) related to sequence A060850.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 7, 10, 9, 12, 16, 11, 14, 15, 20, 18, 24, 32, 13, 22, 21, 28, 25, 30, 40, 27, 36, 48, 64, 17, 26, 33, 44, 35, 42, 56, 50, 45, 60, 80, 54, 72, 96, 128, 19, 34, 39, 52, 55, 66, 88, 49, 70, 63, 84, 112, 75, 100, 90, 120, 160, 81, 108, 144, 192, 256
Offset: 0

Views

Author

Alford Arnold, Mar 31 2007

Keywords

Comments

The tree begins (at height n, n >= 0, nodes represent partitions of n)
0: 1
1: 2
2: 3 4
3: 5 6 8
4: 7 10 9 12 16
5: 11 14 15 20 18 24 32
...
and hence differs from A114622.
Ordering [graded reverse lexicographic order] of partitions (positive integer representation) of nonnegative integers, where part of size i [as summand] is mapped to i-th prime [as multiplicand], where the empty partition for 0 yields the empty product, i.e., 1. Permutation of positive integers, since bijection [1-1 and onto map] between the set of all partitions of nonnegative integers and positive integers. - Daniel Forgues, Aug 07 2018
These are all Heinz numbers of integer partitions in graded reverse-lexicographic order, where The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This is the so-called "Mathematica" order (sum/revlex) of partitions (A080577). Partitions in lexicographic order (sum/lex) are A193073, with Heinz numbers A334434. - Gus Wiseman, May 19 2020

Examples

			The array is a tree structure as described by A128628. If a node value has only one branch the value is twice that of its parent node. If it has two branches one is twice that of its parent node but the other is defined as indicated below:
(1) pick an odd number (e.g., 135)
(2) calculate its prime factorization (135 = 5*3*3*3)
(3) note the least prime factor (LPF(135) = 3)
(4) note the index of the LPF (index(3) = 2)
(5) subtract one from the index (2-1 = 1)
(6) calculate the prime associated with the value in step five (prime(1) = 2)
(7) The parent node of the odd number 135 is (2/3)*135 = 90 = A252461(135).
From _Daniel Forgues_, Aug 07 2018: (Start)
Partitions of 4 in graded reverse lexicographic order:
{4}: p_4 = 7;
{3,1}: p_3 * p_1 = 5 * 2 = 10;
{2,2}: p_2 * p_2 = 3^2 = 9;
{2,1,1}: p_2 * p_1 * p_1 = 3 * 2^2 = 12;
{1,1,1,1}: p_1 * p_1 * p_1 * p_1 = 2^4 = 16. (End)
From _Gus Wiseman_, May 19 2020: (Start)
The sequence together with the corresponding partitions begins:
    1: ()            24: (2,1,1,1)         35: (4,3)
    2: (1)           32: (1,1,1,1,1)       42: (4,2,1)
    3: (2)           13: (6)               56: (4,1,1,1)
    4: (1,1)         22: (5,1)             50: (3,3,1)
    5: (3)           21: (4,2)             45: (3,2,2)
    6: (2,1)         28: (4,1,1)           60: (3,2,1,1)
    8: (1,1,1)       25: (3,3)             80: (3,1,1,1,1)
    7: (4)           30: (3,2,1)           54: (2,2,2,1)
   10: (3,1)         40: (3,1,1,1)         72: (2,2,1,1,1)
    9: (2,2)         27: (2,2,2)           96: (2,1,1,1,1,1)
   12: (2,1,1)       36: (2,2,1,1)        128: (1,1,1,1,1,1,1)
   16: (1,1,1,1)     48: (2,1,1,1,1)       19: (8)
   11: (5)           64: (1,1,1,1,1,1)     34: (7,1)
   14: (4,1)         17: (7)               39: (6,2)
   15: (3,2)         26: (6,1)             52: (6,1,1)
   20: (3,1,1)       33: (5,2)             55: (5,3)
   18: (2,2,1)       44: (5,1,1)           66: (5,2,1)
(End)
		

Crossrefs

Cf. A080577 (the partitions), A252461, A114622, A128628, A215366 (sorted rows).
Row lengths are A000041.
Compositions under the same order are A066099.
The opposite version (sum/lex) is A334434.
The length-sensitive version (sum/length/revlex) is A334438.
The version for reversed (weakly increasing) partitions is A334436.
Lexicographically ordered reversed partitions are A026791.
Reversed partitions in Abramowitz-Stegun order (sum/length/lex) are A036036.
Sum of prime indices is A056239.
Sorting reversed partitions by Heinz number gives A112798.
Partitions in lexicographic order are A193073.
Sorting partitions by Heinz number gives A296150.

Programs

  • Maple
    b:= (n, i)-> `if`(n=0 or i=1, [2^n], [map(x-> x*ithprime(i),
                    b(n-i, min(n-i, i)))[], b(n, i-1)[]]):
    T:= n-> b(n$2)[]:
    seq(T(n), n=0..10);  # Alois P. Heinz, Feb 14 2020
  • Mathematica
    Array[Times @@ # & /@ Prime@ IntegerPartitions@ # &, 9, 0] // Flatten (* Michael De Vlieger, Aug 07 2018 *)
    b[n_, i_] := b[n, i] = If[n == 0 || i == 1, {2^n}, Join[(# Prime[i]&) /@ b[n - i, Min[n - i, i]], b[n, i - 1]]];
    T[n_] := b[n, n];
    T /@ Range[0, 10] // Flatten (* Jean-François Alcover, May 21 2021, after Alois P. Heinz *)

Formula

From Gus Wiseman, May 19 2020: (Start)
A001222(a(n)) = A238966(n).
A001221(a(n)) = A115623(n).
A056239(a(n)) = A036042(n).
A061395(a(n)) = A331581(n).
(End)

A194447 Rank of the n-th region of the set of partitions of j, if 1<=n<=A000041(j).

Original entry on oeis.org

0, 0, 0, 1, -1, 2, -2, 1, 2, 2, -5, 2, 3, 3, -8, 1, 2, 2, 2, 4, 3, -14, 2, 3, 3, 3, 2, 4, 4, -21, 1, 2, 2, 2, 4, 3, 1, 3, 5, 5, 4, -32, 2, 3, 3, 3, 2, 4, 4, 1, 4, 3, 5, 6, 5, -45, 1, 2, 2, 2, 4, 3, 1, 3, 5, 5, 4, -2, 2, 4, 4, 5, 3, 6, 6, 5, -65
Offset: 1

Views

Author

Omar E. Pol, Dec 04 2011

Keywords

Comments

Here the rank of a "region" is defined to be the largest part minus the number of parts (the same idea as the Dyson's rank of a partition).
Also triangle read by rows: T(j,k) = rank of the k-th region of the last section of the set of partitions of j.
The sum of every row is equal to zero.
Note that in some rows there are several negative terms. - Omar E. Pol, Oct 27 2012
For the definition of "region" see A206437. See also A225600 and A225610. - Omar E. Pol, Aug 12 2013

Examples

			In the triangle T(j,k) for j = 6 the number of regions in the last section of the set of partitions of 6 is equal to 4. The first region given by [2] has rank 2-1 = 1. The second region given by [4,2] has rank 4-2 = 2. The third region given by [3] has rank 3-1 = 2. The fourth region given by [6,3,2,2,1,1,1,1,1,1,1] has rank 6-11 = -5 (see below):
From _Omar E. Pol_, Aug 12 2013: (Start)
---------------------------------------------------------
.    Regions       Illustration of ranks of the regions
---------------------------------------------------------
.    For J=6        k=1     k=2      k=3        k=4
.  _ _ _ _ _ _                              _ _ _ _ _ _
. |_ _ _      |                     _ _ _   .          |
. |_ _ _|_    |           _ _ _ _   * * .|    .        |
. |_ _    |   |     _ _   * * .  |              .      |
. |_ _|_ _|_  |     * .|        .|                .    |
.           | |                                     .  |
.           | |                                       .|
.           | |                                       *|
.           | |                                       *|
.           | |                                       *|
.           | |                                       *|
.           |_|                                       *|
.
So row 6 lists:     1       2         2              -5
(End)
Written as a triangle begins:
0;
0;
0;
1,-1;
2,-2;
1,2,2,-5;
2,3,3,-8;
1,2,2,2,4,3,-14;
2,3,3,3,2,4,4,-21;
1,2,2,2,4,3,1,3,5,5,4,-32;
2,3,3,3,2,4,4,1,4,3,5,6,5,-45;
1,2,2,2,4,3,1,3,5,5,4,-2,2,4,4,5,3,6,6,5,-65;
2,3,3,3,2,4,4,1,4,3,5,6,5,-3,3,5,5,4,5,4,7,7,6,-88;
		

Crossrefs

Row j has length A187219(j). The absolute value of the last term of row j is A000094(j+1). Row sums give A000004.

Formula

a(n) = A141285(n) - A194446(n). - Omar E. Pol, Dec 05 2011

A040051 Parity of partition function A000041.

Original entry on oeis.org

1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1
Offset: 0

Views

Author

Keywords

Comments

From M. V. Subbarao (m.v.subbarao(AT)ualberta.ca), Sep 05 2003: (Start)
Essentially this same question was raised by Ramanujan in a letter to P. A. MacMahon around 1920 (see page 1087, MacMahon's Collected Papers). With the help of Jacobi's triple product identity, MacMahon showed that p(1000) is odd (as he says, with five minutes work - there were no computers those days).
Now we know that among the first ten million values of p(n) 5002137 of them are odd. It is conjectured (T. R. Parkin and D. Shanks) that p(n) is equally often even and odd. Lower bound estimates for the number of times p(n) is even among the first N values of p(n) for any given N are known (Scott Ahlgren; and Nicolas, Rusza and Sárközy among others).
Earlier this year a remarkable result was proved by Boylan and Ahlgren (AMS ABSTRACT # 987-11-82) which says that beyond the three eighty-year old Ramanujan congruences - namely, p(5n+4), p(7n+5) and p(11n +6) being divisible respectively by 5,7 and 11 - there are no other simple congruences of this kind.
My 1966 conjecture that in every arithmetic progression r (mod s) for arbitrary integral r and s, there are infinitely many integers n for which p(n) is odd - with a similar statement for p(n) even - was proved for the even case by Ken Ono (1996) and for the odd case for all s up to 10^5 and for all s which are powers of 2 by Bolyan and Ono, 2002.
(End)
a(n) is also the parity of the trace Tr(n) = A183011(n), the numerator of the Bruinier-Ono formula for the partition function, if n >= 1. - Omar E. Pol, Mar 14 2012
Consider the diagram of the regions of n (see A206437). Then, in each odd-indexed region of n, fill each part of size k with k 1's. Then, in each even-indexed region of n, fill each part of size k with k 0's. The successive digits of row 1 of the diagram give the first n elements of this sequence, if n >= 1. - Omar E. Pol, May 02 2012

References

  • H. Gupta, A note on the parity of p(n), J. Indian Math. Soc. (N.S.) 10, (1946). 32-33. MR0020588 (8,566g)
  • K. M. Majumdar, On the parity of the partition function p(n), J. Indian Math. Soc. (N.S.) 13, (1949). 23-24. MR0030553 (11,13d)
  • M. V. Subbarao, A note on the parity of p(n), Indian J. Math. 14 (1972), 147-148. MR0357355 (50 #9823)

Crossrefs

Programs

  • Haskell
    import Data.Bits (xor)
    a040051 n = p 1 n :: Int where
       p _ 0 = 1
       p k m | k <= m = p k (m - k) `xor` p (k+1) m | k > m  = 0
    -- Reinhard Zumkeller, Nov 15 2011
    
  • Mathematica
    Table[ Mod[ PartitionsP@ n, 2], {n, 105}] (* Robert G. Wilson v, Mar 25 2011 *)
  • PARI
    a(n)=if(n<0, 0, numbpart(n)%2)
    
  • PARI
    a(n)=if(n<0, 0, polcoeff(1/eta(x+x*O(x^n)), n)%2)
    
  • PARI
    a(n)=if(n<10^9, return(numbpart(n)%2)); my(r=n%4, u=select(k->k^2%32==8*r+1,[1..31]), st=u[1], m=n\4, s); u=[u[2]-u[1],u[3]-u[2],u[4]-u[3],u[1]+32-u[4]]; forstep(t=[1,3,7,5][r+1], sqrtint(32*m-1), u, k=t^2>>5; if(a(m-k), s++)); s%2 \\ Merca's algorithm, switching to direct computation for n less than 10^9. Very time-consuming but low memory use. - Charles R Greathouse IV, Jan 24 2018
    
  • Python
    from sympy import npartitions
    def a(n): return npartitions(n)%2 # Indranil Ghosh, May 25 2017

Formula

a(n) = pp(n, 1), with Boolean pp(n, k) = if kReinhard Zumkeller, Sep 04 2003
a(n) = Pm(n,1) with Pm(n,k) = if kReinhard Zumkeller, Jun 09 2009
a(n) = A000035(A000041(n)). - Omar E. Pol, Aug 05 2013
a(n) = A000035(A000025(n)). - John M. Campbell, Jun 29 2016

A175804 Square array A(n,k), n>=0, k>=0, read by antidiagonals: A(n,k) is the n-th term of the k-th differences of partition numbers A000041.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, -1, 0, 1, 3, 2, 1, 1, 2, 5, -4, -2, -1, 0, 2, 7, 9, 5, 3, 2, 2, 4, 11, -21, -12, -7, -4, -2, 0, 4, 15, 49, 28, 16, 9, 5, 3, 3, 7, 22, -112, -63, -35, -19, -10, -5, -2, 1, 8, 30, 249, 137, 74, 39, 20, 10, 5, 3, 4, 12, 42, -539, -290, -153, -79, -40, -20, -10, -5, -2, 2, 14, 56
Offset: 0

Views

Author

Alois P. Heinz, Dec 04 2010

Keywords

Comments

Odlyzko showed that the k-th differences of A000041(n) alternate in sign with increasing n up to a certain index n_0(k) and then stay positive.
Are there any zeros after the first four, which all lie in columns k = 1, 2? - Gus Wiseman, Dec 15 2024

Examples

			Square array A(n,k) begins:
   1,  0,  1, -1,  2,  -4,   9,  ...
   1,  1,  0,  1, -2,   5, -12,  ...
   2,  1,  1, -1,  3,  -7,  16,  ...
   3,  2,  0,  2, -4,   9, -19,  ...
   5,  2,  2, -2,  5, -10,  20,  ...
   7,  4,  0,  3, -5,  10, -20,  ...
  11,  4,  3, -2,  5, -10,  22,  ...
		

Crossrefs

Columns k=0-5 give: A000041, A002865, A053445, A072380, A081094, A081095.
Main diagonal gives A379378.
For primes we have A095195 or A376682.
Row n = 0 is A281425.
Row n = 1 is A320590 except first term.
For composites we have A377033.
For squarefree numbers we have A377038.
For nonsquarefree numbers we have A377046.
For prime powers we have A377051.
Antidiagonal sums are A377056, absolute value version A378621.
The version for strict partitions is A378622, first column A293467.
A000009 counts strict integer partitions, differences A087897, A378972.

Programs

  • Maple
    A41:= combinat[numbpart]:
    DD:= proc(p) proc(n) option remember; p(n+1) -p(n) end end:
    A:= (n,k)-> (DD@@k)(A41)(n):
    seq(seq(A(n, d-n), n=0..d), d=0..11);
  • Mathematica
    max = 11; a41 = Array[PartitionsP, max+1, 0]; a[n_, k_] := Differences[a41, k][[n+1]]; Table[a[n, k-n], {k, 0, max}, {n, 0, k}] // Flatten (* Jean-François Alcover, Aug 29 2014 *)
    nn=5;Table[Table[Sum[(-1)^(k-i)*Binomial[k,i]*PartitionsP[n+i],{i,0,k}],{k,0,nn}],{n,0,nn}] (* Gus Wiseman, Dec 15 2024 *)

Formula

A(n,k) = (Delta^(k) A000041)(n).
A(n,k) = Sum_{i=0..k} (-1)^(k-i) * binomial(k,i) * A000041(n+i). In words, row x is the inverse zero-based binomial transform of A000041 shifted left x times. - Gus Wiseman, Dec 15 2024

A186412 Sum of all parts in the n-th region of the set of partitions of j, if 1<=n<=A000041(j).

Original entry on oeis.org

1, 3, 5, 2, 9, 3, 12, 2, 6, 3, 20, 3, 7, 4, 25, 2, 6, 3, 13, 5, 4, 38, 3, 7, 4, 14, 3, 9, 5, 49, 2, 6, 3, 13, 5, 4, 23, 4, 10, 6, 5, 69, 3, 7, 4, 14, 3, 9, 5, 27, 5, 4, 15, 7, 6, 87, 2, 6, 3, 13, 5, 4, 23, 4, 10, 6, 5, 39, 3, 9, 5, 19, 4, 12, 7, 6, 123
Offset: 1

Views

Author

Omar E. Pol, Aug 12 2011

Keywords

Comments

Also triangle read by rows: T(j,k) = sum of all parts in the k-th region of the last section of the set of partitions of j. See Example section. For more information see A135010. - Omar E. Pol, Nov 26 2011
For the definition of "region" see A206437. - Omar E. Pol, Aug 19 2013

Examples

			Contribution from Omar E. Pol, Nov 26 2011 (Start):
Written as a triangle:
1;
3;
5;
2,9;
3,12;
2,6,3,20;
3,7,4,25;
2,6,3,13,5,4,38;
3,7,4,14,3,9,5,49;
2,6,3,13,5,4,23,4,10,6,5,69;
3,7,4,14,3,9,5,27,5,4,15,7,6,87;
2,6,3,13,5,4,23,4,10,6,5,39,3,9,5,19,4,12,7,6,123;
(End)
From _Omar E. Pol_, Aug 18 2013: (Start)
Illustration of initial terms (first seven regions):
.                                             _ _ _ _ _
.                                     _ _ _  |_ _ _ _ _|
.                           _ _ _ _  |_ _ _|       |_ _|
.                     _ _  |_ _ _ _|                 |_|
.             _ _ _  |_ _|     |_ _|                 |_|
.       _ _  |_ _ _|             |_|                 |_|
.   _  |_ _|     |_|             |_|                 |_|
.  |_|   |_|     |_|             |_|                 |_|
.
.   1     3       5     2         9       3          12
.
(End)
		

Crossrefs

Row sums of triangle A186114 and of triangle A193870.
Row j has length A187219(j).
Row sums give A138879.
Right border gives A046746, j >= 1.
Records give A046746, j >= 1.
Partial sums give A182244.

Programs

  • Mathematica
    lex[n_]:=DeleteCases[Sort@PadRight[Reverse /@ IntegerPartitions@n], x_ /; x==0,2];
    A186412 = {}; l = {};
    For[j = 1, j <= 50, j++,
      mx = Max@lex[j][[j]]; AppendTo[l, mx];
      For[i = j, i > 0, i--, If[l[[i]] > mx, Break[]]];
      AppendTo[A186412, Total@Take[Reverse[First /@ lex[mx]], j - i]];
      ];
    A186412  (* Robert Price, Jul 25 2020 *)

Formula

a(A000041(n)) = A046746(n).
Previous Showing 11-20 of 3898 results. Next