cp's OEIS Frontend

The indexing used here follows that given in the classic books by Riordan and Comtet. For two other versions see A173018 and A123125. - N. J. A. Sloane, Nov 21 2010

Coefficients of Eulerian polynomials. Number of permutations of n objects with k-1 rises. Number of increasing rooted trees with n+1 nodes and k leaves.

T(n,k) = number of permutations of [n] with k runs. T(n,k) = number of permutations of [n] requiring k readings (see the Knuth reference). T(n,k) = number of permutations of [n] having k distinct entries in its inversion table. - Emeric Deutsch, Jun 09 2004

T(n,k) = number of ways to write the Coxeter element s_{e1}s_{e1-e2}s_{e2-e3}s_{e3-e4}...s_{e_{n-1}-e_n} of the reflection group of type B_n, using s_{e_k} and as few reflections of the form s_{e_i+e_j}, where i = 1, 2, ..., n and j is not equal to i, as possible. - Pramook Khungurn (pramook(AT)mit.edu), Jul 07 2004

Subtriangle for k>=1 and n>=1 of triangle A123125. - Philippe Deléham, Oct 22 2006

T(n,k)/n! also represents the n-dimensional volume of the portion of the n-dimensional hypercube cut by the (n-1)-dimensional hyperplanes x_1 + x_2 + ... x_n = k, x_1 + x_2 + ... x_n = k-1; or, equivalently, it represents the probability that the sum of n independent random variables with uniform distribution between 0 and 1 is between k-1 and k. - Stefano Zunino, Oct 25 2006

[E(.,t)/(1-t)]^n = n!*Lag[n,-P(.,t)/(1-t)] and [-P(.,t)/(1-t)]^n = n!*Lag[n, E(.,t)/(1-t)] umbrally comprise a combinatorial Laguerre transform pair, where E(n,t) are the Eulerian polynomials and P(n,t) are the polynomials in A131758. - Tom Copeland, Sep 30 2007

From Tom Copeland, Oct 07 2008: (Start)

G(x,t) = 1/(1 + (1-exp(x*t))/t) = 1 + 1*x + (2+t)*x^2/2! + (6+6*t+t^2)*x^3/3! + ... gives row polynomials for A090582, the reverse f-polynomials for the permutohedra (see A019538).

G(x,t-1) = 1 + 1*x + (1+t)*x^2/2! + (1+4*t+t^2)*x^3/3! + ... gives row polynomials for A008292, the h-polynomials for permutohedra (Postnikov et al.).

G((t+1)*x, -1/(t+1)) = 1 + (1+t)*x + (1+3*t+2*t^2)*x^2/2! + ... gives row polynomials for A028246.

(End)

A subexceedant function f on [n] is a map f:[n] -> [n] such that 1 <= f(i) <= i for all i, 1 <= i <= n. T(n,k) equals the number of subexceedant functions f of [n] such that the image of f has cardinality k [Mantaci & Rakotondrajao]. Example T(3,2) = 4: if we identify a subexceedant function f with the word f(1)f(2)...f(n) then the subexceedant functions on [3] are 111, 112, 113, 121, 122 and 123 and four of these functions have an image set of cardinality 2. - Peter Bala, Oct 21 2008

Further to the comments of Tom Copeland above, the n-th row of this triangle is the h-vector of the simplicial complex dual to a permutohedron of type A_(n-1). The corresponding f-vectors are the rows of A019538. For example, 1 + 4*x + x^2 = y^2 + 6*y + 6 and 1 + 11*x + 11*x^2 + x^3 = y^3 + 14*y^2 + 36*y + 24, where x = y + 1, give [1,6,6] and [1,14,36,24] as the third and fourth rows of A019538. The Hilbert transform of this triangle (see A145905 for the definition) is A047969. See A060187 for the triangle of Eulerian numbers of type B (the h-vectors of the simplicial complexes dual to permutohedra of type B). See A066094 for the array of h-vectors of type D. For tables of restricted Eulerian numbers see A144696 - A144699. - Peter Bala, Oct 26 2008

For a natural refinement of A008292 with connections to compositional inversion and iterated derivatives, see A145271. - Tom Copeland, Nov 06 2008

The polynomials E(z,n) = numerator(Sum_{k>=1} (-1)^(n+1)*k^n*z^(k-1)) for n >=1 lead directly to the triangle of Eulerian numbers. - Johannes W. Meijer, May 24 2009

From Walther Janous (walther.janous(AT)tirol.com), Nov 01 2009: (Start)

The (Eulerian) polynomials e(n,x) = Sum_{k=0..n-1} T(n,k+1)*x^k turn out to be also the numerators of the closed-form expressions of the infinite sums:

S(p,x) = Sum_{j>=0} (j+1)^p*x^j, that is

S(p,x) = e(p,x)/(1-x)^(p+1), whenever |x| < 1 and p is a positive integer.

(Note the inconsistent use of T(n,k) in the section listing the formula section. I adhere tacitly to the first one.) (End)

If n is an odd prime, then all numbers of the (n-2)-th and (n-1)-th rows are in the progression k*n+1. - Vladimir Shevelev, Jul 01 2011

