cp's OEIS Frontend

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

The numbers (4, 20, 120, 840, 6720, ...) arise from the divisor values in the general formula a(n) = n*(n+1)*(n+2)*(n+3)* ... *(n+k)*(n*(n+k) + (k-1)*k/6)/((k+3)!/6) (which covers the following sequences: A000578, A000537, A024166, A101094, A101097, A101102). - Alexander R. Povolotsky, May 17 2008

a(n) is also the number of decreasing 3-cycles in the decomposition of permutations as product of disjoint cycles, a(3)=1, a(4)=4, a(5)=20. - Wenjin Woan, Dec 21 2008

Equals eigensequence of triangle A130128 reflected. - Gary W. Adamson, Dec 23 2008

a(n) is the number of n-permutations having 1, 2, and 3 in three distinct cycles. - Geoffrey Critzer, Apr 26 2009

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

The asymptotic expansion of the higher order exponential integral E(x,m=1,n=4) ~ exp(-x)/x*(1 - 4/x + 20/x^2 - 120/x^3 + 840/x^4 - 6720/x^5 + 60480/x^6 - 604800/x^7 + ...) leads to the sequence given above. See A163931 and A130534 for more information.

(End)

The asymptotic expansion of the higher-order exponential integral E(x,m=1,n=5) ~ exp(-x)/x*(1 - 5/x + 30/x^2 - 210/x^3 + 1680/x^4 - 15120/x^5 + 151200/x^6 - 1663200/x^7 + ...) leads to this sequence. See A163931 and A130534 for more information. - Johannes W. Meijer, Oct 20 2009

The (reverse) Bessel polynomials P(n,x):=Sum_{m=0..n} a(n,m)*x^m, the row polynomials, called Theta_n(x) in the Grosswald reference, solve x*(d^2/dx^2)P(n,x) - 2*(x+n)*(d/dx)P(n,x) + 2*n*P(n,x) = 0.

With the related Sheffer associated polynomials defined by Carlitz as

B(0,x) = 1

B(1,x) = x

B(2,x) = x + x^2

B(3,x) = 3 x + 3 x^2 + x^3

B(4,x) = 15 x + 15 x^2 + 6 x^3 + x^4

... (see Mathworld reference), then P(n,x) = 2^n * B(n,x/2) are the Sheffer polynomials described in A119274. - Tom Copeland, Feb 10 2008

Exponential Riordan array [1/sqrt(1-2x), 1-sqrt(1-2x)]. - Paul Barry, Jul 27 2010

From Vladimir Kruchinin, Mar 18 2011: (Start)

For B(n,k){...} the Bell polynomial of the second kind we have

B(n,k){f', f'', f''', ...} = T(n-1,k-1)*(1-2*x)^(k/2-n), where f(x) = 1-sqrt(1-2*x).

The expansions of the first few rows are:

1/sqrt(1-2*x);

1/(1-2*x)^(3/2), 1/(1-2*x);

3/(1-2*x)^(5/2), 3/(1-2*x)^2, 1/(1-2*x)^(3/2);

15/(1-2*x)^(7/2), 15/(1-2*x)^3, 6/(1-2*x)^(5/2), 1/(1-2*x)^2. (End)

Also the Bell transform of A001147 (whithout column 0 which is 1,0,0,...). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 19 2016

Antidiagonals of A099174 are rows of this entry. Dividing each diagonal by its first element generates A054142. - Tom Copeland, Oct 04 2016

The row polynomials p_n(x) of A107102 are (-1)^n B_n(1-x), where B_n(x) are the modified Carlitz-Bessel polynomials above, e.g., (-1)^2 B_2(1-x) = (1-x) + (1-x)^2 = 2 - 3 x + x^2 = p_2(x). - Tom Copeland, Oct 10 2016

a(n-1,m-1) counts rooted unordered binary forests with n labeled leaves and m roots. - David desJardins, Feb 23 2019

From Jianing Song, Nov 29 2021: (Start)

The polynomials P_n(x) = Sum_{k=0..n} T(n,k)*x^k satisfy: P_n(x) - (d/dx)P_n(x) = x*P_{n-1}(x) for n >= 1.

{P(n,x)} are related to the Fourier transform of 1/(1+x^2)^(n+1) and x/(1+x^2)^(n+2):

(i) For n >= 0, real number t, we have Integral_{x=-oo..oo} exp(-i*t*x)/(1+x^2)^(n+1) dx = Pi/(2^n*n!) * P_n(|t|) * exp(-|t|);

(ii) For n >= 0, real number t, we have Integral_{x=-oo..oo} x*exp(-i*t*x)/(1+x^2)^(n+2) dx = Pi/(2^(n+1)*(n+1)!) * ((-t)*P_n(-|t|)) * exp(-|t|). (End)

Suppose that f(x) is an n-times differentiable function defined on (a,b) for 0 <= a < b <= +oo, then for n >= 1, the n-th derivative of f(sqrt(x)) on (a^2,b^2) is Sum_{k=1..n} ((-1)^(n-k)*T(n-1,k-1)*f^(k)(sqrt(x))) / (2^n*x^(n-(k/2))), where f^(k) is the k-th derivative of f. - Jianing Song, Nov 30 2023

Let M(t) = exp(t*T) = lim_{n->oo} (1 + t*T/n)^n.

Pascal matrix = [ binomial(n,k) ] = M(1) = exp(T), truncating the series gives the n X n submatrices.

Inverse Pascal matrix = M(-1) = exp(-T) = matrix for inverse binomial transform.

A(j) = T^j / j! equals the matrix [binomial(n,k) * delta(n-k-j)] where delta(n) = 1 if n=0 and vanishes otherwise (Kronecker delta); i.e., A(j) is a matrix with all the terms 0 except for the j-th lower (or main for j=0) diagonal, which equals that of the Pascal triangle. Hence the A(j)'s form a linearly independent basis for all matrices of the form [binomial(n,k) * d(n-k)] which include as a subset the invertible associated matrices of the list partition transform (LPT) of A133314.

For sequences with b(0) = 1, umbrally,

M[b(.)] = exp(b(.)*T) = [ binomial(n,k) * b(n-k) ] = matrices associated to b by LPT.

[M[b(.)]]^(-1) = exp(c(.)*T) = [ binomial(n,k) * c(n-k) ] = matrices associated to c, where c = LPT(b) . Or,

[M[b(.)]]^(-1) = exp[LPT(b(.))*T] = LPT[M(b(.))] = M[LPT(b(.))]= M[c(.)].

This is related to xDx, the iterated Laguerre transform and the general Euler transformation of a sequence through the comments in A132013 and A132014 and the relation [Sum_{k=0..n} binomial(n,k) * b(n-k) * d(k)] = M(b)*d, (n-th term). See also A132382.

If b(n,x) is a binomial type Sheffer sequence, then M[b(.,x)]*s(y) = s(x+y) when s(y) = (s(0,y),s(1,y),s(2,y),...) is an array for a Sheffer sequence with the same delta operator as b(n,x) and [M[b(.,x)]]^(-1) is given by the formulas above with b(n) replaced by b(n,x) as b(0,x)=1 for a binomial-type Sheffer sequence.

T = I - A132013 and conversely A132013 = I - T, which is the matrix representation for the iterated mixed order Laguerre transform characterized in A132013 (and A132014).

(I-T)^m generates the group [A132013]^m for m = 0,1,2,... discussed in A132014.

The inverse is 1/(I-T) = I + T + T^2 + T^3 + ... = [A132013]^(-1) = A094587 with the associated sequence (0!,1!,2!,3!,...) under the LPT.

And 1/(I-T)^2 = I + 2*T + 3*T^2 + 4*T^3 + ... = [A132013]^(-2) = A132159 with the associated sequence (1!,2!,3!,4!,...) under the LPT.

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

1) b(0) = 0, b(n) = n * a(n-1),

2) B(x) = xDx A(x)

3) B(x) = x * Lag(1,-:xD:) A(x)

4) EB(x) = x * EA(x) where D is the derivative w.r.t. x, (:xD:)^j = x^j*D^j and Lag(n,x) is the Laguerre polynomial.

So the exponentiated operator can be characterized as

5) exp(t*T) A(x) = exp(t*xDx) A(x) = [Sum_{n=0,1,...} (t*x)^n * Lag(n,-:xD:)] A(x) = [exp{[t*u/(1-t*u)]*:xD:} / (1-t*u) ] A(x) (eval. at u=x) = A[x/(1-t*x)]/(1-t*x), a generalized Euler transformation for an o.g.f.,

6) exp(t*T) EA(x) = exp(t*x)*EA(x) = exp[(t+a(.))*x], gen. Euler trf. for an e.g.f.

7) exp(t*T) * a = M(t) * a = [Sum_{k=0..n} binomial(n,k) * t^(n-k) * a(k)].

The umbral extension of formulas 5, 6 and 7 gives formally

