OFFSET
0,4
COMMENTS
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..50
EXAMPLE
Non-isomorphic representatives of the a(2) = 1 through a(4) = 12 strict multiset partitions with empty intersection:
2: {{1},{2}}
3: {{1},{2,2}}
{{1},{2,3}}
{{1},{2},{3}}
4: {{1},{2,2,2}}
{{1},{2,3,3}}
{{1},{2,3,4}}
{{1,1},{2,2}}
{{1,2},{3,3}}
{{1,2},{3,4}}
{{1},{2},{1,2}}
{{1},{2},{2,2}}
{{1},{2},{3,3}}
{{1},{2},{3,4}}
{{1},{3},{2,3}}
{{1},{2},{3},{4}}
PROG
(PARI)
EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
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, k)={EulerT(Vec(sum(j=1, #q, my(g=gcd(t, q[j])); g*x^(q[j]/g)) + O(x*x^k), -k))}
R(q, n)={vector(n, t, subst(x*Ser(K(q, t, n\t)/t), x, x^t))}
a(n)={my(s=0); forpart(q=n, my(f=prod(i=1, #q, 1 - x^q[i]), u=R(q, n)); s+=permcount(q)*sum(k=0, n, my(c=polcoef(f, k)); if(c, c*polcoef(exp(sum(t=1, n\(k+1), x^(t*k)*u[t] - subst(x^(t*k)*u[t] + O(x*x^(n\2)), x, x^2), O(x*x^n) ))*if(k, 1+x^k, 1), n))) ); s/n!} \\ Andrew Howroyd, May 30 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Sep 27 2018
EXTENSIONS
Terms a(11) and beyond from Andrew Howroyd, May 30 2023
STATUS
approved