cp's OEIS Frontend

Total number of permutations of all subsets of an n-set.

Also the number of one-to-one sequences that can be formed from n distinct objects.

Old name "Total number of permutations of a set with n elements", or the same with the word "arrangements", both sound too much like A000142.

Related to number of operations of addition and multiplication to evaluate a determinant of order n by cofactor expansion - see A026243.

a(n) is also the number of paths (without loops) in the complete graph on n+2 vertices starting at one vertex v1 and ending at another v2. Example: when n=2 there are 5 paths in the complete graph with 4 vertices starting at the vertex 1 and ending at the vertex 2: (12),(132),(142),(1342),(1432) so a(2) = 5. - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001; comment corrected by Jonathan Coxhead, Mar 21 2003

Also row sums of Table A008279, which can be generated by taking the derivatives of x^k. For example, for y = 1*x^3, y' = 3x^2, y" = 6x, y'''= 6 so a(4) = 1 + 3 + 6 + 6 = 16. - Alford Arnold, Dec 15 1999

a(n) is the permanent of the n X n matrix with 2s on the diagonal and 1s elsewhere. - Yuval Dekel, Nov 01 2003

(A000166 + this_sequence)/2 = A009179, (A000166 - this_sequence)/2 = A009628.

Stirling transform of A006252(n-1) = [1,1,1,2,4,14,38,...] is a(n-1) = [1,2,5,16,65,...]. - Michael Somos, Mar 04 2004

Number of {12,12*,21*}- and {12,12*,2*1}-avoiding signed permutations in the hyperoctahedral group.

a(n) = b such that Integral_{x=0..1} x^n*exp(-x) dx = a-b*exp(-1). - Sébastien Dumortier, Mar 05 2005

a(n) is the number of permutations on [n+1] whose left-to-right record lows all occur at the start. Example: a(2) counts all permutations on [3] except 231 (the last entry is a record low but its predecessor is not). - David Callan, Jul 20 2005

a(n) is the number of permutations on [n+1] that avoid the (scattered) pattern 1-2-3|. The vertical bar means the "3" must occur at the end of the permutation. For example, 21354 is not counted by a(4): 234 is an offending subpermutation. - David Callan, Nov 02 2005

Number of deco polyominoes of height n+1 having no reentrant corners along the lower contour (i.e., no vertical step that is followed by a horizontal step). In other words, a(n)=A121579(n+1,0). A deco polyomino is a directed column-convex polyomino in which the height, measured along the diagonal, is attained only in the last column. Example: a(1)=2 because the only deco polyominoes of height 2 are the vertical and horizontal dominoes, having no reentrant corners along their lower contours. - Emeric Deutsch, Aug 16 2006

Unreduced numerators of partial sums of the Taylor series for e. - Jonathan Sondow, Aug 18 2006

a(n) is the number of permutations on [n+1] (written in one-line notation) for which the subsequence beginning at 1 is increasing. Example: a(2)=5 counts 123, 213, 231, 312, 321. - David Callan, Oct 06 2006

a(n) is the number of permutations (written in one-line notation) on the set [n + k], k >= 1, for which the subsequence beginning at 1,2,...,k is increasing. Example: n = 2, k = 2. a(2) = 5 counts 1234, 3124, 3412, 4123, 4312. - Peter Bala, Jul 29 2014

a(n) and (1,-2,3,-4,5,-6,7,...) form a reciprocal pair under the list partition transform and associated operations described in A133314. - Tom Copeland, Nov 01 2007

Consider the subsets of the set {1,2,3,...,n} formed by the first n integers. E.g., for n = 3 we have {}, {1}, {2}, {3}, {1,2}, {1,3}, {2,3}, {1,2,3}. Let the variable sbst denote a subset. For each subset sbst we determine its number of parts, that is nprts(sbst). The sum over all possible subsets is written as Sum_{sbst=subsets}. Then a(n) = Sum_{sbst=subsets} nprts(sbst)!. E.g., for n = 3 we have 1!+1!+1!+1!+2!+2!+2!+3!=16. - Thomas Wieder, Jun 17 2006

Equals row sums of triangle A158359(unsigned). - Gary W. Adamson, Mar 17 2009

Equals eigensequence of triangle A158821. - Gary W. Adamson, Mar 30 2009

For positive n, equals 1/BarnesG(n+1) times the determinant of the n X n matrix whose (i,j)-coefficient is the (i+j)th Bell number. - John M. Campbell, Oct 03 2011

a(n) is the number of n X n binary matrices with i) at most one 1 in each row and column and ii) the subset of rows that contain a 1 must also be the columns that contain a 1. Cf. A002720 where restriction ii is removed. - Geoffrey Critzer, Dec 20 2011

Number of restricted growth strings (RGS) [d(1),d(2),...,d(n)] such that d(k) <= k and d(k) <= 1 + (number of nonzero digits in prefix). The positions of nonzero digits determine the subset, and their values (decreased by 1) are the (left) inversion table (a rising factorial number) for the permutation, see example. - Joerg Arndt, Dec 09 2012

Number of a restricted growth strings (RGS) [d(0), d(1), d(2), ..., d(n)] where d(k) >= 0 and d(k) <= 1 + chg([d(0), d(1), d(2), ..., d(k-1)]) and chg(.) gives the number of changes of its argument. Replacing the function chg(.) by a function asc(.) that counts the ascents in the prefix gives A022493 (ascent sequences). - Joerg Arndt, May 10 2013

The sequence t(n) = number of i <= n such that floor(e*i!) is a square is mentioned in the abstract of Luca & Shparlinski. The values are t(n) = 0 for 0 <= n <= 2 and t(n) = 1 for at least 3 <= n <= 300. - R. J. Mathar, Jan 16 2014

a(n) = p(n,1) = q(n,1), where p and q are polynomials defined at A248664 and A248669. - Clark Kimberling, Oct 11 2014

a(n) is the number of ways at most n people can queue up at a (slow) ticket counter when one or more of the people may choose not to queue up. Note that there are C(n,k) sets of k people who quene up and k! ways to queue up. Since k can run from 0 to n, a(n) = Sum_{k=0..n} n!/(n-k)! = Sum_{k=0..n} n!/k!. For example, if n=3 and the people are A(dam), B(eth), and C(arl), a(3)=16 since there are 16 possible lineups: ABC, ACB, BAC, BCA, CAB, CBA, AB, BA, AC, CA, BC, CB, A, B, C, and empty queue. - Dennis P. Walsh, Oct 02 2015

As the row sums of A008279, Motzkin uses the abbreviated notation $n_<^\Sigma$ for a(n).

The piecewise polynomial function f defined by f(x) = a(n)*x^n/n! on each interval [ 1-1/a(n), 1-1/a(n+1) ) is continuous on [0,1) and lim_{x->1} f(x) = e. - Luc Rousseau, Oct 15 2019

a(n) is composite for 3 <= n <= 2015, but a(2016) is prime (or at least a strong pseudoprime): see Johansson link. - Robert Israel, Jul 27 2020 [a(2016) is prime, ECPP certificate generated with CM 0.4.3 and checked by factordb. - Jason H Parker, Jun 15 2025]

In general, sequences of the form a(0)=a, a(n) = n*a(n-1) + k, n>0, will have a closed form of n!*a + floor(n!*(e-1))*k. - Gary Detlefs, Oct 26 2020

From Peter Bala, Apr 03 2022: (Start)

a(2*n) is odd and a(2*n+1) is even. More generally, a(n+k) == a(n) (mod k) for all n and k. It follows that for each positive integer k, the sequence obtained by reducing a(n) modulo k is periodic, with the exact period dividing k. Various divisibility properties of the sequence follow from this; for example, a(5*n+2) == a(5*n+4) == 0 (mod 5), a(25*n+7) == a(25*n+19) == 0 (mod 25) and a(13*n+4) == a(13*n+10)== a(13*n+12) == 0 (mod 13). (End)

Number of possible ranking options on a typical ranked choice voting ballot with n candidates (allowing undervotes). - P. Christopher Staecker, May 05 2024

