A188392 -id:A188392 - cp's OEIS frontend

A318360 Number of set multipartitions (multisets of sets) of a multiset whose multiplicities are the prime indices of n.

1, 1, 1, 2, 1, 2, 1, 5, 3, 2, 1, 6, 1, 2, 3, 15, 1, 9, 1, 6, 3, 2, 1, 21, 4, 2, 16, 6, 1, 10, 1, 52, 3, 2, 4, 35, 1, 2, 3, 22, 1, 10, 1, 6, 19, 2, 1, 83, 5, 13, 3, 6, 1, 66, 4, 22, 3, 2, 1, 41, 1, 2, 20, 203, 4, 10, 1, 6, 3, 14, 1, 153, 1, 2, 26, 6, 5, 10, 1

Offset: 1

Gus Wiseman, Aug 24 2018

nonn

			The a(12) = 6 set multipartitions of {1,1,2,3}:
  {{1},{1,2,3}}
  {{1,2},{1,3}}
  {{1},{1},{2,3}}
  {{1},{2},{1,3}}
  {{1},{3},{1,2}}
  {{1},{1},{2},{3}}

Andrew Howroyd, Table of n, a(n) for n = 1..1000

Cf. A001055, A007716, A049311, A116540, A181821, A188392, A255906.

Cf. A318283, A318284, A318286, A318361, A318362, A318369.

nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
sqfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
Table[Length[sqfacs[Times@@Prime/@nrmptn[n]]],{n,80}]

permcount(v) = {my(m=1,s=0,k=0,t); for(i=1,#v,t=v[i]; k=if(i>1&&t==v[i-1],k+1,1); m*=t*k;s+=t); s!/m}
sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i,2], j, primepi(f[i,1]))))}
count(sig)={my(n=vecsum(sig), s=0); forpart(p=n, my(q=prod(i=1, #p, 1 + x^p[i] + O(x*x^n))); s+=prod(i=1, #sig, polcoef(q,sig[i]))*permcount(p)); s/n!}
a(n)={if(n==1, 1, my(s=sig(n)); if(#s<=2, if(#s==1, 1, min(s[1],s[2])+1), count(sig(n))))} \\ Andrew Howroyd, Dec 10 2018

a(n) = A050320(A181821(n)).
From Andrew Howroyd, Dec 10 2018:(Start)
a(p) = 1 for prime(p).
a(prime(i)*prime(j)) = min(i,j) + 1.
a(prime(n)^k) = A188392(n,k). (End)

A262809 Number A(n,k) of lattice paths from {n}^k to {0}^k using steps that decrement one or more components by one; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 13, 13, 1, 1, 1, 75, 409, 63, 1, 1, 1, 541, 23917, 16081, 321, 1, 1, 1, 4683, 2244361, 10681263, 699121, 1683, 1, 1, 1, 47293, 308682013, 14638956721, 5552351121, 32193253, 8989, 1, 1, 1, 545835, 58514835289, 35941784497263, 117029959485121, 3147728203035, 1538743249, 48639, 1, 1

Offset: 0

Alois P. Heinz, Oct 02 2015

Also, A(n,k) is the number of alignments for k sequences of length n each (Slowinski 1998).
Row r > 0 is asymptotic to sqrt(r*Pi) * (r^(r-1)/(r-1)!)^n * n^(r*n+1/2) / (2^(r/2) * exp(r*n) * (log(2))^(r*n+1)), or equivalently to sqrt(r) * (r^(r-1)/(r-1)!)^n * (n!)^r / (2^r * (Pi*n)^((r-1)/2) * (log(2))^(r*n+1)). - Vaclav Kotesovec, Mar 23 2016
From Vaclav Kotesovec, Mar 23 2016: (Start)
Column k > 0 is asymptotic to sqrt(c(k)) * d(k)^n / (Pi*n)^((k-1)/2), where c(k) and d(k) are roots of polynomial equations of degree k, independent on n.
---------------------------------------------------
k               d(k)
---------------------------------------------------
2              5.8284271247461900976033774484193...
3             56.9476283720414911685286267804411...
4            780.2794068067951456595241495989622...
5          13755.2719024115081712083954421541320...
6         296476.9162644200814909862281498491264...
7        7553550.6198338218721069097516499501996...
8      222082591.6017202421029000117685530884167...
9     7400694480.0494436216324852038000444393262...
10  275651917450.6709238286995776605620357737005...
---------------------------------------------------
d(k) is a root of polynomial:
---------------------------------------------------
k=2, 1 - 6*d + d^2
k=3, -1 + 3*d - 57*d^2 + d^3
k=4, 1 - 12*d - 218*d^2 - 780*d^3 + d^4
k=5, -1 + 5*d - 1260*d^2 - 3740*d^3 - 13755*d^4 + d^5
k=6, 1 - 18*d - 5397*d^2 - 123696*d^3 + 321303*d^4 - 296478*d^5 + d^6
k=7, -1 + 7*d - 24031*d^2 - 374521*d^3 - 24850385*d^4 + 17978709*d^5 - 7553553*d^6 + d^7
k=8, 1 - 24*d - 102692*d^2 - 9298344*d^3 + 536208070*d^4 - 7106080680*d^5 - 1688209700*d^6 - 222082584*d^7 + d^8
(End)
d(k) = (2^(1/k) - 1)^(-k). - David Bevan, Apr 07 2022
d(k) is asymptotic to (k/log(2))^k/sqrt(2). - David Bevan, Apr 07 2022
A(n,k) is the number of binary matrices with k columns and any number of nonzero rows with n ones in every column. - Andrew Howroyd, Jan 23 2020

			A(2,2) = 13: [(2,2),(1,2),(0,2),(0,1),(0,0)], [(2,2),(1,2),(0,1),(0,0)], [(2,2),(1,2),(1,1),(0,1),(0,0)], [(2,2),(1,2),(1,1),(0,0)], [(2,2),(1,2),(1,1),(1,0),(0,0)], [(2,2),(2,1),(1,1),(0,1),(0,0)], [(2,2),(2,1),(1,1),(0,0)], [(2,2),(2,1),(1,1),(1,0),(0,0)], [(2,2),(2,1),(2,0),(0,1),(0,0)], [(2,2),(2,1),(1,0),(0,0)], [(2,2),(1,1),(0,1),(0,0)], [(2,2),(1,1),(0,0)], [(2,2),(1,1),(1,0),(0,0)].
Square array A(n,k) begins:
  1, 1,    1,        1,             1,                   1, ...
  1, 1,    3,       13,            75,                 541, ...
  1, 1,   13,      409,         23917,             2244361, ...
  1, 1,   63,    16081,      10681263,         14638956721, ...
  1, 1,  321,   699121,    5552351121,     117029959485121, ...
  1, 1, 1683, 32193253, 3147728203035, 1050740615666453461, ...

Alois P. Heinz, Antidiagonals n = 0..48, flattened
J. B. Slowinski, The Number of Multiple Alignments, Molecular Phylogenetics and Evolution 10:2 (1998), 264-266. doi:10.1006/mpev.1998.0522

Columns: A000012 (k=0 and k=1), A001850 (k=2), A126086 (k=3), A263064 (k=4), A263065 (k=5), A263066 (k=6), A263067 (k=7), A263068 (k=8), A263069 (k=9), A263070 (k=10).
Rows: A000012 (n=0), A000670 (n=1), A055203 (n=2), A062208 (n=3), A062205 (n=4), A263061 (n=5), A263062 (n=6), A062204 (n=7), A263063 (n=8), A263071 (n=9), A263072 (n=10).
Main diagonal: A262810.
Cf. A210472, A225094, A227578, A227655, A229142, A229345, A263159, A316674.
Cf. A188392, A330942, A331461, A331637.

A:= (n, k)-> add(add((-1)^i*binomial(j, i)*
     binomial(j-i, n)^k, i=0..j), j=0..k*n):
seq(seq(A(n, d-n), n=0..d), d=0..10);

A[, 0] =  1; A[n, k_] := Sum[Sum[(-1)^i*Binomial[j, i]*Binomial[j - i, n]^k, {i, 0, j}], {j, 0, k*n}];
Table[Table[A[n, d - n], {n, 0, d}], {d, 0, 10}] // Flatten (* Jean-François Alcover, Jul 22 2016, after Alois P. Heinz *)

T(n,k) = {my(m=n*k); sum(j=0, m, binomial(j,n)^k*sum(i=j, m, (-1)^(i-j)*binomial(i, j)))} \\ Andrew Howroyd, Jan 23 2020

A(n,k) = Sum_{j=0..k*n} Sum_{i=0..j} (-1)^i*C(j,i)*C(j-i,n)^k.
A(n,k) = Sum_{i >= 0} binomial(i,n)^k/2^(i+1). - Peter Bala, Jan 30 2018
A(n,k) = Sum_{j=0..n*k} binomial(j,n)^k * Sum_{i=j..n*k} (-1)^(i-j) * binomial(i,j). - Andrew Howroyd, Jan 23 2020

A330942 Array read by antidiagonals: A(n,k) is the number of binary matrices with k columns and any number of nonzero rows with n ones in every column and columns in nonincreasing lexicographic order.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 4, 7, 1, 1, 1, 8, 75, 32, 1, 1, 1, 16, 1105, 2712, 161, 1, 1, 1, 32, 20821, 449102, 116681, 842, 1, 1, 1, 64, 478439, 122886128, 231522891, 5366384, 4495, 1, 1, 1, 128, 12977815, 50225389432, 975712562347, 131163390878, 256461703, 24320, 1, 1

Offset: 0

Andrew Howroyd, Jan 13 2020

The condition that the columns be in nonincreasing order is equivalent to considering nonequivalent matrices up to permutation of columns.

A(n,k) is the number of labeled n-uniform hypergraphs with multiple edges allowed and with k edges and no isolated vertices. When n=2 these objects are multigraphs.

			Array begins:
============================================================
n\k | 0 1    2         3              4                5
----+-------------------------------------------------------
  0 | 1 1    1         1              1                1 ...
  1 | 1 1    2         4              8               16 ...
  2 | 1 1    7        75           1105            20821 ...
  3 | 1 1   32      2712         449102        122886128 ...
  4 | 1 1  161    116681      231522891     975712562347 ...
  5 | 1 1  842   5366384   131163390878 8756434117294432 ...
  6 | 1 1 4495 256461703 78650129124911 ...
  ...
The A(2,2) = 7 matrices are:
   [1 0]  [1 0]  [1 0]  [1 1]  [1 0]  [1 0]  [1 1]
   [1 0]  [0 1]  [0 1]  [1 0]  [1 1]  [0 1]  [1 1]
   [0 1]  [1 0]  [0 1]  [0 1]  [0 1]  [1 1]
   [0 1]  [0 1]  [1 0]

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

Rows n=1..3 are A000012, A121316, A136246.
Columns k=0..3 are A000012, A000012, A226994, A137220.
The version with nonnegative integer entries is A331315.
Other variations considering distinct rows and columns and equivalence under different combinations of permutations of rows and columns are:
All solutions: A262809 (all), A331567 (distinct rows).
Up to row permutation: A188392, A188445, A331126, A331039.
Up to column permutation: this sequence, A331571, A331277, A331569.
Nonisomorphic: A331461, A331510, A331508, A331509.
Cf. A331638.

T[n_, k_] := With[{m = n k}, Sum[Binomial[Binomial[j, n] + k - 1, k] Sum[ (-1)^(i - j) Binomial[i, j], {i, j, m}], {j, 0, m}]];
Table[T[n - k, k], {n, 0, 9}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Apr 10 2020, from PARI *)

T(n, k)={my(m=n*k); sum(j=0, m, binomial(binomial(j, n)+k-1, k)*sum(i=j, m, (-1)^(i-j)*binomial(i, j)))}

A(n,k) = Sum_{j=0..n*k} binomial(binomial(j,n)+k-1, k) * (Sum_{i=j..n*k} (-1)^(i-j)*binomial(i,j)).
A(n, k) = Sum_{j=0..k} abs(Stirling1(k, j))*A262809(n, j)/k!.
A(n, k) = Sum_{j=0..k} binomial(k-1, k-j)*A331277(n, j).
A331638(n) = Sum_{d|n} A(n/d, d).

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.

Original entry on oeis.org

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

A188445 T(n,k) is the number of (n*k) X k binary arrays with nonzero rows in decreasing order and n ones in every column.

Original entry on oeis.org

1, 2, 0, 5, 1, 0, 15, 8, 0, 0, 52, 80, 5, 0, 0, 203, 1088, 205, 1, 0, 0, 877, 19232, 11301, 278, 0, 0, 0, 4140, 424400, 904580, 67198, 205, 0, 0, 0, 21147, 11361786, 101173251, 24537905, 250735, 80, 0, 0, 0, 115975, 361058000, 15207243828, 13744869502

Offset: 1

R. H. Hardin, Mar 31 2011

			Array begins:
============================================================================
n\k| 1 2 3   4       5          6             7              8             9
---+------------------------------------------------------------------------
1  | 1 2 5  15      52        203           877           4140         21147
2  | 0 1 8  80    1088      19232        424400       11361786     361058000
3  | 0 0 5 205   11301     904580     101173251    15207243828 2975725761202
4  | 0 0 1 278   67198   24537905   13744869502 11385203921707 ...
5  | 0 0 0 205  250735  425677958 1184910460297 ...
6  | 0 0 0  80  621348 5064948309 ...
7  | 0 0 0  15 1058139 ...
8  | 0 0 0   1 ...
...
Some solutions for 16 X 4:
  1 1 1 0    1 1 1 1    1 1 1 1    1 1 1 0    1 1 1 1
  1 0 1 1    1 1 0 1    1 1 0 0    1 0 1 1    1 1 0 0
  1 0 1 0    1 0 1 1    1 0 1 1    1 0 0 1    1 0 1 1
  1 0 0 1    1 0 0 0    1 0 0 0    1 0 0 0    1 0 0 0
  0 1 1 1    0 1 1 0    0 1 1 1    0 1 1 0    0 1 1 1
  0 1 0 1    0 1 0 0    0 1 0 0    0 1 0 1    0 1 0 0
  0 1 0 0    0 0 1 1    0 0 1 1    0 1 0 0    0 0 1 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 1 1    0 0 0 1
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0

Andrew Howroyd, Table of n, a(n) for n = 1..181 (terms 1..69 from R. H. Hardin)

Rows 1..5 are A000110, A002718, A060486, A188446, A188447.
Columns 5..6 are A331127, A331129.
Column sums are A319190.
Cf. A059443, A060487, A188392, A219585, A318361, A330942, A330964, A331039.

WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, (-1)^(n-1)/n))))-1, -#v)}
D(p, n, k)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); WeighT(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) = 0 for n > 2^(k-1). - Andrew Howroyd, Jan 24 2020

A020554 Number of multigraphs on n labeled edges (without loops).

Original entry on oeis.org

1, 1, 3, 16, 139, 1750, 29388, 624889, 16255738, 504717929, 18353177160, 769917601384, 36803030137203, 1984024379014193, 119571835094300406, 7995677265437541258, 589356399302126773920, 47609742627231823142029, 4193665147256300117666879

Offset: 0

Gilbert Labelle (gilbert(AT)lacim.uqam.ca) and Simon Plouffe

Or, number of bicoverings of an n-set.
Or, number of 2-covers of [1,...,n].
Also the number of set multipartitions (multisets of sets) of {1, 1, 2, 2, 3, 3, ..., n, n}. - Gus Wiseman, Jul 18 2018

			From _Gus Wiseman_, Jul 18 2018: (Start)
The a(3) = 16 set multipartitions of {1, 1, 2, 2, 3, 3}:
  (123)(123)
  (1)(23)(123) (2)(13)(123) (3)(12)(123) (12)(13)(23)
  (1)(1)(23)(23) (1)(2)(3)(123) (1)(2)(13)(23) (1)(3)(12)(23) (2)(2)(13)(13) (2)(3)(12)(13) (3)(3)(12)(12)
  (1)(1)(2)(3)(23) (1)(2)(2)(3)(13) (1)(2)(3)(3)(12)
  (1)(1)(2)(2)(3)(3)
(End)

G. Paquin, Dénombrement de multigraphes enrichis, Mémoire, Math. Dept., Univ. Québec à Montréal, 2004.

Alois P. Heinz, Table of n, a(n) for n = 0..100
Peter Cameron, Thomas Prellberg, and Dudley Stark, Asymptotic enumeration of 2-covers and line graphs, Discrete Math. 310 (2010), no. 2, 230-240 (see s_n).
L. Comtet, Birecouvrements et birevêtements d'un ensemble fini, Studia Sci. Math. Hungar 3 (1968): 137-152. [Annotated scanned copy. Warning: the table of v(n,k) has errors.]
G. Labelle, Counting enriched multigraphs according to the number of their edges (or arcs), Discrete Math., 217 (2000), 237-248.
G. Paquin, Dénombrement de multigraphes enrichis, Mémoire, Math. Dept., Univ. Québec à Montréal, 2004. [Cached copy, with permission]

Row 2 of A188392.

Cf. A002718, A007716, A020555, A050535, A094574, A316974.

Ceiling[ CoefficientList[ Series[ Exp[ -1 + (Exp[ z ] - 1)/2 ]Sum[ Exp[ s(s - 1)z/2 ]/s!, {s, 0, 21} ], {z, 0, 9} ], z ] Table[ n!, {n, 0, 9} ] ] (* Mitch Harris, May 01 2004 *)
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
Table[Length[Select[mps[Ceiling[Range[1/2,n,1/2]]],And@@UnsameQ@@@#&]],{n,5}] (* Gus Wiseman, Jul 18 2018 *)

E.g.f.: exp(-3/2+exp(x)/2) * Sum_{n>=0} exp(binomial(n, 2)*x)/n! [Comtet]. - Vladeta Jovovic, Apr 27 2004
E.g.f. (an equivalent version in Maple format): G:=exp(-1+(exp(z)-1)/2)*sum(exp(s*(s-1)*z/2)/s!, s=0..infinity);
E.g.f.: exp((exp(x)-1)/2)*Sum_{n>=0} A020556(n)*(x/2)^n/n!. - Vladeta Jovovic, May 02 2004
Stirling_2 transform of A014500.
The e.g.f.'s of A020554 (S(x)) and A014500 (U(x)) are related by S(x) = U(e^x-1).

A331461 Array read by antidiagonals: A(n,k) is the number of nonequivalent binary matrices with k columns and any number of nonzero rows with n ones in every column up to permutation of rows and columns.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 3, 1, 1, 1, 5, 8, 4, 1, 1, 1, 7, 23, 16, 5, 1, 1, 1, 11, 66, 93, 30, 6, 1, 1, 1, 15, 212, 652, 332, 50, 7, 1, 1, 1, 22, 686, 6369, 6414, 1062, 80, 8, 1, 1, 1, 30, 2389, 79568, 226041, 56712, 3117, 120, 9, 1, 1, 1, 42, 8682, 1256425, 12848128, 7295812, 441881, 8399, 175, 10, 1, 1

Offset: 0

Andrew Howroyd, Jan 18 2020

A(n,k) is the number of non-isomorphic set multipartitions (multiset of sets) with k parts each part has size n.

			Array begins:
===========================================================
n\k | 0 1 2   3    4       5          6              7
----+-----------------------------------------------------
  0 | 1 1 1   1    1       1          1              1 ...
  1 | 1 1 2   3    5       7         11             15 ...
  2 | 1 1 3   8   23      66        212            686 ...
  3 | 1 1 4  16   93     652       6369          79568 ...
  4 | 1 1 5  30  332    6414     226041       12848128 ...
  5 | 1 1 6  50 1062   56712    7295812     1817321457 ...
  6 | 1 1 7  80 3117  441881  195486906   200065951078 ...
  7 | 1 1 8 120 8399 3006771 4298181107 17131523059493 ...
  ...
The A(2,3) = 8 matrices are:
  [1 0 0]  [1 1 0]  [1 1 1]  [1 1 0]  [1 1 0]  [1 1 1]  [1 1 0]  [1 1 1]
  [1 0 0]  [1 0 0]  [1 0 0]  [1 1 0]  [1 0 1]  [1 1 0]  [1 0 1]  [1 1 1]
  [0 1 0]  [0 1 0]  [0 1 0]  [0 0 1]  [0 1 0]  [0 0 1]  [0 1 1]
  [0 1 0]  [0 0 1]  [0 0 1]  [0 0 1]  [0 0 1]
  [0 0 1]  [0 0 1]
  [0 0 1]

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

Rows n=0..6 are A000012, A000041, A050535, A050913, A058783, A058784, A058785.
Columns k=0..4 are A000012, A000012, A000027(n+1), A002624, A331720.
Cf. A188392, A262809, A304942, A306018, A330942, A331485, A331508, A331509, A331510.

\\ See A304942 for Blocks
T(n,k)={Blocks(k, n*k, n)}
{ for(n=0, 7, for(k=0, 6, print1(T(n,k), ", ")); print) }

A306018(n) = Sum_{d|n} A(n/d, d).

A331126 Array read by antidiagonals: A(n,k) is the number of T_0 n-regular set multipartitions (multisets of sets) on a k-set.

Original entry on oeis.org

1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 2, 1, 1, 0, 1, 9, 3, 1, 1, 0, 1, 70, 29, 4, 1, 1, 0, 1, 794, 666, 68, 5, 1, 1, 0, 1, 12055, 28344, 3642, 134, 6, 1, 1, 0, 1, 233238, 1935054, 469368, 14951, 237, 7, 1, 1, 0, 1, 5556725, 193926796, 119843417, 5289611, 50985, 388, 8, 1, 1

Offset: 0

Andrew Howroyd, Jan 10 2020

An n-regular set multipartition is a finite multiset of nonempty sets in which each element appears in n blocks.
A set multipartition is T_0 if for every two distinct elements there exists a block containing one but not the other element.
A(n,k) is the number of binary matrices with k distinct columns and any number of nonzero rows with n ones in every column and rows in nonincreasing lexicographic order.

			Array begins:
====================================================================
n\k | 0 1 2   3      4         5             6                 7
----+---------------------------------------------------------------
  0 | 1 1 0   0      0         0             0                 0 ...
  1 | 1 1 1   1      1         1             1                 1 ...
  2 | 1 1 2   9     70       794         12055            233238 ...
  3 | 1 1 3  29    666     28344       1935054         193926796 ...
  4 | 1 1 4  68   3642    469368     119843417       53059346010 ...
  5 | 1 1 5 134  14951   5289611    4681749424     8639480647842 ...
  6 | 1 1 6 237  50985  46241343  134332244907   989821806791367 ...
  7 | 1 1 7 388 151901 333750928 3032595328876 85801167516707734 ...
     ...
The A(2,2) = 2 matrices are:
   [1 1]   [1 0]
   [1 0]   [1 0]
   [0 1]   [0 1]
           [0 1]
The corresponding set multipartitions are:
    {{1,2}, {1}, {2}},
    {{1}, {1}, {2}, {2}}.

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

Rows n=1..3 are A000012, A014500, A331389.
Columns k=0..3 are A000012, A000012, A001477, A331390.
Cf. A188392, A188445, A330942, A331039, A331160, A331161, A331391.

WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, (-1)^(n-1)/n))))-1, -#v)}
D(p, n, k)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); binomial(WeighT(v)[n], k)*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, k<=1, sum(j=0, m, my(s=0); forpart(p=j, s+=D(p, n, k), [1, n]); s*q[#q-j]))}

A(n, k) = Sum_{j=1..k} Stirling1(k, j)*A188392(n, j) for n, k >= 1.

A331391(n) = Sum_{d|n} A(n/d, d).

A165434 Number of tri-coverings of a set.

Original entry on oeis.org

1, 1, 4, 39, 862, 35775, 2406208, 238773109, 32867762616, 6009498859909, 1412846181645855, 416415343791239162, 150747204270574506888, 65905473934553360340713, 34305461329980340135062217, 21003556204331356488142290707, 14967168378184553824642693791437

Offset: 0

Doron Zeilberger, Sep 18 2009

nonn

			For n=2, a(2)=4, since if you have two sets of identical triples the A-brothers and the B-sisters, and you want to arrange them into a multiset of nonempty sets, where no one is allowed to cohabitate with his or her sibling, the following are possible 1.{{AB},{AB},{AB}} 2.{{AB},{AB},{A},{B}} 3.{{AB},{A},{A},{B},{B}} 4.{{A},{A},{A},{B},{B},{B}}.

Andrew Howroyd, Table of n, a(n) for n = 0..100
E. A. Bender, Partitions of multisets, Discrete Mathematics 9 (1974) 301-312.
J. S. Devitt and D. M. Jackson, The enumeration of covers of a finite set, J. London Math. Soc.(2) 25 (1982), 1-6.
Doron Zeilberger, In How Many Ways Can You Reassemble Several Russian Dolls?, has links to more terms and related sequences
Doron Zeilberger, In How Many Ways Can You Reassemble Several Russian Dolls?, arXiv:0909.3453 [math.CO], 2009.
Doron Zeilberger, BABUSHKAS; Local copy

Row 3 of A188392.

Cf. A000110 (unicoverings), A020554 (bicoverings).

Do SeqBrn(3,n); in the Maple package BABUSHKAS (see links) where n+1 is the number of desired terms.

Edited by Charles R Greathouse IV, Oct 28 2009

A165435 Number of quad-coverings of a set.

Original entry on oeis.org

1, 1, 5, 81, 4079, 507549, 127126912, 55643064708, 38715666455777, 40095856807088486, 58901884724160709571, 118283825763578012358080, 315297447553856723754810322, 1089117884184737631305395668913, 4778248300569771278502378511288513

Offset: 0

Doron Zeilberger, Sep 18 2009

nonn

Analogous to A165434. See comments and references there.

Row 4 of A188392.

Cf. A165434.

cp's OEIS Frontend

A318360 Number of set multipartitions (multisets of sets) of a multiset whose multiplicities are the prime indices of n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A262809 Number A(n,k) of lattice paths from {n}^k to {0}^k using steps that decrement one or more components by one; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A330942 Array read by antidiagonals: A(n,k) is the number of binary matrices with k columns and any number of nonzero rows with n ones in every column and columns in nonincreasing lexicographic order.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

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

A188445 T(n,k) is the number of (n*k) X k binary arrays with nonzero rows in decreasing order and n ones in every column.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

PARI

Formula

A020554 Number of multigraphs on n labeled edges (without loops).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Formula

A331461 Array read by antidiagonals: A(n,k) is the number of nonequivalent binary matrices with k columns and any number of nonzero rows with n ones in every column up to permutation of rows and columns.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula