cp's OEIS Frontend

A318402 Number of sets of nonempty sets whose multiset union is a strongly normal multiset of size n.

%I A318402 #9 Dec 30 2020 17:34:16
%S A318402 1,2,6,20,74,311,1401,6913,36376,205421,1228288,7786802,51937607,
%T A318402 364250763,2673314121,20504809133,163844631872,1361874185139,
%U A318402 11748149246269,105029750531640,971403871953460,9282643841237360,91519776792040324,929892817423282068,9725646244888190337
%N A318402 Number of sets of nonempty sets whose multiset union is a strongly normal multiset of size n.
%C A318402 A multiset is strongly normal if it spans an initial interval of positive integers with weakly decreasing multiplicities.
%H A318402 Andrew Howroyd, <a href="/A318402/b318402.txt">Table of n, a(n) for n = 1..50</a>
%e A318402 The a(4) = 20 sets of sets:
%e A318402   {{1,2,3,4}}
%e A318402   {{1},{1,2,3}}
%e A318402   {{1},{2,3,4}}
%e A318402   {{2},{1,3,4}}
%e A318402   {{3},{1,2,4}}
%e A318402   {{4},{1,2,3}}
%e A318402   {{1,2},{1,3}}
%e A318402   {{1,2},{3,4}}
%e A318402   {{1,3},{2,4}}
%e A318402   {{1,4},{2,3}}
%e A318402   {{1},{2},{1,2}}
%e A318402   {{1},{2},{1,3}}
%e A318402   {{1},{2},{3,4}}
%e A318402   {{1},{3},{1,2}}
%e A318402   {{1},{3},{2,4}}
%e A318402   {{1},{4},{2,3}}
%e A318402   {{2},{3},{1,4}}
%e A318402   {{2},{4},{1,3}}
%e A318402   {{3},{4},{1,2}}
%e A318402   {{1},{2},{3},{4}}
%o A318402 (PARI)
%o A318402 WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v,n,(-1)^(n-1)/n))))-1,-#v)}
%o A318402 D(p, n)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); my(u=WeighT(v)); Vec(1/prod(k=1, n, 1 - u[k]*x^k + O(x*x^n))-1,-n)/prod(i=1, #v, i^v[i]*v[i]!)}
%o A318402 seq(n)={my(s); for(k=1, n, forpart(p=k, s+=(-1)^(k+#p)*D(p,n))); s[n]+=1; s/2} \\ _Andrew Howroyd_, Dec 30 2020
%Y A318402 Cf. A007716, A049311, A050326, A116540, A283877, A292432, A292444, A318361, A318369, A318370.
%K A318402 nonn
%O A318402 1,2
%A A318402 _Gus Wiseman_, Aug 25 2018
%E A318402 Terms a(10) and beyond from _Andrew Howroyd_, Dec 30 2020