From Thomas Scheuerle, Dec 28 2024: (Start)

Number of decorated permutations of size n.

Number of Le-diagrams with bounding box semiperimeter n, for n > 0.

By counting over all k = 1..n and n > 0, the number of positroid cells for the totally nonnegative real Grassmannian Gr(k, n), equivalently the number of Grassmann necklaces of type (k, n). (End)

Also, table of Pochhammer sequences read by antidiagonals (see Rudolph-Lilith, 2015). - N. J. A. Sloane, Mar 31 2016

Reverse of A008279. Row sums are A000522. Diagonal sums are A003470. Rows of inverse matrix begin {1}, {-1,1}, {0,-2,1}, {0,0,-3,1}, {0,0,0,-4,1} ... The signed lower triangular matrix (-1)^(n+k)n!/k! has as row sums the signed rencontres numbers Sum_{k=0..n} (-1)^(n+k)n!/k!. (See A000166). It has matrix inverse 1 1,1 0,2,1 0,0,3,1 0,0,0,4,1,...

Exponential Riordan array [1/(1-x),x]; column k has e.g.f. x^k/(1-x). - Paul Barry, Mar 27 2007

From Tom Copeland, Nov 01 2007: (Start)

T is the umbral extension of n!*Lag[n,(.)!*Lag[.,x,-1],0] = (1-D)^(-1) x^n = (-1)^n * n! * Lag(n,x,-1-n) = Sum_{j=0..n} binomial(n,j) * j! * x^(n-j) = Sum_{j=0..n} (n!/j!) x^j. The inverse operator is A132013 with generalizations discussed in A132014.

b = T*a can be characterized several ways in terms of a(n) and b(n) or their o.g.f.'s A(x) and B(x).

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

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

3) b(n) = Sum_{j=0..n} (n!/j!) a(j)

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

5) B(x) = Sum_{j=0,1,...} (xDx)^j A(x)

6) B(x) = Sum_{j=0,1,...} x^j * D^j * x^j A(x)

7) B(x) = Sum_{j=0,1,...} j! * x^j * L(j,-:xD:,0) A(x) 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 operator series at the j = n term gives an o.g.f. for b(0) through b(n).

c = (0!,1!,2!,3!,4!,...) is the sequence associated to T under the list partition transform and the associated operations described in A133314 so T(n,k) = binomial(n,k)*c(n-k). The reciprocal sequence is d = (1,-1,0,0,0,...). (End)

From Peter Bala, Jul 10 2008: (Start)

This array is the particular case P(1,1) of the generalized Pascal triangle P(a,b), a lower unit triangular matrix, shown below:

n\k|0.....................1...............2.......3......4

----------------------------------------------------------

0..|1.....................................................

1..|a....................1................................

2..|a(a+b)...............2a..............1................

3..|a(a+b)(a+2b).........3a(a+b).........3a........1......

4..|a(a+b)(a+2b)(a+3b)...4a(a+b)(a+2b)...6a(a+b)...4a....1

...

The entries A(n,k) of this array satisfy the recursion A(n,k) = (a+b*(n-k-1))*A(n-1,k) + A(n-1,k-1), which reduces to the Pascal formula when a = 1, b = 0.

Various cases are recorded in the database, including: P(1,0) = Pascal's triangle A007318, P(2,0) = A038207, P(3,0) = A027465, P(2,1) = A132159, P(1,3) = A136215 and P(2,3) = A136216.

When b <> 0 the array P(a,b) has e.g.f. exp(x*y)/(1-b*y)^(a/b) = 1 + (a+x)*y + (a*(a+b)+2a*x+x^2)*y^2/2! + (a*(a+b)*(a+2b) + 3a*(a+b)*x + 3a*x^2+x^3)*y^3/3! + ...; the array P(a,0) has e.g.f. exp((x+a)*y).

We have the matrix identities P(a,b)*P(a',b) = P(a+a',b); P(a,b)^-1 = P(-a,b).

An analog of the binomial expansion for the row entries of P(a,b) has been proved by [Echi]. Introduce a (generally noncommutative and nonassociative) product ** on the ring of polynomials in two variables by defining F(x,y)**G(x,y) = F(x,y)G(x,y) + by^2*d/dy(G(x,y)).

