Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #26 Nov 23 2016 16:44:05
%S 1,2,9,88,905,12666,220297,4506440,104707521,2758506850
%N E.g.f. equals the limit of the average of all permutations of the compositions of the functions x*exp(x^k), for k=1..n, as n increases.
%F E.g.f.: A(x) = limit_{n->oo} (1/n!) * Sum(all permutations of compositions of functions x*exp(x^k) for k=1..n).
%e E.g.f.: A(x) = x + 2*x^2/2! + 9*x^3/3! + 88*x^4/4! + 905*x^5/5! + 12666*x^6/6! + 220297*x^7/7! + 4506440*x^8/8! + 104707521*x^9/9! + 2758506850*x^10/10! +...
%e Generating method.
%e Define F(n,x) as the average of the sum over all n! permutations of the compositions of x*exp(x^k) for k=1..n, then the e.g.f. of this sequence is the limit of the functions F(n,x) as n grows.
%e Examples of some initial functions F(n,x) are as follows.
%e At n=1, F(1,x) = x*exp(x).
%e At n=2, F(2,x) = (1/2!)*([x*exp(x) o x*exp(x^2)] + [x*exp(x^2) o x*exp(x)]) = x + 2*x^2/2! + 9*x^3/3! + 64*x^4/4! + 425*x^5/5! + 3486*x^6/6! +...
%e At n=3, F(3,x) = (1/3!)*([x*exp(x) o x*exp(x^2) o x*exp(x^3)] + [x*exp(x^2) o x*exp(x) o x*exp(x^3)] + [x*exp(x) o x*exp(x^3) o x*exp(x^2)] + [x*exp(x^3) o x*exp(x) o x*exp(x^2)] + [x*exp(x^2) o x*exp(x^3) o x*exp(x)] + [x*exp(x^3) o x*exp(x^2) o x*exp(x)]) = x + 2*x^2/2! + 9*x^3/3! + 88*x^4/4! + 785*x^5/5! + 9426*x^6/6! +...
%e etc.
%o (PARI) /* Informal code to generate terms (ran over 100 hrs for 10 terms!) */
%o {F(n) = x*exp(x^n +x*O(x^16))}
%o {P(n) = vector(n!,k, numtoperm(n,k) )}
%o {A(n) = my(H=0); for(k=1,n!, G=x; for(j=1,n, G=subst(G,x, F(P(n)[k][j]) )); H=H+G);H/n!}
%o {Vec(serlaplace(A(6)))}
%Y Cf. A278333(n) = a(n+1)/(n+1) for n>=0.
%Y Cf. A277180, A277181.
%K nonn,more
%O 1,2
%A _Paul D. Hanna_, Nov 18 2016