%I #8 Dec 30 2020 17:11:24
%S 1,6,21,104,452,2335,11992,66810,385101,2336352,14738380,96831730,
%T 659809115,4657075074,33974259046,255781455848,1984239830571,
%U 15839628564349,129951186405574,1094486382191624,9453318070371926,83654146992936350,757769011659766015,7020652591448497490
%N Number of non-isomorphic multiset partitions of multiset partitions of multisets of size n.
%e Non-isomorphic representatives of the a(3) = 21 multiset partitions of multiset partitions:
%e {{{1,1,1}}}
%e {{{1,1,2}}}
%e {{{1,2,3}}}
%e {{{1},{1,1}}}
%e {{{1},{1,2}}}
%e {{{1},{2,3}}}
%e {{{2},{1,1}}}
%e {{{1},{1},{1}}}
%e {{{1},{1},{2}}}
%e {{{1},{2},{3}}}
%e {{{1}},{{1,1}}}
%e {{{1}},{{1,2}}}
%e {{{1}},{{2,3}}}
%e {{{2}},{{1,1}}}
%e {{{1}},{{1},{1}}}
%e {{{1}},{{1},{2}}}
%e {{{1}},{{2},{3}}}
%e {{{2}},{{1},{1}}}
%e {{{1}},{{1}},{{1}}}
%e {{{1}},{{1}},{{2}}}
%e {{{1}},{{2}},{{3}}}
%t sps[{}]:={{}};sps[set:{i_,___}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,___}];
%t mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
%t strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]&/@IntegerPartitions[n];
%t dubnorm[m_]:=First[Union[Table[Map[Sort,m/.Rule@@@Table[{Union[Flatten[m]][[i]],Union[Flatten[m]][[perm[[i]]]]},{i,Length[perm]}],{0,2}],{perm,Permutations[Union[Flatten[m]]]}]]];
%t Table[Length[Union[dubnorm/@Join@@mps/@Join@@mps/@strnorm[n]]],{n,5}]
%o (PARI) \\ See links in A339645 for combinatorial species functions.
%o seq(n)={my(A=sExp(symGroupSeries(n))); NumUnlabeledObjsSeq(sCartProd(A, sExp(A)-1))} \\ _Andrew Howroyd_, Dec 30 2020
%Y Cf. A001970, A007716, A050336, A050338, A255906, A269134, A317533, A317791, A318564, A318565.
%K nonn
%O 1,2
%A _Gus Wiseman_, Aug 29 2018
%E Terms a(8) and beyond from _Andrew Howroyd_, Dec 30 2020