Define the iterated product F^(n)(x,y) of a polynomial F(x,y) by setting F^(1) = F(x,y) and F^(n)(x,y) = F(x,y)**F^(n-1)(x,y) for n >= 2. Then (x+a*y)^(n) = x^n + C(n,1)*a*x^(n-1)*y + C(n,2)*a*(a+b)*x^(n-2)*y^2 + ... + C(n,n)*a*(a+b)*(a+2b)*...*(a+(n-1)b)*y^n. (End)

(n+1) * n-th row = reversal of triangle A068424: (1; 2,2; 6,6,3; ...) - Gary W. Adamson, May 03 2009

Let G(m, k, p) = (-p)^k*Product_{j=0..k-1}(j - m - 1/p) and T(n,k,p) = G(n-1,n-k,p) then T(n, k, 1) is this sequence, T(n, k, 2) = A112292(n, k) and T(n, k, 3) = A136214. - Peter Luschny, Jun 01 2009, revised Jun 18 2019

The higher order exponential integrals E(x,m,n) are defined in A163931. For a discussion of the asymptotic expansions of the 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 - ...) see A130534. The asymptotic expansion of E(x,m=1,n) leads for n >= 1 to the left hand columns of the triangle given above. Triangle A165674 is generated by the asymptotic expansions of E(x,m=2,n). - Johannes W. Meijer, Oct 07 2009

T(n,k) = n!/k! = number of permutations of [n+1] with exactly k+1 cycles and with elements 1,2,...,k+1 in separate cycles. See link and example below. - Dennis P. Walsh, Jan 24 2011

T(n,k) is the number of n permutations that leave some size k subset of {1,2,...,n} fixed. Sum_{k=0..n}(-1)^k*T(n,k) = A000166(n) (the derangements). - Geoffrey Critzer, Dec 11 2011

T(n,k) = A162995(n-1,k-1), 2 <= k <= n; T(n,k) = A173333(n,k), 1 <= k <= n. - Reinhard Zumkeller, Jul 05 2012

The row polynomials form an Appell sequence. The matrix is a special case of a group of general matrices sketched in A132382. - Tom Copeland, Dec 03 2013

For interpretations in terms of colored necklaces, see A213936 and A173333. - Tom Copeland, Aug 18 2016

See A008279 for a relation of this entry to the e.g.f.s enumerating the faces of permutahedra and stellahedra. - Tom Copeland, Nov 14 2016

Also, T(n,k) is the number of ways to arrange n-k nonattacking rooks on the n X (n-k) chessboard. - Andrey Zabolotskiy, Dec 16 2016

The infinitesimal generator of this triangle is the generalized exponential Riordan array [-log(1-x), x] and equals the unsigned version of A238363. - Peter Bala, Feb 13 2017

Formulas for exponential and power series infinitesimal generators for this triangle T are given in Copeland's 2012 and 2014 formulas as T = unsigned exp[(I-A238385)] = 1/(I - A132440), where I is the identity matrix. - Tom Copeland, Jul 03 2017

If A(0) = 1/(1-x), and A(n) = d/dx(A(n-1)), then A(n) = n!/(1-x)^(n+1) = Sum_{k>=0} (n+k)!/k!*x^k = Sum_{k>=0} T(n+k, k)*x^k. - Michael Somos, Sep 19 2021

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)

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

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

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

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) = a(0), b(1) = a(n) - 2 a(0), b(n) = a(n) - 2n a(n-1) + n(n-1) a(n-2) for n > 0.

2) b(n) = n! Lag{n,(.)!*Lag[.,a1(.),0],-1}, umbrally, where a1(n) = n! Lag{n,(.)!*Lag[.,a(.),0],-1}.

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

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

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

6) B(x) = Sum_{j=0..2} {(-1)^j * binomial(2,j)*j!*x^j*Lag(j,-:xD:,0)} A(x)

where D is the derivative w.r.t. x, (:xD:)^j = x^j*D^j and Lag(n,x,m) is the associated Laguerre polynomial of order m.

