[go: up one dir, main page]

login
Number of non-isomorphic strict multiset partitions of weight n.
104

%I #15 Jan 21 2023 17:58:10

%S 1,1,3,8,23,63,197,588,1892,6140,20734,71472,254090,923900,3446572,

%T 13149295,51316445,204556612,832467052,3455533022,14621598811,

%U 63023667027,276559371189,1234802595648,5606647482646,25875459311317,121324797470067,577692044073205

%N Number of non-isomorphic strict multiset partitions of weight n.

%C Also the number of nonnegative integer n X n matrices with sum of elements equal to n, under row and column permutations, with no equal rows (or alternatively, with no equal columns).

%C Also the number of non-isomorphic multiset partitions of weight n with no equivalent vertices. In a multiset partition, two vertices are equivalent if in every block the multiplicity of the first is equal to the multiplicity of the second.

%H Andrew Howroyd, <a href="/A316980/b316980.txt">Table of n, a(n) for n = 0..50</a>

%F Euler transform of A319557. - _Gus Wiseman_, Sep 23 2018

%e Non-isomorphic representatives of the a(3) = 8 multiset partitions with no equivalent vertices (first column) and with no equal blocks (second column):

%e (111) <-> (111)

%e (122) <-> (1)(11)

%e (1)(11) <-> (122)

%e (1)(22) <-> (1)(22)

%e (2)(12) <-> (2)(12)

%e (1)(1)(1) <-> (123)

%e (1)(2)(2) <-> (1)(23)

%e (1)(2)(3) <-> (1)(2)(3)

%o (PARI)

%o EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}

%o 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}

%o K(q, t, k)={EulerT(Vec(sum(j=1, #q, my(g=gcd(t, q[j])); g*x^(q[j]/g)) + O(x*x^k), -k))}

%o a(n)={if(n==0, 1, my(s=0); forpart(q=n, my(p=sum(t=1, n, subst(x*Ser(K(q, t, n\t))/t, x, x^t))); s+=permcount(q)*polcoef(exp(p-subst(p,x,x^2)), n)); s/n!)} \\ _Andrew Howroyd_, Jan 21 2023

%Y Cf. A000009, A001055, A007716, A007717, A020555, A045778, A130091.

%Y Cf. A316974, A316978, A316979, A316981, A316983.

%Y Cf. A049311, A059201, A319558, A319560, A319567.

%K nonn

%O 0,3

%A _Gus Wiseman_, Jul 18 2018

%E a(7)-a(10) from _Gus Wiseman_, Sep 23 2018

%E Terms a(11) and beyond from _Andrew Howroyd_, Jan 19 2023