OFFSET
0,1
COMMENTS
Also the number of inequivalent square {0,1} matrices with distinct rows and exactly n 1's modulo permutation of the rows. a(2) = 5: [11/00], [10/01], [100/010/000], [100/001/000], [010/001/000].
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..200
FORMULA
a(n) = Sum_{k=0..n+1} A360693(k,n).
EXAMPLE
a(0) = 2: {}, {b}.
a(1) = 3: {a}, {ab,bb}, {ba,bb}.
a(3) = 16: {aa,ab}, {aa,ba}, {aab,abb,bbb}, {aab,bab,bbb}, {aab,bba,bbb}, {aba,abb,bbb}, {aba,bab,bbb}, {aba,bba,bbb}, {abb,baa,bbb}, {abb,bab,bba}, {baa,bab,bbb}, {baa,bba,bbb}, {abbb,babb,bbab,bbbb}, {abbb,babb,bbba,bbbb}, {abbb,bbab,bbba,bbbb}, {babb,bbab,bbba,bbbb}.
MAPLE
g:= proc(n, i, j) option remember; expand(`if`(j=0, 1, `if`(i<0, 0, add(
g(n, i-1, j-k)*x^(i*k)*binomial(binomial(n, i) , k), k=0..j))))
end:
a:= n-> add(coeff(g(k$3), x, n), k=0..n+1):
seq(a(n), n=0..23);
MATHEMATICA
g[n_, i_, j_] := g[n, i, j] = Expand[If[j == 0, 1, If[i < 0, 0, Sum[g[n, i - 1, j - k]*x^(i*k)*Binomial[Binomial[n, i] , k], {k, 0, j}]]]];
a[n_] := Sum[Coefficient[g[k, k, k], x, n], {k, 0, n + 1}];
Table[a[n], {n, 0, 23}] (* Jean-François Alcover, Dec 13 2023, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Feb 16 2023
STATUS
approved