7) EB(x) = (1-x)^2 EA(x)

8) T = S^2 = A132013^2 = A094587^(-2) = A132159^(-1).

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

The reciprocal sequence to c is d = (1!,2!,3!,4!,...), so the inverse of T is TI(n,k) = binomial(n,k)*d(n-k) = A132159.

These formulas are easily generalized for m applications of the basic operator n! Lag[n,(.)!*Lag[.,a(.),0],-1] by replacing 2 with m in formulas 3, 4, 5, 6 and 7.

The generalized c are given by the generalized coefficients of 6, i.e.,

c(n) = (-1)^n * binomial(m,n)*n! = (-1)^n * m!/(m-n)!.

The generalized d are given by the array at and below the term SI(m-1,m-1) in SI(n,k) = binomial(n,k) * (n-k)!, the inverse of S; i.e.,

d(n) = SI(m-1+n,m-1) = binomial(m-1+n,m-1) * n! = (m-1+n)!/(m-1)!.

As an aside, this shows that the signed falling factorials and the rising factorials form reciprocal pairs under the list partition transform of A133314.

Row sums of T = [formula 3 with all a(n) = 1] = [binomial transform of c] = [coefficients of B(x) with A(x) = 1/(1-x)] = (1,-1,-1,1,5,11,19,...),

with e.g.f. = [EB(x) with EA(x) = exp(x)] = (1-x)^2 * exp(x) = exp(x)*exp(c(.)*x) = exp[(1+c(.))*x].

Alternating row sums of T = [formula 3 with all a(n) = (-1)^n] = [finite differences of c] = [coefficients of B(x) with A(x) = 1/(1+x)] = (1,-3,7,-13,21,-31,...) = (-1)^n A002061(n+1),

with e.g.f. = [EB(x) with EA(x) = exp(-x)] = (1-x)^2 * exp(-x) = exp(- x)*exp(c(.)*x) = exp[-(1-c(.))*x].

See A132159 for a relation to the Poisson-Charlier polynomials. - Tom Copeland, Jan 15 2016

The associated Laguerre polynomials n!*Lag(n,3-n,-x) are related to the rook polynomials of a rectangular 3 X n chessboard by R(3,n,x) = n!*x^n*Lag(n,3-n,-1/x), which are also the matching polynomials, or generating function of the number of k-edge matchings, of the complete bipartite graph K(m,n), or biclique (cf. Wikipedia for details).

The formulas here and below can be naturally extended with 3 replaced by any positive integer m. For m = 1 and 2, see unsigned A132013 and A132014. The formulas there can be extrapolated to apply to this matrix.

Construct the infinite array of polynomials

a(0,t) = 1

a(1,t) = 1

a(2,t) = -1 + 3*t

a(3,t) = 1 - 8*t + 13*t^2

a(4,t) = 1 + 11*t - 61*t^2 + 73*t^3

a(5,t) = -19 + 66*t + 66*t^2 - 494*t^3 + 501*t^4

a(6,t) = 151 - 993*t + 2102*t^2 - 298*t^3 - 4293*t^4 + 4051*t^5

This array is the reciprocal array of the following array b(n,t) under the list partition transform and its associated operations described in A133314.

b(0,t) = 1 and b(n,t) = -A000262(n)*(t-1)^(n-1) for n > 0.

Then A111884(n) = a(n,0).

Lower triangular matrix A094587 = binomial(n,k)*a(n-k,1).

A084358(n) = a(n,2).

Signed A128229 = matrix inverse of binomial(n,k)*a(n-k,1) = binomial(n,k)*b(n-k,1) = A132013.

As t tends to infinity, a(n,t)/t^(n-1) tends to A000262(n) for n > 0.

The P(n,t) of A131758 can be constructed from T(n,k,t) = binomial(n,k)*a(n-k,t) by letting T(n,k,t) multiply the column vector c(n,t) given by c(0,t) = 0! and c(n,t) = n!*(t-1)^(n-1) for n > 0. The P(n,t) have rich associations to other sequences.