A005380 -id:A005380 - cp's OEIS frontend

A000219 Number of plane partitions (or planar partitions) of n.

1, 1, 3, 6, 13, 24, 48, 86, 160, 282, 500, 859, 1479, 2485, 4167, 6879, 11297, 18334, 29601, 47330, 75278, 118794, 186475, 290783, 451194, 696033, 1068745, 1632658, 2483234, 3759612, 5668963, 8512309, 12733429, 18974973, 28175955, 41691046, 61484961, 90379784, 132441995, 193487501, 281846923

Offset: 0

N. J. A. Sloane

Two-dimensional partitions of n in which no row or column is longer than the one before it (compare A001970). E.g., a(4) = 13:
4.31.3.22.2.211.21..2.1111.111.11.11.1 but not 2
.....1....2.....1...1......1...11.1..1........ 11
....................1.............1..1
.....................................1
In the above, one also must require that rows & columns are nondecreasing, e.g., [1,1; 2] is also forbidden (which implies that row and column lengths are nondecreasing, if empty cells are identified with cells filled with 0's). - M. F. Hasler, Sep 22 2018
Can also be regarded as number of "safe pilings" of cubes in the corner of a room: the height should not increase away from the corner. - Wouter Meeussen
Also number of partitions of n objects of 2 colors, each part containing at least one black object; see example. - Christian G. Bower, Jan 08 2004
Number of partitions of n into 1 type of part 1, 2 types of part 2, ..., k types of part k. E.g., n=3 gives 111, 12, 12', 3, 3', 3''. - Jon Perry, May 27 2004
The bijection between the partitions in the two preceding comments goes by identifying a part with k black objects with a part of type k. - David Scambler and Joerg Arndt, May 01 2013
Can also be regarded as the number of Jordan canonical forms for an n X n matrix. (I.e., a 5 X 5 matrix has 24 distinct Jordan canonical forms, dependent on the algebraic and geometric multiplicity of each eigenvalue.) - Aaron Gable (agable(AT)hmc.edu), May 26 2009
(1/n) * convolution product of n terms * A001157 (sum of squares of divisors of n): (1, 5, 10, 21, 26, 50, 50, 85, ...) = a(n). As shown by [Bressoud, p. 12]: 1/6 * [1*24 + 5*13 + 10*6 + 21*3 + 26*1 + 50*1] = 288/6 = 48. - Gary W. Adamson, Jun 13 2009
Convolved with the aerated version (1, 0, 1, 0, 3, 0, 6, 0, 13, ...) = A026007: (1, 1, 2, 5, 8, 16, 28, 49, 83, ...). - Gary W. Adamson, Jun 13 2009
Starting with offset 1 = row sums of triangle A162453. - Gary W. Adamson, Jul 03 2009
Unfortunately, Wright's formula is also incomplete in the paper by G. Almkvist: "Asymptotic formulas and generalized Dedekind sums", p. 344, (the denominator should have sqrt(3*Pi) not sqrt(Pi)). This error was already corrected in the paper by Steven Finch: "Integer Partitions". - Vaclav Kotesovec, Aug 17 2015
Also the number of non-isomorphic weight-n chains of multisets whose dual is also a chain of multisets. The dual of a multiset partition has, for each vertex, one block consisting of the indices (or positions) of the blocks containing that vertex, counted with multiplicity. The weight of a multiset partition is the sum of sizes of its parts. - Gus Wiseman, Sep 25 2018

			A planar partition of 13:
  4 3 1 1
  2 1
  1
a(5) = (1/5!)*(sigma_2(1)^5+10*sigma_2(2)*sigma_2(1)^3+20*sigma_2(3)*sigma_2(1)^2+ 15*sigma_2(1)*sigma_2(2)^2+30*sigma_2(4)*sigma_2(1)+20*sigma_2(2)*sigma_2(3)+24*sigma_2(5)) = 24. - _Vladeta Jovovic_, Jan 10 2003
From _David Scambler_ and _Joerg Arndt_, May 01 2013: (Start)
There are a(4) = 13 partitions of 4 objects of 2 colors ('b' and 'w'), each part containing at least one black object:
1 black part:
  [ bwww ]
2 black parts:
  [ bbww ]
  [ bww, b ]
  [ bw, bw ]
3 black parts:
  [ bbbw ]
  [ bbw, b ]
  [ bb, bw ]
(but not: [bw, bb ] )
  [ bw, b, b ]
4 black parts:
  [ bbbb ]
  [ bbb, b ]
  [ bb, bb ]
  [ bb, b, b ]
  [ b, b, b, b ]
(End)
From _Geoffrey Critzer_, Nov 29 2014: (Start)
The corresponding partitions of the integer 4 are:
  4'''
  4''
  3'' + 1
  2' + 2'
  4'
  3' + 1
  2 + 2'
  2' + 1 + 1
  4
  3 + 1
  2 + 2
  2 + 1 + 1
  1 + 1 + 1 + 1.
(End)
From _Gus Wiseman_, Sep 25 2018: (Start)
Non-isomorphic representatives of the a(4) = 13 chains of multisets whose dual is also a chain of multisets:
  {{1,1,1,1}}
  {{1,1,2,2}}
  {{1,2,2,2}}
  {{1,2,3,3}}
  {{1,2,3,4}}
  {{1},{1,1,1}}
  {{2},{1,2,2}}
  {{3},{1,2,3}}
  {{1,1},{1,1}}
  {{1,2},{1,2}}
  {{1},{1},{1,1}}
  {{2},{2},{1,2}}
  {{1},{1},{1},{1}}
(End)
G.f. = 1 + x + 3*x^2 + 6*x^3 + 13*x^4 + 24*x^5 + 48*x^6 + 86*x^7 + 160*x^8 + ...

G. Almkvist, The differences of the number of plane partitions, Manuscript, circa 1991.
G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976, p. 241.
D. M. Bressoud, Proofs and Confirmations, Camb. Univ. Press, 1999; pp(n) on p. 10.
Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 575.
L. Carlitz, Generating functions and partition problems, pp. 144-169 of A. L. Whiteman, ed., Theory of Numbers, Proc. Sympos. Pure Math., 8 (1965). Amer. Math. Soc., see p. 145, eq. (1.6).
I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, Wiley, N.Y., 1983, (5.4.5).
P. A. MacMahon, Memoir on the theory of partitions of numbers - Part VI, Phil. Trans. Royal Soc., 211 (1912), 345-373.
P. A. MacMahon, Combinatory Analysis. Cambridge Univ. Press, London and New York, Vol. 1, 1915 and Vol. 2, 1916; see vol. 2, p 332.
P. A. MacMahon, The connexion between the sum of the squares of the divisors and the number of partitions of a given number, Messenger Math., 54 (1924), 113-116. Collected Papers, MIT Press, 1978, Vol. I, pp. 1364-1367. See Table II. - N. J. A. Sloane, May 21 2014
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).

Suresh Govindarajan, Table of n, a(n) for n = 0..6500 (first 401 terms from T. D. Noe)
G. Almkvist, Asymptotic formulas and generalized Dedekind sums, Exper. Math., 7 (No. 4, 1998), pp. 343-359.
G. E. Andrews and P. Paule, MacMahon's partition analysis XII: Plane Partitions, J. Lond. Math. Soc., 76 (2007), 647-666.
A. O. L. Atkin, P. Bratley, I. G. McDonald and J. K. S. McKay, Some computations for m-dimensional partitions, Proc. Camb. Phil. Soc., 63 (1967), 1097-1100.
A. O. L. Atkin, P. Bratley, I. G. McDonald and J. K. S. McKay, Some computations for m-dimensional partitions, Proc. Camb. Phil. Soc., 63 (1967), 1097-1100. [Annotated scanned copy]
Michael Beeler, R. William Gosper and Richard C. Schroeppel, HAKMEM, ITEM 18, Memo AIM-239, Artificial Intelligence Laboratory, Massachusetts Institute of Technology, Cambridge, Mass., 1972.
Edward A. Bender, Asymptotic methods in enumeration, SIAM Review 16 (1974), no. 4, p. 509.
E. A. Bender and D. E. Knuth, Enumeration of Plane Partitions, J. Combin. Theory A. 13, 40-54, 1972.
S. Benvenuti, B. Feng, A. Hanany and Y. H. He, Counting BPS operators in gauge theories: Quivers, syzygies and plethystics, arXiv:hep-th/0608050, p. 41-42.
Henry Bottomley, Illustration of initial terms
D. M. Bressoud and J. Propp, How the alternating sign matrix conjecture was solved, Notices Amer. Math. Soc., 46 (No. 6, 1999), 637-646.
Shouvik Datta, M. R. Gaberdiel, W. Li, and C. Peng, Twisted sectors from plane partitions, arXiv preprint arXiv:1606.07070 [hep-th], 2016.
Wenjie Fang, Hsien-Kuei Hwang, and Mihyun Kang, Phase transitions from exp(n^(1/2)) to exp(n^(2/3)) in the asymptotics of banded plane partitions, arXiv:2004.08901 [math.CO], 2020.
Steven Finch, Integer Partitions, September 22, 2004. [Cached copy, with permission of the author]
P. Flajolet and R. Sedgewick, Analytic Combinatorics, 2009; see page 580.
Bernhard Heim, Markus Neuhauser and Robert Tröger, Inequalities for Plane Partitions, arXiv:2109.15145 [math.CO], 2021.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 141
Vaclav Kotesovec, A method of finding the asymptotics of q-series based on the convolution of generating functions, arXiv:1509.08708 [math.CO], 2015-2016, p. 18.
Vaclav Kotesovec, Graphs - The asymptotic ratio (250000 terms)
D. E. Knuth, A Note on Solid Partitions, Math. Comp. 24, 955-961, 1970.
Oleg Lazarev, Matt Mizuhara and Ben Reid, Some Results in Partitions, Plane Partitions, and Multipartitions, 13 August 2010.
P. A. MacMahon, Combinatory analysis.
J. Mangual, McMahon's Formula via Free Fermions, arXiv preprint arXiv:1210.7109 [math.CO], 2012. - From _N. J. A. Sloane_, Jan 01 2013
Ville Mustonen and R. Rajesh, Numerical Estimation of the Asymptotic Behaviour of Solid Partitions ..., arXiv:cond-mat/0303607 [cond-mat.stat-mech], 2003.
L. Mutafchiev and E. Kamenov, On The Asymptotic Formula for the Number of Plane Partitions..., arXiv:math/0601253 [math.CO], 2006; C. R. Acad. Bulgare Sci. 59(2006), No. 4, 361-366.
Ken Ono, Sudhir Pujahari and Larry Rolen, Turán inequalities for the plane partition function, arXiv:2201.01352 [math.NT], 2022.
I. Pak, Partition bijections, a survey, Ramanujan J. 12 (2006) 5-75.
A. Rovenchak, Enumeration of plane partitions with a restricted number of parts, arXiv preprint arXiv:1401.4367 [math-ph], 2014.
Raphael Schumacher, The self-counting identity, Fib. Quart., 55 (No. 2 2017), 157-167.
N. J. A. Sloane, Transforms
J. Stienstra, Mahler measure, Eisenstein series and dimers, arXiv:math/0502197 [math.NT], 2005.
Balázs Szendrői, Non-commutative Donaldson-Thomas invariants and the conifold, Geometry & Topology 12.2 (2008): 1171-1202.
Eric Weisstein's World of Mathematics, Plane Partition
E. M. Wright, Rotatable partitions, J. London Math. Soc., 43 (1968), 501-505.
Index entries for "core" sequences

Cf. A000784, A000785, A000786, A005380, A005987, A048141, A048142, A089300.
Cf. A023871-A023878, A026007, A001157, A162453, A285216.
Differences: A191659, A191660, A191661.
Row sums of A089353 and A091438 and A091298.
Column k=1 of A144048. - Alois P. Heinz, Nov 02 2012
Sequences "number of r-line partitions": A000041 (r=1), A000990 (r=2), A000991 (r=3), A002799 (r=4), A001452 (r=5), A225196 (r=6), A225197 (r=7), A225198 (r=8), A225199 (r=9).
Cf. A249386, A249387.
Cf. A161870, A255610, A255611, A255612, A255613, A255614, A193427.

using Nemo, Memoize
@memoize function a(n)
    if n == 0 return 1 end
    s = sum(a(n - j) * divisor_sigma(j, 2) for j in 1:n)
    return div(s, n)
end
[a(n) for n in 0:20] # Peter Luschny, May 03 2020

series(mul((1-x^k)^(-k),k=1..64),x,63);
# second Maple program:
a:= proc(n) option remember; `if`(n=0, 1, add(
      a(n-j)*numtheory[sigma][2](j), j=1..n)/n)
    end:
seq(a(n), n=0..50);  # Alois P. Heinz, Aug 17 2015

CoefficientList[Series[Product[(1 - x^k)^-k, {k, 64}], {x, 0, 64}], x]
Zeta[3]^(7/36)/2^(11/36)/Sqrt[3 Pi]/Glaisher E^(3 Zeta[3]^(1/3) (n/2)^(2/3) + 1/12)/n^(25/36) (* asymptotic formula after Wright; Vaclav Kotesovec, Jun 23 2014 *)
a[0] = 1; a[n_] := a[n] = Sum[a[n - j] DivisorSigma[2, j], {j, n}]/n; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Sep 21 2015, after Alois P. Heinz *)
CoefficientList[Series[Exp[Sum[DivisorSigma[2, n] x^n/n, {n, 50}]], {x, 0, 50}], x] (* Eric W. Weisstein, Feb 01 2018 *)

{a(n) = if( n<0, 0, polcoeff( exp( sum( k=1, n, x^k / (1 - x^k)^2 / k, x * O(x^n))), n))}; /* Michael Somos, Jan 29 2005 */

{a(n) = if( n<0, 0, polcoeff( prod( k=1, n, (1 - x^k + x * O(x^n))^-k), n))}; /* Michael Somos, Jan 29 2005 */

my(N=66, x='x+O('x^N)); Vec( prod(n=1,N, (1-x^n)^-n) ) \\ Joerg Arndt, Mar 25 2014

A000219(n)=#PlanePartitions(n) \\ See A091298 for PlanePartitions(). For illustrative use: much slower than the above. - M. F. Hasler, Sep 24 2018

from sympy import cacheit
from sympy.ntheory import divisor_sigma
@cacheit
def A000219(n):
    if n <= 1:
        return 1
    return sum(A000219(n - k) * divisor_sigma(k, 2) for k in range(1, n + 1)) // n
print([A000219(n) for n in range(20)])
# R. J. Mathar, Oct 18 2009

# uses[EulerTransform from A166861]
b = EulerTransform(lambda n: n)
print([b(n) for n in range(37)]) # Peter Luschny, Nov 11 2020

G.f.: Product_{k >= 1} 1/(1 - x^k)^k. - MacMahon, 1912.
Euler transform of sequence [1, 2, 3, ...].
a(n) ~ (c_2 / n^(25/36)) * exp( c_1 * n^(2/3) ), where c_1 = A249387 = 2.00945... and c_2 = A249386 = 0.23151... - Wright, 1931. Corrected Jun 01 2010 by Rod Canfield - see Mutafchiev and Kamenov. The exact value of c_2 is e^(2c)*2^(-11/36)*zeta(3)^(7/36)*(3*Pi)^(-1/2), where c = Integral_{y=0..inf} (y*log(y)/(e^(2*Pi*y)-1))dy = (1/2)*zeta'(-1).
The exact value of c_1 is 3*2^(-2/3)*Zeta(3)^(1/3) = 2.0094456608770137530649... - Vaclav Kotesovec, Sep 14 2014
a(n) = (1/n) * Sum_{k=1..n} a(n-k)*sigma_2(k), n > 0, a(0)=1, where sigma_2(n) = A001157(n) = sum of squares of divisors of n. - Vladeta Jovovic, Jan 20 2002
G.f.: exp(Sum_{n>0} sigma_2(n)*x^n/n). a(n) = Sum_{pi} Product_{i=1..n} binomial(k(i)+i-1, k(i)) where pi runs through all nonnegative solutions of k(1)+2*k(2)+..+n*k(n)=n. - Vladeta Jovovic, Jan 10 2003
From Vaclav Kotesovec, Nov 07 2016: (Start)
More precise asymptotics: a(n) ~ Zeta(3)^(7/36) * exp(3 * Zeta(3)^(1/3) * (n/2)^(2/3) + 1/12) / (A * sqrt(3*Pi) * 2^(11/36) * n^(25/36))
* (1 + c1/n^(2/3) + c2/n^(4/3) + c3/n^2), where
c1 = -0.23994424421250649114273759... = -277/(864*(2*Zeta(3))^(1/3)) - Zeta(3)^(2/3)/(1440*2^(1/3))
c2 = -0.02576771365117401620018082... = 353*Zeta(3)^(1/3)/(248832*2^(2/3)) - 17*Zeta(3)^(4/3)/(3225600*2^(2/3)) - 71575/(1492992*(2*Zeta(3))^(2/3))
c3 = -0.00533195302658826100834286... = -629557/859963392 - 42944125/(7739670528*Zeta(3)) + 14977*Zeta(3)/1114767360 - 22567*Zeta(3)^2/250822656000
and A = A074962 is the Glaisher-Kinkelin constant.
(End)

Corrected by N. J. A. Sloane, Jul 29 2006

Minor edits by Vaclav Kotesovec, Oct 27 2014

A052847 G.f.: 1 / Product_{k>=1} (1-x^k)^(k-1).

Original entry on oeis.org

1, 0, 1, 2, 4, 6, 12, 18, 33, 52, 88, 138, 229, 354, 568, 880, 1378, 2110, 3260, 4942, 7527, 11320, 17031, 25394, 37842, 55956, 82630, 121300, 177677, 258980, 376626, 545352, 787784, 1133764, 1627657, 2329020, 3324559, 4731396, 6717774, 9512060

Offset: 0

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Euler transform of sequence [0,1,2,3,...]. - Michael Somos, Jul 02 2004
Number of partitions of n objects of 2 colors, where each part must contain at least one of each color. - Franklin T. Adams-Watters, Jan 23 2006
Number of partitions of n without 1s, one kind of 2s, two kinds of 3s, etc. - Joerg Arndt, Jul 31 2011
From Vaclav Kotesovec, Oct 17 2015: (Start)
In general, if v>=0 and g.f. = Product_{k>=1} 1/(1-x^(k+v))^k, then a(n) ~ d1(v) * n^(v^2/6 - 25/36) * exp(-Pi^4 * v^2 / (432*Zeta(3)) + 3*Zeta(3)^(1/3) * n^(2/3)/2^(2/3) - v * Pi^2 * n^(1/3) / (3 * 2^(4/3) * Zeta(3)^(1/3))) / (sqrt(3*Pi) * 2^(v^2/6 + 11/36) * Zeta(3)^(v^2/6 - 7/36)), where Zeta(3) = A002117.
d1(v) = exp(Integral_{x=0..infinity} (1/(x*exp((v-1)*x) * (exp(x)-1)^2) - (6*v^2-1) / (12*x*exp(x)) + v/x^2 - 1/x^3) dx).
d1(v) = (exp(Zeta'(-1) - v*Zeta'(0))) / Product_{j=0..v-1} j!, where Zeta'(0) = -A075700, Zeta'(-1) = A084448 and Product_{j=0..v-1} j! = A000178(v-1).
d1(v) = exp(1/12) * (2*Pi)^(v/2) / (A * G(v+1)), where A = A074962 is the Glaisher-Kinkelin constant and G is the Barnes G-function.
(End)

			1 + x^2 + 2*x^3 + 4*x^4 + 6*x^5 + 12*x^6 + 18*x^7 + 33*x^8 + 52*x^9 + ...
From _Gus Wiseman_, Jan 22 2019: (Start)
The partitions described in Franklin T. Adams-Watters's comment are (n = 2 through 6):
  {{12}}  {{112}}  {{1112}}    {{11112}}    {{111112}}
          {{122}}  {{1122}}    {{11122}}    {{111122}}
                   {{1222}}    {{11222}}    {{111222}}
                   {{12}{12}}  {{12222}}    {{112222}}
                               {{12}{112}}  {{122222}}
                               {{12}{122}}  {{112}{112}}
                                            {{112}{122}}
                                            {{12}{1112}}
                                            {{12}{1122}}
                                            {{12}{1222}}
                                            {{122}{122}}
                                            {{12}{12}{12}}
(End)

Vaclav Kotesovec, Table of n, a(n) for n = 0..1000
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 815
Vaclav Kotesovec, Graph - The asymptotic ratio

Cf. A005380, A000203, A001157, A052812.
Cf. A000219 (v=0), A052847 (v=1), A263358 (v=2), A263359 (v=3), A263360 (v=4), A263361 (v=5), A263362 (v=6), A263363 (v=7), A263364 (v=8).
Cf. A054974, A321407, A321760, A323654, A323655, A323656.

spec := [S,{B=Sequence(Z,1 <= card),C=Prod(B,B),S= Set(C)},unlabeled]: seq(combstruct[count](spec, size=n), n=0..20);
with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d, j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: a:=etr(n-> n-1): seq(a(n), n=0..50); # Vaclav Kotesovec, Mar 04 2015 after Alois P. Heinz

Clear[a]; a[n_]:= a[n] = 1/n*Sum[(DivisorSigma[2,k]-DivisorSigma[1,k])*a[n-k],{k,1,n}]; a[0]=1; Table[a[n],{n,0,100}] (* Vaclav Kotesovec, Mar 04 2015 *)
nmax = 40; CoefficientList[Series[Product[1/(1-x^(k+1))^k, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Oct 16 2015 *)

{a(n) = if( n<0, 0, polcoeff( 1 / prod(k=1, n, (1 - x^k + x*O(x^n))^(k-1)), n))}

a(n) = 1/n*Sum_{k=1..n} (sigma[2](k)-sigma[1](k))*a(n-k).
G.f.: exp( Sum_{k>0} ( x^k / (1 - x^k) )^2 / k ).
G.f.: exp( sum(n>=0, (sigma[2](n)-sigma[1](n)) *x^n/n ) ). - Joerg Arndt, Jul 31 2011
a(n) ~ 2^(1/36) * Zeta(3)^(1/36) * exp(1/12 - Pi^4/(432*Zeta(3)) - Pi^2 * n^(1/3) / (3 * 2^(4/3) * Zeta(3)^(1/3)) + 3 * Zeta(3)^(1/3) * n^(2/3) / 2^(2/3)) / (A * 3^(1/2) * n^(19/36)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Mar 07 2015

Edited by Vladeta Jovovic, Sep 10 2002

A054225 Triangle read by rows: row n (n>=0) gives the number of partitions of (n,0), (n-1,1), (n-2,2), ..., (0,n) respectively into sums of pairs.

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 3, 4, 4, 3, 5, 7, 9, 7, 5, 7, 12, 16, 16, 12, 7, 11, 19, 29, 31, 29, 19, 11, 15, 30, 47, 57, 57, 47, 30, 15, 22, 45, 77, 97, 109, 97, 77, 45, 22, 30, 67, 118, 162, 189, 189, 162, 118, 67, 30, 42, 97, 181, 257, 323, 339, 323, 257, 181, 97, 42, 56, 139, 267, 401, 522, 589, 589, 522, 401, 267, 139, 56

Offset: 0

Marc LeBrun, Feb 04 2000

By analogy with ordinary partitions (A000041). The empty partition gives T(0,0)=1 by definition. A054225 and A201377 give partitions of pairs into sums of distinct pairs. Parts (i,j) are "positive" in the sense that min {i,j} >= 0 and max {i,j} >0. The empty partition of (0,0) is counted as 1.
Or, triangle T(n,k) of bipartite partitions of n objects, k of which are black.
Or, number of ways to factor p^(n-k)*q^k where p and q are distinct primes.
In the paper by F. C. Auluck: "On partitions of bipartite numbers", p.74, in the formula for fixed m there should be factor 1/m!. The correct asymptotic formula is p(m, n) ~ (sqrt(6*n)/Pi)^m * exp(Pi*sqrt(2*n/3)) / (4*sqrt(3)*m!*n). - Vaclav Kotesovec, Feb 01 2016
T(n,k)=T(n,k-n) is the number of multiset partitions of the multiset {1^k, 2^(n-k)}, see example link. - Joerg Arndt, Jan 01 2024
Let R be the ring of power series in two countably infinite sets of variables x_1,y_1,x_2,y_2,... that are invariant under the diagonal action (i.e, the group S of permutations of positive integers acts by w(x_i)=x_{w(i)} and w(y_i)=y_{w(i)}).  Then T(n,k) is the dimension of the (n,k)-bigraded piece of R, i.e., the bihomogeneous power series of degree n in the x-variables and k in the y-variables that are S-invariant. - Jeremy L. Martin, Nov 27 2024

			The second row (n=1) is 1,1 since (1,0) and (0,1) each have a single partition.
The third row (n=2) is 2, 2, 2 from (2,0) = (1,0)+(1,0), (1,1) = (1,0)+(0,1), (0,2) = (0,1)+(0,1).
In the fourth row (n=3), T(2,1)=4 from (2,1) = (2,0)+(0,1) = (1,0)+(1,1) = (1,0)+(1,0)+(0,1).
The triangle begins:
   1;
   1,  1;
   2,  2,  2;
   3,  4,  4,  3;
   5,  7,  9,  7,   5;
   7, 12, 16, 16,  12,  7;
  11, 19, 29, 31,  29, 19, 11;
  15, 30, 47, 57,  57, 47, 30, 15;
  22, 45, 77, 97, 109, 97, 77, 45, 22;
  ...
A further example: T(2,2) = 9:
[(2,2)],
[(2,1),(0,1)],
[(2,0),(0,2)],
[(2,0),(0,1),(0,1)],
[(1,2),(1,0)],
[(1,1),(1,1)],
[(1,1),(1,0),(0,1)],
[(1,0),(1,0),(0,2)],
[(1,0),(1,0),(0,1),(0,1)].

M. S. Cheema, Tables of partitions of Gaussian integers, National Institute of Sciences of India, New Delhi, 1956.
D. E. Knuth, The Art of Computer Programming, Vol. 4A, Table A-1, page 778. - N. J. A. Sloane, Dec 30 2018

Alois P. Heinz, Rows n = 0..200, flattened
Joerg Arndt, Multiset partitions of the multisets {1^r,2^s} for r+s=5.
F. C. Auluck, On partitions of bipartite numbers, Proc. Cambridge Philos. Soc. 49, (1953). 72-83.
F. C. Auluck, On partitions of bipartite numbers, annotated scan of a few pages.
F. C. Auluck, On partitions of bipartite numbers, Mathematical Proceedings of the Cambridge Philosophical Society, Volume 49, Issue 01, January 1953, pp. 72-83. (full article)
P. A. MacMahon, Memoir on symmetric functions of the roots of systems of equations, Phil. Trans. Royal Soc. London, 181 (1890), 481-536; Coll. Papers II, 32-87.
Reinhard Zumkeller, Haskell programs for A201376, A054225, A201377, A054242

See A201376 for the same triangle formatted in a different way.
Columns 0-10: A000041, A000070, A000291, A000412, A000465, A000491, A002755, A002756, A002757, A002758, A002759.
Row sums: A005380. a(2n, n): A002774. a(n, [n/2]): A091437. Cf. A060244.
The outer edges are T(n,0) = T(0,n) = A000041(n).
A054242 gives partitions into sums of distinct pairs.

Haskell
```
see Zumkeller link.
```

read transforms; t1 := mul( mul( 1/(1-x^(i-j)*y^j), j=0..i), i=1..11): SERIES2(t1,x,y,6);

rows = 11; se = Series[ Product[ 1/(1-x^(n-k)*y^k), {n, 1, rows}, {k, 0, n}], {x, 0, rows}, {y, 0, rows}]; coes = CoefficientList[ se, {x, y}]; Flatten[ Table[ coes[[n-k+1, k]], {n, 1, rows+1}, {k, 1, n}]] (* Jean-François Alcover, Nov 21 2011, after g.f. *)
p = 2; q = 3; b[n_, k_] := b[n, k] = If[n>k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d>k, 0, b[n/d, d]], {d, DeleteCases[Divisors[n], 1|n]}]]; t[n_, k_] := b[p^(n-k)*q^k, p^(n-k)*q^k]; Table[t[n, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 13 2014, after Alois P. Heinz *)

{T(n, k) = if( n<0 || k<0, 0, polcoeff( polcoeff( prod( i=1, n, prod( j=0, i, 1 / (1 - x^i * y^j), 1 + x * O(x^n))),n),k))} /* Michael Somos, Apr 19 2005 */

G.f.: Product_{i>=1, j=0..i} 1/(1-x^(i-j)*y^j).

Series ends ... + 7*x^5 + 12*x^4*y + 16*x^3*y^2 + 16*x^2*y^3 + 12*x*y^4 + 7*y^5 + 5*x^4 + 7*x^3*y + 9*x^2*y^2 + 7*x*y^3 + 5*y^4 + 3*x^3 + 4*x^2*y + 4*x*y^2 + 3*y^3 + 2*x^2 + 2*x*y + 2*y^2 + x + y + 1.

Entry revised by N. J. A. Sloane, Nov 30 2011, to incorporate corrections provided by Reinhard Zumkeller, who also contributed the alternative version A201376. Once the errors were corrected, this sequence coincided with A060243, due to N. J. A. Sloane, Mar 22 2001, which included edits by Vladeta Jovovic, Mar 23 2001, and Christian G. Bower, Jan 08 2004. The two entries have now been merged.

A002755 Number of bipartite partitions of n white objects and 6 black ones.

Original entry on oeis.org

11, 30, 77, 162, 323, 589, 1043, 1752, 2876, 4571, 7128, 10860, 16306, 24051, 35040, 50355, 71609, 100697, 140349, 193784, 265505, 360889, 487214, 653243, 870613, 1153322, 1519658, 1991689, 2597762, 3372107, 4358198, 5608418, 7188632

Offset: 0

N. J. A. Sloane

nonn

Number of ways to factor p^n*q^6 where p and q are distinct primes.

a(n) is the number of multiset partitions of the multiset {r^n, s^6}. - Joerg Arndt, Jan 01 2024

M. S. Cheema and H. Gupta, Tables of Partitions of Gaussian Integers. National Institute of Sciences of India, Mathematical Tables, Vol. 1, New Delhi, 1956, p. 1.
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).

Alois P. Heinz, Table of n, a(n) for n = 0..1000
M. S. Cheema and H. Gupta, Tables of Partitions of Gaussian Integers. National Institute of Sciences of India, Mathematical Tables, Vol. 1, New Delhi, 1956 (Annotated scanned pages from, plus a review)

Column 6 of A054225.

Cf. A005380.

b[n_, k_] := b[n, k] = If[n>k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d>k, 0, b[n/d, d]], {d, DeleteCases[Divisors[n], 1|n]}]]; a[n_] := b[3^6*2^n, 3^6*2^n]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Mar 13 2014, after Alois P. Heinz *)
nmax = 50; CoefficientList[Series[(11 + 8*x + 6*x^2 - 7*x^4 - 13*x^5 - 19*x^6 - 10*x^7 - 3*x^8 + 7*x^9 + 11*x^10 + 15*x^11 + 6*x^12 - 2*x^14 - 7*x^15 - 4*x^16 - 2*x^17 + 3*x^18 + 2*x^19 - x^20)/((1-x) * (1-x^2) * (1-x^3) * (1-x^4) * (1-x^5) * (1-x^6)) * Product[1/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Feb 01 2016 *)

from sympy import divisors, isprime
from functools import cache
@cache
def T(n, m): # after Indranil Ghosh in A001055
    if isprime(n): return 1 if n <= m else 0
    s = sum(T(n//d, d) for d in divisors(n)[1:-1] if d <= m)
    return s + 1 if n <= m else s
def a(n): return (lambda x: T(x, x))(2**n * 3**6)
print([a(n) for n in range(33)]) # Michael S. Branicky, Nov 30 2021

a(n) = if n <= 6 then A054225(6,n) else A054225(n,6). - Reinhard Zumkeller, Nov 30 2011

a(n) ~ sqrt(3) * n^2 * exp(Pi*sqrt(2*n/3)) / (40*Pi^6). - Vaclav Kotesovec, Feb 01 2016

Edited by Christian G. Bower, Jan 08 2004

A002774 Number of bipartite partitions of n white objects and n black ones.

Original entry on oeis.org

1, 2, 9, 31, 109, 339, 1043, 2998, 8406, 22652, 59521, 151958, 379693, 927622, 2224235, 5236586, 12130780, 27669593, 62229990, 138095696, 302673029, 655627975, 1404599867, 2977831389, 6251060785, 12999299705, 26791990052, 54750235190, 110977389012

Offset: 0

N. J. A. Sloane

nonn

Number of ways to factor p^n * q^n where p and q are distinct primes.

M. S. Cheema and H. Gupta, Tables of Partitions of Gaussian Integers. National Institute of Sciences of India, Mathematical Tables, Vol. 1, New Delhi, 1956, p. 1.
D. E. Knuth, The Art of Computer Programming, Vol. 4A, Table A-1, see p(n,n), page 778. - N. J. A. Sloane, Dec 30 2018
A. Murthy, Generalization of partition function, introducing Smarandache factor partitions. Smarandache Notions Journal, Vol. 11, No. 1-2-3, Spring 2000.
A. Murthy, Program for finding out the number of Smarandache factor partitions. (To be published in Smarandache Notions Journal).
Amarnath Murthy and Charles Ashbacher, Generalized Partitions and Some New Ideas on Number Theory and Smarandache Sequences, Hexis, Phoenix; USA 2005. See Section 1.4, 1.14.
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).

Alois P. Heinz and Vaclav Kotesovec, Table of n, a(n) for n = 0..400 (terms 0..100 from Alois P. Heinz)
F. C. Auluck, On partitions of bipartite numbers, Proc. Cambridge Philos. Soc. 49, (1953). 72-83.
F. C. Auluck, On partitions of bipartite numbers, annotated scan of a few pages.
F. C. Auluck, On partitions of bipartite numbers, Mathematical Proceedings of the Cambridge Philosophical Society, Volume 49, Issue 01, January 1953, pp. 72-83. (full article)
Katherine Ormeño Bastías, Paul Martin, and Steen Ryom-Hansen, On the spherical partition algebra, arXiv:2402.01890 [math.RT], 2024.
M. S. Cheema and H. Gupta, Tables of Partitions of Gaussian Integers. National Institute of Sciences of India, Mathematical Tables, Vol. 1, New Delhi, 1956 (Annotated scanned pages from, plus a review)
Edward Pearce-Crump, Permutation Equivariant Neural Networks for Symmetric Tensors, arXiv:2503.11276 [cs.LG], 2025. See p. 20.

Cf. A005380.
Cf. A219554. Column k=2 of A219727. - Alois P. Heinz, Nov 26 2012
Cf. also A000041, A000070, A000291, A000412, A000465, A000491, A002755, A002756, A002757, A002758, A002759, A277239.
Main diagonal of A054225 if that entry is drawn as a square array. - N. J. A. Sloane, Dec 30 2018

with(numtheory):
b:= proc(n, k) option remember; `if`(n>k, 0, 1) +`if`(isprime(n), 0,
      add(`if`(d>k, 0, b(n/d, d)), d=divisors(n) minus {1, n}))
    end:
a:= n-> b(6^n$2):
seq(a(n), n=0..30);  # Alois P. Heinz, Jun 27 2013

max = 26; se = Series[ Sum[ Log[1 - x^(n-k)*y^k], {n, 1, 2max}, {k, 0, n}], {x, 0, 2max}, {y, 0, 2max}]; coes = CoefficientList[ Series[ Exp[-se], {x, 0, 2max}, {y, 0, 2max}], {x, y}]; a[n_] := coes[[n+1, n+1]]; Table[a[n], {n, 0, max} ](* Jean-François Alcover, Dec 06 2011 *)

a(n) = A054225(2n, n) = A091437(2n).
a(n) ~ Zeta(3)^(19/36) * exp(3*Zeta(3)^(1/3) * n^(2/3) + Pi^2 * n^(1/3) / (6*Zeta(3)^(1/3)) + Zeta'(-1) - Pi^4/(432*Zeta(3))) / (sqrt(3) * (2*Pi)^(3/2) * n^(55/36)). - Vaclav Kotesovec, Jan 30 2016
Formula (25) in the article by Auluck is incorrect. The correct formula is: p(n,n) ~ c^(19/12) * exp(3*c*n^(2/3) + 3*d*n^(1/3) + Zeta'(-1) - 3*d^2/(4*c)) / (sqrt(3) * (2*Pi)^(3/2) * n^(55/36)), where c = Zeta(3)^(1/3), d = Zeta(2)/(3*c). Also formula (24) is incorrect. - Vaclav Kotesovec, Jan 30 2016
From Vaclav Kotesovec, Feb 04 2016: (Start)
The correct formula (24) is p(m,n) ~ c^(7/4)/(2*Pi*sqrt(3)) * exp(3*c*(m*n)^(1/3) + 3*d*(m+n)/(2*(m*n)^(1/3)) - 19*log(m*n)/24 - ((m/n - 2*n/m)*log(m) + (n/m - 2*m/n)*log(n))/36 - (m/n + n/m)*(log(c)/12 + Zeta'(-1) - 1/12 + 3*d^2/(4*c)) + 3*d^2/(4*c) - 3*log(2*Pi)/4 + fi((n/m)^(1/2))),
where m and n are of the same order, c = Zeta(3)^(1/3), d = Zeta(2)/(3*c) and fi(alfa) = Integral_{t=0..infinity} (1/t)*(1/(exp(alfa*t)-1)/(exp(t/alfa)-1) - (alfa/t)/(exp(alfa*t)-1) - ((1/alfa)/t)/(exp(t/alfa)-1) + 1/t^2 + (1/4)/(exp(alfa*t)-1) + (1/4)/(exp(t/alfa)-1) - (alfa/4)/t - ((1/4)/alfa)/t).
If m = n then alfa = 1 and fi(1) = 3*Zeta'(-1) + log(2*Pi)/4 - 1/6.
For the asymptotic formula for fixed m see A054225.
(End)

Corrected using A000491.

Edited by Christian G. Bower, Jan 08 2004

A219555 Number of bipartite partitions of (i,j) with i+j = n into distinct pairs.

Original entry on oeis.org

1, 2, 4, 10, 19, 38, 73, 134, 242, 430, 749, 1282, 2171, 3622, 5979, 9770, 15802, 25334, 40288, 63560, 99554, 154884, 239397, 367800, 561846, 853584, 1290107, 1940304, 2904447, 4328184, 6422164, 9489940, 13967783, 20480534, 29920277, 43557272, 63194864

Offset: 0

Alois P. Heinz, Nov 22 2012

nonn

			a(2) = 4: [(2,0)], [(1,1)], [(1,0),(0,1)], [(0,2)].

Alois P. Heinz, Table of n, a(n) for n = 0..8000 (terms n=101..1000 from Vaclav Kotesovec)

Row sums of A054242.

Cf. A026007, A052812, A005380, A255834, A255836.

b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0, add(
       b(n-i*j, min(n-i*j, i-1))*binomial(i+1, j), j=0..n/i)))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..42);  # Alois P. Heinz, Sep 19 2019

nmax=50; CoefficientList[Series[Product[(1+x^k)^(k+1),{k,1,nmax}],{x,0,nmax}],x] (* Vaclav Kotesovec, Mar 07 2015 *)

a(n) = Sum_{i+j=n} [x^i*y^j] 1/2 * Product_{k,m>=0} (1+x^k*y^m).
G.f.: Product_{k>=1} (1+x^k)^(k+1). - Vaclav Kotesovec, Mar 07 2015
a(n) ~ Zeta(3)^(1/6) * exp(-Pi^4 / (1296*Zeta(3)) + Pi^2 * n^(1/3) / (2^(5/3) * 3^(4/3) * Zeta(3)^(1/3)) + (3/2)^(4/3) * Zeta(3)^(1/3) * n^(2/3)) / (2^(5/4) * 3^(1/3) * sqrt(Pi) * n^(2/3)), where Zeta(3) = A002117. - Vaclav Kotesovec, Mar 07 2015
G.f.: exp(Sum_{k>=1} (-1)^(k+1)*x^k*(2 - x^k)/(k*(1 - x^k)^2)). - Ilya Gutkovskiy, Aug 11 2018

A217093 Number of partitions of n objects of 3 colors.

Original entry on oeis.org

1, 3, 12, 38, 117, 330, 906, 2367, 6027, 14873, 35892, 84657, 196018, 445746, 997962, 2201438, 4792005, 10300950, 21889368, 46012119, 95746284, 197344937, 403121547, 816501180, 1640549317, 3271188702, 6475456896, 12730032791, 24861111315, 48246729411, 93065426256

Offset: 0

Geoffrey Critzer, Sep 26 2012

a(n) is also the number of unlabeled simple graphs with n nodes of 3 colors whose components are complete graphs.
Number of (integer) partitions of n into 3 sorts of part 1, 6 sorts of part 2, 10 sorts of part 3, ..., (k+2)*(k+1)/2 sorts of part k. - Joerg Arndt, Dec 07 2014
In general the g.f. 1 / prod(n>=1, (1-x^k)^m(k) ) gives the number of (integer) partitions where there are m(k) sorts of part k. - Joerg Arndt, Mar 10 2015

			We represent each summand, k, in a partition of n as k identical objects. Then we color each object. We have no regard for the order of the colored objects.
a(2) = 12 because we have: ww; wg; wb; gg; gb; bb; w + w; w + g; w + b; g + g; g + b; b + b, where the 3 colors are white w, gray g, and black b.

Alois P. Heinz, Table of n, a(n) for n = 0..1000
S. Benvenuti, B. Feng, A. Hanany and Y. H. He, Counting BPS operators in gauge theories: Quivers, syzygies and plethystics, arXiv:hep-th/0608050, Aug 2006, p.42.
Carlos A. A. Florentino, Plethystic exponential calculus and permutation polynomials, arXiv:2105.13049 [math.CO], 2021. Mentions this sequence.
Vaclav Kotesovec, Graph - The asymptotic ratio

Cf. A005380, A120844, A253289, A255050, A255052, A255802, A255803, A255835, A255836, A363601, A363602.

Column k=3 of A075196.

with(numtheory):
a:= proc(n) option remember; `if`(n=0, 1, add(add(
      d*binomial(d+2, 2), d=divisors(j))*a(n-j), j=1..n)/n)
    end:
seq(a(n), n=0..30);  # Alois P. Heinz, Sep 26 2012
with(numtheory):
series(exp(add(((1/2)*sigma[3](k) + (3/2)*sigma[2](k) + sigma[1](k))*x^k/k, k = 1..30)), x, 31):
seq(coeftayl(%, x = 0, n), n = 0..30); # Peter Bala, Jan 16 2025

nn=30; p=Product[1/(1- x^i)^Binomial[i+2,2],{i,1,nn}]; CoefficientList[Series[p,{x,0,nn}],x]

from functools import lru_cache
from sympy import divisors
@lru_cache(maxsize=None)
def A217093_aux(n): return sum(d*(d+1)*(d+2)>>1 for d in divisors(n,generator=True))
@lru_cache(maxsize=None)
def A217093(n): return 1 if n == 0 else (A217093_aux(n)+sum(A217093_aux(k)*A217093(n-k) for k in range(1,n)))//n # Chai Wah Wu, Mar 19 2025

G.f.: Product_{i>=1} 1/(1-x^i)^binomial(i+2,2).
EULER transform of 3, 6, 10, 15, ... .
Generally for the number of partitions of k colors the generating function is Product_{i>=1} 1/(1-x^i)^binomial(i+k-1,k-1).
a(n) ~ Pi^(1/8) * exp(1/8 + 3^4 * 5^2 * Zeta(3)^3 / (2*Pi^8) - 31*Zeta(3) / (8*Pi^2) + 5^(1/4) * Pi * n^(1/4) / 6^(3/4) - 3^(13/4) * 5^(5/4) * Zeta(3)^2 * n^(1/4) / (2^(7/4) * Pi^5) + 3^(3/2) * 5^(1/2) * Zeta(3) * n^(1/2) / (2^(1/2) * Pi^2) + 2^(7/4) * Pi * n^(3/4) / (3^(5/4) * 5^(1/4))) / (A^(3/2) * 2^(73/32) * 15^(9/32) * n^(25/32)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant and Zeta(3) = A002117 = 1.202056903... . - Vaclav Kotesovec, Mar 08 2015
G.f.: exp(Sum_{k >= 1} ((1/2)*sigma_3(k) + (3/2)*sigma_2(k) + sigma_1(k))*x^k/k) = 1 + 3*x + 12*x^2 + 38*x^3 + 117*x^4 + .... - Peter Bala, Jan 16 2025

cp's OEIS Frontend

A000219 Number of plane partitions (or planar partitions) of n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Julia

Maple

Mathematica

PARI

PARI

PARI

PARI

Python

SageMath

Formula

Extensions

A052847 G.f.: 1 / Product_{k>=1} (1-x^k)^(k-1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A054225 Triangle read by rows: row n (n>=0) gives the number of partitions of (n,0), (n-1,1), (n-2,2), ..., (0,n) respectively into sums of pairs.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

A002755 Number of bipartite partitions of n white objects and 6 black ones.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

Python

Formula

Extensions

A002774 Number of bipartite partitions of n white objects and n black ones.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A219555 Number of bipartite partitions of (i,j) with i+j = n into distinct pairs.

Views