The Eulerian triangle is an element of the formula for the r-th successive summation of Sum_{k=1..n} k^j which appears to be Sum_{k=1..n} T(j,k-1) * binomial(j-k+n+r, j+r). - Gary Detlefs, Nov 11 2011

Li and Wong show that T(n,k) counts the combinatorially inequivalent star polygons with n+1 vertices and sum of angles (2*k-n-1)*Pi. An equivalent formulation is: define the total sign change S(p) of a permutation p in the symmetric group S_n to be equal to Sum_{i=1..n} sign(p(i)-p(i+1)), where we take p(n+1) = p(1). T(n,k) gives the number of permutations q in S_(n+1) with q(1) = 1 and S(q) = 2*k-n-1. For example, T(3,2) = 4 since in S_4 the permutations (1243), (1324), (1342) and (1423) have total sign change 0. - Peter Bala, Dec 27 2011

Xiong, Hall and Tsao refer to Riordan and mention that a traditional Eulerian number A(n,k) is the number of permutations of (1,2...n) with k weak exceedances. - Susanne Wienand, Aug 25 2014

Connections to algebraic geometry/topology and characteristic classes are discussed in the Buchstaber and Bunkova, the Copeland, the Hirzebruch, the Lenart and Zainoulline, the Losev and Manin, and the Sheppeard links; to the Grassmannian, in the Copeland, the Farber and Postnikov, the Sheppeard, and the Williams links; and to compositional inversion and differential operators, in the Copeland and the Parker links. - Tom Copeland, Oct 20 2015

The bivariate e.g.f. noted in the formulas is related to multiplying edges in certain graphs discussed in the Aluffi-Marcolli link. See p. 42. - Tom Copeland, Dec 18 2016

Distribution of left children in treeshelves is given by a shift of the Eulerian numbers. Treeshelves are ordered binary (0-1-2) increasing trees where every child is connected to its parent by a left or a right link. See A278677, A278678 or A278679 for more definitions and examples. - Sergey Kirgizov, Dec 24 2016

The row polynomial P(n, x) = Sum_{k=1..n} T(n, k)*x^k appears in the numerator of the o.g.f. G(n, x) = Sum_{m>=0} S(n, m)*x^m with S(n, m) = Sum_{j=0..m} j^n for n >= 1 as G(n, x) = Sum_{k=1..n} P(n, x)/(1 - x)^(n+2) for n >= 0 (with 0^0=1). See also triangle A131689 with a Mar 31 2017 comment for a rewritten form. For the e.g.f see A028246 with a Mar 13 2017 comment. - Wolfdieter Lang, Mar 31 2017

For relations to Ehrhart polynomials, volumes of polytopes, polylogarithms, the Todd operator, and other special functions, polynomials, and sequences, see A131758 and the references therein. - Tom Copeland, Jun 20 2017

For relations to values of the Riemann zeta function at integral arguments, see A131758 and the Dupont reference. - Tom Copeland, Mar 19 2018

Normalized volumes of the hypersimplices, attributed to Laplace. (Cf. the De Loera et al. reference, p. 327.) - Tom Copeland, Jun 25 2018

This triangle is an unsigned version of the triangle of Stirling numbers of the first kind, A008275, which is the main entry for these numbers. - N. J. A. Sloane, Jan 25 2011

Or, triangle T(n,k), 0 <= k <= n, read by rows given by [1,1,2,2,3,3,4,4,5,5,6,6,...] DELTA [1,0,1,0,1,0,1,0,1,0,1,0,...] where DELTA is the operator defined in A084938.

Reversal of A094638.

Equals A132393*A007318, as infinite lower triangular matrices. - Philippe Deléham, Nov 13 2007

From Johannes W. Meijer, Oct 07 2009: (Start)

The higher order exponential integrals E(x,m,n) are defined in A163931. The asymptotic expansion of the exponential integrals E(x,m=1,n) ~ (exp(-x)/x)*(1 - n/x + n*(n+1)/x^2 - n*(n+1)*(n+2)/x^3 + ...), see Abramowitz and Stegun. This formula follows from the general formula for the asymptotic expansion, see A163932. We rewrite E(x,m=1,n) ~ (exp(-x)/x)*(1 - n/x + (n^2+n)/x^2 - (2*n+3*n^2+n^3)/x^3 + (6*n+11*n^2+6*n^3+n^4)/x^3 - ...) and observe that the T(n,m) are the polynomials coefficients in the denominators. Looking at the a(n,m) formula of A028421, A163932 and A163934, and shifting the offset given above to 1, we can write T(n-1,m-1) = a(n,m) = (-1)^(n+m)*Stirling1(n,m), see the Maple program.

The asymptotic expansion leads for values of n from one to eleven to known sequences, see the cross-references. With these sequences one can form the triangles A008279 (right-hand columns) and A094587 (left-hand columns).

See A163936 for information about the o.g.f.s. of the right-hand columns of this triangle.

(End)

The number of elements greater than i to the left of i in a permutation gives the i-th element of the inversion vector. (Skiena-Pemmaraju 2003, p. 69.) T(n,k) is the number of n-permutations that have exactly k 0's in their inversion vector. See evidence in Mathematica code below. - Geoffrey Critzer, May 07 2010

T(n,k) counts the rooted trees with k+1 trunks in forests of "naturally grown" rooted trees with n+2 nodes. This corresponds to sums of coefficients of iterated derivatives representing vectors, Lie derivatives, or infinitesimal generators for flow fields and formal group laws. Cf. links in A139605. - Tom Copeland, Mar 23 2014

A refinement is A036039. - Tom Copeland, Mar 30 2014

From Tom Copeland, Apr 05 2014: (Start)

With initial n=1 and row polynomials of T as p(n,x)=x(x+1)...(x+n-1), the powers of x correspond to the number of trunks of the rooted trees of the "naturally-grown" forest referred to above. With each trunk allowed m colors, p(n,m) gives the number of such non-plane colored trees for the forest with each tree having n+1 vertices.

p(2,m) = m + m^2 = A002378(m) = 2*A000217(m) = 2*(first subdiag of |A238363|).

p(3,m) = 2m + 3m^2 + m^3 = A007531(m+2) = 3*A007290(m+2) = 3*(second subdiag A238363).

p(4,m) = 6m + 11m^2 + 6m^3 + m^4 = A052762(m+3) = 4*A033487(m) = 4*(third subdiag).

From the Joni et al. link, p(n,m) also represents the disposition of n distinguishable flags on m distinguishable flagpoles.

The chromatic polynomial for the complete graph K_n is the falling factorial, which encodes the colorings of the n vertices of K_n and gives a shifted version of p(n,m).

E.g.f. for the row polynomials: (1-y)^(-x).

(End)

A relation to derivatives of the determinant |V(n)| of the n X n Vandermonde matrix V(n) in the indeterminates c(1) thru c(n):

|V(n)| = Product_{1<=jTom Copeland, Apr 10 2014

From Peter Bala, Jul 21 2014: (Start)

Let M denote the lower unit triangular array A094587 and for k = 0,1,2,... define M(k) to be the lower unit triangular block array

/I_k 0\

\ 0 M/

having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. Then the present triangle equals the infinite matrix product M(0)*M(1)*M(2)*... (which is clearly well defined). See the Example section. (End)

For the relation of this rising factorial to the moments of Viennot's Laguerre stories, see the Hetyei link, p. 4. - Tom Copeland, Oct 01 2015

Can also be seen as the Bell transform of n! without column 0 (and shifted enumeration). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 27 2016

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

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

From Tom Copeland, Oct 14 2011: (Start)

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:

W(s,x) = exp(s g(x)d/dx)x = F(s+FI(x)) ,

W(s,0) = F(s) ,

W(0,x) = x ,

dW(0,x)/ds = g(x) = F'[FI(x)] , implying

dW(0,F(x))/ds = g(F(x)) = F'(x) , and

W(s,W(r,x)) = F(s+FI(F(r+FI(x)))) = F(s+r+FI(x)) = W(s+r,x) . (See MF link below.) (End)

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) . - Tom Copeland, Oct 26 2011

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

Related to formal group laws for elliptic curves (see Hoffman). - Tom Copeland, Feb 24 2012

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

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

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

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

Coefficients are listed in Abramowitz and Stegun order (A036036).

Given an invertible function f(t) analytic about t=0 (or a formal power series) with f(0)=0 and Df(0) not equal to 0, form h(t) = t / f(t) and let h_n denote the coefficient of t^n in h(t).

