cp's OEIS Frontend

a(n) = A000110(n) - A000108(n).

a(n) = Sum_{k=0..n} S2(n,k) - binomial(2*n,n)/(n+1); S2(n,k) Stirling numbers of the second kind.

E.g.f.: exp(exp(x)-1) - (BesselI(0,2*x) - BesselI(1,2*x))*exp(2*x). - Ilya Gutkovskiy, Aug 31 2016

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

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

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

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

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

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

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')

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')

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.

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

Confer A190015 for another way to compute numbers for the array for each partition. - Tom Copeland, Oct 17 2014

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.)

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

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

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

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

G.f. G = G(t, z) satisfies (1+t)*G^2 - z*(1-z-2*t*z)*G + t*z^4 = 0.

T(n, k) = binomial(n-3, k)*binomial(n+k-1, k)/(k+1) for n >= 3, 0 <= k <= n-3.

From Tom Copeland, Nov 03 2008: (Start)

Two g.f.s (f1 and f2) for A033282 and their inverses (x1 and x2) can be derived from the Drake and Barry references.

1. a: f1(x,t) = y = {1 - (2t+1) x - sqrt[1 - (2t+1) 2x + x^2]}/[2x (t+1)] = t x + (t + 2 t^2) x^2 + (t + 5 t^2 + 5 t^3) x^3 + ...

b: x1 = y/[t + (2t+1)y + (t+1)y^2] = y {1/[t/(t+1) + y] - 1/(1+y)} = (y/t) - (1+2t)(y/t)^2 + (1+ 3t + 3t^2)(y/t)^3 +...

2. a: f2(x,t) = y = {1 - x - sqrt[(1-x)^2 - 4xt]}/[2(t+1)] = (t/(t+1)) x + t x^2 + (t + 2 t^2) x^3 + (t + 5 t^2 + 5 t^3) x^4 + ...

b: x2 = y(t+1) [1- y(t+1)]/[t + y(t+1)] = (t+1) (y/t) - (t+1)^3 (y/t)^2 + (t+1)^4 (y/t)^3 + ...

c: y/x2(y,t) = [t/(t+1) + y] / [1- y(t+1)] = t/(t+1) + (1+t) y + (1+t)^2 y^2 + (1+t)^3 y^3 + ...

x2(y,t) can be used along with the Lagrange inversion for an o.g.f. (A133437) to generate A033282 and show that A133437 is a refinement of A033282, i.e., a refinement of the f-polynomials of the associahedra, the Stasheff polytopes.

y/x2(y,t) can be used along with the indirect Lagrange inversion (A134264) to generate A033282 and show that A134264 is a refinement of A001263, i.e., a refinement of the h-polynomials of the associahedra.

f1[x,t](t+1) gives a generator for A088617.

f1[xt,1/t](t+1) gives a generator for A060693, with inverse y/[1 + t + (2+t) y + y^2].

f1[x(t-1),1/(t-1)]t gives a generator for A001263, with inverse y/[t + (1+t) y + y^2].

The unsigned coefficients of x1(y t,t) are A074909, reverse rows of A135278. (End)

