A360638 - OEIS

%I #16 Dec 09 2023 09:00:03

%S 1,3,16,100,593,3497,20316,116378,658214,3679450,20350028,111459648,

%T 605060633,3257784589,17408647968,92378535290,487031130699,

%U 2552197485757,13298890952222,68930923717598,355507581655752,1824924721216084,9326440815314046,47464093855706540

%N Number of sets of nonempty words over binary alphabet where each letter occurs n times.

%H Alois P. Heinz, <a href="/A360638/b360638.txt">Table of n, a(n) for n = 0..380</a>

%F a(n) = A360634(2n,n).

%F a(n) mod 2 = 1 <=> n in { A080277 } U {0}.

%e a(0) = 1: {}.

%e a(1) = 3: {ab}, {ba}, {a,b}.

%e a(2) = 16: {aabb}, {abab}, {abba}, {baab}, {baba}, {bbaa}, {a,abb}, {a,bab}, {a,bba}, {aa,bb}, {aab,b}, {ab,ba}, {aba,b}, {b,baa}, {a,ab,b}, {a,b,ba}.

%p g:= proc(n, i, j) option remember; expand(`if`(j=0, 1, `if`(i<0, 0, add(

%p       g(n, i-1, j-k)*x^(i*k)*binomial(binomial(n, i), k), k=0..j))))

%p     end:

%p b:= proc(n, i) option remember; expand(`if`(n=0, 1,

%p      `if`(i<1, 0, add(b(n-i*j, i-1)*g(i$2, j), j=0..n/i))))

%p     end:

%p a:= n-> coeff(b(2*n$2), x, n):

%p seq(a(n), n=0..31);

%t g[n_, i_, j_] := g[n, i, j] = Expand[If[j == 0, 1, If[i < 0, 0, Sum[

%t     g[n, i - 1, j - k]*x^(i*k)*Binomial[Binomial[n, i], k], {k, 0, j}]]]];

%t b[n_, i_] := b[n, i] = Expand[If[n == 0, 1,

%t     If[i < 1, 0, Sum[b[n - i*j, i - 1]*g[i, i, j], {j, 0, n/i}]]]];

%t a[n_] := Coefficient[b[2n, 2n], x, n];

%t Table[a[n], {n, 0, 31}] (* _Jean-François Alcover_, Dec 09 2023, after _Alois P. Heinz_ *)

%Y Cf. A080277, A360626 (the same for multisets), A360634.

%K nonn

%O 0,2

%A _Alois P. Heinz_, Feb 14 2023