OFFSET
0,3
COMMENTS
The ordering of the functions f_k is defined in A215703: f_1, f_2, ... = x, x^x, x^(x^2), x^(x^x), x^(x^3), x^(x^x*x), x^(x^(x^2)), x^(x^(x^x)), x^(x^4), x^(x^x*x^2), ... .
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..400
FORMULA
a(n) = A215703(n,n+1).
EXAMPLE
a(0) = x_{x=1} = 1.
a(1) = (d/dx x^x)_{x=1} = (x^x*(log(x)+1))_{x=1} = log(1)+1 = 1.
a(2) = (d^2/dx^2 x^(x^2))_{x=1} = (x^(x^2) * (2*x*log(x)+x)^2 + x^(x^2) * (2*log(x)+3))_{x=1} = (2*log(1)+1)^2 + 2*log(1)+3 = 4.
a(3) = (d^3/dx^3 x^(x^x))_{x=1} = 9.
a(4) = (d^4/dx^4 x^(x^3))_{x=1} = 156.
MAPLE
T:= proc(n) T(n):=`if`(n=1, [x], map(h-> x^h, g(n-1$2))) end:
g:= proc(n, i) option remember; `if`(i=1, [x^n], [seq(seq(
seq(mul(T(i)[w[t]-t+1], t=1..j)*v, v=g(n-i*j, i-1)), w=
combinat[choose]([$1..nops(T(i))+j-1], j)), j=0..n/i)])
end:
f:= proc() local i, l; i, l:= 0, []; proc(n) while n>
nops(l) do i:= i+1; l:= [l[], T(i)[]] od; l[n] end
end():
a:= n-> n!*coeff(series(subs(x=x+1, f(n+1)), x, n+1), x, n):
seq(a(n), n=0..23);
MATHEMATICA
T[n_] := If[n == 1, {x}, Map[x^# &, g[n - 1, n - 1]]];
g[n_, i_] := g[n, i] = If[i == 1, {x^n}, Flatten@ Table[ Table[ Table[ Product[T[i][[w[[t]] - t + 1]], {t, 1, j}]*v, {v, g[n - i*j, i - 1]}], {w, Subsets[Range[Length[T[i]] + j - 1], {j}]}], {j, 0, n/i}]];
f[n_] := Module[{i = 0, l = {}}, While[n > Length[l], i++; l = Join[l, T[i]]]; l[[n]]];
a[n_] := n!*SeriesCoefficient[f[n+1] /. x -> x+1, {x, 0, n}];
Table[a[n], {n, 0, 23}] (* Jean-François Alcover, Sep 15 2022, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Mar 06 2019
STATUS
approved