G.f.: 1/(1-x*y-(x+x*y)/(1-x*y/(1-(x+x*y)/(1-x*y/(1-(x+x*y)/(1-x*y/(1-.... (continued fraction). - Paul Barry, Feb 06 2009

Let h(t) = (1-t)^2/(1+(u-1)*(1-t)^2) = 1/(u + 2*t + 3*t^2 + 4*t^3 + ...), then a signed (n-1)-th row polynomial of A033282 is given by u^(2n-1)*(1/n!)*((h(t)*d/dt)^n) t, evaluated at t=0, with initial n=2. The power series expansion of h(t) is related to A181289 (cf. A086810). - Tom Copeland, Sep 06 2011

With a different offset, the row polynomials equal 1/(1 + x)*Integral_{0..x} R(n,t) dt, where R(n,t) = Sum_{k = 0..n} binomial(n,k)*binomial(n+k,k)*t^k are the row polynomials of A063007. - Peter Bala, Jun 23 2016

n-th row polynomial = ( LegendreP(n-1,2*x + 1) - LegendreP(n-3,2*x + 1) )/((4*n - 6)*x*(x + 1)), n >= 3. - Peter Bala, Feb 22 2017

n*T(n+1, k) = (4n-6)*T(n, k-1) + (2n-3)*T(n, k) - (n-3)*T(n-1, k) for n >= 4. - Fang Lixing, May 07 2019

T(n, k) = [binomial(n+1, k)/(n+1)]*Sum_{j=1..floor((n-k)/2)} binomial(n+1-k, j)*binomial(n-k-j-1, j-1) for kn. G.f.=G=G(t, z) satisfies z(1+z-tz)G^2-(1+z-tz)G+1=0. T(n, k)=r(n-k)*binomial(n, k), where r(n)=A005043(n) are the Riordan numbers.

G.f.: 1/(1-xy-x^2/(1-x-xy-x^2/(1-x-xy-x^2/(1-x-xy-x^2/(1-... (continued fraction). - Paul Barry, Aug 03 2009

Sum_{k=0..n} T(n,k)*x^k = A126930(n), A005043(n), A000108(n), A007317(n), A064613(n), A104455(n) for x = -1,0,1,2,3,4 respectively. - Philippe Deléham, Dec 03 2009

Sum_{k=0..n} (-1)^(n-k)*T(n,k)*x^k = A168491(n), A099323(n+1), A001405(n), A005773(n+1), A001700(n), A026378(n+1), A005573(n), A122898(n) for x = -1, 0, 1, 2, 3, 4, 5, 6 respectively. - Philippe Deléham, Dec 03 2009

E.g.f.: e^(x+xy)*(Bessel_I(0,2x)-Bessel_I(1,2x)). - Paul Barry, Mar 10 2010

From Tom Copeland, Nov 06 2014: (Start)

O.g.f.: H(x,t) = {1-sqrt[1-4x/(1-(t-1)x)]}/2 (shifted index, as given in Copeland's comment, see comp. inverse there).

H(x,t)= x / [1-(C.+(t-1))x] = Sum_{n>=1} (C.+ (t-1))^(n-1)*x^n umbrally, e.g., (a.+b.)^2 = a_0*b_2 + 2 a_1*b1_+ a_0*b_2, where (C.)^n = C_n are the Catalan numbers (1,1,2,5,14,..) of A000108.

This shows directly that the lowering operator for the polynomials is D=d/dt, i.e., D p(n,t)= D(C. + (t-1))^n = n * (C. + (t-1))^(n-1) = n*p(n-1,t), so that the polynomials form an Appell sequence, and that p(n,0) gives a Motzkin sum, or Riordan, number A005043.

(End)

T(n,k) = (-1)^(n+k)*binomial(n,k)*hypergeom([k-n,1/2],[2],4). - Peter Luschny, Jul 27 2016

E.g.f.: (y-1)*exp(x*y)/(y-exp((y-1)*x)). - Vladeta Jovovic, Sep 20 2003

p(t,x) = (1 - x)*exp(t)/(1 - x*exp(t*(1 - x))). - Roger L. Bagula, Nov 21 2009

With offset=0, T(n,0)=1 otherwise T(n,k) = sum_{i=0..k-1} C(n,i)((i-k)^i*(k-i+1)^(n-i) - (i-k+1)^i*(k-i)^(n-i)) (cf. Williams). - Tom Copeland, Oct 10 2014

With offset 0, T = A007318 * A123125. Second column is A000225; 3rd, appears to be A066810. - Tom Copeland, Jan 23 2015

A raising operator (with D = d/dx) associated with this entry's row polynomials is R = x + t + (1-t) / [1-t e^{(1-t)D}] = x + t + 1 + t D + (t+t^2) D^2/2! + (t+4t^2+t^3) D^3/3! + ... , containing the e.g.f. for the Eulerian polynomials of A123125. Then R^n 1 = (p.(0;t)+x)^n = p_n(x;t) are the Appell polynomials with this entry's row polynomials p_n(0;t) as the base sequence. Examples of this formalism are given in A028246 and A248727. - Tom Copeland, Jan 24 2015

With offset 0, T = A007318*(padded A090582)*(inverse of A097805) = A007318*(padded A090582)*(padded A130595) = A007318*A123125 = A007318*(padded A090582)*A007318*A097808*A130595, where padded matrices are of the form of padded A007318, which is A097805. Inverses of padded matrices are just the padded versions of inverses of the unpadded matrices. This relates the face vectors, or f-vectors, and h-vectors of the permutahedra / permutohedra to those of the stellahedra / stellohedra. - Tom Copeland, Nov 13 2016

Umbrally, the row polynomials (offset 0) are r_n(x) = (1 + q.(x))^n, where (q.(x))^k = q_k(x) are the row polynomials of A123125. - Tom Copeland, Nov 16 2016

From the previous umbral statement, OP(x,d/dy) y^n = (y + q.(x))^n, where OP(x,y) = exp[y * q.(x)] = (1-x)/(1-x*exp((1-x)y)), the e.g.f. of A123125, so OP(x,d/dy) y^n evaluated at y = 1 is r_n(x), the n-th row polynomial of this entry, with offset 0. - Tom Copeland, Jun 25 2018

Consolidating some formulas in this entry and A248727, in umbral notation for concision, with all offsets 0: Let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125. Then the row polynomials of this entry (A046802, the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of A248727 (the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020

From Peter Luschny, Apr 30 2021: (Start)

Sum_{k=0..n} (-1)^k*T(n, k) = A122045(n).

Sum_{k=0..n} 2^(n-k)*T(n,k) = A007047(n).

Sum_{k=0..n} T(n, n-k) = A000522(n).

Sum_{k=0..n} T(n-k, k) = Sum_{k=0..n} (n - k)^k = A026898(n-1) for n >= 1.

Sum_{k=0..n} k*T(n, k) = A036919(n) = floor(n*n!*e/2).

(End)

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))! ].

From Tom Copeland, Sep 06 2011: (Start)

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)

From Tom Copeland, Oct 20 2011: (Start)

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

From Tom Copeland, Sep 20 2016: (Start)

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)

From Tom Copeland, Sep 22 2016: (Start)

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

T(n,k) = 2^(n - 2*k - 1)*binomial(n-1,2*k)*binomial(2*k,k)/(k + 1), T(0,0) = 1, for 0 <= k <= floor((n-1)/2).

G.f.: G = G(t,z) satisfies: t*z*G^2 - (1 - 2*z + 2*t*z)*G + 1 - z + t*z = 0.

With first row zero, the o.g.f. is g(x,t) = (1 - 2*x - sqrt((1 - 2*x)^2 - 4*t*x^2)) / (2*t*x) with the inverse ginv(x,t) = x / (1 + 2*x + t*x^2), an o.g.f. for shifted A207538 and A133156 mod signs, so A134264 and A125181 can be used to interpret the polynomials of this entry. Cf. A097610. - Tom Copeland, Feb 08 2016

If we delete the first 1 from the data these are the coefficients of the polynomials p(n) = 2^n*hypergeom([(1 - n)/2, - n/2], [2], x). - Peter Luschny, Jan 23 2018

G.f.: [1-tz-sqrt(1-2tz+t^2*z^2-4z^2)]/(2z^2).

T(n, k) = n!/[k!((n-k)/2)!((n-k)/2-1)! ] = A055151(n, (n-k)/2) if n-k is a nonnegative even number; otherwise T(n, k) = 0.

T(n, k) = C(n, k)*C((n-k)/2)*(1+(-1)^(n-k))/2 if k <= n, 0 otherwise. - Paul Barry, May 18 2005

T(n,k) = A121448(n,k)/2^k. - Philippe Deléham, Aug 17 2006

Sum_{k=0..n} T(n,k)*2^k = A000108(n+1). - Philippe Deléham, Aug 22 2006

Sum_{k=0..n} T(n,k)*3^k = A002212(n+1). - Philippe Deléham, Oct 03 2007

G.f.: 1/(1-x*y-x^2/(1-x*y-x^2/(1-x*y-x^2/.... (continued fraction). - Paul Barry, Dec 15 2008

Sum_{k=0..n} T(n,k)*4^k = A005572(n). - Philippe Deléham, Dec 03 2009

T(n,k) = A007318(n,k)*A126120(n-k). - Philippe Deléham, Dec 12 2009

From Tom Copeland, Jan 23 2016: (Start)

E.g.f.: M(x,t) = e^(xt) AC(t) = e^(xt) I_1(2t)/t = e(xt) * e.g.f.(A126120(t)) = e^(xt) Sum_{n>=0} t^(2n)/(n!(n+1)!) = exp[t P(.,x)].

The e.g.f. of this Appell sequence of polynomials P(n,x) is e^(xt) times the e.g.f. AC(t) of the aerated Catalan numbers A126120. AC(t) = I_1(2t)/t, where I_n(x) = T_n(d/dx) I_0(x) are the modified Bessel functions of the first kind and T_n, the Chebyshev polynomials of the first kind.

P(n,x) has the lowering and raising operators L = d/dx = D and R = d/dD log{M(x,D)} = x + d/dD log{AC(D)} = x + Sum_{n>=0} c(n) D^(2n+1)/(2n+1)! with c(n) = (-1)^n A180874(n+1), i.e., L P(n,x) = n P(n-1,x) and R P(n,x) = P(n+1,x).

(P(.,x) + y)^n = P(n,x+y) = Sum_{k=0..n} binomial(n,k) P(k,x) y^(n-k) = (b.+x+y)^n, where (b.)^k = b_k = A126120(k).

Exp(b.D) e^(xt) = exp[(x+b.)t] = exp[P(.,x)t] = e^(b.t) e^(xt) = e^(xt) AC(t).

See p. 12 of the Alexeev et al. link and A055151 for a refinement.

Shifted o.g.f: G(x,t) = [1-tx-sqrt[(1-tx)^2-4x^2]] / 2x = x + t x^2 + (1+t) x^3 + ... has the compositional inverse Ginv(x,t) = x / [1 + tx + x^2] = x - t x^2 +(-1+t^2) x^3 + (2t-t^3) x^4 + (1-3t^2+t^4) x^5 + ..., a shifted o.g.f. for the signed Chebyshev polynomials of the second kind of A049310 (cf. also the Fibonacci polynomials of A011973). Then the inversion formula of A134264, involving non-crossing partitions and free probability with their multitude of interpretations (cf. A125181 also), can be used with h_0 = 1, h_1 = t, and h_2 = 1 to interpret the coefficients of the Motzkin polynomials combinatorially.

(End)

From Tom Copeland, Jan 29 2016: (Start)

Provides coefficients of the inverse of f(x) = x / [1 + h1 x + h2 x^2], a bivariate generating function of A049310 (mod signs).

finv(x) = [(1-h1*x) - sqrt[(1-h1*x)^2-4h2*x^2]]/(2*h2*x) = x + h1 x^2 + (h2 + h1^2) x^3 + (3 h1 h2 + h1^3) x^4 + ... is a bivariate o.g.f. for this entry.

The infinitesimal generator for finv(x) is g(x) d/dx with g(x) = 1 /[df(x)/dx] = x^2 / [(f(x))^2 (1 - h2 x^2)] = (1 + h1 x + h2 x^2)^2 / (1 - h2 x^2) so that [g(x)d/dx]^n/n! x evaluated at x = 0 gives the row polynomials FI(n,h1,h2) of the compositional inverse of f(x), i.e., exp[x g(u)d/du] u |_(u=0) = finv(x) = 1 / [1 -x FI(.,h1,h2)]. Cf. A145271. E.g.,

FI(0,h1,h2) = 0

FI(1,h1,h2) = 1

FI(2,h1,h2) = 1 h1

FI(3,h1,h2) = 1 h2 + 1 h1^2

FI(4,h1,h2) = 3 h2 h1 + 1 h1^3

FI(5,h1,h2) = 2 h2^2 + 6 h2 h1^2 + 1 h1^4

FI(6,h1,h2) = 10 h2^2 h1 + 10 h2 h1^3 + 1 h1^5.

And with D = d/dh1, FI(n+1, h1,h2) = MT(n,h1,h2) = (b.y + h1)^n = Sum_{k=0..n} binomial(n,k) b(k) y^k h1^(n-k) = exp[(b.y D] (h1)^n = AC(y D) (h1)^n, where b(k) = A126120(k), y = sqrt(h2), and AC(t) is defined in my Jan 23 formulas above. Equivalently, AC(y D) e^(x h1) = exp[x MT(.,h1,h2)].

The MT polynomials comprise an Appell sequence in h1 with e.g.f. e^(h1*x) AC(xy) = exp[x MT(.,h1,h2)] with lowering operator L = d/dh1 = D, i.e., L MT(n,h1,h2) = dMT(n,h1,h2)/dh1 = n MT(n-1,h1,h2) and raising operator R = h1 + dlog{AC(y L)}/dL = h1 + Sum_{n>=0} c(n) h2^(n+1) D^(2n+1)/(2n+1)! = h1 + h2 d/dh1 - h2^2 (d/dh1)^3/3! + 5 h2^3 (d/dh1)^5/5! - ... with c(n) = (-1)^n A180874(n+1) (consistent with the raising operator in the Jan 23 formulas).

The compositional inverse finv(x) is also obtained from the non-crossing partitions of A134264 (or A125181) with h_0 = 1, h_1 = h1, h_2 = h2, and h_n = 0 for all other n.

See A238390 for the umbral compositional inverse in h1 of MT(n,h1,h2) and inverse matrix.

(End)

From Tom Copeland, Feb 13 2016: (Start)

z1(x,h1,h2) = finv(x), the bivariate o.g.f. above for this entry, is the zero that vanishes for x=0 for the quadratic polynomial Q(z;z1(x,h1,h2),z2(x,h1,h2)) = (z-z1)(z-z2) = z^2 - (z1+z2) z + (z1*z2) = z^2 - e1 z + e2 = z^2 - [(1-h1*x)/(h2*x)] z + 1/h2, where e1 and e2 are the elementary symmetric polynomials for two indeterminates.

The other zero is given by z2(x,h1,h2) = (1 - h1*x)/(h2*x) - z1(x,h1,h2) = [1 - h1*x + sqrt[(1-h1*x)^2 - 4 h2*x^2]] / (2h2*x).

The two are the nontrivial zeros of the elliptic curve in Legendre normal form y^2 = z (z-z1)(z-z2), (see Landweber et al., p. 14, Ellingsrud, and A121448), and the zeros for the Riccati equation z' = (z - z1)(z - z2), associated to soliton solutions of the KdV equation (see Copeland link).

(End)

Comparing the shifted o.g.f. S(x) = x / (1 - h_1 x + h_2 x^2) for the bivariate Chebyshev polynomials S_n(h_1,h_2) of A049310 with the shifted o.g.f. H(x) = x / ((1 - a x)(1 - b x)) for the complete homogeneous symmetric polynomials H_n(a,b) = (a^(n+1)-b^(n+1)) / (a - b) shows that S_n(h_1,h_2) = H_n(a,b) for h_1 = a + b and h_2 = ab and, conversely, a = (h_1 + sqrt(h_1^2 - 4 h_2)) / 2 and b = (h_1 - sqrt(h_1^2 - 4 h_2)) / 2. The compositional inverse about the origin of S(x) gives a bivariate o.g.f. for signed Motzkin polynomials M_n(h_1,h_2) of this entry, and that of H(x) gives one for signed Narayana polynomials N_n(a,b) of A001263, thereby relating the bivariate Motzkin and Narayana polynomials by the indeterminate transformations. E.g., M_2(h_1,h_2) = h_2 + h_1^2 = ab + (a + b)^2 = a^2 + 3 ab + b^2 = N_2(a,b). - Tom Copeland, Jan 27 2024

T(n,k) = n!/((n-2k)! k! (k+1)!) = A007318(n, 2k)*A000108(k). - Henry Bottomley, Jan 31 2003

E.g.f. row polynomials R(n,y): exp(x)*BesselI(1, 2*x*sqrt(y))/(x*sqrt(y)). - Vladeta Jovovic, Aug 20 2003

G.f. row polynomials R(n,y): 2 / (1 - x + sqrt((1 - x)^2 - 4 *y * x^2)).

From Peter Bala, Oct 28 2008: (Start)

The rows of this triangle are the gamma vectors of the n-dimensional (type A) associahedra (Postnikov et al., p. 38). Cf. A089627 and A101280.

The row polynomials R(n,x) = Sum_{k = 0..n} T(n,k)*x^k begin R(0,x) = 1, R(1,x) = 1, R(2,x) = 1 + x, R(3,x) = 1 + 3*x. They are related to the Narayana polynomials N(n,x) := Sum_{k = 1..n} (1/n)*C(n,k)*C(n,k-1)*x^k through x*(1 + x)^n*R(n, x/(1 + x)^2) = N(n+1,x). For example, for n = 3, x*(1 + x)^3*(1 + 3*x/(1 + x)^2) = x + 6*x^2 + 6*x^3 + x^4, the 4th Narayana polynomial.

Recursion relation: (n + 2)*R(n,x) = (2*n + 1)*R(n-1,x) - (n - 1)*(1 - 4*x)*R(n-2,x), R(0,x) = 1, R(1,x) = 1. (End)

G.f.: M(x,y) satisfies: M(x,y)= 1 + x M(x,y) + y*x^2*M(x,y)^2. - Geoffrey Critzer, Feb 05 2014

T(n,k) = A161642(n,k)*A258820(n,k) = (binomial(n,k)/A003989(n+1, k+1))* A258820(n,k). - Tom Copeland, Jun 18 2015

Let T(n,k;q) = n!*(1+k)/((n-2*k)!*(1+k)!^2)*hypergeom([k,2*k-n],[k+2],q) then T(n,k;0) = A055151(n,k), T(n,k;1) = A008315(n,k) and T(n,k;-1) = A091156(n,k). - Peter Luschny, Oct 16 2015

From Tom Copeland, Jan 21 2016: (Start)

Reversed rows of A107131 are rows of this entry, and the diagonals of A107131 are the columns of this entry. The diagonals of this entry are the rows of A088617. The antidiagonals (bottom to top) of A088617 are the rows of this entry.

O.g.f.: [1-x-sqrt[(1-x)^2-4tx^2]]/(2tx^2), from the relation to A107131.

Re-indexed and signed, this triangle gives the row polynomials of the compositional inverse of the shifted o.g.f. for the Fibonacci polynomials of A011973, x / [1-x-tx^2] = x + x^2 + (1+t) x^3 + (1+2t) x^4 + ... . (End)

Row polynomials are P(n,x) = (1 + b.y)^n = Sum{k=0 to n} binomial(n,k) b(k) y^k = y^n M(n,1/y), where b(k) = A126120(k), y = sqrt(x), and M(n,y) are the Motzkin polynomials of A097610. - Tom Copeland, Jan 29 2016

Coefficients of the polynomials p(n,x) = hypergeom([(1-n)/2, -n/2], [2], 4x). - Peter Luschny, Jan 23 2018

Sum_{k=1..floor(n/2)} k * T(n,k) = A014531(n-1) for n>1. - Alois P. Heinz, Mar 29 2020

T(n,k) = binomial(n+1,2k-n).

G.f.: 1/(1 - 2*t*z - t*(1-t)*z^2).

T(n,k) = A034867(n,n-k)

From Tom Copeland, Sep 30 2011: (Start)

With K(x,t) = 1/{d/dx{x/[t-1+1/(1-x)]}} = [t-1+1/(1-x)]^2/{t-[x/(1-x)]^2}, the g.f. of A119900 = K(x*t,t)-t+1.

From formulas in A134264: K(x,t)d/dx is a generator for A001263. A refinement of A119900 to partition polynomials is given by umbralizing

K(x,t) roughly as K(h.x,h_0) and precisely as in A134264 as

W(x)= 1/{d/dx[f(x)]}=1/{d/dx[x/h(x)]}. (End)

T(n,k) = 2*T(n-1,k-1) + T(n-2,k-1) - T(n-2,k-2). - Philippe Deléham, Oct 02 2011

From Tom Copeland, Dec 07 2015: (Start)

An alternate o.g.f. is (1/(x*t)) {-1 + 1 / [1 - (1/t)[x*t/(1-x*t)]^2]} = Sum_{n>0} x^(2(n-1)+1) t^(n-1) / (1-t*x)^(2n) = x + 2t x^2 + (t+3t^2) x^3 + ... .

The n-th diagonal has elements binomial(2n+1+k,k), starting with k=0 for the first non-vanishing element, with o.g.f. (1-x)^(-2(n+1)). The first few subdiagonals are shifted versions of A000292, A000389, and A000580. Cf. A049310.

See A034867 for the matrix representation for the infinitesimal generator K(x,t) d/dx for the Narayana polynomials. (End)

From Peter Bala, Aug 17 2016: (Start)

Let S(k,n) = Sum_{i = 1..n} i^k. Calculations in Zielinski 2016 suggest the following identity holds involving the p-th row elements of this triangle:

Sum_{k = 0..p} T(p,k)*S(2*k + 1,n) = (n*(n + 1)/2)^(p+1).

For example, for row 6 we find S(7,n) + 21*S(9,n) + 35*S(11,n) + 7*S(13,n) = (n*(n + 1)/2)^7.

There appears to be a similar result for the even power sums S(2*k,n) involving A207543. (End)