login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A145271 Coefficients for expansion of (g(x)d/dx)^n g(x); refined Eulerian numbers for calculating compositional inverse of h(x) = (d/dx)^(-1) 1/g(x); iterated derivatives as infinitesimal generators of flows. 41

%I #259 Nov 22 2022 22:55:41

%S 1,1,1,1,1,4,1,1,11,4,7,1,1,26,34,32,15,11,1,1,57,180,122,34,192,76,

%T 15,26,16,1,1,120,768,423,496,1494,426,294,267,474,156,56,42,22,1,1,

%U 247,2904,1389,4288,9204,2127,496,5946,2829,5142,1206,855,768,1344,1038,288,56,98,64,29,1

%N Coefficients for expansion of (g(x)d/dx)^n g(x); refined Eulerian numbers for calculating compositional inverse of h(x) = (d/dx)^(-1) 1/g(x); iterated derivatives as infinitesimal generators of flows.

%C For more detail, including connections to Legendre transformations, rooted trees, A139605, A139002 and A074060, see Mathemagical Forests p. 9.

%C For connections to the h-polynomials associated to the refined f-polynomials of permutohedra see my comments in A008292 and A049019.

%C From _Tom Copeland_, Oct 14 2011: (Start)

%C Given analytic functions F(x) and FI(x) such that F(FI(x))=FI(F(x))=x about 0, i.e., they are compositional inverses of each other, then, with g(x) = 1/dFI(x)/dx, a flow function W(s,x) can be defined with the following relations:

%C W(s,x) = exp(s g(x)d/dx)x = F(s+FI(x)) <flow fct.>,

%C W(s,0) = F(s) <orbit of the flow>,

%C W(0,x) = x <identity property>,

%C dW(0,x)/ds = g(x) = F'[FI(x)] <infinitesimal generator>, implying

%C dW(0,F(x))/ds = g(F(x)) = F'(x) <autonomous diff. eqn.>, and

%C W(s,W(r,x)) = F(s+FI(F(r+FI(x)))) = F(s+r+FI(x)) = W(s+r,x) <group property>. (See MF link below.) (End)

%C dW(s,x)/ds - g(x)dW(s,x)/dx = 0, so (1,-g(x)) are the components of a vector orthogonal to the gradient of W and, therefore, tangent to the contour of W, at (s,x) <tangency property>. - _Tom Copeland_, Oct 26 2011

%C Though A139605 contains A145271, the op. of A145271 contains that of A139605 in the sense that exp(s g(x)d/dx) w(x) = w(F(s+FI(x))) = exp((exp(s g(x)d/dx)x)d/du)w(u) evaluated at u=0. This is reflected in the fact that the forest of rooted trees assoc. to (g(x)d/dx)^n, FOR_n, can be generated by removing the single trunk of the planted rooted trees of FOR_(n+1). - _Tom Copeland_, Nov 29 2011

%C Related to formal group laws for elliptic curves (see Hoffman). - _Tom Copeland_, Feb 24 2012

%C The functional equation W(s,x) = F(s+FI(x)), or a restriction of it, is sometimes called the Abel equation or Abel's functional equation (see Houzel and Wikipedia) and is related to Schröder's functional equation and Koenigs functions for compositional iterates (Alexander, Goryainov and Kudryavtseva). - _Tom Copeland_, Apr 04 2012

%C g(W(s,x)) = F'(s + FI(x)) = dW(s,x)/ds = g(x) dW(s,x)/dx, connecting the operators here to presentations of the Koenigs / Königs function and Loewner / Löwner evolution equations of the Contreras et al. papers. - _Tom Copeland_, Jun 03 2018

%C The autonomous differential equation above also appears with a change in variable of the form x = log(u) in the renormalization group equation, or Beta function. See Wikipedia, Zinn-Justin equations 2.10 and 3.11, and Krajewski and Martinetti equation 21. - _Tom Copeland_, Jul 23 2020

%C A variant of these partition polynomials appears on p. 83 of Petreolle et al. with the indeterminates e_n there related to those given in the examples below by e_n = n!*(n'). The coefficients are interpreted as enumerating certain types of trees. See also A190015. - _Tom Copeland_, Oct 03 2022

%D D. S. Alexander, A History of Complex Dynamics: From Schröder to Fatou to Julia, Friedrich Vieweg & Sohn, 1994.

%D T. Mansour and M. Schork, Commutation Relations, Normal Ordering, and Stirling Numbers, Chapman and Hall/CRC, 2015.

%H Luc Rousseau, <a href="/A145271/b145271.txt">Table of n, a(n) for n = 0..9295</a> (rows 0 to 25, flattened).

%H M. Abramowitz and I. A. Stegun, eds., <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of Mathematical Functions</a>, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972.

%H V. E. Adler, <a href="http://arxiv.org/abs/1510.02900">Set partitions and integrable hierarchies</a>, arXiv:1510.02900 [nlin.SI], 2015.

%H F. Bracci, M. Contreras, and S. Díaz-Madrigal, <a href="https://arxiv.org/abs/1804.10465">On the Königs function of semigroups of holomorphic self-maps of the unit disc</a>, arXiv:1804.10465 [math.CV], p. 15, 2018.

%H F. Bracci, M. Contreras, S. Díaz-Madrigal, and A. Vasil'ev, <a href="https://arxiv.org/abs/1309.6423">Classical and stochastic Löwner-Kufarev equations </a>, arXiv:1309.6423 [math.CV], (cf., e.g., p. 23), 2013.

%H C. Brouder, <a href="http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.180.7535&amp;rep=rep1&amp;type=pdf">Trees, renormalization, and differential equations</a>, BIT Numerical Mathematics, 44: 425-438, 2004, (Flow / autonomous differential equation, p. 429).

%H M. Contreras, S. Díaz-Madrigal, and P. Gumenyuk, <a href="https://arxiv.org/abs/0902.3116">Loewner chains in the unit disk</a>, arXiv:0902.3116 [math.CV], p. 29, 2009.

%H Tom Copeland, <a href="https://tcjpn.wordpress.com/2015/10/12/the-elliptic-lie-triad-kdv-and-ricattt-equations-infinigens-and-elliptic-genera/">The Elliptic Lie Triad: KdV and Ricatti Equations, Infinigens, and Elliptic Genera</a>

%H Tom Copeland, <a href="https://tcjpn.wordpress.com/2008/10/09/flipping-functions-with-permutohedra/">Flipping Functions with Permutohedra</a>, Posted Oct 2008.

%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/2010/12/28/14/">Addendum to Mathemagical Forests</a>, 2010.

%H Tom Copeland, <a href="https://mathoverflow.net/questions/214927/important-formulas-in-combinatorics/215203#215203">Important formulas in combinatorics</a>, MathOverflow answer, 2015.

%H Tom Copeland, <a href="https://tcjpn.wordpress.com/2018/01/23/formal-group-laws-and-binomial-sheffer-sequences/">Formal group laws and binomial Sheffer sequences</a>, 2018.

%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 Tom Copeland, <a href="https://tcjpn.wordpress.com/2022/10/08/compilation-of-oeis-partition-polynomials-a133314-a134685-a145271-a356144-and-a356145/">Compilation of OEIS Partition Polynomials A133314, A134685, A145271, A356144, and A356145</a>, 2022.

%H P. Feinsilver, <a href="http://chanoir.math.siu.edu/MATH/Merida/PDF/Merida.pdf">Lie algebras, representations, and analytic semigroups through dual vector fields</a>, Cimpa-Unesco-Venezuela School, Mérida, Venezuela, Jan-Feb 2006.

%H P. Feinsilver, R. Schott, <a href="http://dx.doi.org/10.1007/BF01046735">Appell systems on Lie groups</a>, J. Theor. Probab. 5 (2) (1992) 251-281.

%H P. Flajolet and R. Sedgewick, <a href="http://algo.inria.fr/flajolet/Publications/books.html">Analytic Combinatorics</a>, 2009, p. 526 and 527.

%H V. Goryainov and O. Kudryavtseva, <a href="http://dx.doi.org/10.1070/SM2011v202n07ABEH004173">One-parameter semigroups of analytic functions, fixed points and the Koenigs function</a>, Sbornik: Mathematics, 202:7 (2011), 971-1000.

%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 J. Hoffman, <a href="http://www.math.lsu.edu/~hoffman/papers/elmod.pdf">Topics in Elliptic Curves and Modular Forms</a>, p. 10.

%H C. Houzel, <a href="http://www.abelprisen.no/c53052/binfil/download.php?tid=53200">The Work of Niels Henrik Abel</a>, The Legacy of Niels Henrik Abel-The Abel Bicentennial, Oslo 2002 (Editors O. Laudal and R Piene), Springer-Verlag (2004), pp. 24-25.

%H T. Krajewski and T. Martinetti, <a href="https://arxiv.org/abs/0806.4309">Wilsonian renormalization, differential equations and Hopf algebras</a>, arXiv:0806.4309 [hep-th], 2008.

%H Peter Luschny, <a href="/A145271/a145271.pdf">Expansion A145271</a> (added Jul 21 2016)

%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/291888/characterizing-positivity-of-formal-group-laws">Characterizing positivity of formal group laws</a>, a MO question posed by Jair Taylor, 2018.

%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/415617/a-leibniz-like-formula-for-fx-fracddxn-fx">A Leibniz-like formula for (f(x)D_x)^n f(x)?</a>, a MO question posed by the user M.G. and answered by Tom Copeland, 2022.

%H MathStackExchange, <a href="http://math.stackexchange.com/questions/859883/closed-form-for-sequence-a145271">Closed form for sequence A145271</a>, posed 2014, response by Tom Copeland in 2016.

%H M. Mendez, <a href="https://arxiv.org/abs/1610.03602">Combinatorial differential operators in: Faà di Bruno formula, enumeration of ballot paths, enriched rooted trees and increasing rooted trees</a>, arXiv:1610.03602 [math.CO], p. 28 Example 6, 2016.

%H M. Petreolle, A. Sokal, and B. Zhu, <a href="https://arxiv.org/abs/1807.03271">Lattice paths and branched continued fractions: An infinite sequence of generalizations of the Stieltjes--Rogers and Thron--Rogers polynomials, with coefficientwise Hankel-total positivity</a>, arXiv:1807.03271 [math.CO], 2020.

%H J. Taylor, <a href="https://digital.lib.washington.edu/researchworks/handle/1773/36757">Formal group laws and hypergraph colorings</a>, doctoral thesis, Univ. of Wash., 2016, p. 66, eqn. 9.3.

%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Abel_equation">Abel equation</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Renormalization_group#Exact_renormalization_group_equations">Renormalization Group</a>

%H B. Zhu, <a href="https://arxiv.org/abs/2202.03793">Coefficientwise Hankel-total positivity of row-generating polynomials for the m-Jacobi-Rogers triangle</a>, arXiv:2202.03793 [math.CO], 2022.

%H J. Zinn-Justin, <a href="http://www.bourbaphy.fr/Zinn.pdf">Phase Transitions and Renormalization Group: from Theory to Numbers</a>, Séminaire Poincaré 2, pp. 55-74, 2002.

%F Let R = g(x)d/dx; then

%F R^0 g(x) = 1 (0')^1

%F R^1 g(x) = 1 (0')^1 (1')^1

%F R^2 g(x) = 1 (0')^1 (1')^2 + 1 (0')^2 (2')^1

%F R^3 g(x) = 1 (0')^1 (1')^3 + 4 (0')^2 (1')^1 (2')^1 + 1 (0')^3 (3')^1

%F R^4 g(x) = 1 (0')^1 (1')^4 + 11 (0')^2 (1')^2 (2')^1 + 4 (0')^3 (2')^2 + 7 (0')^3 (1')^1 (3')^1 + 1 (0')^4 (4')^1

%F R^5 g(x) = 1 (0') (1')^5 + 26 (0')^2 (1')^3 (2') + (0')^3 [34 (1') (2')^2 + 32 (1')^2 (3')] + (0')^4 [ 15 (2') (3') + 11 (1') (4')] + (0')^5 (5')

%F R^6 g(x) = 1 (0') (1')^6 + 57 (0')^2 (1')^4 (2') + (0')^3 [180 (1')^2 (2')^2 + 122 (1')^3 (3')] + (0')^4 [ 34 (2')^3 + 192 (1') (2') (3') + 76 (1')^2 (4')] + (0')^5 [15 (3')^2 + 26 (2') (4') + 16 (1') (5')] + (0')^6 (6')

%F where (j')^k = ((d/dx)^j g(x))^k. And R^(n-1) g(x) evaluated at x=0 is the n-th Taylor series coefficient of the compositional inverse of h(x) = (d/dx)^(-1) 1/g(x), with the integral from 0 to x.

%F The partitions are in reverse order to those in Abramowitz and Stegun p. 831. Summing over coefficients with like powers of (0') gives A008292.

%F Confer A190015 for another way to compute numbers for the array for each partition. - _Tom Copeland_, Oct 17 2014

%F Equivalent matrix computation: Multiply the n-th diagonal (with n=0 the main diagonal) of the lower triangular Pascal matrix by g_n = (d/dx)^n g(x) to obtain the matrix VP with VP(n,k) = binomial(n,k) g_(n-k). Then R^n g(x) = (1, 0, 0, 0, ...) [VP * S]^n (g_0, g_1, g_2, ...)^T, where S is the shift matrix A129185, representing differentiation in the divided powers basis x^n/n!. - _Tom Copeland_, Feb 10 2016 (An evaluation removed by author on Jul 19 2016. Cf. A139605 and A134685.)

%F Also, R^n g(x) = (1, 0, 0, 0, ...) [VP * S]^(n+1) (0, 1, 0, ...)^T in agreement with A139605. - _Tom Copeland_, Jul 21 2016

%F A recursion relation for computing each partition polynomial of this entry from the lower order polynomials and the coefficients of the cycle index polynomials of A036039 is presented in the blog entry "Formal group laws and binomial Sheffer sequences". - _Tom Copeland_, Feb 06 2018

%F A formula for computing the polynomials of each row of this matrix is presented as T_{n,1} on p. 196 of the Ihara reference in A139605. - _Tom Copeland_, Mar 25 2020

%F Indeterminate substitutions as illustrated in A356145 lead to [E] = [L][P] = [P][E]^(-1)[P] = [P][RT] and [E]^(-1) = [P][L] = [P][E][P] = [RT][P], where [E] contains the refined Eulerian partition polynomials of this entry; [E]^(-1), A356145, the inverse set to [E]; [P], the permutahedra polynomials of A133314; [L], the classic Lagrange inversion polynomials of A134685; and [RT], the reciprocal tangent polynomials of A356144. Since [L]^2 = [P]^2 = [RT]^2 = [I], the substitutional identity, [L] = [E][P] = [P][E]^(-1) = [RT][P], [RT] = [E]^(-1)[P] = [P][L][P] = [P][E], and [P] = [L][E] = [E][RT] = [E]^(-1)[L] = [RT][E]^(-1). - _Tom Copeland_, Oct 05 2022

%e From _Tom Copeland_, Sep 19 2014: (Start)

%e Let h(x) = log((1+a*x)/(1+b*x))/(a-b); then, g(x) = 1/(dh(x)/dx) = (1+ax)(1+bx), so (0')=1, (1')=a+b, (2')=2ab, evaluated at x=0, and higher order derivatives of g(x) vanish. Therefore, evaluated at x=0,

%e R^0 g(x) = 1

%e R^1 g(x) = a+b

%e R^2 g(x) = (a+b)^2 + 2ab = a^2 + 4 ab + b^2

%e R^3 g(x) = (a+b)^3 + 4*(a+b)*2ab = a^3 + 11 a^2*b + 11 ab^2 + b^3

%e R^4 g(x) = (a+b)^4 + 11*(a+b)^2*2ab + 4*(2ab)^2

%e = a^4 + 26 a^3*b + 66 a^2*b^2 + 26 ab^3 + b^4,

%e etc., and these bivariate Eulerian polynomials (A008292) are the first few coefficients of h^(-1)(x) = (e^(ax) - e^(bx))/(a*e^(bx) - b*e^(ax)), the inverse of h(x). (End)

%e Triangle starts:

%e 1;

%e 1;

%e 1, 1;

%e 1, 4, 1;

%e 1, 11, 4, 7, 1;

%e 1, 26, 34, 32, 15, 11, 1;

%e 1, 57, 180, 122, 34, 192, 76, 15, 26, 16, 1;

%e 1, 120, 768, 423, 496, 1494, 426, 294, 267, 474, 156, 56, 42, 22, 1;

%e 1, 247, 2904, 1389, 4288, 9204, 2127, 496, 5946, 2829, 5142, 1206, 855, 768, 1344, 1038, 288, 56, 98, 64, 29, 1;

%p with(LinearAlgebra): with(ListTools):

%p A145271_row := proc(n) local b, M, V, U, G, R, T;

%p if n < 2 then return 1 fi;

%p b := (n,k) -> `if`(k=1 or k>n+1,0,binomial(n-1,k-2)*g[n-k+1]);

%p M := n -> Matrix(n, b):

%p V := n -> Vector[row]([1, seq(0,i=2..n)]):

%p U := n -> VectorMatrixMultiply(V(n), M(n)^(n-1)):

%p G := n -> Vector([seq(g[i], i=0..n-1)]);

%p R := n -> VectorMatrixMultiply(U(n), G(n)):

%p T := Reverse([op(sort(expand(R(n+1))))]);

%p seq(subs({seq(g[i]=1, i=0..n)},T[j]),j=1..nops(T)) end:

%p for n from 0 to 9 do A145271_row(n) od; # _Peter Luschny_, Jul 20 2016

%Y Cf. (A133437, A086810, A181289) = (LIF, reduced LIF, associated g(x)), where LIF is a Lagrange inversion formula. Similarly for (A134264, A001263, A119900), (A134685, A134991, A019538), (A133932, A111999, A007318).

%Y Cf. A008292, A049019, A074060, A129185, A139002, A139605, A190015.

%Y Second column is A000295, subdiagonal is A000124, row sums are A000142, row lengths are A000041. - _Peter Luschny_, Jul 21 2016

%Y Cf. A036039, A133314, A356144, A356145.

%K easy,nonn,tabf,look

%O 0,6

%A _Tom Copeland_, Oct 06 2008

%E Title amplified by _Tom Copeland_, Mar 17 2014

%E R^5 and R^6 formulas and terms a(19)-a(29) added by _Tom Copeland_, Jul 11 2016

%E More terms from _Peter Luschny_, Jul 20 2016

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 12:26 EDT 2024. Contains 371254 sequences. (Running on oeis4.)