Lagrange inversion gives the compositional inverse about t=0 as g(t) = Sum_{j>=1} ( t^j * (1/j) * Sum_{permutations s with s(1) + s(2) + ... + s(j) = j - 1} h_s(1) * h_s(2) * ... * h_s(j) ) = t * T(1,1) * h_0 + Sum_{j>=2} ( t^j * Sum_{k=1..(# of partitions for j-1)} T(j,k) * H(j-1,k ; h_0,h_1,...) ), where H(j-1,k ; h_0,h_1,...) is the k-th partition for h_1 through h_(j-1) corresponding to n=j-1 on page 831 of Abramowitz and Stegun (ordered as in A&S) with (h_0)^(j-m)=(h_0)^(n+1-m) appended to each partition subsumed under n and m of A&S.

Denoting h_n by (n') for brevity, to 8th order in t,

g(t) = t * (0')

+ t^2 * [ (0') (1') ]

+ t^3 * [ (0')^2 (2') + (0') (1')^2 ]

+ t^4 * [ (0')^3 (3') + 3 (0')^2 (1') (2') + (0') (1')^3 ]

+ t^5 * [ (0')^4 (4') + 4 (0')^3 (1') (3') + 2 (0')^3 (2')^2 + 6 (0')^2 (1')^2 (2') + (0') (1')^4 ]

+ t^6 * [ (0')^5 (5') + 5 (0')^4 (1') (4') + 5 (0')^4 (2') (3') + 10 (0')^3 (1')^2 (3') + 10 (0')^3 (1') (2')^2 + 10 (0')^2 (1')^3 (2') + (0') (1')^5 ]

+ t^7 * [ (0')^6 (6') + 6 (0')^5 (1') (5') + 6 (0')^5 (2') (4') + 3 (0')^5 (3')^2 + 15 (0')^4 (1')^2 (4') + 30 (0')^4 (1') (2') (3') + 5 (0')^4 (2')^3 + 20 (0')^3 (1')^3 (3') + 30 (0')^3 (1')^2 (2')^2 + 15 (0')^2 (1')^4 (2') + (0') (1')^6]

+ t^8 * [ (0')^7 (7') + 7 (0')^6 (1') (6') + 7 (0')^6 (2') (5') + 7 (0')^6 (3') (4') + 21 (0')^5 (1')^2* (5') + 42 (0')^5 (1') (2') (4') + 21 (0')^5 (1') (3')^2 + 21 (0')^5 (2')^2 (3') + 35 (0')^4 (1')^3 (4') + 105 (0)^4 (1')^2 (2') (3') + 35 (0')^4 (1') (2')^3 + 35 (0')^3 (1')^4 (3') + 70 (0')^3 (1')^3 (2')^2 + 21 (0')^2 (1')^5 (2') + (0') (1')^7 ]

+ ..., where from the formula section, for example, T(8,1',2',...,7') = 7! / ((8 - (1'+ 2' + ... + 7'))! * 1'! * 2'! * ... * 7'!) are the coefficients of the integer partitions (1')^1' (2')^2' ... (7')^7' in the t^8 term.

A125181 is an extended, reordered version of the above sequence, omitting the leading 1, with alternate interpretations.

If the coefficients of partitions with the same exponent for h_0 are summed within rows, A001263 is obtained, omitting the leading 1.

From identification of the elements of the inversion with those on page 25 of the Ardila et al. link, the coefficients of the irregular table enumerate non-crossing partitions on [n]. - Tom Copeland, Oct 13 2014

From Tom Copeland, Oct 28-29 2014: (Start)

Operating with d/d(1') = d/d(h_1) on the n-th partition polynomial Prt(n;h_0,h_1,..,h_n) in square brackets above associated with t^(n+1) generates n * Prt(n-1;h_0,h_1,..,h_(n-1)); therefore, the polynomials are an Appell sequence of polynomials in the indeterminate h_1 when h_0=1 (a special type of Sheffer sequence).

Consequently, umbrally, [Prt(.;1,x,h_2,..) + y]^n = Prt(n;1,x+y,h_2,..); that is, Sum_{k=0..n} binomial(n,k) * Prt(k;1,x,h_2,..) * y^(n-k) = Prt(n;1,x+y,h_2,..).

Or, e^(x*z) * exp[Prt(.;1,0,h_2,..) * z] = exp[Prt(.;1,x,h_2,..) * z]. Then with x = h_1 = -(1/2) * d^2[f(t)]/dt^2 evaluated at t=0, the formal Laplace transform from z to 1/t of this expression generates g(t), the comp. inverse of f(t), when h_0 = 1 = df(t)/dt eval. at t=0.

I.e., t / (1 - t*(x + Prt(.;1,0,h_2,..))) = t / (1 - t*Prt(.;1,x,h_2,..)) = g(t), interpreted umbrally, when h_0 = 1.

(End)

Connections to and between arrays associated to the Catalan (A000108 and A007317), Riordan (A005043), Fibonacci (A000045), and Fine (A000957) numbers and to lattice paths, e.g., the Motzkin, Dyck, and Łukasiewicz, can be made explicit by considering the inverse in x of the o.g.f. of A104597(x,-t), i.e., f(x) = P(Cinv(x),t-1) = Cinv(x) / (1 + (t-1)*Cinv(x)) = x*(1-x) / (1 + (t-1)*x*(1-x)) = (x-x^2) / (1 + (t-1)*(x-x^2)), where Cinv(x) = x*(1-x) is the inverse of C(x) = (1 - sqrt(1-4*x)) / 2, a shifted o.g.f. for the Catalan numbers, and P(x,t) = x / (1+t*x) with inverse Pinv(x,t) = -P(-x,t) = x / (1-t*x). Then h(x,t) = x / f(x,t) = x * (1+(t-1)Cinv(x)) / Cinv(x) = 1 + t*x + x^2 + x^3 + ..., i.e., h_1=t and all other coefficients are 1, so the inverse of f(x,t) in x, which is explicitly in closed form finv(x,t) = C(Pinv(x,t-1)), is given by A091867, whose coefficients are sums of the refined Narayana numbers above obtained by setting h_1=(1')=t in the partition polynomials and all other coefficients to one. The group generators C(x) and P(x,t) and their inverses allow associations to be easily made between these classic number arrays. - Tom Copeland, Nov 03 2014

From Tom Copeland, Nov 10 2014: (Start)

Inverting in x with t a parameter, let F(x;t,n) = x - t*x^(n+1). Then h(x) = x / F(x;t,n) = 1 / (1-t*x^n) = 1 + t*x^n + t^2*x^(2n) + t^3*x^(3n) + ..., so h_k vanishes unless k = m*n with m an integer in which case h_k = t^m.

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). [Added braces to disambiguate the formula. - N. J. A. Sloane, Oct 20 2015]

This relation reveals properties of the partitions and sums of the coefficients of the array. For n=1, h_k = t^k for all k, implying that the row sums are the Catalan numbers. For n = 2, h_k for k odd vanishes, implying that there are no blocks with only even-indexed h_k on the even-numbered rows and that only the blocks containing only even-sized bins contribute to the odd-row sums giving the Fuss-Catalan numbers for n=2. And so on, for n > 2.

These relations are reflected in any combinatorial structures enumerated by this array and the partitions, such as the noncrossing partitions depicted for a five-element set (a pentagon) in Wikipedia.

(End)

From Tom Copeland, Nov 12 2014: (Start)

An Appell sequence possesses an umbral inverse sequence (cf. A249548). The partition polynomials here, Prt(n;1,h_1,...), are an Appell sequence in the indeterminate h_1=u, so have an e.g.f. exp[Prt(.;1,u,h_2...)*t] = e^(u*t) * exp[Prt(.;1,0,h2,...)*t] with umbral inverses with an e.g.f e^(-u*t) / exp[Prt(.;1,0,h2,...)*t]. This makes contact with the formalism of A133314 (cf. also A049019 and A019538) and the signed, refined face partition polynomials of the permutahedra (or their duals), which determine the reciprocal of exp[Prt(.,0,u,h2...)*t] (cf. A249548) or exp[Prt(.;1,u,h2,...)*t], forming connections among the combinatorics of permutahedra and the noncrossing partitions, Dyck paths and trees (cf. A125181), and many other important structures isomorphic to the partitions of this entry, as well as to formal cumulants through A127671 and algebraic structures of Lie algebras. (Cf. relationship of permutahedra with the Eulerians A008292.)

(End)

From Tom Copeland, Nov 24 2014: (Start)

The n-th row multiplied by n gives the number of terms in the homogeneous symmetric monomials generated by [x(1) + x(2) + ... + x(n+1)]^n under the umbral mapping x(m)^j = h_j, for any m. E.g., [a + b + c]^2 = [a^2 + b^2 + c^2] + 2 * [a*b + a*c + b*c] is mapped to [3 * h_2] + 2 * [3 * h_1^2], and 3 * A134264(3) = 3 *(1,1)= (3,3) the number of summands in the two homogeneous polynomials in the square brackets. For n=3, [a + b + c + d]^3 = [a^3 + b^3 + ...] + 3 [a*b^2 + a*c^2 + ...] + 6 [a*b*c + a*c*d + ...] maps to [4 * h_3] + 3 [12 * h_1 * h_2] + 6 [4 * (h_1)^3], and the number of terms in the brackets is given by 4 * A134264(4) = 4 * (1,3,1) = (4,12,4).

The further reduced expression is 4 h_3 + 36 h_1 h_2 + 24 (h_1)^3 = A248120(4) with h_0 = 1. The general relation is n * A134264(n) = A248120(n) / A036038(n-1) where the arithmetic is performed on the coefficients of matching partitions in each row n.

Abramowitz and Stegun give combinatorial interpretations of A036038 and relations to other number arrays.

This can also be related to repeated umbral composition of Appell sequences and topology with the Bernoulli numbers playing a special role. See the Todd class link.

(End)

These partition polynomials are dubbed the Voiculescu polynomials on page 11 of the He and Jejjala link. - Tom Copeland, Jan 16 2015

See page 5 of the Josuat-Verges et al. reference for a refinement of these partition polynomials into a noncommutative version composed of nondecreasing parking functions. - Tom Copeland, Oct 05 2016

(Per Copeland's Oct 13 2014 comment.) The number of non-crossing set partitions whose block sizes are the parts of the n-th integer partition, where the ordering of integer partitions is first by total, then by length, then lexicographically by the reversed sequence of parts. - Gus Wiseman, Feb 15 2019

With h_0 = 1 and the other h_n replaced by suitably signed partition polynomials of A263633, the refined face partition polynomials for the associahedra of normalized A133437 with a shift in indices are obtained (cf. In the Realm of Shadows). - Tom Copeland, Sep 09 2019

Number of primitive parking functions associated to each partition of n. See Lemma 3.8 on p. 28 of Rattan. - Tom Copeland, Sep 10 2019

With h_n = n + 1, the d_k (A006013) of Table 2, p. 18, of Jong et al. are obtained, counting the n-point correlation functions in a quantum field theory. - Tom Copeland, Dec 25 2019

By inspection of the diagrams on Robert Dickau's website, one can see the relationship between the monomials of this entry and the connectivity of the line segments of the noncrossing partitions. - Tom Copeland, Dec 25 2019

Speicher has examples of the first four inversion partition polynomials on pp. 22 and 23 with his k_n equivalent to h_n = (n') here with h_0 = 1. Identifying z = t, C(z) = t/f(t) = h(t), and M(z) = f^(-1)(t)/t, then statement (3), on p. 43, of Theorem 3.26, C(z M(z)) = M(z), is equivalent to substituting f^(-1)(t) for t in t/f(t), and statement (4), M(z/C(z)) = C(z), to substituting f(t) for t in f^(-1)(t)/t. - Tom Copeland, Dec 08 2021

Given a Laurent series of the form f(z) = 1/z + h_1 + h_2 z + h_3 z^2 + ..., the compositional inverse is f^(-1)(z) = 1/z + Prt(1;1,h_1)/z^2 + Prt(2;1,h_1,h_2)/z^3 + ... = 1/z + h_1/z^2 + (h_1^2 + h_2)/z^3 + (h_1^3 + 3 h_1 h_2 + h_3)/z^4 + (h_1^4 + 6 h_1^2 h_2 + 4 h_1 h_3 + 2 h_2^2 + h_4)/z^5 + ... for which the polynomials in the numerators are the partition polynomials of this entry. For example, this formula applied to the q-expansion of Klein's j-invariant / function with coefficients A000521, related to monstrous moonshine, gives the compositional inverse with the coefficients A091406 (see He and Jejjala). - Tom Copeland, Dec 18 2021

The partition polynomials of A350499 'invert' the polynomials of this entry giving the indeterminates h_n. A multinomial formula for the coefficients of the partition polynomials of this entry, equivalent to the multinomial formula presented in the first four sentences of the formula section below, is presented in the MathOverflow question referenced in A350499. - Tom Copeland, Feb 19 2022

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!

...

See A134685 for more information.

A111785 is obtained from A133437 by dividing through the bracketed terms of the P(n,t) by n! and unsigned A111785 is a refinement of A033282 and A126216. - Tom Copeland, Sep 28 2008

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

From Tom Copeland, Mar 24 2020: (Start)

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)

Let f(t) = u(t) - u(0) = Ev[exp(u.* t) - u(0)] = log{Ev[(exp(z.* t))/z_0]} = Ev[-log(1- a.* t)], where the operator Ev denotes umbral evaluation of the umbral variables u., z. or a., e.g., Ev[a.^n + a.^m] = a_n + a_m . The relation between z_n and u_n is given in reference in A127671 and u_n = (n-1)! * a_n .

If u_1 is not equal to 0, then the compositional inverse for these expressions is given by g(t) = Sum_{j>=1} P(j,t) where, with u_n denoted by (n') for brevity,

P(1,t) = (1')^(-1) * [ 1 ] * t

P(2,t) = (1')^(-3) * [ -(2') ] * t^2 / 2!

P(3,t) = (1')^(-5) * [ 3 (2')^2 - (1')(3') ] * t^3 / 3!

P(4,t) = (1')^(-7) * [ -15 (2')^3 + 10 (1')(2')(3') - (1')^2 (4') ] * t^4 / 4!

P(5,t) = (1')^(-9) * [ 105 (2')^4 - 105 (1') (2')^2 (3') + 10 (1')^2 (3')^2 + 15 (1')^2 (2') (4') - (1')^3 (5') ] * t^5 / 5!

P(6,t) = (1')^(-11) * [ -945 (2')^5 + 1260 (1') (2')^3 (3') - 280 (1')^2 (2') (3')^2 - 210 (1')^2 (2')^2 (4') + 35 (1')^3 (3')(4') + 21 (1')^3 (2')(5') - (1')^4 (6') ] * t^6 / 6!

P(7,t) = (1')^(-13) * [ 10395 (2')^6 - 17325 (1') (2')^4 (3') + (1')^2 [ 6300 (2')^2 (3')^2 + 3150 (2')^3 (4')] - (1')^3 [280 (3')^3 + 1260 (2')(3')(4') + 378 (2')^2(5')] + (1')^4 [35 (4')^2 + 56 (3')(5') + 28 (2')(6')] - (1')^5 (7') ] * t^7 / 7!

P(8,t) = (1')^(-15) * [ -135135 (2')^7 + 270270 (1') (2')^5 (3') - (1')^2 [ 138600 (2')^3 (3')^2 + 51975 (2')^4 (4')] + (1')^3 [15400 (2')(3')^3 + 34650 (2')^2(3')(4') + 6930 (2')^3(5')] - (1')^4 [2100 (3')^2(4') + 1575 (2')(4')^2 + 2520 (2')(3')(5') + 630 (2')^2(6') ] + (1')^5 [126 (4')(5') + 84 (3')(6') + 36 (2')(7')] - (1')^6 (8') ] * t^8 / 8!

...

Substituting ((m-1)') for (m') in each partition and ignoring the (0') factors, the partitions in the brackets of P(n,t) become those of n-1 listed in Abramowitz and Stegun on page 831 (in the reversed order) and the number of partitions in P(n,t) is given by A000041(n-1).

Combinatorial interpretations are given in the link.

From Tom Copeland, Jul 10 2018: (Start)

Coefficients occurring in prolongation for the special Euclidean group SE(2) and special affine group SA(2) in the Olver presentation on moving frames (MFP) in slides 33 and 42. These are a result of applying an iterated derivative of the form h(x)d/dx = d/dy as in this entry (more generally as g(x) d/dx as discussed in A145271). See also p. 6 of Olver's paper on contact forms, but note that the 12 should be a 15 in the formula for the compositional inverse of S(t).

Change variables in the MFP to obtain connections to the partition polynomials Prt_n = n! * P(n,1) above. Let delta and beta in the formulas for the equi-affine curves in MFP be L and B, respectively, and D_y = (1/(L-B*u_x)) d/dx = (1/w_x) d/dx. Then v_(yy) = (1/B) [-w_(xx)/(w_x)^3] in MFP (there is an overall sign error in MFP for v_(yy) and higher derivatives w.r.t. y), and (d/dy)^n v = v_n = (1/B)* [(1/w_1)*(d/dx)]^(n-2) [-w_2/(w_1)^3] for n > 1, with w_n = (d/dx)^n w. Consequently, in the partition polynomials Prt_n for n > 1 here substitute (n') = -B*u_n = w_n for n > 1 and (1') = L-B*u_1 = w_1, where u_n = (d/dx)^n u, and then divide by B. For example, v_4 = (1/B)*Prt_4 = (1/B)*4!*P(4,1) = (1/B) (L-B*u_n)^(-7) [-15*(-B*u_2)^3 + 10 (L-B*u_1)(-B*u_2)(-B*u_3) - (L-B*u_1)^2 (-B*u_4)], agreeing with v_4 in MFP except for the overall sign.

For the SE(2) transformation formulas in MFP, let w_x = cos(phi) + sin(phi)*u_x, and then the same transformations apply as above with cos(phi) and sin(phi) substituted for L and -B, respectively. (End)

A permutation w has a double fall at k if w(k) > w(k+1) > w(k+2) and has an initial fall if w(1) > w(2).

exp(x*(1-y+y^2)*D_y)*f(y)|_{y=0} = f(1-E(-x)) for any function f with a Taylor series. D_y means differentiation with respect to y and E(x) is the e.g.f. given below. For a proof of exp(x*g(y)*D_y)*f(y) = f(F^{-1}(x+F(y))) with the compositional inverse F^{-1} of F(y)=int(1/g(y),y) with F(0)=0 see, e.g., the Datolli et al. reference.

Number of increasing ordered trees on vertex set {1,2,...,n}, rooted at 1, in which all outdegrees are <= 2. - David Callan, Mar 30 2007

Number of increasing colored 1-2 trees of order n with choice of two colors for the right branches of the vertices of outdegree 2. - Wenjin Woan, May 21 2011

Let f(t) = -log(1 - u(.)*t) = Sum_{n>=1} (u_n / 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) * [ -1 (2') ] * t^2 / 2!

P(3,t) = (1')^(-5) * [ 3 (2')^2 - 2 (1')(3') ] * t^3 / 3!

P(4,t) = (1')^(-7) * [ -15 (2')^3 + 20 (1')(2')(3') - 6 (1')^2 (4') ] * t^4 / 4!

P(5,t) = (1')^(-9) * [ 105 (2')^4 - 210 (1') (2')^2 (3') + 40 (1')^2 (3')^2 + 90 (1')^2 (2') (4') - 24 (1')^3 (5') ] * t^5 / 5!

P(6,t) = (1')^(-11) * [ -945 (2')^5 + 2520 (1') (2')^3 (3') - 1120 (1')^2 (2') (3')^2 - 1260 (1')^2 (2')^2 (4') + 420 (1')^3 (3')(4') + 504 (1')^3 (2')(5') - 120 (1')^4 (6') ] * t^6 / 6!

See A134685 for more information.

From Tom Copeland, Sep 28 2016: (Start)

P(7,t) = (1')^(-13) * [ 10395 (2')^6 - 34650 (1')(2')^4(3') + (1')^2 [25200 (2')^2(3')^2 + 18900 (2')^3(4')] - (1')^3 [2240 (3')^3 + 15120 (2')(3')(4') + 9072 (2')^2(5')] + (1')^4 [1260 (4')^2 + 2688 (3')(5') + 3360 (2')(6')] - 720 (1')^5(7')] * t^7 / 7!

P(8,t) = (1')^(-15) * [ -135135 (2')^7 + 540540 (1')(2')^5(3') - (1')^2 [554400 (2')^3(3')^2 + 311850 (2')^4(4')] + (1')^3 [123200 (2')(3')^3 + 415800 (2')^2(3')(4') + 166320 (2')^3(5')] - (1')^4 [50400 (3')^2(4') + 56700 (2')(4')^2 + 120960 (2')(3')(5') + 75600 (2')^2(6')] + (1')^5 [18144 (4')(5') + 20160 (3')(6') + 25920 (2')(7')] - 5040 (1')^6(8')] * t^8 / 8! (End)

Coefficients are listed in reverse graded colexicographic order (A228100). This is the reverse of Abramowitz and Stegun order (A036036).

Coefficients for Lagrange (compositional) inversion of a function in terms of the Taylor series expansion of its shifted reciprocal. Complementary to A134264 for formal power series and a scaled version of A248927. A refinement of A055302, which enumerates the number of labeled rooted trees with n nodes and k leaves, with row sums A000169.

Given an invertible function f(t) analytic about t=0 with f(0)=0 and df(0)/dt not 0, form h(t) = t / f(t) and denote h_n = (n') as the coefficient of t^n/n! in h(t). Then the compositional inverse of f(t), g(t), as a formal Taylor series, or e.g.f., is given up to the first few orders by

g(t) = [ 1 (0') ] * t

+ [ 2 (0') (1') ] * t^2/2!

+ [ 6 (0') (1')^2 + 3 (0')^2 (2') ] * t^3/3!

+ [24 (0') (1')^3 + 36 (0')^2 (1') (2') + 4 (0')^3 (3')] * t^4/4!

+ [120 (0') (1')^4 + 360 (0')^2 (1')^2 (2') + (0')^3 [60 (2')^2

+ 80 (1') (3')] + 5 (0')^4 (4')] * t^5/5!

+ [720 (0')(1')^5 + 3600 (0')^2 (1')^3(2') + (0')^3 [1800 (1')(2')^2 + 1200 ( 1')^2(3')] + (0')^4 [300 (2')(3') + 150 (1')(4')] + 6 (0')^5 (5')] * t^6/6! + ... .

Operating with [1/(n*(n-1))] d/d(1') = [1/(n*(n-1))] d/d(h_1) on the n-th partition polynomial in square brackets above associated with t^n/n! generates the (n-1)-th partition polynomial.

Each n-th partition polynomial here is n times the (n-1)-th partition polynomial of A248927.

From Tom Copeland, Nov 24 2014: (Start)

The n-th row is a mapping of the homogeneous symmetric monomials generated by [x(1) + x(2) + ... + x(n)]^(n-1) under the umbral mapping x(m)^j = h_j, for any m. E.g., [a + b + c]^2 = [a^2 + b^2 + c^2] + 2 * [a*b + a*c + b*c] is mapped to [3 * h_2] + 2 * [3 * h_1 * h_1] = 3 * h_2 + 6 * h_1^2 = A248120(3) with h_0 = 1. (Example corrected Jul 14 2015.)

For another example and relations to A134264 and A036038, see A134264. The general relation is n * A134264(n) = A248120(n) / A036038(n-1) where the arithmetic is performed on the coefficients of matching partitions in each row n.

The Abramowitz and Stegun reference in A036038 gives combinatorial interpretations of A036038 and relations to other number arrays.

This can also be related to repeated umbral composition of Appell sequences and topology with the Bernoulli numbers playing a special role. See the Todd class link. (End)

As presented above and in the Copeland link, this entry is related to exponentiation of e.g.f.s and, therefore, to discussions in the Scott and Sokal preprint (see eqn. 3.1 on p. 10 and eqn. 3.62 p. 24). - Tom Copeland, Jan 17 2017

Coefficients are listed in reverse graded colexicographic order (A228100). This is the reverse of Abramowitz and Stegun order (A036036).

Coefficients for Lagrange (compositional) inversion of a function in terms of the Taylor series expansion of its shifted reciprocal. Complementary to A134264 for formal power series. A refinement of A141618 with row sums A000272.

Given an invertible function f(t) analytic about t=0 with f(0)=0 and df(0)/dt not 0, form h(t) = t / f(t) and denote h_n = (n') as the coefficient of t^n/n! in h(t). Then the compositional inverse of f(t), g(t), as a formal Taylor series, or e.g.f., is given up to the first few orders by

g(t)/t = [ 1 (0') ]

+ [ 1 (0') (1') ] * t

+ [ 2 (0') (1')^2 + 1 (0')^2 (2') ] * t^2/2!

+ [ 6 (0') (1')^3 + 9 (0')^2 (1') (2') + 1 (0')^3 (3') ] * t^3/3!

+ [24 (0') (1')^4 + 72 (0')^2 (1')^2 (2') + (0')^3 [12 (2')^2

+ 16 (1') (3')] + (0')^4 (4')] * t^4/4!

+ [120 (0')(1')^5 + 600 (0')^2 (1')^3(2') + (0')^3 [300 (1')(2')^2 + 200 ( 1')^2(3')] + (0')^4 [50 (2')(3') + 25 (1')(4')] + (0')^5 (5')] * t^5/5! + [720 (0')(1')^6 + (0')^2 (1')^4(2')+(0')^3 [5400 (1')^2(2')^2 + 2400 (1')^3(3')] + (0')^4 [450 (2')^3+ 1800 (1')(2')(3') + 450( 1')^2(4')]+ (0')^5 [60 (3')^2 + 90 (2')(4') + 36 (1')(5')] + (0')^6 (6')] * t^6/6! + ...

..........

From Tom Copeland, Oct 28 2014: (Start)

Expressing g(t) as a Taylor series or formal e.g.f. in the indeterminates h_n generates a refinement of A055302, which enumerates the number of labeled root trees with n nodes and k leaves, with row sum A000169.

Operating with (1/n^2) d/d(1') = (1/n^2) d/d(h_1) on the n-th partition polynomial in square brackets above associated with t^n/n! generates the (n-1)-th partition polynomial.

Multiplying the n-th partition polynomial here by (n + 1) gives the (n + 1)-th partition polynomial of A248120. (End)

These are also the coefficients in the expansion of a series related to the Lagrange reversion theorem presented in Wikipedia of which the Lagrange inversion formula about the origin is a special case. Cf. Copeland link. - Tom Copeland, Nov 01 2016