|
|
A133437
|
|
Irregular triangle of coefficients of a partition transform for direct Lagrange inversion of an o.g.f., complementary to A134685 for an e.g.f.; normalized by the factorials, these are signed, refined face polynomials of the associahedra.
|
|
25
|
|
|
1, -2, 12, -6, -120, 120, -24, 1680, -2520, 360, 720, -120, -30240, 60480, -20160, -20160, 5040, 5040, -720, 665280, -1663200, 907200, 604800, -60480, -362880, -181440, 20160, 40320, 40320, -5040, -17297280, 51891840, -39916800, -19958400, 6652800, 19958400, 6652800, -1814400, -1814400, -3628800, -1814400, 362880, 362880, 362880, -40320
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Let f(t) = u(t) - u(0) = Sum_{n>=1} u_n * t^n.
If u_1 is not equal to 0, then the compositional inverse for f(t) is given by g(t) = Sum_{j>=1) P(n,t) where, with u_n denoted by (n'),
P(1,t) = (1')^(-1) * [ 1 ] * t
P(2,t) = (1')^(-3) * [ -2 (2') ] * t^2 / 2!
P(3,t) = (1')^(-5) * [ 12 (2')^2 - 6 (1')(3') ] * t^3 / 3!
P(4,t) = (1')^(-7) * [ -120 (2')^3 + 120 (1')(2')(3') - 24 (1')^2 (4') ] * t^4 / 4!
P(5,t) = (1')^(-9) * [ 1680 (2')^4 - 2520 (1') (2')^2 (3') + 360 (1')^2 (3')^2 + 720 (1')^2 (2') (4') - 120 (1')^3 (5') ] * t^5 / 5!
P(6,t) = (1')^(-11) * [ -30240 (2')^5 + 60480 (1') (2')^3 (3') - 20160 (1')^2 (2') (3')^2 - 20160 (1')^2 (2')^2 (4') + 5040 (1')^3 (3')(4') + 5040 (1')^3 (2')(5') - 720 (1')^4 (6') ] * t^6 / 6!
P(7,t) = (1')^(-13) * [ 665280 (2')^6 - 1663200 (1')(2')^4(3') + (1')^2 [907200 (2')^2(3')^2 + 604800 (2')^3(4')] - (1')^3 [60480 (3')^3 + 362880 (2')(3')(4') + 181440 (2')^2(5')] + (1')^4 [20160 (4')^2 + 40320 (3')(5') + 40320 (2')(6')] - 5040 (1')^5(7')] * t^7 / 7!
P(8,t) = (1')^(-15) * [ -17297280 (2')^7 + 51891840 (1')(2')^5(3') - (1')^2 [39916800 (2')^3(3')^2 + 19958400 (2')^4(4')] + (1')^3 [6652800 (2')(3')^3 + 19958400 (2')^2(3')(4') + 6652800 (2')^3(5')] - (1')^4 [1814400 (3')^2(4') + 1814400 (2')(4')^2 + 3628800 (2')(3')(5') + 1814400 (2')^2(6')] + (1')^5 [362880 (4')(5') + 362880 (3')(6') + 362880 (2')(7')] - 40320 (1')^6(8')] * t^8 / 8!
...
For relation to the geometry of associahedra or Stasheff polytopes (and other combinatorial objects) see the Loday and McCammond links. E.g., P(5,t) = (1')^(-9) * [ 14 (2')^4 - 21 (1') (2')^2 (3') + 6 (1')^2 (2') (4')+ 3 (1')^2 (3')^2 - 1 (1')^3 (5') ] * t^5 is related to the 3-D associahedron with 14 vertices (0-D faces), 21 edges (1-D faces), 6 pentagons (2-D faces), 3 rectangles (2-D faces), 1 3-D polytope (3-D faces). Summing over faces of the same dimension gives A033282 or A126216. - Tom Copeland, Sep 29 2008
A relation between this Lagrange inversion for an o.g.f. and partition polynomials formed from the "refined Lah numbers" A130561 is presented in the link "Lagrange a la Lah" along with umbral binary tree representations.
With f(x,t) = x + t*Sum_{n>=2} u_n*x^n, the compositional inverse in x is related to the velocity profile of particles governed by the inviscid Burgers's, or Hopf, eqn. See A001764 and A086810. - Tom Copeland, Feb 15 2014
Newton was aware of this power series expansion for series reversion. See the Ferraro reference p. 75 eqn. 52. - Tom Copeland, Jan 22 2017
The coefficients of the partition polynomials divided by the associated factorial enumerate the faces of the convex, bounded polytopes called the associahedra, and the absolute value of the sum of the renormalized coefficients gives the Euler characteristic of unity for each polytope; i.e., the absolute value of the sum of each row of the array is either n! (unnormalized) or unity (normalized). In addition, the signs of the faces alternate with dimension, and the coefficients of faces with the same dimension for each polytope have the same sign. - Tom Copeland, Nov 13 2019
With u_1 = 1 and the other u_n replaced by suitably signed partition polynomials of A263633, the partition polynomials enumerating the refined noncrossing partitions of A134264 with a shift in indices are obtained (cf. In the Realm of Shadows). - Tom Copeland, Nov 16 2019
Relations between associahedra and oriented n-simplices are presented by Halvorson and by Street. - Tom Copeland, Dec 08 2019
Let f(x;t,n) = x - t * x^(n+1), giving u_1 = (1') = 1 and u_(n+1) = (n+1) = -t. Then inverting in x with t a parameter gives finv(x;t,n) = Sum_{j>=0} {binomial((n+1)*j,j) / (n*j + 1)} * t^j * x^(n*j + 1), which gives the Catalan numbers for n=1, and the Fuss-Catalan sequences for n>1 (see A001764, n=2). Comparing this with the same result in A134264 gives relations between the faces of associahedra and noncrossing partitions (and other combinatorial constructs related to these inversion formulas and those listed in A145271). - Tom Copeland, Jan 27 2020
There is a mapping between the faces of K_n, the associahedron of dimension (n-1), and polygon dissections. The dissecting noncrossing diagonals (i.e., nonintersecting in the interior) form subpolygons. Assign the indeterminate x_k to a subpolygon where k = (number of vertices of the subpolygon) - 1. Multiply the x_k together to form the monomials for the inversion formula.
For the 3-dimensional associahedron K_4, the fundamental polygon is the hexagon, which can be dissected into pentagons, associated to x_4; tetragons, to x_3; and triangles, to x_2; for example, there are six distinguished partitions of the hexagon into one triangle and one pentagon, sharing two vertices, associated to the monomial 6 * x_2 * x_4 since the unshared vertex of the triangle can be moved consecutively from one vertex of the hexagon to the next. This term corresponds to 720 (1')^2 (2') (4') / 5! in P(5,t) above, denumerating the six pentagonal facets of K_4. (End)
|
|
REFERENCES
|
G. Ferraro, The Rise and Development of the Theory of Series up to the Early 1820s, Springer Science and Business Media, 2007.
H. Halvorson (editor), Deep Beauty: Understanding the Quantum World Through Innovation, Cambridge Univ. Press, 2011.
H. Turnbull (editor), The Correspondence of Isaac Newton Vol. II 1676-1687, Cambridge Univ. Press, 1960, p. 147.
|
|
LINKS
|
S. He, Scattering from Geometries, presentation at the Interdisciplinary Center for Theoretical Study at the University of Science and Technology of China, 2017.
D. Herceg, L. Petkovic, and M. Petkovic, On Schröder’s families of root-finding methods, Journal of Computational and Applied Mathematics, Vol. 233, Issue 8, 15 February 2010, Pages 1755-1762, (see p. 1760, below eqn. 19).
J. Stasheff, What is ... an Operad?, Notices of the American Mathematical Society 51 (6), June-July 2004, 630-631.
|
|
FORMULA
|
The bracketed partitions of P(n,t) are of the form (u_1)^e(1) (u_2)^e(2) ... (u_n)^e(n) with coefficients given by (-1)^(n-1+e(1)) * [2*(n-1)-e(1)]! / [ (e(2))! * (e(3))! * ... * (e(n))! ].
Let h(t) = 1/(df(t)/dt)
= 1/Ev[u./(1-u.t)^2]
= 1/((u_1) + 2*(u_2)*t + 3*(u_3)*t^2 + 4*(u_4)*t^3 + ...),
where Ev denotes umbral evaluation.
Then for the partition polynomials of A133437,
n!*P(n,t) = ((t*h(y)*d/dy)^n) y evaluated at y=0,
and the compositional inverse of f(t) is
g(t) = exp(t*h(y)*d/dy) y evaluated at y=0.
Also, dg(t)/dt = h(g(t)). (End)
With exp[x* PS(.,t)] = exp[t*g(x)] = exp[x*h(y)d/dy] exp(t*y) eval. at y=0, the raising/creation and lowering/annihilation operators defined by R PS(n,t)=PS(n+1,t) and L PS(n,t) = n*PS(n-1,t) are
R = t*h(d/dt) = t* 1/[(u_1) + 2*(u_2)*d/dt + 3*(u_3)*(d/dt)^2 + ...] and
L = f(d/dt) = (u_1)*d/dt + (u_2)*(d/dt)^2 + (u_3)*(d/dt)^3 + ....
Then P(n,t) = (t^n/n!) dPS(n,z)/dz eval. at z=0. (Cf. A139605, A145271, and link therein to Mathemagical Forests for relation to planted trees on p. 13.) (End)
The bracketed partition polynomials of P(n,t) are also given by (d/dx)^(n-1) 1/[u_1 + u_2 * x + u_3 * x^2 + ... + u_n * x^(n-1)]^n evaluated at x=0. - Tom Copeland, Jul 07 2015
Let PS(n,u1,u2,...,un) = P(n,t) / t^n, i.e., the square-bracketed part of the partition polynomials in the expansion for the inverse in the comment section, with u_k = uk.
Also let PS(n,u1=1,u2,...,un) = PB(n,b1,b2,...,bK,...) where each bK represents the partitions of PS, with u1 = 1, that have K components or blocks, e.g., PS(5,1,u2,...,u5) = PB(5,b1,b2,b3,b4) = b1 + b2 + b3 + b4 with b1 = -u5, b2 = 6 u2 u4 + 3 u3^2, b3 = -21 u2^2 u3, and b4 = 14 u2^4.
The relation between solutions of the inviscid Burgers' equation and compositional inverse pairs (cf. A086810) implies that, for n > 2, PB(n, 0 * b1, 1 * b2, ..., (K-1) * bK, ...) = [(n+1)/2] * Sum_{k = 2..n-1} PS(n-k+1,u_1=1,u_2,...,u_(n-k+1)) * PS(k,u_1=1,u_2,...,u_k).
For example, PB(5,0 * b1, 1 * b2, 2 * b3, 3 * b4) = 3 * 14 u2^4 - 2 * 21 u2^2 u3 + 1 * 6 u2 u4 + 1 * 3 u3^2 - 0 * u5 = 42 u2^4 - 42 u2^2 u3 + 6 u2 u4 + 3 u3^2 = 3 * [2 * PS(2,1,u2) * PS(4,1,u2,...,u4) + PS(3,1,u2,u3)^2] = 3 * [ 2 * (-u2) (-5 u2^3 + 5 u2 u3 - u4) + (2 u2^2 - u3)^2].
Also, PB(n,0*b1,1*b2,...,(K-1)*bK,...) = d/dt t^(n-2)*PS(n,u1=1/t,u2,...,un)|_{t=1} = d/dt (1/t)*PS(n,u1=1,t*u2,...,t*un)|_{t=1}.
(End)
Equivalent matrix computation: Multiply the m-th diagonal (with m=1 the index of the main diagonal) of the lower triangular Pascal matrix A007318 by f_m = m!*u_m = (d/dx)^m f(x) evaluated at x=0 to obtain the matrix UP with UP(n,k) = binomial(n,k) f_{n+1-k}, or equivalently multiply the diagonals of A132159 by u_m. Then P(n,t) = (1, 0, 0, 0, ...) [UP^(-1) * S]^(n-1) FC * t^n/n!, where S is the shift matrix A129185, representing differentiation in the basis x^n//n!, and FC is the first column of UP^(-1), the inverse matrix of UP. These results follow from A145271 and A133314.
Also, P(n,t) = (1, 0, 0, 0, ...) [UP^(-1) * S]^n (0, 1, 0, ...)^T * t^n/n! in agreement with A139605. (End)
A recursion relation for computing each partition polynomial of this entry from the lower order polynomials and the coefficients of the refined Lah polynomials of A130561 is presented in the blog entry "Formal group laws and binomial Sheffer sequences." - Tom Copeland, Feb 06 2018
The derivative of the partition polynomials of A350499 with respect to a distinguished indeterminate give polynomials proportional to those of this entry. The connection of this derivative relation to the inviscid Burgers-Hopf evolution equation is given in a reference for that entry. - Tom Copeland, Feb 19 2022
|
|
MATHEMATICA
|
rows[nn_] := {{1}}~Join~With[{s = InverseSeries[t (1 + Sum[u[k] t^k, {k, nn}] + O[t]^(nn+1))]}, Table[(n+1)! Coefficient[s, t^(n+1) Product[u[w], {w, p}]], {n, nn}, {p, Reverse[Sort[Sort /@ IntegerPartitions[n]]]}]];
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign,tabf,changed
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Missing coefficient in P(6,t) replaced by Tom Copeland, Nov 06 2008
Terms ordered according to the reversed Abramowitz-Stegun ordering of partitions (with every k' replaced by (k-1)') by Andrey Zabolotskiy, Mar 07 2024
|
|
STATUS
|
approved
|
|
|
|