%I #161 May 01 2024 10:04:24
%S 1,1,1,1,1,3,1,1,4,1,4,7,6,1,1,7,4,11,1,5,30,15,10,25,10,1,1,11,15,32,
%T 34,26,1,6,57,34,146,31,15,120,90,20,65,15,1,1,16,26,15,76,192,34,122,
%U 180,57,1,7,98,140,406,462,588,63,21,252,154,896,301
%N Weights for expansion of iterated derivatives, powers of a Lie derivative, or infinitesimal generator in vector form, (f(x)D_x)^n; coefficients of A-polynomials of Comtet; Scherk partition polynomials.
%C This entry and the references differ slightly among themselves in the order of coefficients for higher order terms. Table on p. 167 of Comtet has at least one index error.
%C Let F[FI(x)] = FI[F(x)] = x (i.e., F and FI are a compositional inverse pair) about x=0 with F(0)=FI(0)=0. Define f(x) = 1/[dFI(x)/dx], then for w(x) analytic about the origin, exp[t f(x)d/dx] w(x) = w{F[t+FI(x)]} = q(t,x) with q{t,F[s+FI(x)]} = q(t+s,x). See A145271 for w(x)=x and note that A145271 is embedded in A139605. E.g.f. of the binomial Sheffer sequence associated to F(x) is exp[x f(z)d/dz] exp(t*z)= exp{t*F[x+FI(z)]} evaluated at z=0. - _Tom Copeland_, Oct 19 2011
%C dq(t,x)/dt - f(x)dq(t,x)/dx = 0, so (1,-f(x)) gives the components of a vector orthogonal to the gradient of q and therefore tangent to the contour of q at (t,x). - _Tom Copeland_, Oct 26 2011
%C The formula exp[t f(x)d/dx] w(x)= w{F[t+FI(x)]} above is implicit in the chain rule formulas on pages 10 and 12 of Mathemagical Forests. Another derivation is alluded to in the Dattoli reference in A080635 (repeated below). - _Tom Copeland_, Nov 28 2011
%C Let f(x) and g(x) be two infinitely differential functions. Denote f_0 = f(x), f_1 = df_0/dx, f_2 = df_1/dx, and so on. Same with g_0 = g(x). Define the linear operator L(u(x)) = g(x) * du(x)/dx. Denote F_1 = L(f(x)), F_2 = L(F_1), and so on. When n>0, F_n is a linear combination of f_1, ..., f_n where each f_k is multiplied by a homogeneous polynomial (P(n,k)) of degree n in g_0, ..., g_{n-1}. The triangle of the sum of the coefficients of P(n,k) is A130534. - _Michael Somos_, Mar 23 2014
%C Triangle with row n length A000070(n+1) and row n consists of the coefficients: P(n,1), ..., P(n,n). The order of coefficients in P(n,k) is Abramowitz and Stegun order for partitions of n-k with parts g_1, ..., g_{n-k}. - _Michael Somos_, Mar 23 2014
%C A130534(n,k) gives the number of rooted trees with (k+1) trunks that are associated with D^(k+1) in the forest of "naturally grown" rooted trees with (n+2) nodes, or vertices, that are associated with R^(n+1) in the example below. Cf. MF link. - _Tom Copeland_, Mar 23 2014
%C These partition polynomials appeared in 1823 in a dissertation by Heinrich Scherk. See p. 76 of Blasiak and Flajolet. - _Tom Copeland_, Jul 14 2021
%C Schröder made use of iterated derivatives, or iterated infinitesimal generators (IGs), ((1/f') D)^n in his investigations of functional iteration, or iterated functional composition, related to extensions of Newton's method of finding zeros of an equation. He constructs the series, in terms of the IGs, for finv[t+f(z)] evaluated at t = -f(z), giving z_1 = finv(0) although he doesn't present his analysis this way. - _Tom Copeland_, Jul 19 2021
%D F. Bergeron, G. Labelle and P. Leroux, Combinatorial Species and Tree-like Structures, (1997), Cambridge University Press, p. 386.
%D H. Davis, The Theory of Linear Operators, (1936), The Principia Press, p. 13.
%D T. Mansour and M. Schork, Commutation Relations, Normal Ordering, and Stirling Numbers, Chapman and Hall/CRC, 2015.
%H P. Blasiak and P. Flajolet, <a href="https://arxiv.org/abs/1010.0354">Combinatorial Models of Creation-Annihilation</a>, arXiv:1010.0354 [math.CO], 2011.
%H Louis Comtet, <a href="https://gallica.bnf.fr/ark:/12148/bpt6k6217213f/f179.item">Une formule explicite pour les puissances successives de l'opérateur de dérivation de Lie</a>, Comptes Rendus Acad. Sc. Paris, Serie A tome 276 (1973), pp. 165 - 168.
%H Tom Copeland, <a href="https://tcjpn.wordpress.com/2008/06/12/mathemagical-forests/">Mathemagical Forests v2</a>, 2008.
%H Tom Copeland, <a href="http://tcjpn.wordpress.com/2012/11/29/infinigens-the-pascal-pyramid-and-the-witt-and-virasoro-algebras/">Infinigens, the Pascal Triangle, and the Witt and Virasoro Algebras</a>, 2012.
%H Tom Copeland, <a href="https://tcjpn.wordpress.com/2018/07/10/pre-lie-algebra-and-cayley-trees//">Pre-Lie algebras, Cayley's analytic trees, and mathemagical forests</a>, 2018.
%H G. Datolli, P. L. Ottaviani, A. Torre and L. Vazquez, <a href="http://dx.doi.org/10.1007/BF02907529">Evolution operator equations: integration with algebraic and finite differences methods.[...]</a>, La Rivista del Nuovo Cimento 20,2 (1997) 1-133. eq. (I.2.18).
%H D. Grinberg, <a href="https://arxiv.org/abs/1908.09179">Commutators, matrices, and an identity of Copeland</a>, arXiv:1908.09179 [math.RA], 2019.
%H Guo-Niu Han and Shi-mei Ma, <a href="https://irma.math.unistra.fr/~guoniu/papers/p118gindex.pdf">Eulerian polynomials and the g-indices of Young tableaux</a>, Proc. Amer. Math. Soc., (2023).
%H Kentaro Ihara, <a href="http://dx.doi.org/10.1016/j.jpaa.2011.06.004">Derivations and automorphisms on non-commutative power series,</a> Journal of Pure and Applied Algebra, Volume 216, Issue 1, January 2012, Pages 192-201.
%H T. Mansour and M. Shork, <a href="http://dx.doi.org/10.1016/j.amc.2013.04.010">The generalized Touchard polynomials revisited</a>, Journal of Applied Mathematics and Computation, Volume 219, Issue 19, June 2013, Pages 9978-9991.
%H MathOverflow, <a href="https://mathoverflow.net/questions/214927/important-formulas-in-combinatorics/215203#215203">Important formulas in combinatorics: The combinatorics underlying iterated derivatives (infinitesimal Lie generators) for compositional inversion and flow maps for vector fields</a>, answer by Tom Copeland to an MO question posed by Gil Kalai, 2015.
%H MathOverflow, <a href="https://mathoverflow.net/questions/337766/expansions-of-iterated-or-nested-derivatives-or-vectors-conjectured-matrix-c">Expansions of iterated, or nested, derivatives, or vectors--conjectured matrix computation</a>, a MO question posed by Tom Copeland, answered by Darij Grinberg, 2019.
%H Mathoverflow, <a href="https://mathoverflow.net/questions/41039/formula-for-n-th-iteration-of-dx-dt-bx">Formula for n-th iteration of dx/dt=B(x)</a>, a question on MathOverflow posed by the user resolvent and answered by Tom Copeland, 2021.
%H H. Scherk, <a href="https://gdz.sub.uni-goettingen.de/id/PPN271057963">De evolvenda functione (yd.yd.yd...ydX/dxn) disquisitiones nonnullae analyticae</a>, PhD thesis, Berlin, 1823. Scanned copy at Göttinger Digitalisierungszentrum (GDZ).
%H E. Schröder, <a href="https://www.semanticscholar.org/paper/Ueber-unendlich-viele-Algorithmen-zur-Aufl%C3%B6sung-der-Schr%C3%B6der/7d37b4dbf960770e926575456b9504f5e785b048">Ueber unendlich viele Algorithmen zur Auflösung der Gleichungen</a>, Mathematische Annalen vol. 2, 317-365, 1870.
%H G. Stewart, <a href="https://drum.lib.umd.edu/handle/1903/577">On infinitely many algorithms for solving equations</a>, 1993, (translation into English of Schröder's paper above).
%F Equivalent matrix computation: Multiply the n-th diagonal (with n=0 the main diagonal) of the lower triangular Pascal matrix by f_n = (d/dx)^n f(x) to obtain the matrix VP with VP(n,k) = binomial(n,k) f_(n-k). Then R^n = (1, 0, 0, 0,..) [VP * S]^n (1, D, D^2, ..)^T, where S is the shift matrix A129185, representing differentiation in the basis x^n/n!. Cf. A145271. - _Tom Copeland_, Jul 17 2016
%F A formula for the coefficients of this matrix is presented in Ihara, obtained from Comtet. - _Tom Copeland_, Mar 25 2020
%F Elaborating on my 2011 comments: Let exp[x F(t)] = exp[t p.(x)] be the e.g.f. for the binomial Sheffer sequence of polynomials (p.(x))^n = p_n(x). Then, evaluated at x = t = 0, the coefficient p_(n,k) = (D_x^k/k!) p_n(x) = D_t^n [F(t)]^k/k! = (f(x)D_x)^n x^k/k! = R^n x^k/k!, and so p_(n,k) is the coefficient of D^k of the operator R^n below evaluated at x=0. - _Tom Copeland_, May 14 2020
%F Per earlier comments, the action of the differentials of this entry on an exponential is exp(x g(u)D_u) e^(ut) = e^(t H^{(-1)}(H(u)+x)) with g(x) = 1/D(H(x)) and H^{(-1)}(x) the compositional inverse of H(x). With H^{(-1)}(x) = -log(1-x), the inverse about x=0 is H(x) = 1-e^(-x), giving g(x) = e^x and the resulting action e^(-t log(1-x)) = (1-x)^(-t) for u=0, an e.g.f. for the unsigned Stirling numbers of the first kind (see A008275, A048994, and A130534). Consequently, summing the coefficients of this entry over each associated derivative gives these Stirling numbers. E.g., the fifth row in the examples reduces to (1+4+1) D + (7+4) D^2 + 6 D^3 + D^4 = 6 D + 11 D^2 + 6 D^3 + D^4. The Connes-Moscovici weights A139002 are a refinement of this entry. - _Tom Copeland_, Jul 14 2021
%e Let R = f(x)d/dx = f(x)D and (j,k) = [(d/dx)^j f(x)]^k ; then
%e R^0 = 1.
%e R^1 = (0,1)D.
%e R^2 = (0,1)(1,1)D + (0,2)D^2.
%e R^3 = [(0,1)(1,2) + (0,2)(2,1)]D + 3 (0,2)(1,1)D^2 + (0,3)D^3.
%e R^4 = [(0,1)(1,3) + 4 (0,2)(1,1)(2,1) + (0,3)(3,1)]D +
%e [7 (0,2)(1,2) + 4 (0,3)(2,1)]D^2 + 6 (0,3)(1,1)D^3 + (0,4)D^4. - _Tom Copeland_, Jun 12 2008
%e R^5 = [(0,1)(1,4) + 11 (0,2)(1,2)(2,2) + 4 (0,3)(2,2) + (0,4)(4,1) + 7 (0,3)(1,1)(3,1)]D + [15 (0,2)(1,3) + 30 (0,3)(1,1)(2,1) + 5 (0,4)(1,3)]D^2 + [25 (0,3)(1,2) + 10 (0,4)(2,1) + 25 (0,3)(1,2)]D^3 + 10 (0,4)(1,1)D^4 + (0,5)D^5. - _Tom Copeland_, Jul 17 2016
%e R^6 = [(0,1)(1,5) + 26 (0,2)(1,3)(2,1) + 34 (0,3)(1,1)(2,2) + 32 (0,3)(1,2)(3,1) + 11 (0,4)(1,1)(4,1) + 15 (0,4)(2,1)(3,1) + (0,5)(1,5)]D + [31 (0,2)(1,4) + 146 (0,3)(1,2)(2,1) + 57 (0,4)(1,1)(3,1) + 34 (0,4)(2,2) + 6 (0,5)(4,1)]D^2 + [90 (0,3)(1,3) + 120 (0,4)(1,1)(2,1) + 15 (0,5)(3,1)]D^3 + [65 (0,4)(1,2) + 20 (0,5)(2,1)]D^4 + 15 (0,5)(1,1)D^5 + (0,6)D^6. - _Tom Copeland_, Jul 17 2016
%e ------------
%e F_1 = (1*g_0) * f_1, F_2 = (1*g_0*g_1) * f_1 + (1*g_0^2) * f_2, F_3 = (1*g_0*g_1^2 + 1*g_0^2*g_2) * f_1 + (3*g_0^2*g_1) * f_2 + (1*g_0^3) * f_3. - _Michael Somos_, Mar 23 2014
%e P(4,2) = 4*g0^3*g2 + 7*g0^2*g1^2. P(5,2) = 5*g0^4*g3 + 30*g0^3*g1*g2 + 15*g0^2*g1^3. - _Michael Somos_, Mar 23 2014
%e 1
%e 1 , 1
%e 1 1 , 3 , 1
%e 1 4 1 , 4 7 , 6 , 1
%e 1 7 4 11 1, 5 30 15 , 10 25 , 10 , 1
%e 1 11 15 32 34 26 1 , 6 57 34 146 31 , 15 120 90 , 20 65 , 15 , 1
%t row[n_] := With[{pn = CoefficientRules[Nest[g[x] D[#, x] &, f[x], n], Derivative[#][f][x] & /@ Range[n]][[;; , 2]] /. {Derivative[k_][g][x] :> h[k], g[x] -> 1}}, Table[Coefficient[pn[[k]], Product[h[x], {x, p}]], {k, n - 1}, {p, Sort[Sort /@ IntegerPartitions[n - k]]}]~Join~{{1}}];
%t Table[row[n], {n, 7}] // Flatten (* _Andrey Zabolotskiy_, Mar 08 2024 *)
%Y Cf. A000070 (number of distinct terms for each order).
%Y Cf. A130534 (sum of numerical coefficients of the derivatives).
%Y Cf. A008275, A048994, A080635, A124796, A129185, A139002, A145271.
%K nonn,tabf
%O 1,6
%A _Tom Copeland_, Jun 12 2008
%E Title expanded by _Tom Copeland_, Mar 17 2014
%E Sequence terms rearranged in Abramowitz and Stegun order by _Michael Somos_, Mar 23 2014
%E Title expanded by _Tom Copeland_, Jul 14 2021