8) exp[c(.)*T] A(x) = exp(c(.)*xDx) A(x) = [Sum_{n>=0} (c(.)*x)^n * Lag(n,-:xD:)] A(x) = [exp{[c(.)*u/(1-c(.)*u)]*:xD:} / (1-c(.)*u) ] A(x) (eval. at u=x) = A[x/(1-c(.)*x)]/(1-c(.)*x), where the umbral evaluation should be applied only after a power series in c is obtained,

9) exp[c(.)*T] EA(x) = exp(c(.)*x)*EA(x) = exp[(c(.)+a(.))*x]

10) exp[c(.)*T] * a = M[c(.)] * a = [Sum_{k=0..n} binomial(n,k) * c(n-k) * a(k)] .

The n X n principal submatrix of T is nilpotent, in particular, [Tsub_n]^(n+1) = 0, n=0,1,2,3,....

Note (xDx)^n = x^n D^n x^n = x^n n! (:Dx:)^n/n! = x^n n! Lag(n,-:xD:).

The operator xDx is an important, classical operator explored by among others Dattoli, Al-Salam, Carlitz and Stokes and even earlier investigators.

For a recent treatment of xDx, DxD and more general operators see the paper "Laguerre-type derivatives: Dobinski relations and combinatorial identities". - Karol A. Penson, Sep 15 2009

See Copeland's link for generalized Laguerre functions and connection to fractional differ-integrals in exercises through (:Dx:)^a/a!=(D^a x^a)/a!. - Tom Copeland, Nov 17 2011

From Tom Copeland, Apr 25 2014: (Start)

Conjugation or "similarity" transformations of [dP]=A132440 have an operator interpretation (cf. A074909 and A238363):

In general, select two operators A and B such that A^n = F1(n,B) and B^n = F2(n,A); then A^n = F1(n,F2(.,A)) and B^n = F2(n,F1(.,B)), evaluated umbrally, i.e., F1(n,F2(.,x))=F2(n,F1(.,x))=x^n, implying the polynomials F1 and F2 are an umbral compositional inverse pair.

One such pair are the Bell polynomials Bell(n,x) and falling factorials (x)_n with Bell(n,:xD:)=(xD)^n and (xD)_n=:xD:^n (cf. A074909). Another are the Laguerre polynomials LN(n,x)= n!*Lag(n,x) (A021009), which are umbrally self-inverse, with LN(n,-:xD:)=:Dx:^n and LN(n,:Dx:)= (-:xD:)^n with :Dx:^n=D^n*x^n.

Evaluating, for n>=0, the operator derivative d(B^n)/dA = d(F2(n,A))/dA in the basis B^n, i.e., with A^n finally replaced by F1(n,B), or A^n=F1(.,B)^n=F1(n,B), is equivalent to the matrix conjugation

A) [F2]*[dP]*[F1]

B) = [F2]*[dP]*[F2]^(-1)

C) = [F1]^(-1)*[dP]*[F1],

where [F1] is the lower triangular matrix with the n-th row the coefficients of F1(n,x) and analogously for [F2].

So, given the row vector Rv=(c0 c1 c2 c3 ...) and the column vector Cv(x)=(1 x x^2 x^3 ...)^Transpose, form the power series V(x)=Rv*Cv(x).

D) dV(B)/dA = Rv * [F2]*[dP]*[F1] * Cv(B).

E) With A=D and B=D, F1(n,x)=F2(n,x)=x^n and [F1]=[F2]=I. Then d(B^n)/dA = d(D^n)/dD = n * D^(n-1); therefore, consistently [F2]*[dP]*[F1] = [dP] and dV(D)/dD = Rv * [dP] * Cv(D). (End)

Stirling transform of (-1)^n*a(n-1) = [0, 1, -1, 5, -14, 94, ...] is A000629(n-2) = [0, 1, 2, 6, 26, ...]. - Michael Somos, Mar 04 2004

Stirling transform of a(n) = [1, 1, 5, 14, 94, ...] is A052882(n) = [1, 2, 9, 52, 375, ...]. - Michael Somos, Mar 04 2004

a(n) is the number of n-permutations that have a cycle with length greater than n/2. - Geoffrey Critzer, May 28 2009

From Jens Voß, May 07 2010: (Start)

a(4n) is divisible by 6*n + 1 for all n >= 1; the quotient of a(4*n) and 6*n+1 is A177188(n).

a(4*n+3) is divisible by 6*n + 5 for all n >= 0; the quotient of a(4*n+3) and 6*n + 5 is A177174(n). (End)

From Wolfdieter Lang, Jun 27 2012: (Start)

