login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A334550
Triangle read by rows: T(n,k) is the number of binary matrices with n ones, k columns and no zero rows or columns, up to permutations of rows and columns.
4
1, 1, 2, 1, 2, 3, 1, 5, 5, 5, 1, 5, 12, 9, 7, 1, 9, 23, 29, 17, 11, 1, 9, 39, 62, 57, 28, 15, 1, 14, 63, 147, 154, 110, 47, 22, 1, 14, 102, 278, 409, 329, 194, 73, 30, 1, 20, 150, 568, 991, 1023, 664, 335, 114, 42, 1, 20, 221, 1020, 2334, 2844, 2267, 1243, 549, 170, 56
OFFSET
1,3
COMMENTS
T(n,k) is also the number of non-isomorphic set multipartitions (multisets of sets) of weight n with k parts.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..1275 (rows 1..50)
EXAMPLE
Triangle begins:
1;
1, 2;
1, 2, 3;
1, 5, 5, 5;
1, 5, 12, 9, 7;
1, 9, 23, 29, 17, 11;
1, 9, 39, 62, 57, 28, 15;
1, 14, 63, 147, 154, 110, 47, 22;
...
The T(4,3) = 5 matrices are:
[1 0 0] [1 0 0] [1 1 0] [1 1 1] [1 1 0]
[1 0 0] [1 0 0] [1 0 0] [1 0 0] [1 0 1]
[0 1 0] [0 1 1] [0 0 1]
[0 0 1]
The T(4,3) = 5 the set multipartitions are:
{{1,2}, {3}, {4}},
{{1,2}, {3}, {3}},
{{1,2}, {1}, {3}},
{{1,2}, {1}, {1}},
{{1,2}, {1}, {2}}.
PROG
(PARI) \\ See A321609 for definition of M.
T(n, k)={M(k, n, n) - M(k-1, n, n)}
for(n=1, 10, for(k=1, n, print1(T(n, k), ", ")); print)
(PARI) \\ Faster version.
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}
K(q, t, n)={prod(j=1, #q, (1+x^lcm(t, q[j]) + O(x*x^n))^gcd(t, q[j]))}
G(m, n)={my(s=0); forpart(q=m, s+=permcount(q)*exp(sum(t=1, n, (K(q, t, n)-1)/t) + O(x*x^n))); s/m!}
A(n, m=n)={my(p=sum(k=0, m, G(k, n)*y^k)*(1-y)); matrix(n, m, n, k, polcoef(polcoef(p, n, x), k, y))}
{ my(T=A(10)); for(n=1, #T, print(T[n, 1..n])) }
CROSSREFS
Row sums are A049311.
Main diagonal is A000041.
Sequence in context: A078032 A162453 A008313 * A232177 A111377 A014046
KEYWORD
nonn,tabl
AUTHOR
Andrew Howroyd, Jul 03 2020
STATUS
approved