# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a317653
Showing 1-1 of 1

%I A317653 #11 Sep 15 2018 02:05:27
%S A317653 1,3,34,602,14872,472138,18323359,840503724,44489123726,2668985463839,
%T A317653 178960530393633,13263068003965046,1076580864432281157,
%U A317653 94987639225399100006,9051397653144246683937,926407121115738135640677,101357200280211387377806719,11804887470887800839909147484
%N A317653 Number of free pure symmetric multifunctions whose leaves are a normal multiset of size n.
%C A317653 A multiset is normal if it spans an initial interval of positive integers. A free pure symmetric multifunction f in EPSM is either (case 1) a positive integer, or (case 2) an expression of the form h[g_1, ..., g_k] where k > 0, h is in EPSM, each of the g_i for i = 1, ..., k is in EPSM, and for i < j we have g_i <= g_j under a canonical total ordering of EPSM, such as the Mathematica ordering of expressions.
%H A317653 Andrew Howroyd, <a href="/A317653/b317653.txt">Table of n, a(n) for n = 1..200</a>
%e A317653 The a(3) = 34 free pure symmetric multifunctions:
%e A317653 1[1[1]], 1[1,1], 1[1][1],
%e A317653 1[2[2]], 1[2,2], 2[1[2]], 2[2[1]], 2[1,2], 1[2][2], 2[1][2], 2[2][1],
%e A317653 1[1[2]], 1[2[1]], 1[1,2], 2[1[1]], 2[1,1], 1[1][2], 1[2][1], 2[1][1],
%e A317653 1[2[3]], 1[3[2]], 1[2,3], 2[1[3]], 2[3[1]], 2[1,3], 3[1[2]], 3[2[1]], 3[1,2], 1[2][3], 2[1][3], 1[3][2], 3[1][2], 2[3][1], 3[2][1].
%t A317653 sps[{}]:={{}};sps[set:{i_,___}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,___}];
%t A317653 mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
%t A317653 exprUsing[m_]:=exprUsing[m]=If[Length[m]==0,{},If[Length[m]==1,{First[m]},Join@@Cases[Union[Table[PR[m[[s]],m[[Complement[Range[Length[m]],s]]]],{s,Take[Subsets[Range[Length[m]]],{2,-2}]}]],PR[h_,g_]:>Join@@Table[Apply@@@Tuples[{exprUsing[h],Union[Sort/@Tuples[exprUsing/@p]]}],{p,mps[g]}]]]];
%t A317653 got[y_]:=Join@@Table[Table[i,{y[[i]]}],{i,Range[Length[y]]}];
%t A317653 Table[Sum[Length[exprUsing[got[y]]],{y,Join@@Permutations/@IntegerPartitions[n]}],{n,6}]
%o A317653 (PARI) \\ here R(n,1) is A052893.
%o A317653 EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
%o A317653 R(n,k)={my(v=[k]); for(n=2, n, my(t=EulerT(v)); v=concat(v, sum(k=1, n-1, v[k]*t[n-k]))); v}
%o A317653 seq(n)={sum(k=1, n, R(n,k)*sum(r=k, n, binomial(r,k)*(-1)^(r-k)) )} \\ _Andrew Howroyd_, Sep 14 2018
%Y A317653 Cf. A001003, A052893, A053492, A255906, A277996, A279944, A280000.
%Y A317653 Cf. A317652, A317654, A317655, A317656, A317658.
%K A317653 nonn
%O A317653 1,2
%A A317653 _Gus Wiseman_, Aug 03 2018
%E A317653 Terms a(8) and beyond from _Andrew Howroyd_, Sep 14 2018

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE