cp's OEIS Frontend

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)

This is the triangle of associated Stirling numbers of the second kind, A008299, read along the diagonals.

This is also a row-reversed version of A181996 (with an additional leading 1) - see the table on p. 92 in the Ward reference. A134685 is a refinement of the Ward table.

The first and second diagonals are A001147 and A000457 and appear in the diagonals of several OEIS entries. The polynomials also appear in Carlitz (p. 85), Drake et al. (p. 8) and Smiley (p. 7).

First few polynomials (with a different offset) are

P(0,t) = 0

P(1,t) = 1

P(2,t) = t

P(3,t) = t + 3*t^2

P(4,t) = t + 10*t^2 + 15*t^3

P(5,t) = t + 25*t^2 + 105*t^3 + 105*t^4

These are the "face" numbers of the tropical Grassmannian G(2,n),related to phylogenetic trees (with offset 0 beginning with P(2,t)). Corresponding h-vectors are A008517. - Tom Copeland, Oct 03 2011

A133314 applied to the derivative of A(x,t) implies (a.+b.)^n = 0^n, for (b_n)=P(n+1,t) and (a_0)=1, (a_1)=-t, and (a_n)=-(1+t) P(n,t) otherwise. E.g., umbrally, (a.+b.)^2 = a_2*b_0 + 2 a_1*b_1 + a_0*b_2 = 0. - Tom Copeland, Oct 08 2011

Beginning with the second column, the rows give the faces of the Whitehouse simplicial complex with the fourth-order complex being three isolated vertices and the fifth-order being the Petersen graph with 10 vertices and 15 edges (cf. Readdy). - Tom Copeland, Oct 03 2014

Stratifications of smooth projective varieties which are fine moduli spaces for stable n-pointed rational curves. Cf. pages 20 and 30 of the Kock and Vainsencher reference and references in A134685. - Tom Copeland, May 18 2017

Named after the American mathematician Morgan Ward (1901-1963). - Amiram Eldar, Jun 26 2021

Let b(n) = LPT[ A001147 ] = -A001147(n-1) for n > 0 and 1 for n=0, where LPT represents the action of the list partition transform described in A133314.

Then T(n,k) = binomial(n,k) * b(n-k) .

Form the matrix of polynomials TB(n,k,t) = T(n,k) * t^(n-k) = binomial(n,k) * b(n-k) * t^(n-k) = binomial(n,k) * Pb(n-k,t),

beginning as

1;

-1, 1;

-1*t, -2, 1;

-3*t^2, -3*t, -3, 1;

-15*t^3, -12*t^2, -6*t, -4, 1;

-105*t^4, -75*t^3, -30*t^2, -10*t, -5, 1;

Let Pc(n,t) = LPT(Pb(.,t)).

Then [TB(t)]^(-1) = TC(t) = [ binomial(n,k) * Pc(n-k,t) ] = LPT(TB),

whose first column is

Pc(0,t) = 1

Pc(1,t) = 1

Pc(2,t) = 2 + t

Pc(3,t) = 6 + 6*t + 3*t^2

Pc(4,t) = 24 + 36*t + 30*t^2 + 15*t^3

Pc(5,t) = 120 + 240*t + 270*t^2 + 210*t^3 + 105*t^4.

The coefficients of these polynomials are given by the reverse of A102625 with the highest order coefficients given by A001147 with an additional leading 1.

Note this is not the complete matrix TC. The complete matrix is formed by multiplying along the diagonal of the lower triangular Pascal matrix by these polynomials, embedding trees of coefficients in the matrix.

exp[Pb(.,t)*x] = 1 + [(1-2t*x)^(1/2) - 1] / (t-0) = [1 + a finite diff. of [(1-2t*x)^(1/2)] with step t] = e.g.f. of the first column of TB.

exp[Pc(.,t)*x] = 1 / { 1 + [(1-2t*x)^(1/2) - 1] / t } = 1 / exp[Pb(.,t)*x) = e.g.f. of the first column of TC.

TB(t) and TC(t), being inverse to each other, are the generators of an Abelian group.

TB(0) and TC(0) are generators for a subgroup representing the iterated Laguerre operator described in A132013 and A132014.

Let sb(t,m) and sc(t,m) be the associated sequences under the LPT to TB(t)^m = B(t,m) and TC(t)^m = C(t,m).

Let Esb(t,m) and Esc(t,m) be e.g.f.'s for sb(t,m) and sc(t,m), rB(t,m) and rC(t,m) be the row sums of B(t,m) and C(t,m) and aB(t,m) and aC(t,m) be the alternating row sums.

Then B(t,m) is the inverse of C(t,m), Esb(t,m) is the reciprocal of Esc(t,m) and sb(t,m) and sc(t,m) form a reciprocal pair under the LPT. Similar relations hold among the row sums and the alternating sign row sums and associated quantities.

All the group members have the form B(t,m) * C(u,p) = TB(t)^m * TC(u)^p = [ binomial(n,k) * s(n-k) ]

with associated e.g.f. Es(x) = exp[m * Pb(.,t) * x] * exp[p * Pc(.,u) * x] for the first column of the matrix, with terms s(n), so group multiplication is isomorphic to matrix multiplication and to multiplication of the e.g.f.'s for the associated sequences (see examples).

These results can be extended to other groups of integer-valued arrays by replacing the 2 by any natural number in the expression for exp[Pb(.,t)*x].

More generally,

[ G.f. for M = Product_{i=0..j} B[s(i),m(i)] * C[t(i),n(i)] ]

= exp(u*x) * Product_{i=0..j} { exp[m(i) * Pb(.,s(i)) * x] * exp[n(i) * Pc(.,t(i)) * x] }

= exp(u*x) * Product_{i=0..j} { 1 + [ (1 - 2*s(i)*x)^(1/2) - 1 ] / s(i) }^m(i) / { 1 + [ (1 - 2*t(i)*x)^(1/2) - 1 ] / t(i) }^n(i)

= exp(u*x) * H(x)

[ E.g.f. for M ] = I_o[2*(u*x)^(1/2)] * H(x).

M is an integer-valued matrix for m(i) and n(i) positive integers and s(i) and t(i) integers. To invert M, change B to C in Product for M.

H(x) is the e.g.f. for the first column of M and diagonally multiplying the Pascal matrix by the terms of this column generates M. See examples.

The G.f. for M, i.e., the e.g.f. for the row polynomials of M, implies that the row polynomials form an Appell sequence (see Wikipedia and Mathworld). - Tom Copeland, Dec 03 2013

Row sums have e.g.f. exp(x)*sech(x) (signed version of A009006). Inverse of masked Pascal triangle A119467. Transforms the sequence with e.g.f. g(x) to the sequence with e.g.f. g(x)*sech(x).

Coefficients of the Swiss-Knife polynomials for the computation of Euler, tangent and Bernoulli number (triangle read by rows). Another version in A153641. - Philippe Deléham, Oct 26 2013

Relations to Green functions and raising/creation and lowering/annihilation/destruction operators are presented in Hodges and Sukumar and in Copeland's discussion of this sequence and 2020 pdf. - Tom Copeland, Jul 24 2020

Connected objects from general (disconnected) objects.

The row lengths of this array is p(n):=A000041(n) (partition numbers).

In row n the partitions of n are taken in the Abramowitz-Stegun order.

One could call the unsigned numbers |a(n,k)| M_5 (similar to M_0, M_1, M_2, M_3 and M_4 given in A111786, A036038, A036039, A036040 and A117506, resp.).

The inverse relation (disconnected from connected objects) is found in A036040.

(d/da(1))p_n[a(1),a(2),...,a(n)] = n b_(n-1)[a(1),a(2),...,a(n-1)], where p_n are the partition polynomials of the cumulant generator A127671 and b_n are the partition polynomials for A133314. - Tom Copeland, Oct 13 2012

See notes on relation to Appell sequences in a differently ordered version of this array A263634. - Tom Copeland, Sep 13 2016

Given a binomial Sheffer polynomial sequence defined by the e.g.f. exp[t * f(x)] = Sum_{n >= 0} p_n(t) * x^n/n!, the cumulants formed from these polynomials are the Taylor series coefficients of f(x) multipied by t. An example is the sequence of the Stirling polynomials of the first kind A008275 with f(x) = log(1+x), so the n-th cumulant is (-1)^(n-1) * (n-1)! * t. - Tom Copeland, Jul 25 2019

From Tom Copeland, Oct 15 2020: (Start)

With a_n = n! * b_n = (n-1)! * c_n for n > 0, represent a function with f(0) = a_0 = b_0 = 1 as an

A) exponential generating function (e.g.f), or formal Taylor series: f(x) = e^{a.x} = 1 + Sum_{n > 0} a_n * x^n/n!

B) ordinary generating function (o.g.f.), or formal power series: f(x) = 1/(1-b.x) = 1 + Sum_{n > 0} b_n * x^n

C) logarithmic generating function (l.g.f): f(x) = 1 - log(1 - c.x) = 1 + Sum_{n > 0} c_n * x^n /n.

Expansions of log(f(x)) are given in

I) A127671 and A263634 for the e.g.f: log[ e^{a.*x} ] = e^{L.(a_1,a_2,...)x} = Sum_{n > 0} L_n(a_1,...,a_n) * x^n/n!, the logarithmic polynomials, cumulant expansion polynomials

II) A263916 for the o.g.f.: log[ 1/(1-b.x) ] = log[ 1 - F.(b_1,b_2,...)x ] = -Sum_{n > 0} F_n(b_1,...,b_n) * x^n/n, the Faber polynomials.

Expansions of exp(f(x)-1) are given in

III) A036040 for an e.g.f: exp[ e^{a.x} - 1 ] = e^{BELL.(a_1,...)x}, the Bell/Touchard/exponential partition polynomials, a.k.a. the Stirling partition polynomials of the second kind

IV) A130561 for an o.g.f.: exp[ b.x/(1-b.x) ] = e^{LAH.(b.,...)x}, the Lah partition polynomials

V) A036039 for an l.g.f.: exp[ -log(1-c.x) ] = e^{CIP.(c_1,...)x}, the cycle index polynomials of the symmetric groups S_n, a.k.a. the Stirling partition polynomials of the first kind.

Since exp and log are a compositional inverse pair, one can extract the indeterminates of the log set of partition polynomials from the exp set and vice versa. For a discussion of the relations among these polynomials and the combinatorics of connected and disconnected graphs/maps, see Novak and LaCroix on classical moments and cumulants and the two books on statistical mechanics referenced in A036040. (End)

Ignoring signs, these polynomials appear in Schröder in the set of equations (II) on p. 343 and in Stewart's translation on p. 31. - Tom Copeland, Aug 25 2021

This is a refinement of A019538 with row sums in A000670.

From Tom Copeland, Sep 29 2008: (Start)

This array is related to the reciprocal of an e.g.f. as sketched in A133314. For example, the coefficient of the fourth-order term in the Taylor series expansion of 1/(a(0) + a(1) x + a(2) x^2/2! + a(3) x^3/3! + ...) is a(0)^(-5) * {24 a(1)^4 - 36 a(1)^2 a(2) a(0) + [8 a(1) a(3) + 6 a(2)^2] a(0)^2 - a(4) a(0)^3}.

The unsigned coefficients characterize the P3 permutohedron depicted on page 10 in the Loday link with 24 vertices (0-D faces), 36 edges (1-D faces), 6 squares (2-D faces), 8 hexagons (2-D faces) and 1 3-D permutohedron. Summing coefficients over like dimensions gives A019538 and A090582. Compare to A133437 for the associahedron.

Given the n X n lower triangular matrix M = [ binomial(j,k) u(j-k) ], the first column of the inverse matrix M^(-1) contains the (n-1) rows of A049019 as the coefficients of the multinomials formed from the u(j). M^(-1) can be computed as (1/u(0)){I - [I- M/u(0)]^n} / {I - [I- M/u(0)]} = - u(0)^(-n) {sum(j=1 to n)(-1)^j bin(n,j) u(0)^(n-j) M^(j-1)} where I is the identity matrix.

Another method for computing the coefficients and partitions up to (n-1) rows is to use (1-x^n)/ (1-x) = 1+x^2+x^3+ ... + x^(n-1) with x replaced either by [I- M/a(0)] or [1- g(x)/a(0)] with the n X n matrix M = [bin(j,k) a(j-k)] and g(x)= a(0) + a(1)x + a(2)x^2/2! + ... + a(n) x^n/n!. The first n terms (rows of the first column) of the resulting series (matrix) divided by a(0) contain the (n-1) rows of signed coefficients and associated partitions for A049019.

To obtain unsigned coefficients, change a(j) to -a(j) for j>0. A133314 contains other matrices and recursion formulas that could be used. The Faa di Bruno formula gives the coefficients as n! [e(1)+e(2)+...+e(n)]! / [1!^e(1) e(1)! 2!^e(2) e(2)!... n!^e(n) e(n)! ] for the partition of form [a(1)^e(1)...a(n)^e(n)] with [e(1)+2e(2)+...+ n e(n)] = n (see Abramowitz and Stegun pages 823 and 831) in agreement with Arnold's formula. (End)

The matrix operation b = T*a can be characterized several ways in terms of the coefficients a(n) and b(n), their o.g.f.'s A(x) and B(x), or their e.g.f.'s EA(x) and EB(x).

1) b(n) = n! Lag[n,(.)!*Lag[.,a1(.),-1],0], umbrally,

where a1(n) = n! Lag[n,(.)!*Lag[.,a(.),-1],0]

2) b(n) = (-1)^n * n! * Lag(n,a(.),-2-n)

3) b(n) = Sum_{j=0..n} (-1)^j * binomial(n,j) * binomial(-2,j) * j! * a(n-j)

4) b(n) = Sum_{j=0..n} binomial(n,j) * (j+1)! * a(n-j)

5) B(x) = (1-xDx))^(-2) A(x), formally

6) B(x) = Sum_{j>=0} (-1)^j * binomial(-2,j) * (xDx)^j A(x)

= Sum_{j>=0} (j+1) * (xDx)^j A(x)

7) B(x) = Sum_{j>=0} (j+1) * x^j * D^j * x^j A(x)

8) B(x) = Sum_{j>=0} (j+1)! * x^j * Lag(j,-:xD:,0) A(x)

9) EB(x) = Sum_{j>=0} x^j * Lag[j,(.)! * Lag[.,a1(.),-1],0]

10) EB(x) = Sum_{j>=0} Lag[j,a1(.),-1] * (-x)^j / (1-x)^(j+1)

11) EB(x) = Sum_{j>=0} x^n * Sum_{j=0..n} (j+1)!/j! * a(n-j) / (n-j)!

12) EB(x) = Sum_{j>=0} (-x)^j * Lag[j,a(.),-2-j]

13) EB(x) = exp(a(.)*x) / (1-x)^2 = (1-x)^(-2) * EA(x)

14) T = A094587^2 = A132013^(-2) = A132014^(-1)

where Lag(n,x,m) are the Laguerre polynomials of order m, D the derivative w.r.t. x and (:xD:)^j = x^j * D^j. Truncating the D operator series at the j = n term gives an o.g.f. for b(0) through b(n).

c = (1!,2!,3!,4!,...) is the sequence associated to T under the list partition transform and associated operations described in A133314. Thus T(n,k) = binomial(n,k)*c(n-k) . c are also the coefficients in formulas 4 and 8.

The reciprocal sequence to c is d = (1,-2,2,0,0,0,...), so the inverse of T is TI(n,k) = binomial(n,k)*d(n-k) = A132014. (A121757 is the reverse of T.)

These formulas are easily generalized for m applications of the basic operator n! Lag[n,(.)!*Lag[.,a(.),-1],0] by replacing 2 by m in formulas 2, 3, 5, 6, 12, 13 and 14, or (j+1)! by (m-1+j)!/(m-1)! in 4, 8 and 11. For further discussion of repeated applications of T, see A132014.

The row sums of T = [formula 4 with a(n) all 1] = [binomial transform of c] = [coefficients of B(x) with A(x) = 1/(1-x)] = A001339. Therefore the e.g.f. of A001339 = [formula 13 with a(n) all 1] = exp(x)*(1-x)^(-2) = exp(x)*exp[c(.)*x)] = exp[(1+c(.))*x].

Note the reciprocal is 1/{exp[(1+c(.))*x]} = exp(-x)*(1-x)^2 = e.g.f. of signed A002061 with leading 1 removed], which makes A001339 and the signed, shifted A002061 reciprocal arrays under the list partition transform of A133314.

The e.g.f. for the row polynomials (see A132382) implies they form an Appell sequence (see Wikipedia). - Tom Copeland, Dec 03 2013

As noted in item 12 above and reiterated in the Bala formula below, the e.g.f. is e^(x*t)/(1-x)^2, and the Poisson-Charlier polynomials P_n(t,y) have the e.g.f. (1+x)^y e^(-xt) (Feinsilver, p. 5), so the row polynomials R_n(t) of this entry are (-1)^n P_n(t,-2). The associated Appell sequence IR_n(t) that is the umbral compositional inverse of this entry's polynomials has the e.g.f. (1-x)^2 e^(xt), i.e., the e.g.f. of A132014 (noted above), and, therefore, the row polynomials (-1)^n PC(t,2). As umbral compositional inverses, R_n(IR.(t)) = t^n = IR_n(R.(t)), where, by definition, P.(t)^n = P_n(t), is the umbral evaluation. - Tom Copeland, Jan 15 2016

T(n,k) is the number of ways to place (n-k) rooks in a 2 x (n-1) Ferrers board (or diagram) under the Goldman-Haglund i-row creation rook mode for i=2. Triangular recurrence relation is given by T(n,k) = T(n-1,k-1) + (n+1-k)*T(n-1,k). - Ken Joffaniel M. Gonzales, Jan 21 2016

"Standard order" here means as produced by Maple's "sort" command.

From Petros Hadjicostas, May 27 2020: (Start)

According to the Maple help files for the "sort" command, polynomials in multiple variables are "sorted in total degree with ties broken by lexicographic order (this is called graded lexicographic order)."

Thus for example, x_1^2*x_3 = x_1*x_1*x_3 > x_1*x_2*x_2 = x_1*x_2^2, while x_1^2*x_4 = x_1*x_1*x_4 > x_1*x_2*x_3. (End)

Row sums are 0 (for n > 1). Numbers of terms in rows are partition numbers A000041.

From Tom Copeland, Nov 06 2015: (Start)

With the formal Taylor series f(x) = 1 + x[1] x + x[2] x^2/2! + ... , the partition polynomials of this entry give d[log(f(x))]/dx = L_1(x[1]) + L_2(x[1], x[2]) x + L_3(...) x^2/2! + ..., and the coefficients of the reduced polynomials with x[n] = t are signed A028246.

The raising operator R = x + d[log(f(D)]/dD = x + L_1(x[1]) + L_2[x[1], x[2]) D + L_3(x[1], x[2], x[3]) D^2/2! + ... with D = d/dx generates an Appell sequence of polynomials, given umbrally by P_n(x[1], ..., x[n]; x) = (x[.] + x)^n = Sum_{k=0..n} binomial(n,k) x[k] * x^(n-k) = R^n 1 with the e.g.f. f(t)*e^(x*t) = exp[t P.(x[1], ..., x[.]; x)]. P_0 = x[0] = 1.

The umbral compositional inverse Appell sequence is generated by R = x - d[log(f(D))]/dD with e.g.f. e^(x*t)/f(t) = exp[t IP.(x[1], ..., x[.]; x)], so umbrally IP_n(x[1], ..., x[n]; P.(x[1], ..., x[n]; x)) = x^n = P_n(x[1], ..., x[n]; IP.(x[1], ..., x[n]; x)). An unsigned array for the reduced IP_n(x[1], ..., x[n]; x) polynomials with IP_0 = x[0] = 1 and x[n] = -1 for n > 0 is A154921, for which f(t) = 2 - e^t. (End)

From Tom Copeland, Sep 08 2016: (Start)

The Appell formalism allows a matrix representation in the power basis x^n of the raising operator R that incorporates this array's partition polynomials L_n(x[1], ..., x[n]):

VP_(n+1) = VP_n * R = VP_n * XPS^(-1) * MX * XPS, where XPS is the matrix formed from multiplying the n-th diagonal of the Pascal matrix PS of A007318 by the indeterminate x[n], with x[0] = 1 for the main diagonal of ones, i.e., XPS[n,k] = PS[n,k] * x[n-k]; the matrix MX is A129185; the matrix XPS^(-1) is the inverse of XPS, which can be formed by multiplying the diagonals of the Pascal matrix by the partition polynomials IPT(n, x[1], ..., x[n]) of A133314, i.e., XPS^(-1)[n,k] = PS[n,k] * IPT(n-k, x[1], ...); and VP_n is the row vector in the power basis representing the Appell polynomial P_n(x) formed from the basic sequence of moments 1, x[1], x[2], ..., i.e., umbrally P_n(x) = (x[.] + x)^n = Sum_{k=0..n} binomial(n,k) * x[k] * x^(n-k).

Then R = XPS^(-1) * MX * XPS is the Pascal matrix PS with an additional first superdiagonal of ones and the other lower diagonals multiplied by the partition polynomials of this array, i.e., R[n,k] = PS[n,k] * L_{n+1-k}(x[1], ..., x[n+1-k]) except for the first superdiagonal of ones.

Consistently, VP_n = (1, 0, 0, ...) * R^n = (1, 0, 0, ...) * XPS^(-1) * MX^n * XPS = (1, 0, 0, ...) * MX^n * XPS = the n-th row vector of XPS, which is the vector representation of P_n(x) = (x[.] + x)^n with x[0] = 1.

See the Copeland link for the umbral representation R = exp[g.*D] * x * exp[h.*D] that reflects the matrix representations.

The Stirling partition polynomials of the first kind St1_n(a[1], a[2], ..., a[n]) of A036039, the Stirling partition polynomials of the second kind St2_n(b[1], b[2], ..., b[n]) of A036040, and the refined Lah polynomials Lah_n[c[1], c[2], ..., c[n]) of A130561 are Appell sequences in the respective distinguished indeterminates a[1], b[1], and c[1]. Comparing the formulas for their raising operators with that in this entry, L_n(x[1], x[2], ..., x[n]) evaluates to

A) (n-1)! * a[n] for x[n] = St1_n(a[1], a[2], ..., a[n]);

B) b[n] for x[n] = St2_n(b[1], b[2], ..., b[n]);

C) n! * c[n] for x[n] = Lah_n(c[1], c[2], ..., c[n]).

Conversely, from the respective e.g.f.s (added Sep 12 2016)

D) x[n] = St1_n(L_1(x[1])/0!, ..., L_n(x[1], ..., x[n])/(n-1)!);

E) x[n] = St2_n(L_1(x[1]), ..., L_n(x[1], ..., x[n]));

F) x[n] = Lah_n(L_1(x[1])/1!, ..., L_n(x[1], ..., x[n])/n!).

Given only the Appell sequence with no closed form for the e.g.f., the raising operator can be generated using this formalism, as has been partially done for A134264. (End)

For the Appell sequences above, the raising operator is related to the recursion P_(n+1)(x) = x * P_n(x) + Sum_{k=0..n} binomial(n,k) * L_(n-k+1)(x[1], ..., x[n+k-1]) * P_k(x). For a derivation and connections to formal cumulants (c_n = L_n(x[1], ...)) and moments (m_n = x[n]), see the Copeland link on noncrossing partitions. With x = 0, the recursion reduces to x[n+1] = Sum_{k = 0..n} binomial(n,k) * L_(n-k+1)(x[1], ..., x[n+k-1]) * x[k] with x[0] = 1. This array is a differently ordered version of A127671. - Tom Copeland, Sep 13 2016

With x[n] = x^(n-1), a signed version of A130850 is obtained. - Tom Copeland, Nov 14 2016

See p. 2 of Getzler for a relation to stable graphs called necklaces used in computations for Deligne-Mumford-Knudsen moduli spaces of stable curves of genus 1. - Tom Copeland, Nov 15 2019

For a relation to a combinatorial Faa di Bruno Hopf algebra related to functional composition, as presented by Connes and Moscovici, see Figueroa et al. - Tom Copeland, Jan 17 2020

From Tom Copeland, May 17 2020: (Start)

The e.g.f. of an Appell sequence is f(t) e^(x*t) with f(0) = 1. Given the Laguerre-Polya class function f(t) = e^(-a*t^2 + b*t) Product_m (1 - t/z_m) e^(t/z_m) with a = 0 for simplicity (more generally a >= 0) and b real and where the product runs over all the zeros z_m of f(t) with all zeros real and Sum_m 1/(z_m)^2 convergent, the raising operator of the Appell polynomials is R = x + b - Sum_{k > 0} c_(k+1) D^k with c_k = Sum_m (1/(z_m)^k), i.e., traces of powers of the reciprocals of the zeros. From R in earlier comments, b = L_1(x_1) and otherwise c_k = -L_k(x_1, ..., x_k).

The Laguerre / Turan / de Gua inequalities (Csordas and Williamson and Skovgaard) imply that all the zeros of each Appell polynomial are real and simple and its extrema are local maxima above the x-axis and local minima below and are located above or below the zeros of the next lower degree Appell polynomial. (End)

From Tom Copeland, Oct 15 2020: (Start)

With a_n = n! * b_n = (n-1)! * c_n for n > 0, represent a function with f(0) = a_0 = b_0 = 1 as an

A) exponential generating function (e.g.f), or formal Taylor series: f(x) = e^{a.x} = 1 + Sum_{n > 0} a_n * x^n/n!

B) ordinary generating function (o.g.f.), or formal power series: f(x) = 1/(1-b.x) = 1 + Sum_{n > 0} b_n * x^n

C) logarithmic generating function (l.g.f): f(x) = 1 - log(1 - c.x) = 1 + Sum_{n > 0} c_n * x^n /n.

Expansions of log(f(x)) are given in

I) A127671 and A263634 for the e.g.f: log[ e^{a.*x} ] = e^{L.(a_1,a_2,...)x} = Sum_{n > 0} L_n(a_1,...,a_n) * x^n/n!, the logarithmic polynomials, cumulant expansion polynomials

II) A263916 for the o.g.f.: log[ 1/(1-b.x) ] = log[ 1 - F.(b_1,b_2,...)x ] = -Sum_{n > 0} F_n(b_1,...,b_n) * x^n/n, the Faber polynomials.

Expansions of exp(f(x)-1) are given in

III) A036040 for an e.g.f: exp[ e^{a.x} - 1 ] = e^{BELL.(a_1,...)x}, the Bell/Touchard/exponential partition polynomials, a.k.a. the Stirling partition polynomials of the second kind

IV) A130561 for an o.g.f.: exp[ b.x/(1-b.x) ] = e^{LAH.(b.,...)x}, the Lah partition polynomials

V) A036039 for an l.g.f.: exp[ -log(1-c.x) ] = e^{CIP.(c_1,...)x}, the cycle index polynomials of the symmetric groups S_n, a.k.a. the Stirling partition polynomials of the first kind.

Since exp and log are a compositional inverse pair, one can extract the indeterminates of the log set of partition polynomials from the exp set and vice versa. For a discussion of the relations among these polynomials and the combinatorics of connected and disconnected graphs/maps, see Novak and LaCroix on classical moments and cumulants and the two books on statistical mechanics referenced in A036040. (End)

Ignoring signs, these polynomials appear in Schröder in the set of equations (II) on p. 343 and in Stewart's translation on p. 31. - Tom Copeland, Aug 25 2021

The k-th diagonal of |A008275| appears as the k-th column in |A008276| with k-1 leading zeros.

The recurrence, given below, is derived from (d/dx)g1(k,x) - g1(k,x)= x*(d/dx)g1(k-1,x) + g1(k-1,x), k >= 1, with input g(-1,x):=0 and initial condition g1(k,0)=1, k >= 0. This differential recurrence for the e.g.f. g1(k,x) follows from the one for unsigned Stirling1 numbers.

The column sequences start with A000142 (factorials), A001705, A112487- A112491, for m=0,...,5.

The main diagonal gives (2*k-1)!! = A001147(k), k >= 1.

This computation was inspired by the Bender article (see links), where the Stirling polynomials are discussed.

The e.g.f. for the k-th diagonal, k >= 1, of the unsigned Stirling1 triangle |A008275| with k-1 leading zeros is g1(k-1,x) = exp(x)*Sum_{m=0..k-1} a(k,m)*(x^(k-1+m))/(k-1+m)!.

a(k,n) = number of lists with entries from [n] such that (i) each element of [n] occurs at least once and at most twice, (ii) for each i that occurs twice, all entries between the two occurrences of i are > i, and (iii) exactly k elements of [n] occur twice. Example: a(1,2)=5 counts 112, 121, 122, 211, 221, and a(2,2)=3 counts 1122,1221,2211. - David Callan, Nov 21 2011