T(n-1,k), k=1,...,n-1, gives the number of representative necklaces with n beads (C_N symmetry) of n+1-k distinct colors, say c[1],c[2],...,c[n-k+1], corresponding to the color signature determined by the partition k,1^(n-k) of n. The representative necklaces have k beads of color c[1]. E.g., n=4, k=2: partition 2,1,1, color signature (parts as exponents) c[1]c[1]c[2]c[3], 3=T(3,2) necklaces (write j for color c[j]): cyclic(1123), cyclic(1132) and cyclic(1213). See A212359 for the numbers for general partitions or color signatures. (End)

The coefficients of the Faber polynomials F(n,b(1),b(2),...,b(n)) (Bouali, p. 52) in the order of the partitions of Abramowitz and Stegun. Compare with A115131 and A210258.

These polynomials occur in discussions of the Virasoro algebra, univalent function spaces and the Schwarzian derivative, symmetric functions, and free probability theory. They are intimately related to symmetric functions, free probability, and Appell sequences through the raising operator R = x - d log(H(D))/dD for the Appell sequence inverse pair associated to the e.g.f.s H(t)e^(xt) (cf. A094587) and (1/H(t))e^(xt) with H(0)=1.

Instances of the Faber polynomials occur in discussions of modular invariants and modular functions in the papers by Asai, Kaneko, and Ninomiya, by Ono and Rolen, and by Zagier. - Tom Copeland, Aug 13 2019

The Faber polynomials, denoted by s_n(a(t)) where a(t) is a formal power series defined by a product formula, are implicitly defined by equation 13.4 on p. 62 of Hazewinkel so as to extract the power sums of the reciprocals of the zeros of a(t). This is the Newton identity expressing the power sum symmetric polynomials in terms of the elementary symmetric polynomials/functions. - Tom Copeland, Jun 06 2020

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)

Triangle in which the n-th row begins with n and the k-th term is obtained by multiplying the (k-1)-th term by (n-k+1) until n-k+1 = 1. - Amarnath Murthy, Nov 11 2002

Has many diagonals in common with A105725. - Zerinvary Lajos, Apr 14 2006

Also: the square array of rising factorials A(n,k) = n*(n+1)*(n+2)*...*(n+k-1) read by antidiagonals downwards. There are no perfect squares in T(n,k) for k >= 2 [see Rigge]. T(n,k) is divisible by a prime exceeding k, if n >= 2*k [see Saradha and Shorey]. - R. J. Mathar, May 02 2007

T(n,k) is the number of injective functions f from {1,...,k} into {1,...,n}, since there are n choices for f(1), then (n-1) choices for f(2), ... and (n-k+1) choices for f(k). E.g., T(3,2)=6 because there are exactly 6 injective functions f:{1,2}->{1,2,3}, namely, f1={(1,1),(2,2)}, f2={(1,1),(2,3)}, f3={(1,2),(2,1)}, f4={(1,2),(2,3)}, f5={(1,3),(2,1)} and f6={(1,3),(2,2)}. - Dennis P. Walsh, Oct 18 2007

Permuted words defined by the connectivity of regular simplices are related to T by T = A135278 * (1!, 2!, 3!, 4!, ...). E.g., for T(4,k) with k-1 = simplex number, label the vertices of a tetrahedron with a, b, c, d, then the 0-simplex, the points, a,b,c,d gives 4 * 1 = 4 words; the 1-simplex, the edges: (ab or ba), (ac or ca), (ad or da), (bc or cb), (bd or db), (cd or dc) gives 6 * 2 = 12 words; the 2-simplex, the faces: (abc or ...), (acd or ...), (adb or ...), (bcd or ...) gives 4 * 6 = 24 words; the 3-simplex, (abcd or ....) gives 1 * 24 = 24 words. - Tom Copeland, Dec 08 2007

Reversal of the triangle by rows = (n+1) * n-th row of triangle A094587. - Gary W. Adamson, May 03 2009

From Geoffrey Critzer, May 06 2009: (Start)

The rectangular array R(n,k), read by diagonals is the number of ways n people can queue up in k (possibly empty) distinct queues. R(n,k) = (n+k-1)!/(k-1)!; R(n,k) = (n+k-1)*R(n-1,k).

Northwest corner:

1, 2, 3, 4, 5, ...;

2, 6, 12, 20, 30, ...;

6, 24, 60, 120, 210, ...;

24, 120, 360, 840, 1680, ...;

120, 720, 2520, 6720, 15120, ...;

...

Example: R(2,2)=6 because there are six ways that two people can get in line at a fast food restaurant that has two order windows open. Let 1 and 2 represent the two people and a | will separate the lines. 12|; 21|; |12; |21; 1|2; 2|1. (End)

Cf. [Hardy and Wright], Theorem 34.

The e.g.f. of the Norlund generalized Bernoulli (Appell) polynomials of order m, NB(n,x;m), is given by exponentiation of the e.g.f. of the Bernoulli numbers, i.e., multiple binomial self-convolutions of the Bernoulli numbers, through the e.g.f. exp[NB(.,x;m)t] = [t/(e^t-1)]^(m+1) * e^(xt). Norlund gave the relation to the factorials (x-1)!/(x-1-k)! = (x-1) ... (x-k) = NB(k,x;k), so T(n,k) = NB(k,n+1;k). - Tom Copeland, Oct 01 2015

T(n,k) is the number of sequences without repetition (partial permutations) of k elements taken from an n-set. For sequences with repetition (k-tuples) cf. A075363. - Manfred Boergens, Jun 18 2023

Let D=d/dx and [A,B]=A·B-B·A. Then each row corresponds to the coefficients of the operators :xD:^k = x^k D^k in the expansion of the commutator [log(D),:xD:^n]=[-log(x),:xD:^n]=sum(k=0 to n-1, a(n,k) :xD:^k). The e.g.f. is derived from [log(D), exp(t:xD:)]=[-log(x), exp(t:xD:)]= log(1+t)exp(t:xD:), using the shift property exp(t:xD:)f(x)=f((1+t)x).

The reversed unsigned array is A111492.

See the mathoverflow link and link therein to an associated mathstackexchange question for other formulas for log(D). In addition, R_x = log(D) = -log(x) + c - sum[n=1 to infnty, (-1)^n 1/n :xD:^n/n!]=

-log(x) + Psi(1+xD) = -log(x) + c + Ein(:xD:), where c is the Euler-Mascheroni constant, Psi(x), the digamma function, and Ein(x), a breed of the exponential integrals (cf. Wikipedia). The :xD:^k ops. commute; therefore, the commutator reduces to the -log(x) term.

Also the n-th row corresponds to the expansion of d[(xD)!/(xD-n)!]/d(xD) = d[:xD:^n]/d(xD) in the operators :xD:^k, or, equivalently, the coefficients of x in d[z!/(z-n)!]/dz=d[St1(n,z)]]/dz evaluated umbrally with z=St2(.,x), i.e., z^n replaced by St2(n,x), where St1(n,x) and St2(n,x) are the signed and unsigned Stirling polynomials of the first (A008275) and second (A008277) kinds. The derivatives of the unsigned St1 are A028421. See examples. This formalism follows from the relations between the raising and lowering operators presented in the MathOverflow link and the Pincherle derivative. The results can be generalized through the operator relations in A094638, which are related to the celebrated Witt Lie algebra and pseudodifferential operators / symbols, to encompass other integral arrays.

A002741(n)*(-1)^(n+1) (row sums), A002104(n)*(-1)^(n+1) (alternating row sums). Column sequences: A133942(n-1), A001048(n-1), A238474, ... - Wolfdieter Lang, Mar 01 2014

Add an additional head row of zeros to the lower triangular array and denote it as T (with initial indexing in columns and rows being 0). Let dP = A132440, the infinitesimal generator for the Pascal matrix, and I, the identity matrix, then exp(T)=I+dP, i.e., T=log(I+dP). Also, (T_n)^n=0, where T_n denotes the n X n submatrix, i.e., T_n is nilpotent of order n. - Tom Copeland, Mar 01 2014

Any pair of lowering and raising ops. L p(n,x) = n·p(n-1,x) and R p(n,x) = p(n+1,x) satisfy [L,R]=1 which implies (RL)^n = St2(n,:RL:), and since (St2(·,u))!/(St2(·,u)-n)!= u^n, when evaluated umbrally, d[(RL)!/(RL-n)!]/d(RL) = d[:RL:^n]/d(RL) is well-defined and gives A238363 when the LHS is reduced to a sum of :RL:^k terms, exactly as for L=d/dx and R=x above. (Note that R_x above is a raising op. different from x, with associated L_x=-xD.) - Tom Copeland, Mar 02 2014

For relations to colored forests, disposition of flags on flagpoles, and the colorings of the vertices of the complete graphs K_n, encoded in their chromatic polynomials, see A130534. - Tom Copeland, Apr 05 2014

The unsigned triangle, omitting the main diagonal, gives A211603. See also A092271. Related to the infinitesimal generator of A008290. - Peter Bala, Feb 13 2017