A219561 -id:A219561 - cp's OEIS frontend

A219585 Number A(n,k) of k-partite partitions of {n}^k into distinct k-tuples; square array A(n,k), n>=0, k>=0, read by antidiagonals.

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 5, 5, 2, 1, 1, 15, 40, 17, 2, 1, 1, 52, 457, 364, 46, 3, 1, 1, 203, 6995, 14595, 2897, 123, 4, 1, 1, 877, 136771, 937776, 407287, 21369, 323, 5, 1, 1, 4140, 3299218, 88507276, 107652681, 10200931, 148257, 809, 6, 1

Offset: 0

Alois P. Heinz, Nov 23 2012

A(n,k) is the number of factorizations of m^n into distinct factors where m is a product of k distinct primes. A(2,2) = 5: (2*3)^2 = 36 has 5 factorizations into distinct factors: 36, 3*12, 4*9, 2*18, 2*3*6.

			A(1,3) = 5: [(1,1,1)], [(1,1,0),(0,0,1)], [(1,0,1),(0,1,0)], [(1,0,0),(0,1,0),(0,0,1)], [(0,1,1),(1,0,0)].
A(3,2) = 17: [(3,3)], [(3,0),(0,3)], [(3,2),(0,1)], [(2,3),(1,0)], [(3,1),(0,2)], [(2,2),(1,1)], [(1,3),(2,0)], [(2,1),(1,2)], [(2,1),(1,1),(0,1)], [(3,0),(0,2),(0,1)], [(2,2),(1,0),(0,1)], [(2,1),(0,2),(1,0)], [(1,2),(2,0),(0,1)], [(1,2),(1,1),(1,0)], [(0,3),(2,0),(1,0)], [(2,0),(1,1),(0,2)], [(2,0),(0,2),(1,0),(0,1)].
Square array A(n,k) begins:
  1,  1,   1,      1,          1,            1,         1, ...
  1,  1,   2,      5,         15,           52,       203, ...
  1,  1,   5,     40,        457,         6995,    136771, ...
  1,  2,  17,    364,      14595,       937776,  88507276, ...
  1,  2,  46,   2897,     407287,    107652681,  ...
  1,  3, 123,  21369,   10200931,  10781201973,  ...
  1,  4, 323, 148257,  233051939,  ...
  1,  5, 809, 970246, 4909342744,  ...

Andrew Howroyd, Table of n, a(n) for n = 0..209

Columns k=0..5 give: A000012, A000009, A219554, A219560, A219561, A219565.
Rows n=0..3 give: A000012, A000110, A094574, A319591.
Cf. A188445, A219727 (partitions of {n}^k into k-tuples), A318286.

f[n_, k_] := f[n, k] = 1/2 Product[Sum[O[x[j]]^(n+1), {j, 1, k}]+1+ Product[x[j]^i[j], {j, 1, k}], Evaluate[Sequence @@ Table[{i[j], 0, n}, {j, 1, k}]]];
a[0, ] = a[, 0] = 1; a[n_, k_] := SeriesCoefficient[f[n, k], Sequence @@ Table[{x[j], 0, n}, {j, 1, k}]];
Table[Print[a[n-k, k]]; a[n-k, k], {n, 0, 9}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Dec 11 2013, updated Sep 16 2019 *)

EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
D(p, n, k)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); EulerT(v)[n]^k/prod(i=1, #v, i^v[i]*v[i]!)}
T(n, k)={my(m=n*k+1, q=Vec(exp(intformal(O(x^m) - x^n/(1-x)))/(1+x))); if(n==0, 1, (-1)^m*sum(j=0, m, my(s=0); forpart(p=j, s+=(-1)^#p*D(p, n, k), [1, n]); s*q[#q-j])/2)} \\ Andrew Howroyd, Dec 16 2018

A(n,k) = [(Product_{j=1..k} x_j)^n] 1/2 * Product_{i_1,...,i_k>=0} (1+Product_{j=1..k} x_j^i_j).

A219554 Number of bipartite partitions of (n,n) into distinct pairs.

Original entry on oeis.org

1, 2, 5, 17, 46, 123, 323, 809, 1966, 4660, 10792, 24447, 54344, 118681, 254991, 539852, 1127279, 2323849, 4733680, 9535079, 19005282, 37507802, 73333494, 142112402, 273092320, 520612305, 984944052, 1849920722, 3450476080, 6393203741, 11770416313, 21538246251

Offset: 0

Alois P. Heinz, Nov 22 2012

Number of factorizations of p^n*q^n into distinct factors where p, q are distinct primes.
From Vaclav Kotesovec, Feb 05 2016: (Start)
Formula (15) in the article by S. M. Luthra: "Partitions of bipartite numbers when the summands are unequal", p. 376, is incorrect. The similar error is also in the article by F. C. Auluck: "On partitions of bipartite numbers" (see A002774).
The correct formula (15) is q(m, n) ~ c/(2*sqrt(3)*Pi) * exp(3*c*(m*n)^(1/3) + 3*d*(m^(2/3)/n^(1/3) + n^(2/3)/m^(1/3)) - 3*log(2)/4 + (m/n + n/m)*log(2)/12 + 3*d^2/c - 3*d^2*(m/n + n/m)/c - 2*log(m*n)/3), where m and n are of the same order, c = (3/4*Zeta(3))^(1/3), d = Zeta(2)/(12*c).
If m = n then q(m,n) = a(n).
For the asymptotic formula for fixed m see A054242.
(End)

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

Alois P. Heinz and Vaclav Kotesovec, Table of n, a(n) for n = 0..400 (terms 0..100 from Alois P. Heinz)
S. M. Luthra, Partitions of bipartite numbers when the summands are unequal, Proceedings of the Indian National Science Academy, vol. 23, 1957, issue 5A, p. 370-376. [broken link]

Column k=2 of A219585.
Cf. A002774, A219560, A219561, A219565.
Cf. A054242, A000009, A036469, A268345, A268346, A268347, A268348.

(* This program is not convenient for a large number of terms *)
a[n_] := If[n == 0, 1, (1/2) Coefficient[Product[O[x]^(n+1) + O[y]^(n+1) + (1 + x^i y^j ), {i, 0, n}, {j, 0, n}] // Normal, (x y)^n]];
a /@ Range[0, 31] (* Jean-François Alcover, Jun 26 2013, updated Sep 16 2019 *)
nmax = 20; p = 1; Do[Do[p = Expand[p*(1 + x^i*y^j)]; If[i*j != 0, p = Select[p, (Exponent[#, x] <= nmax) && (Exponent[#, y] <= nmax) &]], {i, 0, nmax}], {j, 0, nmax}]; p = Select[p, Exponent[#, x] == Exponent[#, y] &]; Flatten[{1, Table[Coefficient[p, x^n*y^n]/2, {n, 1, nmax}]}] (* Vaclav Kotesovec, Jan 15 2016 *)

a(n) = [x^n*y^n] 1/2 * Product_{i,j>=0} (1+x^i*y^j).
a(n) = A054242(2*n,n) = A201377(n,n).
a(n) ~ Zeta(3)^(1/3) * exp(3^(4/3) * Zeta(3)^(1/3) * n^(2/3) / 2^(2/3) + Pi^2 * n^(1/3) / (6^(4/3) * Zeta(3)^(1/3)) - Pi^4/(1296*Zeta(3))) / (2^(9/4) * 3^(1/6) * Pi * n^(4/3)). - Vaclav Kotesovec, Jan 31 2016

A280486 G.f.: Product_{i>=1, j>=1, k>=1, l>=1} (1 + x^(ijk*l)).

Original entry on oeis.org

1, 1, 4, 8, 20, 36, 86, 150, 314, 564, 1088, 1902, 3557, 6085, 10902, 18506, 32124, 53584, 91133, 149749, 249315, 405121, 662582, 1063152, 1714580, 2719842, 4327302, 6797316, 10686005, 16622003, 25861855, 39866017, 61422891, 93910783, 143406552, 217537696

Offset: 0

Vaclav Kotesovec, Jan 04 2017

nonn

Vaclav Kotesovec, Table of n, a(n) for n = 0..10000
Lida Ahmadi, Ricardo Gómez Aíza, and Mark Daniel Ward, A unified treatment of families of partition functions, La Matematica (2024). Preprint available as arXiv:2303.02240 [math.CO], 2023.

Cf. A000009, A007426, A107742, A280473, A280487, A219561.

nmax = 50; CoefficientList[Series[Product[(1+x^(i*j*k*l)), {i, 1, nmax}, {j, 1, nmax/i}, {k, 1, nmax/i/j}, {l, 1, nmax/i/j/k}], {x, 0, nmax}], x]
nmax = 50; tau4 = Table[DivisorSum[n, DivisorSigma[0, n/#] * DivisorSigma[0, #] &], {n, 1, nmax}]; s = 1 + x; Do[s *= Sum[Binomial[tau4[[k]], j]*x^(j*k), {j, 0, nmax/k}]; s = Expand[s]; s = Take[s, Min[nmax + 1, Exponent[s, x] + 1, Length[s]]];, {k, 2, nmax}]; Take[CoefficientList[s, x], nmax + 1] (* Vaclav Kotesovec, Sep 08 2018 *)

G.f.: Product_{k>=1} (1 + x^k)^tau_4(k), where tau_4() = A007426. - Ilya Gutkovskiy, May 22 2018

A219560 Number of tripartite partitions of (n,n,n) into distinct triples.

Original entry on oeis.org

1, 5, 40, 364, 2897, 21369, 148257, 970246, 6032341, 35850410, 204646488, 1126463948, 5999145787, 30999381232, 155798366059, 763194776551, 3650648583934, 17079277343463, 78262895082681, 351708874155894, 1551843168854346

Offset: 0

Alois P. Heinz, Nov 23 2012

Number of factorizations of (p*q*r)^n into distinct factors where p, q, r are distinct primes.

			a(0) = 1: [].
a(1) = 5: [(1,1,1)], [(1,1,0),(0,0,1)], [(1,0,1),(0,1,0)], [(0,1,1),(1,0,0)], [(0,0,1),(0,1,0),(1,0,0)].

Column k=3 of A219585.

Cf. A002774, A219554, A219561, A219565, A219678.

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-1)), d=divisors(n) minus {1, n}))
    end:
a:= n-> b(30^n$2):
seq(a(n), n=0..10);  # Alois P. Heinz, May 26 2013

b[n_, k_] := b[n, k] = If[n > k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d > k, 0, b[n/d, d - 1]], {d, Divisors[n][[2 ;; -2]]}]]; a[0] = 1; a[n_] := b[30^n, 30^n];  Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 0, 20}] (* Jean-François Alcover, Jan 15 2016, after Alois P. Heinz *)

a(n) = [(x*y*z)^n] 1/2 * Product_{i,j,k>=0} (1+x^i*y^j*z^k).

a(16) from Alois P. Heinz, May 26 2013
a(17) from Alois P. Heinz, Sep 24 2014
More terms from Jean-François Alcover, Jan 15 2016

A219565 Number of 5-partite partitions of (n,n,n,n,n) into distinct quintuples.

Original entry on oeis.org

1, 52, 6995, 937776, 107652681, 10781201973, 958919976957, 76861542428397, 5620227129073491, 378709513816248475, 23713852762539359688, 1389561695379881634055, 76647024053735036288641, 3999799865715906390697377, 198328846122797866982616805, 9379277765981012067789260214

Offset: 0

Alois P. Heinz, Nov 23 2012

nonn

Number of factorizations of (p*q*r*s*t)^n into distinct factors where p, q, r, s, t are distinct primes.

Column k=5 of A219585.

Cf. A002774, A219554, A219560, A219561.

a[n_] := a[n] = If[n == 0, 1, (1/2) Coefficient[Product[O[v]^(n+1) + O[w]^(n+1) + O[x]^(n+1) + O[y]^(n+1) + O[z]^(n+1) + (1 + v^i w^j x^k y^l z^m), {i, 0, n}, {j, 0, n}, {k, 0, n}, {l, 0, n}, {m, 0, n}] // Normal, (v w x y z)^n]];
Table[Print[n, " ", a[n]]; a[n], {n, 0, 7}] (* Jean-François Alcover, Sep 24 2019 *)

a(n) = [(v*w*x*y*z)^n] 1/2 * Product_{h,i,j,k,m>=0} (1+v^h*w^i*x^j*y^k*z^m).

a(6) from Alois P. Heinz, Sep 25 2014

a(7)-a(15) from Andrew Howroyd, Dec 16 2018

cp's OEIS Frontend

A219585 Number A(n,k) of k-partite partitions of {n}^k into distinct k-tuples; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A219554 Number of bipartite partitions of (n,n) into distinct pairs.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A280486 G.f.: Product_{i>=1, j>=1, k>=1, l>=1} (1 + x^(i*j*k*l)).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

A219560 Number of tripartite partitions of (n,n,n) into distinct triples.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A219565 Number of 5-partite partitions of (n,n,n,n,n) into distinct quintuples.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

Formula

Extensions

A280486 G.f.: Product_{i>=1, j>=1, k>=1, l>=1} (1 + x^(ijk*l)).