A219678 -id:A219678 - cp's OEIS frontend

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

1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 5, 9, 3, 1, 1, 15, 66, 31, 5, 1, 1, 52, 712, 686, 109, 7, 1, 1, 203, 10457, 27036, 6721, 339, 11, 1, 1, 877, 198091, 1688360, 911838, 58616, 1043, 15, 1, 1, 4140, 4659138, 154703688, 231575143, 26908756, 476781, 2998, 22, 1

Offset: 0

Alois P. Heinz, Nov 26 2012

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

A(n,k) is the number of (n*k) X k matrices with nonnegative integer entries and column sums n up to permutation of rows. - Andrew Howroyd, Dec 10 2018

			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(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)].
Square array A(n,k) begins:
  1,   1,    1,      1,        1,         1,         1,       1, ...
  1,   1,    2,      5,       15,        52,       203,     877, ...
  1,   2,    9,     66,      712,     10457,    198091, 4659138, ...
  1,   3,   31,    686,    27036,   1688360, 154703688, ...
  1,   5,  109,   6721,   911838, 231575143, ...
  1,   7,  339,  58616, 26908756, ...
  1,  11, 1043, 476781, ...
  1,  15, 2998, ...

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

Columns k=0..3 give: A000012, A000041, A002774, A219678.
Rows n=0..4 give: A000012, A000110, A020555, A322487, A358781.
Main diagonal gives A322488.
Cf. A188392, A219585 (partitions of {n}^k into distinct k-tuples), A256384, A318284, A318951.

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, q=Vec(exp(O(x*x^m) + intformal((x^n-1)/(1-x)))/(1-x))); if(n==0, 1, sum(j=0, m, my(s=0); forpart(p=j, s+=D(p,n,k), [1,n]); s*q[#q-j]))} \\ Andrew Howroyd, Dec 11 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

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

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

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions