A105278 -id:A105278 - cp's OEIS frontend

A133314 Coefficients of list partition transform: reciprocal of an exponential generating function (e.g.f.).

1, -1, -1, 2, -1, 6, -6, -1, 8, 6, -36, 24, -1, 10, 20, -60, -90, 240, -120, -1, 12, 30, -90, 20, -360, 480, -90, 1080, -1800, 720, -1, 14, 42, -126, 70, -630, 840, -420, -630, 5040, -4200, 2520, -12600, 15120, -5040, -1, 16, 56, -168, 112, -1008, 1344, 70

Offset: 0

Tom Copeland, Oct 18 2007, Oct 29 2007, Nov 16 2007

The list partition transform of a sequence a(n) for which a(0)=1 is illustrated by:
b_0 = 1
b_1 = -a_1
b_2 = -a_2 + 2 a_1^2
b_3 = -a_3 + 6 a_2 a_1 - 6 a_1^3
b_4 = -a_4 + 8 a_3 a_1 + 6 a_2^2 - 36 a_2 a_1^2 + 24 a_1^4
... .
The unsigned coefficients are A049019 with a leading 1. The sign is dependent on the partition as evident from inspection (replace a_n's by -1).
Expressed umbrally, i.e., with the umbral operation (a.)^n := a_n,
  exp(a.x) exp(b.x) = exp[(a.+b.)x] = 1; i.e., (a.+b.)^n = 1 for n=0 and 0 for all other values of n.
Expressed recursively,
b_0 = 1, b_n = -Sum_{j=1..n} binomial(n,j) a_j b_{n-j}; which is conditionally self-inverse, i.e., the roles of a_k and b_k may be reversed with a_0 = b_0 = 1.
Expressed in matrix form, b_n form the first column of B = matrix inverse of A .
A = Pascal matrix diagonally multiplied by a_n, i.e., A_{n,k} = binomial(n,k)* a_{n-k}.
Some examples of reciprocal pairs of sequences under these operations are:
1) A084358 and -A000262 with the first term set to 1.
2) (1,-1,0,0,...) and (0!,1!,2!,3!,...) with the unsigned associated matrices A128229 and A094587.
3) (1,-1,-1,-1,...) and A000670.
4) A000110 and A000587.
5) (1,-2,-2,0,0,0,...) and (0! c_1,1! c_2,2! c_3,3! c_4,...) where c_n = A000129(n) with the associated matrices A110327 and A110330.
6) (1,-2,2,0,0,0,...) and (1!,2!,3!,4!,...).
7) Sequences of rising and signed lowering factorials form reciprocal pairs where a_n = (-1)^n m!/(m-n)! and b_n = (m-1+n)!/(m-1)! for m=0,1,2,... .
Denote the action of the list partition transform on the sequence a. or an invertible matrix M by LPT(a.) = b. or LPT(M)= M^(-1).
If the matrix equation M = exp(T) also holds, then exp[a.*T]*exp[b.*T] = exp[(a.+b.)*T] = I, the identity matrix, because (a.+b.)^n = delta_n, the Kronecker delta with delta_n = 1 and delta_n = 0 otherwise, i.e., (0)^n = delta_n.
Therefore, [exp(a.*T)]^(-1) = exp[b.*T] = exp[LPT(a.)*T] = LPT[exp(a.*T)].
The fundamental Pascal (A007318), unsigned Lah (A105278) and associated Laguerre matrices can be generated by exponentiation of special infinitesimal matrices (see A132440, A132710 and A132681) such that finding LPT(a.) amounts to multiplying the k'th diagonal of the fundamental matrices by a_k for every diagonal followed by matrix inversion and then extraction of the b_n factors from the first column (simplest for the Pascal formulas above).
Conversely, the inverses of matrices formed by diagonally multiplying the three fundamental matrices by a_k are given by diagonally multiplying the fundamental matrices by b_k.
If LPT(M) is defined differently as application of the top formula to a_n = M^n, then b_n = (-M)^n and the formalism could even be applied to more general sequences of matrices M., providing the reciprocal of exp[t*M.].
The group of fundamental lower triangular matrices M = exp(T) such that LPT[exp(a.*T)] = exp[LPT(a.)*T] = [exp[a.*T]]^(-1) are obtained by infinitesimal generator matrices of the form T =
  0;
  t(0), 0;
  0, t(1), 0;
  0, 0, t(2), 0;
  0, 0, 0, t(3), 0;
  ... .
T^m has trivially vanishing terms except along the m'th subdiagonal, which is a sequence of generalized factorials:
[ t(0)*t(1)...t(m-2)*t(m-1), t(1)*t(2)...t(m-1)*t(m), t(2)*t(3)...t(m)*t(m+1), ... ].
Therefore the principal submatrices of T (given by setting t(j) = 0 for j > n-1) are nilpotent with at least [Tsub_n]^(n+1) = 0.
The general group of matrices GM[a.] = exp[a.*T] can also be obtained through diagonal multiplication of M = exp(T) by the sequence a_n, as in the Pascal matrix example above and their inverses by diagonal multiplication by b. = LPT(a.).
Weighted-mappings interpretation for the top partition equation:
Given n pre-nodes (Pre) and k post-nodes (Post), each Pre is connected to only one Post and each Post has at least one Pre connected to it (surjections or onto functions/maps). Weight each Post by -a_m where m is the number of connections to the Post.
Weight each map by the product of the Post weights and multiply by the number of maps that share the same connectivity. Sum over the possible mappings for n Pre. The result is b_n.
E.g., b_3 = [ 3 Pre to 1 Post ] + [ 3 Pre to 2 Post ] + [ 3 Pre to 3 Post ]
= [1 map with 1 Post with 3 connections] + [ 6 maps with 1 Post with 2 connections and 1 Post with 1 connection] + [6 maps with 3 Post with 1 connection each]
= -a_3 + 6 * [-a_2*(-a_1)] + 6 * [-a_1*(-a_1)*(-a_1)].
See A263633 for the complementary formulation for the reciprocal of o.g.f.s rather than e.g.f.s and computations of these partition polynomials as Gram determinants. - Tom Copeland, Dec 04 2016
The coefficients of the partition polynomials enumerate the faces of the convex, bounded polytopes called permutohedra, and the absolute value of the sum of the 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 unity. 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 the fundamental matrix chosen to be the lower triangular Pascal matrix M, the matrix MA whose n-th diagonals are multiplied by a_n (i.e., MA_{i,j} = PM_{i,j} * a_{i-j}) gives a matrix representation of the e.g.f. associated to the Appell polynomial sequence defined by e^{a.t}e^{xt}= e^{(a.+x)t} = e^{A.(x)t} where umbrally (A.(x))^n = A_n(x) = (a. + x)^n = sum_{k=0..n} binomial(n,k) a_k x^{n-k} are the associated Appell polynomials. Left multiplication of the column vector (1,x,x^2,..) by MA gives the Appell polynomial sequence, and multiplication of the two e.g.f.s e^{a.t} and e^{b.t} corresponds to multiplication of their respective matrix representations MA and MB. Forming the reciprocal of an e.g.f. corresponds to taking the matrix inverse of its matrix representation as noted above. A263634 gives an associated modified Pascal matrix representation of the raising operator for the Appell sequence. - Tom Copeland, Nov 13 2019
The diagonal of MA consists of all ones. Let MAN be the truncated square submatrix of MA containing the coefficients of the first N Appell polynomials A_k=(a.+x)^k = Sum(j=0 to k) MAN(k,j) x^j. Then by the Cayley-Hamilton theorem (I-MAN)^N = 0; therefore, MAN^(-1) = Sum(k=1 to N) binomial(N,k) (-MAN)^{k-1} = MBN, the inverse of MAN, containing the coefficients of the first N rows of the Appell polynomials B_k(x) = (b. + x)^k = Sum(j=0 to k) MBN(k,j) x^j, which are the umbral compositional inverses of the Appell row polynomials A_k(x) of MAN; that is, A_k(B.(x)) = x^k = B_k(A.(x)), where, e.g., (A.(x))^k = A_k(x). - Tom Copeland, May 13 2020
The use of the term 'list partition transform' resulted from one of my first uses of these partition polynomials in relating A000262 to A084358 with their simple e.g.f.s. Other appropriate names would be the permutohedra polynomials since they are refined Euler characteristics of the permutohedra or the reciprocal polynomials since they give the multiplicative inverses of e.g.f.s with a constant of 1. - Tom Copeland, Oct 09 2022

			Table starts:
[0] [ 1]
[1] [-1]
[2] [-1,  2]
[3] [-1,  6, -6]
[4] [-1,  8,  6, -36,  24]
[5] [-1, 10, 20, -60, -90,  240, -120]
[6] [-1, 12, 30, -90,  20, -360,  480, -90, 1080, -1800, 720]

Vincenzo Librandi, Table of n, a(n) for n = 0..250
M. Aguiar, Math 7410: Lie Combinatorics and Hyperplane Arrangements, notes made by D. Mehrle for course taught by Aguiar at Cornell, 2016, p. 49.
M. Aguiar and F. Ardila, The algebraic and combinatorial structure of generalized permutahedra [sic], MSRI Summer School July 19, 2017.
M. Aguiar and F. Ardila, Hopf monoids and generalized permutahedra [sic], arXiv:1709.07504 [math.CO], p. 5, 2017.
N. Arkani-Hamed, Y. Bai, S. He, and G. Yan, Scattering forms and the positive geometry of kinematics, color, and the worldsheet , arXiv:1711.09102 [hep-th], 2017.
J. Bergström and S. Minabe, On the cohomology of the Losev-Manin moduli space, arXiv:1108.0338 [math.AG], (cf. p. 1), 2013.
Z. Bern, J.Carrasco, M. Chiodaroli, H. Johansson, and R. Roiban , The Duality Between Color and Kinematics and its Applications, arXiv preprint arXiv:1909.01358 [hep-th], 2019.
L. Berry, S. Forcey, M. Ronco, and P. Showers, Polytopes and Hopf algebras of painted trees: Fan graphs and Stellohedra, arXiv:1608.08546 [math.CO], 2018.
Tom Copeland, Bijective mapping between face polytopes of permutohedra and partitions of integers, Math StackExchange question, 2016.
Tom Copeland, Compilation of OEIS Partition Polynomials A133314, A134685, A145271, A356144, and A356145, 2022.
Tom Copeland, In the Realm of Shadows: Umbral inverses and associahedra, noncrossing partitions, symmetric polynomials, and similarity transforms, 2019.
Karl-Dieter Crisman, The Borda Count, the Kemeny Rule, and the permutahedron [sic], preprint, 2014.
Karl-Dieter Crisman, The Borda Count, the Kemeny Rule, and the permutahedron [sic], in: Karl-Dieter Crisman and Michael A. Jones (eds.), The Mathematics of Decisions, Elections, and Games, Contemporary Mathematics, AMS, Vol. 624, 2014, pp. 101-134.
R. Da Rosa, D. Jensen, and D. Ranganathan, Toric graph associahedra and compactifications of M_(0,n), arXiv:1411.0537 [math.AG], 2015.
Nick Early, Honeycomb tessellations and canonical bases for permutohedral blades, arXiv:1810.03246 [math.CO], 2018.
S. Forcey, The Hedra Zoo
S. Franco and A. Hasan, Graded Quivers, Generalized Dimer Models and Toric Geometry , arXiv preprint arXiv:1904.07954 [hep-th], 2019.
M. Futaki and K. Ueda, Tropical coamoeba and torus-equivariant homological mirror symmetry for the projective space , arXiv preprint arXiv:1001.4858 [math.SG], 2010-2014.
X. Gao, S. He, Y. Zhan, Labelled tree graphs, Feynman diagrams and disk integrals , arXiv:1708.08701 [hep-th], 2017.
A. Hasan, Physics and Mathematics of Graded Quivers, dissertation, Graduate Center, City University of New York, 2019.
D. Karp, D. Ranganathan, P. Riggins, and U. Whitcher, Toric Symmetry of CP^3, arXiv:1109.5157 [math.AG], 2011.
R. Kaufmann and Y. Zhang, Permutohedral structures on E2-operads, arXiv preprint arXiv:1602.08247 [math.AT], 2016.
M. Lin, Graph Cohomology, 2016.
J. Loday, The Multiple Facets of the Associahedron
MathOverflow, Analogue of conic sections for the permutohedra, associahedra, and noncrossing partitions, an MO question posed by T. Copeland, 2017.
W. Norledge and A. Ocneanu, Hopf monoids, permutohedral cones, and generalized retarded functions, arXiv preprint arXiv:1911.11736 [math.CO], 2020.
P. Olver, The canonical contact form, p. 9.
V. Pilaud, The Associahedron and its Friends, presentation for Séminaire Lotharingien de Combinatoire, April 4 - 6, 2016.
J . Pitman and R. Stanley, A polytope related to empirical distributions, plane trees, parking functions, and the associahedron, arxiv: 9908029 [math.CO], 1999.
A. Postnikov, Positive Grassmannian and Polyhedral Subdivisions, arXiv:1806.05307 [math.CO], (cf. p. 17), 2018.
D. Ranganathan, Gromov-Witten Theory of Blowups of Toric Threefolds, senior thesis, Harvey Mudd College, 2012.
D. Ranganathan, Gromov-Witten Theory of Blowups of Toric Threefolds (poster), poster for senior thesis, Harvey Mudd College, 2012.
E. Schröder, Ueber unendlich viele Algorithmen zur Auflösung der Gleichungen, Mathematische Annalen vol. 2, 317-365, 1870, (bottom of p. 342, see Stewart below for a translation).
G. Stewart, On infinitely many algorithms for solving equations, 1993, (translation into English of the Schröder paper above, see top of p. 31 for differently normalized partition polynomials of this entry).

Row lengths in A000041.
Cf. A000110, A000129, A000262, A000587, A000670, A001286, A007318, A019538, A049019, A084358, A094587, A105278, A110327, A110330, A128229, A132440, A132681, A132710, A263633.
Cf. A134685, A145271, A356144, A356145.

b[0] = 1; b[n_] := b[n] = -Sum[Binomial[n, j]*a[j]*b[n-j], {j, 1, n}];
row[0] = {1}; row[n_] := Coefficient[b[n], #]& /@ (Times @@ (a /@ #)&) /@ IntegerPartitions[n];
Table[row[n], {n, 0, 8}] // Flatten (* Jean-François Alcover, Apr 23 2014 *)

def A133314_row(n): return [(-1)^len(s)*factorial(len(s))*SetPartitions(sum(s), s).cardinality() for s in Partitions(n)]
for n in (0..10): print(A133314_row(n)) # Peter Luschny, Sep 18 2015

b_{n-1} = (1/n)(d/da(1))p_n[a_1, a_2, ..., a_n] where p_n are the row partition polynomials of the cumulant generator A127671. - Tom Copeland, Oct 13 2012
(E.g.f. of matrix B) = (e.g.f. of b)·exp(xt) = exp(b.t)·exp(xt) = exp(xt)/exp(a.t) = (e.g.f. of A^(-1)) and (e.g.f. of matrix A) = exp(a.t)·exp(xt) = exp(xt)/exp(b.t) = (e.g.f. of B^(-1)), where the umbral evaluation of exp(b.t) = Sum{n >= 0} (b.t)^n / n! = Sum_{n >= 0} b_n t^n / n! is understood in the denominator. These e.g.f.s define Appell sequences of polynomials. - Tom Copeland, Mar 22 2014
Sum of the n-th row is (-1)^n. - Peter Luschny, Sep 18 2015
The unsigned coefficients for the partitions a_2*a_1^n for n >= 0 are the Lah numbers A001286. - Tom Copeland, Aug 06 2016
G.f.: 1 / (1 + Sum_{n > 0} a_n x^n/n!) = 1 / exp(a.x). - Tom Copeland, Oct 18 2016
Let a_1 = 1 + x + B_1 = x + 1/2 and a_n = B_n = (B.)^n, where B_n are the Bernoulli numbers defined by e^(B.t) = t / (e^t-1), then t / e^(a.t) = t / [(x + 1) * t + exp(B.t)] =  (e^t - 1) /[ 1 + (x + 1) (e^t - 1)] = exp(p.(x)t), where (p.(x))^n = p_n(x) are the shifted signed polynomials of A019538: p_0(x) = 0, p_1(x) = 1, p_2(x) = -(1 + 2 x), p_3(x) = 1 + 6 x + 6 x^2, ... , p_n(x) = n * b_{n-1}. - Tom Copeland, Oct 18 2016
With a_n = 1/(n+1), b_n = B_n, the Bernoulli numbers. - Tom Copeland, Nov 08 2016
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 A145271; [E]^(-1), A356145, the inverse set to [E]; [P], the permutohedra polynomials of this entry; [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

More terms from Jean-François Alcover, Apr 23 2014

A035342 The convolution matrix of the double factorial of odd numbers (A001147).

Original entry on oeis.org

1, 3, 1, 15, 9, 1, 105, 87, 18, 1, 945, 975, 285, 30, 1, 10395, 12645, 4680, 705, 45, 1, 135135, 187425, 82845, 15960, 1470, 63, 1, 2027025, 3133935, 1595790, 370125, 43890, 2730, 84, 1, 34459425, 58437855, 33453945, 8998290

Offset: 1

Wolfdieter Lang

Previous name was: A triangle of numbers related to the triangle A035324; generalization of Stirling numbers of second kind A008277 and Lah numbers A008297.
If one replaces in the recurrence the '2' by '0', resp. '1', one obtains the Lah-number, resp. Stirling-number of 2nd kind, triangle A008297, resp. A008277.
The product of two lower triangular Jabotinsky matrices (see A039692 for the Knuth 1992 reference) is again such a Jabotinsky matrix: J(n,m) = Sum_{j=m..n} J1(n,j)*J2(j,m). The e.g.f.s of the first columns of these triangular matrices are composed in the reversed order: f(x)=f2(f1(x)). With f1(x)=-(log(1-2*x))/2 for J1(n,m)=|A039683(n,m)| and f2(x)=exp(x)-1 for J2(n,m)=A008277(n,m) one has therefore f2(f1(x))=1/sqrt(1-2*x) - 1 = f(x) for J(n,m)=a(n,m). This proves the matrix product given below. The m-th column of a Jabotinsky matrix J(n,m) has e.g.f. (f(x)^m)/m!, m>=1.
a(n,m) gives the number of forests with m rooted ordered trees with n non-root vertices labeled in an organic way. Organic labeling means that the vertex labels along the (unique) path from the root with label 0 to any leaf (non-root vertex of degree 1) is increasing. Proof: first for m=1 then for m>=2 using the recurrence relation for a(n,m) given below. - Wolfdieter Lang, Aug 07 2007
Also the Bell transform of A001147(n+1) (adding 1,0,0,... as column 0). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 19 2016

			Matrix begins:
    1;
    3,   1;
   15,   9,   1;
  105,  87,  18,   1;
  945, 975, 285,  30,   1;
  ...
Combinatoric meaning of a(3,2)=9: The nine increasing path sequences for the three rooted ordered trees with leaves labeled with 1,2,3 and the root labels 0 are: {(0,3),[(0,1),(0,2)]}; {(0,3),[(0,2),(0,1)]}; {(0,3),(0,1,2)}; {(0,1),[(0,3),(0,2)]}; [(0,1),[(0,2),(0,3)]]; [(0,2),[(0,1),(0,3)]]; {(0,2),[(0,3),(0,1)]}; {(0,1),(0,2,3)}; {(0,2),(0,1,3)}.

Reinhard Zumkeller, Rows n = 1..125 of triangle, flattened
Peter Bala, Generalized Dobinski formulas
J. Fernando Barbero G., Jesús Salas, and Eduardo J. S. Villaseñor, Bivariate Generating Functions for a Class of Linear Recurrences. I. General Structure, arXiv:1307.2010 [math.CO], 2013-2014.
Pawel Blasiak, Karol A. Penson, and Allan I. Solomon, The Boson Normal Ordering Problem and Generalized Bell Numbers, arXiv:quant-ph/0212072, 2002.
Pawel Blasiak, Karol A. Penson, and Allan I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.
Richell O. Celeste, Roberto B. Corcino and Ken Joffaniel M. Gonzales. Two Approaches to Normal Order Coefficients. Journal of Integer Sequences, Vol. 20 (2017), Article 17.3.5.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers, 2015.
Tom Copeland, Addendum to Mathemagical Forests, 2010.
Tom Copeland, Mathemagical Forests, 2008.
Askar Dzhumadildaev and Damir Yeliussizov, Path decompositions of digraphs and their applications to Weyl algebra, arXiv preprint arXiv:1408.6764v1 [math.CO], 2014. [Version 1 contained many references to the OEIS, which were removed in Version 2. - _N. J. A. Sloane_, Mar 28 2015]
Askar Dzhumadil'daev and Damir Yeliussizov, Walks, partitions, and normal ordering, Electronic Journal of Combinatorics, 22(4) (2015), #P4.10.
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) #09.8.3.
Donald E. Knuth, Convolution polynomials, arXiv:math/9207221 [math.CA], 1992; Mathematica J. 2.1 (1992), no. 4, 67-78.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
Wolfdieter Lang, First 10 rows.
Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104v2 [math.CO], 2012.
Toufik Mansour, Matthias Schork and Mark Shattuck, The Generalized Stirling and Bell Numbers Revisited, Journal of Integer Sequences, Vol. 15 (2012), #12.8.3.
Erich Neuwirth, Recursively defined combinatorial functions: Extending Galton's board, Discrete Math. (2001) Vol. 239, 33-51.
Mathias Pétréolle and Alan D. Sokal, Lattice paths and branched continued fractions. II. Multivariate Lah polynomials and Lah symmetric functions, arXiv:1907.02645 [math.CO], 2019.
Umesh Shankar, Log-concavity of rows of triangular arrays satisfying a certain super-recurrence, arXiv:2508.12467 [math.CO], 2025. See p. 4.

The column sequences are A001147, A035101, A035119, ...
Row sums: A049118(n), n >= 1.
Cf. A000108, A035324, A008277, A008297, A094638.

a035342 n k = a035342_tabl !! (n-1) !! (k-1)
a035342_row n = a035342_tabl !! (n-1)
a035342_tabl = map fst $ iterate (\(xs, i) -> (zipWith (+)
   ([0] ++ xs) $ zipWith (*) [i..] (xs ++ [0]), i + 2)) ([1], 3)
-- Reinhard Zumkeller, Mar 12 2014

T := (n,k) -> 2^(k-n)*hypergeom([k-n,k+1],[k-2*n+1],2)*GAMMA(2*n-k)/
(GAMMA(k)*GAMMA(n-k+1)); for n from 1 to 9 do seq(simplify(T(n,k)),k=1..n) od; # Peter Luschny, Mar 31 2015
T := (n, k) -> local j; 2^n*add((-1)^(k-j)*binomial(k, j)*pochhammer(j/2, n), j = 1..k)/k!: for n from 1 to 6 do seq(T(n, k), k=1..n) od;  # Peter Luschny, Mar 04 2024

a[n_, k_] := 2^(n+k)*n!/(4^n*n*k!)*Sum[(j+k)*2^(j)*Binomial[j + k - 1, k-1]*Binomial[2*n - j - k - 1, n-1], {j, 0, n-k}]; Flatten[Table[a[n,k], {n, 1, 9}, {k, 1, n}] ] [[1 ;; 40]] (* Jean-François Alcover, Jun 01 2011, after Vladimir Kruchinin *)

a(n,k):=2^(n+k)*n!/(4^n*n*k!)*sum((j+k)*2^(j)*binomial(j+k-1,k-1)*binomial(2*n-j-k-1,n-1),j,0,n-k); /* Vladimir Kruchinin, Mar 30 2011 */

# uses[bell_matrix from A264428]
# Adds a column 1,0,0,0, ... at the left side of the triangle.
print(bell_matrix(lambda n: A001147(n+1), 9)) # Peter Luschny, Jan 19 2016

a(n, m) = Sum_{j=m..n} |A039683(n, j)|*S2(j, m) (matrix product), with S2(j, m) := A008277(j, m) (Stirling2 triangle). Priv. comm. to Wolfdieter Lang by E. Neuwirth, Feb 15 2001; see also the 2001 Neuwirth reference. See the comment on products of Jabotinsky matrices.

a(n, m) = n!*A035324(n, m)/(m!*2^(n-m)), n >= m >= 1; a(n+1, m)= (2*n+m)*a(n, m)+a(n, m-1); a(n, m) := 0, n

E.g.f. of m-th column: ((x*c(x/2)/sqrt(1-2*x))^m)/m!, where c(x) = g.f. for Catalan numbers A000108.

From Vladimir Kruchinin, Mar 30 2011: (Start)

G.f. (1/sqrt(1-2*x) - 1)^k = Sum_{n>=k} (k!/n!)*a(n,k)*x^n.

a(n,k) = 2^(n+k) * n! / (4^n*n*k!) * Sum_{j=0..n-k} (j+k) * 2^(j) * binomial(j+k-1,k-1) * binomial(2*n-j-k-1,n-1). (End)

From Peter Bala, Nov 25 2011: (Start)

E.g.f.: G(x,t) = exp(t*A(x)) = 1 + t*x + (3*t + t^2)*x^2/2! + (15*t + 9*t^2 + t^3)*x^3/3! + ..., where A(x) = -1 + 1/sqrt(1-2*x) satisfies the autonomous differential equation A'(x) = (1+A(x))^3.

The generating function G(x,t) satisfies the partial differential equation t*(dG/dt+G) = (1-2*x)*dG/dx, from which follows the recurrence given above.

The row polynomials are given by D^n(exp(x*t)) evaluated at x = 0, where D is the operator (1+x)^3*d/dx. Cf. A008277 (D = (1+x)*d/dx), A105278 (D = (1+x)^2*d/dx), A035469 (D = (1+x)^4*d/dx) and A049029 (D = (1+x)^5*d/dx). (End)

The n-th row polynomial R(n,x) is given by the Dobinski-type formula R(n,x) = exp(-x)*Sum_{k>=1} k*(k+2)*...*(k+2*n-2)*x^k/k!. - Peter Bala, Jun 22 2014

T(n,k) = 2^(k-n)*hypergeom([k-n,k+1],[k-2*n+1],2)*Gamma(2*n-k)/(Gamma(k)*Gamma(n-k+1)). - Peter Luschny, Mar 31 2015

T(n,k) = 2^n*Sum_{j=1..k} ((-1)^(k-j)*binomial(k, j)*Pochhammer(j/2, n)) / k!. - Peter Luschny, Mar 04 2024

Simpler name from Peter Luschny, Mar 31 2015

A021009 Triangle of coefficients of Laguerre polynomials n!*L_n(x) (rising powers of x).

Original entry on oeis.org

1, 1, -1, 2, -4, 1, 6, -18, 9, -1, 24, -96, 72, -16, 1, 120, -600, 600, -200, 25, -1, 720, -4320, 5400, -2400, 450, -36, 1, 5040, -35280, 52920, -29400, 7350, -882, 49, -1, 40320, -322560, 564480, -376320, 117600, -18816, 1568, -64, 1, 362880, -3265920

Offset: 0

N. J. A. Sloane

In absolute values, this sequence also gives the lower triangular readout of the exponential of a matrix whose entry {j+1,j} equals (j-1)^2 (and all other entries are zero). - Joseph Biberstine (jrbibers(AT)indiana.edu), May 26 2006
A partial permutation on a set X is a bijection between two subsets of X. |T(n,n-k)| equals the numbers of partial permutations of an n-set having domain cardinality equal to k. Let E denote the operator D*x*D, where D is the derivative operator d/dx. Then E^n = Sum_{k = 0..n} |T(n,k)|*x^k*D^(n+k). - Peter Bala, Oct 28 2008
The unsigned triangle is the generalized Riordan array (exp(x), x) with respect to the sequence n!^2 as defined by Wang and Wang (the generalized Riordan array (exp(x), x) with respect to the sequence n! is Pascal's triangle A007318, and with respect to the sequence n!*(n+1)! is A105278). - Peter Bala, Aug 15 2013
The unsigned triangle appears on page 83 of Ser (1933). - N. J. A. Sloane, Jan 16 2020

			The triangle a(n,m) starts:
n\m   0       1      2       3      4      5    6  7  8
0:    1
1:    1      -1
2:    2      -4      1
3:    6     -18      9      -1
4:   24     -96     72     -16      1
5:  120    -600    600    -200     25     -1
6:  720   -4320   5400   -2400    450    -36    1
7: 5040  -35280  52920  -29400   7350   -882   49  -1
8:40320 -322560 564480 -376320 117600 -18816 1568 -64 1
...
From _Wolfdieter Lang_, Jan 31 2013 (Start)
Recurrence (usual one): a(4,1) = 7*(-18) - 6 - 3^2*(-4) = -96.
Recurrence (simplified version): a(4,1) = 5*(-18) - 6 = -96.
Recurrence (Sage program): |a(4,1)| = 6 + 3*18 + 4*9 = 96. (End)
Embedded recurrence (Maple program): a(4,1) = -4!*(1 + 3) = -96.

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 799.
G. Rota, Finite Operator Calculus, Academic Press, New York, 1975.
J. Ser, Les Calculs Formels des Séries de Factorielles. Gauthier-Villars, Paris, 1933, p. 83.

T. D. Noe, Rows n = 0..50 of triangle, flattened
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
W. A. Al-Salam, Operational representations for the Laguerre and other polynomials, Duke Math. Jour., vol 31 (1964), pp. 127-142.
Paul Barry, The Restricted Toda Chain, Exponential Riordan Arrays, and Hankel Transforms, J. Int. Seq. 13 (2010) # 10.8.4, example 5.
Paul Barry, Exponential Riordan Arrays and Permutation Enumeration, J. Int. Seq. 13 (2010) # 10.9.1, example 7.
Paul Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, example 21.
Paul Barry, Combinatorial polynomials as moments, Hankel transforms and exponential Riordan arrays, arXiv preprint arXiv:1105.3044 [math.CO], 2011, also J. Int. Seq. 14 (2011) # 11.6.7.
Paul Barry, On a transformation of Riordan moment sequences, arXiv:1802.03443 [math.CO], 2018.
A. Belov-Kanel and M. Kontsevich, Automorphisms of the Weyl algebra, arXiv preprint arXiv:0512169 [math.QA], 2005.
A. Belov-Kanel and M. Kontsevich, The Jacobian Conjecture is stably equivalent to the Dixmier Conjecture, arXiv preprint arXiv:0512171 [math.RA], 2005.
I. Gessel, Applications of the classical umbral calculus
G. Hetyei, Meixner polynomials of the second kind and quantum algebras representing su(1,1), arXiv preprint arXiv:0909.4352 [math.QA], 2009, p. 4.
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3.
Robert S. Maier, Boson Operator Ordering Identities from Generalized Stirling and Eulerian Numbers, arXiv:2308.10332 [math.CO], 2023. See. p. 19.
Massimo Nocentini, An algebraic and combinatorial study of some infinite sequences of numbers supported by symbolic and logic computation, PhD Thesis, University of Florence, 2019. See p. 31.
J. Ser, Les Calculs Formels des Séries de Factorielles, Gauthier-Villars, Paris, 1933 [Local copy].
J. Ser, Les Calculs Formels des Séries de Factorielles (Annotated scans of some selected pages)
M. Z. Spivey, On Solutions to a General Combinatorial Recurrence, J. Int. Seq. 14 (2011) # 11.9.7.
W. Wang and T. Wang, Generalized Riordan arrays, Discrete Mathematics, Vol. 308, No. 24, 6466-6500.
Eric Weisstein's World of Mathematics, Laguerre Polynomial
Index entries for sequences related to Laguerre polynomials

Row sums give A009940, alternating row sums are A002720.
Column sequences (unsigned): A000142, A001563, A001809-A001812 for m=0..5.
Central terms: A295383.
For generators and generalizations see A132440.
Cf. A021010, A025166, A025167, A062137, A062138, A062139, A062140, A066667, A321620.

/* As triangle: */ [[((-1)^k)*Factorial(n)*Binomial(n, k)/Factorial(k): k in [0..n]]: n in [0.. 10]]; // Vincenzo Librandi, Jan 18 2020

A021009 := proc(n,k) local S; S := proc(n,k) option remember; `if`(k = 0, 1, `if`( k > n, 0, S(n-1,k-1)/k + S(n-1,k))) end: (-1)^k*n!*S(n,k) end: seq(seq(A021009(n,k), k=0..n), n=0..8); # Peter Luschny, Jun 21 2017
# Alternative for the unsigned case (function RiordanSquare defined in A321620):
RiordanSquare(add(x^m, m=0..10), 10, true); # Peter Luschny, Dec 06 2018

Flatten[ Table[ CoefficientList[ n!*LaguerreL[n, x], x], {n, 0, 9}]] (* Jean-François Alcover, Dec 13 2011 *)

p(n) = denominator(bestapprPade(Ser(vector(2*n, k, (k-1)!))));
concat(1, concat(vector(9, n, Vec(-p(n)))))  \\ Gheorghe Coserea, Dec 01 2016

{T(n, k) = if( n<0, 0, n! * polcoeff( sum(i=0, n, binomial(n, n-i) * (-x)^i / i!), k))}; /* Michael Somos, Dec 01 2016 */

row(n) = Vecrev(n!*pollaguerre(n)); \\ Michel Marcus, Feb 06 2021

def A021009_triangle(dim): # computes unsigned T(n,k).
    M = matrix(ZZ,dim,dim)
    for n in (0..dim-1): M[n,n] = 1
    for n in (1..dim-1):
        for k in (0..n-1):
            M[n,k] = M[n-1,k-1]+(2*k+1)*M[n-1,k]+(k+1)^2*M[n-1,k+1]
    return M
A021009_triangle(9) # Peter Luschny, Sep 19 2012

a(n, m) = ((-1)^m)*n!*binomial(n, m)/m! = ((-1)^m)*((n!/m!)^2)/(n-m)! if n >= m, otherwise 0.
E.g.f. for m-th column: (-x/(1-x))^m /((1-x)*m!), m >= 0.
Representation (of unsigned a(n, m)) as special values of Gauss hypergeometric function 2F1, in Maple notation: n!*(-1)^m*hypergeom([ -m, n+1 ], [ 1 ], 1)/m!. - Karol A. Penson, Oct 02 2003
Sum_{m>=0} (-1)^m*a(n, m) = A002720(n). - Philippe Deléham, Mar 10 2004
E.g.f.: (1/(1-x))*exp(x*y/(x-1)). - Vladeta Jovovic, Apr 07 2005
Sum_{n>=0, m>=0} a(n, m)*(x^n/n!^2)*y^m = exp(x)*BesselJ(0, 2*sqrt(x*y)). - Vladeta Jovovic, Apr 07 2005
Matrix square yields the identity matrix: L^2 = I. - Paul D. Hanna, Nov 22 2008
From Tom Copeland, Oct 20 2012: (Start)
Symbolically, with D=d/dx and LN(n,x)=n!L_n(x), define :Dx:^j = D^j x^j, :xD:^j = x^j D^j, and LN(.,x)^j = LN(j,x) = row polynomials of A021009.
Then some useful relations are
1) (:Dx:)^n = LN(n,-:xD:)    [Rodriguez formula]
2) (xDx)^n = x^n D^n x^n = x^n LN(n,-:xD:)  [See Al-Salam ref./A132440]
3) (DxD)^n = D^n x^n D^n = LN(n,-:xD:) D^n  [See ref. in A132440]
4) umbral composition LN(n,LN(.,x))= x^n   [See Rota ref.]
5) umbral comp. LN(n,-:Dx:) = LN(n,-LN(.,-:xD:)) = 2^n LN(n,-:xD:/2)= n! * (n-th row e.g.f.(x) of A038207 with x replaced by :xD:).
An example for 2) is the operator (xDx)^2 = (xDx)(xDx) = xD(x^2 + x^3D)= 2x^2 + 4x^3 D + x^4 D^2 = x^2 (2 + 4x D + x^2 D^2) = x^2 (2 + 4 :xD: + :xD:^2) = x^2 LN(2,-:xD:) = x^2 2! L_2(-:xD:).
An example of the umbral composition in 5) is given in A038207.
The op. xDx is related to the Euler/binomial transformation for power series/o.g.f.s. through exp(t*xDx) f(x) = f[x/(1-t*x)]/(1-t*x) and to the special Moebius/linear fractional/projective transformation z exp(-t*zDz)(1/z)f(z) = f(z/(1+t*z)).
For a general discussion of umbral calculus see the Gessel link. (End)
From Wolfdieter Lang, Jan 31 2013: (Start)
Standard recurrence derived from the three term recurrence of the orthogonal polynomials system {n!*L(n,x)}: L(n,x) = (2*n  - 1 - x)*L(n-1,x) - (n-1)^2*L(n-2,x), n>=1, L(-1,x) = 0, L(0,x) = 1.
  a(n,m) = (2*n-1)*a(n-1,m) - a(n-1,m-1) - (n-1)^2*a(n-2,m),
  n >=1, with a(n,-1) = 0, a(0,0) = 1, a(n,m) = 0 if n < m. (compare this with Peter Luschny's program for the unsigned case |a(n,m)| = (-1)^m*a(n,m)).
Simplified recurrence (using column recurrence from explicit form for a(n,m) given above):
a(n,m) = (n+m)*a(n-1,m) - a(n-1,m-1), n >= 1, a(0,0) = 1, a(n,-1) = 0, a(n,m) = 0 if n < m. (End)
|T(n,k)| = [x^k] (-1)^n*U(-n,1,-x), where U(a,b,x) is Kummer's hypergeometric U function. - Peter Luschny, Apr 11 2015
T(n,k) = (-1)^k*n!*S(n,k) where S(n,k) is recursively defined by: "if k = 0 then 1 else if k > n then 0 else S(n-1,k-1)/k + S(n-1,k)". - Peter Luschny, Jun 21 2017
The unsigned case is the exponential Riordan square (see A321620) of the factorial numbers. - Peter Luschny, Dec 06 2018
Omitting the diagonal and signs, this array is generated by the commutator [D^n,x^n] = D^n x^n - x^n D^n = Sum_{i=0..n-1} ((n!/i!)^2/(n-i)!) x^i D^i on p. 9 of both papers by Belov-Kanel and Kontsevich. - Tom Copeland, Jan 23 2020

Name changed and table given by Wolfdieter Lang, Nov 28 2011

A049029 Triangle read by rows, the Bell transform of the quartic factorial numbers A007696(n+1) without column 0.

Original entry on oeis.org

1, 5, 1, 45, 15, 1, 585, 255, 30, 1, 9945, 5175, 825, 50, 1, 208845, 123795, 24150, 2025, 75, 1, 5221125, 3427515, 775845, 80850, 4200, 105, 1, 151412625, 108046575, 27478710, 3363045, 219450, 7770, 140, 1, 4996616625, 3824996175, 1069801425

Offset: 1

Wolfdieter Lang

Previous name was: Triangle of numbers related to triangle A048882; generalization of Stirling numbers of second kind A008277, Lah-numbers A008297, ...
a(n,m) enumerates unordered n-vertex m-forests composed of m plane increasing quintic (5-ary) trees. Proof based on the a(n,m) recurrence. See also the F. Bergeron et al. reference, especially Table 1, first row and Example 1 for the e.g.f. for m=1. - Wolfdieter Lang, Sep 14 2007
Also the Bell transform of A007696(n+1). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 28 2016

			Triangle starts:
{1};
{5,1};
{45,15,1};
{585,255,30,1};
{9945,5175,825,50,1};
...

F. Bergeron, Ph. Flajolet and B. Salvy, Varieties of Increasing Trees, in Lecture Notes in Computer Science vol. 581, ed. J.-C. Raoult, Springer 1992, pp. 24-48.
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.
T. Copeland, Mathemagical Forests
T. Copeland, Addendum to Mathemagical Forests
T. Copeland, A Class of Differential Operators and the Stirling Numbers
M. Janjic, Some classes of numbers and derivatives, JIS 12 (2009) #09.8.3.
W. Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
W. Lang, First 10 rows.
Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104v2 [math.CO]. - From N. J. A. Sloane, Aug 21 2012
E. Neuwirth, Recursively defined combinatorial Functions: Extending Galton's board, Discr. Maths. 239 (2001) 33-51.
Mathias Pétréolle, Alan D. Sokal, Lattice paths and branched continued fractions. II. Multivariate Lah polynomials and Lah symmetric functions, arXiv:1907.02645 [math.CO], 2019.

a(n, m) := S2(5, n, m) is the fifth triangle of numbers in the sequence S2(1, n, m) := A008277(n, m) (Stirling 2nd kind), S2(2, n, m) := A008297(n, m) (Lah), S2(3, n, m) := A035342(n, m), S2(4, n, m) := A035469(n, m). a(n, 1)= A007696(n). A007559(n).
Cf. A048882, A007696. Row sums: A049120(n), n >= 1.
Cf. A094638

# The function BellMatrix is defined in A264428.
# Adds (1,0,0,0, ..) as column 0.
BellMatrix(n -> mul(4*k+1, k=0..n), 9); # Peter Luschny, Jan 28 2016

a[n_, m_] /; n >= m >= 1 := a[n, m] = (4(n-1) + m)*a[n-1, m] + a[n-1, m-1]; a[n_, m_] /; n < m = 0; a[, 0] = 0; a[1, 1] = 1; Flatten[Table[a[n, m], {n, 1, 9}, {m, 1, n}]] (* _Jean-François Alcover, Jul 22 2011 *)
rows = 9;
a[n_, m_] := BellY[n, m, Table[Product[4k+1, {k, 0, j}], {j, 0, rows}]];
Table[a[n, m], {n, 1, rows}, {m, 1, n}] // Flatten (* Jean-François Alcover, Jun 22 2018 *)

a(n, m) = n!*A048882(n, m)/(m!*4^(n-m)); a(n+1, m) = (4*n+m)*a(n, m)+ a(n, m-1), n >= m >= 1; a(n, m) := 0, n

a(n, m) = sum(|A051142(n, j)|*S2(j, m), j=m..n) (matrix product), with S2(j, m) := A008277(j, m) (Stirling2 triangle). Priv. comm. to W. Lang by E. Neuwirth, Feb 15 2001; see also the 2001 Neuwirth reference. See the general comment on products of Jabotinsky matrices given under A035342.

From Peter Bala, Nov 25 2011: (Start)

E.g.f.: G(x,t) = exp(t*A(x)) = 1+t*x+(5*t+t^2)*x^2/2!+(45*t+15*t^2+t^3)*x^3/3!+..., where A(x) = -1+(1-4*x)^(-1/4) satisfies the autonomous differential equation A'(x) = (1+A(x))^5.

The generating function G(x,t) satisfies the partial differential equation t*(dG/dt+G) = (1-4*x)*dG/dx, from which follows the recurrence given above.

The row polynomials are given by D^n(exp(x*t)) evaluated at x = 0, where D is the operator (1+x)^5*d/dx. Cf. A008277 (D = (1+x)*d/dx), A105278 (D = (1+x)^2*d/dx), A035342 (D = (1+x)^3*d/dx) and A035469 (D = (1+x)^4*d/dx).

(End)

New name from Peter Luschny, Jan 30 2016

A111596 The matrix inverse of the unsigned Lah numbers A271703.

Original entry on oeis.org

1, 0, 1, 0, -2, 1, 0, 6, -6, 1, 0, -24, 36, -12, 1, 0, 120, -240, 120, -20, 1, 0, -720, 1800, -1200, 300, -30, 1, 0, 5040, -15120, 12600, -4200, 630, -42, 1, 0, -40320, 141120, -141120, 58800, -11760, 1176, -56, 1, 0, 362880, -1451520, 1693440, -846720, 211680, -28224, 2016, -72, 1

Offset: 0

Wolfdieter Lang, Aug 23 2005

Also the associated Sheffer triangle to Sheffer triangle A111595.
Coefficients of Laguerre polynomials (-1)^n * n! * L(n,-1,x), which equals (-1)^n * Lag(n,x,-1) below. Lag(n,Lag(.,x,-1),-1) = x^n evaluated umbrally, i.e., with (Lag(.,x,-1))^k = Lag(k,x,-1). - Tom Copeland, Apr 26 2014
Without row n=0 and column m=0 this is, up to signs, the Lah triangle A008297.
The unsigned column sequences are (with leading zeros): A000142, A001286, A001754, A001755, A001777, A001778, A111597-A111600 for m=1..10.
The row polynomials p(n,x) := Sum_{m=0..n} a(n,m)*x^m, together with the row polynomials s(n,x) of A111595 satisfy the exponential (or binomial) convolution identity s(n,x+y) = Sum_{k=0..n} binomial(n,k)*s(k,x)*p(n-k,y), n>=0.
Exponential Riordan array [1,x/(1+x)]. Inverse of the exponential Riordan array [1,x/(1-x)], which is the unsigned version of A111596. - Paul Barry, Apr 12 2007
For the unsigned subtriangle without column number m=0 and row number n=0, see A105278.
Unsigned triangle also matrix product |S1|*S2 of Stirling number matrices.
The unsigned row polynomials are Lag(n,-x,-1), the associated Laguerre polynomials of order -1 with negated argument. See Gradshteyn and Ryzhik, Abramowitz and Stegun and Rota (Finite Operator Calculus) for extensive formulas. - Tom Copeland, Nov 17 2007, Sep 09 2008
An infinitesimal matrix generator for unsigned A111596 is given by A132792. - Tom Copeland, Nov 22 2007
From the formalism of A132792 and A133314 for n > k, unsigned A111596(n,k) = a(k) * a(k+1)...a(n-1) / (n-k)! = a generalized factorial, where a(n) = A002378(n) = n-th term of first subdiagonal of unsigned A111596. Hence Deutsch's remark in A002378 provides an interpretation of A111596(n,k) in terms of combinations of certain circular binary words. - Tom Copeland, Nov 22 2007
Given T(n,k)= A111596(n,k) and matrices A and B with A(n,k) = T(n,k)*a(n-k) and B(n,k) = T(n,k)*b(n-k), then A*B = C where C(n,k) = T(n,k)*[a(.)+b(.)]^(n-k), umbrally. - Tom Copeland, Aug 27 2008
Operationally, the unsigned row polynomials may be expressed as p_n(:xD:) = x*:Dx:^n*x^{-1}=x*D^nx^n*x^{-1}= n!*binomial(xD+n-1,n) = (-1)^n n! binomial(-xD,n) = n!L(n,-1,-:xD:), where, by definition, :AB:^n = A^nB^n for any two operators A and B, D = d/dx, and L(n,-1,x) is the Laguerre polynomial of order -1. A similarity transformation of the operators :Dx:^n generates the higher order Laguerre polynomials, which can also be expressed in terms of rising or falling factorials or Kummer's confluent hypergeometric functions (cf. the Mathoverflow post). -  Tom Copeland, Sep 21 2019

			Binomial convolution of row polynomials: p(3,x) = 6*x-6*x^2+x^3; p(2,x) = -2*x+x^2, p(1,x) = x, p(0,x) = 1,
together with those from A111595: s(3,x) = 9*x-6*x^2+x^3; s(2,x) = 1-2*x+x^2, s(1,x) = x, s(0,x) = 1; therefore
9*(x+y)-6*(x+y)^2+(x+y)^3 = s(3,x+y) = 1*s(0,x)*p(3,y) + 3*s(1,x)*p(2,y) + 3*s(2,x)*p(1,y) +1*s(3,x)*p(0,y) = (6*y-6*y^2+y^3) + 3*x*(-2*y+y^2) + 3*(1-2*x+x^2)*y + 9*x-6*x^2+x^3.
From _Wolfdieter Lang_, Apr 28 2014: (Start)
The triangle a(n,m) begins:
n\m  0     1       2       3      4     5   6  7
0:   1
1:   0     1
2:   0    -2       1
3:   0     6      -6       1
4:   0   -24      36     -12      1
5:   0   120    -240     120    -20     1
6:   0  -720    1800   -1200    300   -30   1
7:   0  5040  -15120   12600  -4200   630 -42  1
...
For more rows see the link.
(End)

G. C. Greubel, Rows n=0..100 of triangle, flattened
Wolfdieter Lang, The first 11 rows of the triangle.
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Paul Barry, The Restricted Toda Chain, Exponential Riordan Arrays, and Hankel Transforms, J. Int. Seq. 13 (2010) # 10.8.4, example 4.
Paul Barry, Exponential Riordan Arrays and Permutation Enumeration, J. Int. Seq. 13 (2010) # 10.9.1, example 6.
Paul Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, example 20.
Paul Barry, Combinatorial polynomials as moments, Hankel transforms and exponential Riordan arrays, arXiv preprint arXiv:1105.3044 [math.CO], 2011, also J. Int. Seq. 14 (2011) 11.6.7.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers; Generators, Inversion, and Matrix, Binomial, and Integral Transforms; Lagrange a la Lah
A. Hennessy and P. Barry, Generalized Stirling Numbers, Exponential Riordan Arrays, and Orthogonal Polynomials, J. Int. Seq. 14 (2011) # 11.8.2.
M. Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3.
Mathoverflow, Pochhammer symbol of a differential, and hypergeometric polynomials, a question posed by Emilio Pisanty and answered by Tom Copeland, 2012.
J. Taylor, Counting words with Laguerre polynomials, DMTCS Proc., Vol. AS, 2013, p. 1131-1142. [_Tom Copeland_, Jan 08 2016] [Broken link]
J. Taylor, Formal group laws and hypergraph colorings, doctoral thesis, Univ. of Wash., 2016, p. 96. [_Tom Copeland_, Dec 20 2018]
Jian Zhou, On Some Mathematics Related to the Interpolating Statistics, arXiv:2108.10514 [math-ph], 2021.

Row sums: A111884. Unsigned row sums: A000262.
A002868 gives maximal element (in magnitude) in each row.
Cf. A130561 for a natural refinement.
Cf. A264428, A264429, A271703 (unsigned).
Cf. A008297, A089231, A105278 (variants).
Cf. A002378, A008277, A132792, A133314, A132792, A001263.

# The function BellMatrix is defined in A264428.
BellMatrix(n -> `if`(n::odd, -(n+1)!, (n+1)!), 9); # Peter Luschny, Jan 27 2016

a[0, 0] = 1; a[n_, m_] := ((-1)^(n-m))*(n!/m!)*Binomial[n-1, m-1]; Table[a[n, m], {n, 0, 10}, {m, 0, n}] // Flatten (* Jean-François Alcover, Jul 05 2013 *)
T[ n_, k_] := (-1)^n n! Coefficient[ LaguerreL[ n, -1, x], x, k]; (* Michael Somos, Dec 15 2014 *)
rows = 9;
t = Table[(-1)^(n+1) n!, {n, 1, rows}];
T[n_, k_] := BellY[n, k, t];
Table[T[n, k], {n, 0, rows}, {k, 0, n}]  // Flatten (* Jean-François Alcover, Jun 22 2018, after Peter Luschny *)

{T(n, k) = if( n<1 || k<1, n==0 && k==0, (-1)^n * n! * polcoeff( sum(k=1, n, binomial( n-1, k-1) * (-x)^k / k!), k))}; /* Michael Somos, Dec 15 2014 */

lah_number = lambda n, k: factorial(n-k)*binomial(n,n-k)*binomial(n-1,n-k)
A111596_row = lambda n: [(-1)^(n-k)*lah_number(n, k) for k in (0..n)]
for n in range(10): print(A111596_row(n)) # Peter Luschny, Oct 05 2014

# uses[inverse_bell_transform from A264429]
def A111596_matrix(dim):
    fact = [factorial(n) for n in (1..dim)]
    return inverse_bell_transform(dim, fact)
A111596_matrix(10) # Peter Luschny, Dec 20 2015

E.g.f. m-th column: ((x/(1+x))^m)/m!, m>=0.
E.g.f. for row polynomials p(n, x) is exp(x*y/(1+y)).
a(n, m) = ((-1)^(n-m))*|A008297(n, m)| = ((-1)^(n-m))*(n!/m!)*binomial(n-1, m-1), n>=m>=1; a(0, 0)=1; else 0.
a(n, m) = -(n-1+m)*a(n-1, m) + a(n-1, m-1), n>=m>=0, a(n, -1):=0, a(0, 0)=1; a(n, m)=0 if n|a(n,m)| = Sum_{k=m..n} |S1(n,k)|*S2(k,m), n>=0. S2(n,m):=A048993. S1(n,m):=A048994. - Wolfdieter Lang, May 04 2007
From Tom Copeland, Nov 21 2011: (Start)
For this Lah triangle, the n-th row polynomial is given umbrally by
  (-1)^n n! binomial(-Bell.(-x),n), where Bell_n(-x)= exp(x)(xd/dx)^n exp(-x), the n-th Bell / Touchard / exponential polynomial with neg. arg., (cf. A008277). E.g., 2! binomial(-Bell.(-x),2) = -Bell.(-x)*(-Bell.(-x)-1) = Bell_2(-x)+Bell_1(-x) = -2x+x^2.
A Dobinski relation is (-1)^n n! binomial(-Bell.(-x),n)= (-1)^n n! e^x Sum_{j>=0} (-1)^j binomial(-j,n)x^j/j!= n! e^x Sum_{j>=0} (-1)^j binomial(j-1+n,n)x^j/j!. See the Copeland link for the relation to inverse Mellin transform. (End)
The n-th row polynomial is (-1/x)^n e^x (x^2*D_x)^n e^(-x). - Tom Copeland, Oct 29 2012
Let f(.,x)^n = f(n,x) = x!/(x-n)!, the falling factorial,and r(.,x)^n = r(n,x) = (x-1+n)!/(x-1)!, the rising factorial, then the Lah polynomials, Lah(n,t)= n!*Sum{k=1..n} binomial(n-1,k-1)(-t)^k/k! (extra sign factor on odd rows), give the transform Lah(n,-f(.,x))= r(n,x), and Lah(n,r(.,x))= (-1)^n * f(n,x). - Tom Copeland, Oct 04 2014
|T(n,k)| = Sum_{j=0..2*(n-k)} A254881(n-k,j)*k^j/(n-k)!. Note that A254883 is constructed analogously from A254882. - Peter Luschny, Feb 10 2015
The T(n,k) are the inverse Bell transform of [1!,2!,3!,...] and |T(n,k)| are the Bell transform of [1!,2!,3!,...]. See A264428 for the definition of the Bell transform and A264429 for the definition of the inverse Bell transform. - Peter Luschny, Dec 20 2015
Dividing each n-th diagonal by n!, where the main diagonal is n=1, generates a shifted, signed Narayana matrix A001263. - Tom Copeland, Sep 23 2020

New name using a comment from Wolfdieter Lang by Peter Luschny, May 10 2021

A035469 Triangle read by rows, the Bell transform of the triple factorial numbers A007559(n+1) without column 0.

Original entry on oeis.org

1, 4, 1, 28, 12, 1, 280, 160, 24, 1, 3640, 2520, 520, 40, 1, 58240, 46480, 11880, 1280, 60, 1, 1106560, 987840, 295960, 40040, 2660, 84, 1, 24344320, 23826880, 8090880, 1296960, 109200, 4928, 112, 1, 608608000, 643843200

Offset: 1

Wolfdieter Lang

Previous name was: Triangle of numbers related to triangle A035529; generalization of Stirling numbers of second kind A008277, Lah-numbers A008297 and A035342.
a(n,m) enumerates unordered n-vertex m-forests composed of m plane increasing quartic (4-ary) trees. Proof based on the a(n,m) recurrence. See a D. Callan comment on the m=1 case A007559. See also the F. Bergeron et al. reference, especially Table 1, first row and Example 1 for the e.g.f. for m=1. - Wolfdieter Lang, Sep 14 2007
For the definition of the Bell transform see A264428. - Peter Luschny, Jan 19 2016

			Triangle starts:
     {1}
     {4,    1}
    {28,   12,    1}
   {280,  160,   24,    1}
  {3640, 2520,  520,   40,    1}

F. Bergeron, Ph. Flajolet and B. Salvy, Varieties of Increasing Trees, in Lecture Notes in Computer Science vol. 581, ed. J.-C. Raoult, Springer 1922, pp. 24-48.

Peter Bala, Generalized Dobinski formulas
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.
Richell O. Celeste, Roberto B. Corcino and Ken Joffaniel M. Gonzales. Two Approaches to Normal Order Coefficients, Journal of Integer Sequences, Vol. 20 (2017), Article 17.3.5.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers
Tom Copeland, Mathemagical Forests
Tom Copeland, Addendum to Mathemagical Forests
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
Wolfdieter Lang, First 10 rows.
Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104 [math.CO], 2012. - From _N. J. A. Sloane_, Aug 21 2012
E. Neuwirth, Recursively defined combinatorial functions: Extending Galton's board, Discrete Math. 239 (2001) 33-51.
Mathias Pétréolle and Alan D. Sokal, Lattice paths and branched continued fractions. II. Multivariate Lah polynomials and Lah symmetric functions, arXiv:1907.02645 [math.CO], 2019.

a(n, m)=: S2(4, n, m) is the fourth triangle of numbers in the sequence S2(1, n, m) := A008277(n, m) (Stirling 2nd kind), S2(2, n, m) := A008297(n, m) (Lah), S2(3, n, m) := A035342(n, m). a(n, 1)= A007559(n).
Row sums: A049119(n), n >= 1.
Cf. A094638.

a[n_, m_] /; n >= m >= 1 := a[n, m] = (3(n-1) + m)*a[n-1, m] + a[n-1, m-1]; a[n_, m_] /; n < m = 0; a[, 0] = 0; a[1, 1] = 1; Flatten[Table[a[n, m], {n, 1, 9}, {m, 1, n}]] (* _Jean-François Alcover, Jul 22 2011 *)
rows = 9;
a[n_, m_] := BellY[n, m, Table[Product[3k+1, {k, 0, j}], {j, 0, rows}]];
Table[a[n, m], {n, 1, rows}, {m, 1, n}] // Flatten (* Jean-François Alcover, Jun 22 2018 *)

# uses[bell_matrix from A264428]
# Adds a column 1,0,0,0, ... at the left side of the triangle.
bell_matrix(lambda n: A007559(n+1) , 9) # Peter Luschny, Jan 19 2016

a(n, m) = Sum_{j=m..n} |A051141(n, j)|*S2(j, m) (matrix product), with S2(j, m):=A008277(j, m) (Stirling2 triangle). Priv. comm. to Wolfdieter Lang by E. Neuwirth, Feb 15 2001; see also the 2001 Neuwirth reference. See the general comment on products of Jabotinsky matrices given under A035342.
a(n, m) = n!*A035529(n, m)/(m!*3^(n-m)); a(n+1, m) = (3*n+m)*a(n, m) + a(n, m-1), n >= m >= 1; a(n, m) := 0, n < m; a(n, 0) := 0, a(1, 1)=1;
E.g.f. of m-th column: ((-1+(1-3*x)^(-1/3))^m)/m!.
From Peter Bala, Nov 25 2011: (Start)
E.g.f.: G(x,t) = exp(t*A(x)) = 1 + t*x + (4*t+t^2)*x^2/2! + (28*t + 12*t^2 + t^3)*x^3/3! + ..., where A(x) = -1 + (1-3*x)^(-1/3) satisfies the autonomous differential equation A'(x) = (1+A(x))^4.
The generating function G(x,t) satisfies the partial differential equation t*(dG/dt+G) = (1-3*x)*dG/dx, from which follows the recurrence given above.
The row polynomials are given by D^n(exp(x*t)) evaluated at x = 0, where D is the operator (1+x)^4*d/dx. Cf. A008277 (D = (1+x)*d/dx), A105278 (D = (1+x)^2*d/dx), A035342 (D = (1+x)^3*d/dx) and A049029 (D = (1+x)^5*d/dx).
(End)
Dobinski-type formula for the row polynomials: R(n,x) = exp(-x)*Sum_{k>=0} k*(k+3)*(k+6)*...*(k+3*(n-1))*x^k/k!. - Peter Bala, Jun 23 2014

New name from Peter Luschny, Jan 19 2016

A271703 Triangle read by rows: the unsigned Lah numbers T(n, k) = binomial(n-1, k-1)*n!/k! if n > 0 and k > 0, T(n, 0) = 0^n and otherwise 0, for n >= 0 and 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 6, 6, 1, 0, 24, 36, 12, 1, 0, 120, 240, 120, 20, 1, 0, 720, 1800, 1200, 300, 30, 1, 0, 5040, 15120, 12600, 4200, 630, 42, 1, 0, 40320, 141120, 141120, 58800, 11760, 1176, 56, 1, 0, 362880, 1451520, 1693440, 846720, 211680, 28224, 2016, 72, 1

Offset: 0

Peter Luschny, Apr 14 2016

The Lah numbers can be seen as the case m=1 of the family of triangles T_{m}(n,k) = T_{m}(n-1,k-1)+(k^m+(n-1)^m)*T_{m}(n-1,k) (see the link 'Partition transform').

This is the Sheffer triangle (lower triangular infinite matrix) (1, x/(1-x)), an element of the Jabotinsky subgroup of the Sheffer group. - Wolfdieter Lang, Jun 12 2017

			As a rectangular array (diagonals of the triangle):
  1,      1,       1,       1,       1,       1,       ... A000012
  0,      2,       6,       12,      20,      30,      ... A002378
  0,      6,       36,      120,     300,     630,     ... A083374
  0,      24,      240,     1200,    4200,    11760,   ... A253285
  0,      120,     1800,    12600,   58800,   211680,  ...
  0,      720,     15120,   141120,  846720,  3810240, ...
A000007, A000142, A001286, A001754, A001755,  A001777.
The triangle T(n,k) begins:
n\k 0       1        2        3        4       5      6     7    8  9 10 ...
0:  1
1:  0       1
2:  0       2        1
3:  0       6        6        1
4:  0      24       36       12        1
5:  0     120      240      120       20       1
6:  0     720     1800     1200      300      30      1
7:  0    5040    15120    12600     4200     630     42     1
8:  0   40320   141120   141120    58800   11760   1176    56    1
9:  0  362880  1451520  1693440   846720  211680  28224  2016   72  1
10: 0 3628800 16329600 21772800 12700800 3810240 635040 60480 3240 90  1
...  - _Wolfdieter Lang_, Jun 12 2017

R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, 2nd ed., pp. 312, 552.
I. Lah, Eine neue Art von Zahlen, ihre Eigenschaften und Anwendung in der mathematischen Statistik, Mitt.-Bl. Math. Statistik, 7:203-213, 1955.
T. Mansour, M. Schork, Commutation Relations, Normal Ordering, and Stirling Numbers, CRC Press, 2016

Michael De Vlieger, Table of n, a(n) for n = 0..11475 (rows 0 <= n <= 150, flattened)
Richell O. Celeste, Roberto B. Corcino, and Ken Joffaniel M. Gonzales, Two Approaches to Normal Order Coefficients, Journal of Integer Sequences, Vol. 20 (2017), Article 17.3.5.
M. F. Hasler and Peter Luschny, Formulas for A271703, OEIS Wiki, Aug. 2017.
S. A. Joni, G.-C. Rota, and B. Sagan, From sets to functions: Three elementary examples, Discrete Mathematics, Volume 37, Issues 2-3, 1981, 193-202.
Marin Knežević, Vedran Krčadinac, and Lucija Relić, Matrix products of binomial coefficients and unsigned Stirling numbers, arXiv:2012.15307 [math.CO], 2020.
D. E. Knuth, Convolution polynomials, Mathematica J. 2.1 (1992), no. 4, 67-78.
Peter Luschny, Lah numbers
Peter Luschny, Partition transform
Robert S. Maier, Boson Operator Ordering Identities from Generalized Stirling and Eulerian Numbers, arXiv:2308.10332 [math.CO], 2023. See. p. 18.
Piotr Miska and Maciej Ulas, On some properties of the number of permutations being products of pairwise disjoint d-cycles, arXiv:1904.03395 [math.NT], 2019.
Emanuele Munarini, Combinatorial identities involving the central coefficients of a Sheffer matrix, Applicable Analysis and Discrete Mathematics (2019) Vol. 13, 495-517.
Elena L. Wang and Guoce Xin, On Ward Numbers and Increasing Schröder Trees, arXiv:2507.15654 [math.CO], 2025. See p. 12.

Variants: A008297 the main entry for these numbers, A105278, A111596 (signed).
A000262 (row sums). Largest number of the n-th row in A002868.
Cf. A097805, A354795, A360177.

T := (n, k) -> `if`(n=k, 1, binomial(n-1,k-1)*n!/k!):
seq(seq(T(n, k), k=0..n), n=0..9);

T[n_, k_] := Binomial[n, k]*FactorialPower[n-1, n-k];
Table[T[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 20 2017 *)

@cached_function
def T(n,k):
    if k<0 : return 0
    if k==n: return 1
    return T(n-1,k-1) + (k+n-1)*T(n-1,k)
for n in (0..8): print([T(n,k) for k in (0..n)])

For a collection of formulas see the 'Lah numbers' link.
T(n, k) = A097805(n, k)*n!/k! = (-1)^k*P_{n, k}(1,1,1,...) where P_{n, k}(s) is the partition transform of s.
T(n, k) = coeff(n! * P(n), x, k) with P(n) = (1/n)*(Sum_{k=0..n-1}(x(n-k)*P(k))), for n >= 1 and P(n=0) = 1, with x(n) = n*x. See A036039. - Johannes W. Meijer, Jul 08 2016
From Wolfdieter Lang, Jun 12 2017: (Start)
E.g.f. of row polynomials R(n, x) = Sum_{k=0..n} T(n, k)*x^k (that is egf of the triangle) is exp(x*t/(1-t)) (a Sheffer triangle of the Jabotinsky type).
E.g.f. column k: (t/(1-t))^k/k!.
Three term recurrence: T(n, k) = T(n-1, k-1) + (n-1+k)*T(n, k-1), n >= 1, k = 0..n, with T(0, 0) =1, T(n, -1) = 0, T(n, k) = 0 if n < k.
T(n, k) = binomial(n, k)*fallfac(x=n-1, n-k), with fallfac(x, n) = Product_{j=0..(n-1)} (x - j), for n >= 1, and 0 for n = 0.
risefac(x, n) = Sum_{k=0..n} T(n, k)*fallfac(k), with risefac(x, n) = Product_{j=0..(n-1)} (x + j), for n >= 1, and 0 for n = 0.
See Graham et al., exercise 31, p. 312, solution p. 552. (End)
Formally, let f_n(x) = Sum_{k>n} (k-1)!*x^k; then f_n(x) = Sum_{k=0..n} T(n, k)* x^(n+k)*f_0^((k))(x), where ^((k)) stands for the k-th derivative. - Luc Rousseau, Dec 27 2020
T(n, k) = Sum_{j=k..n} A354795(n, j)*A360177(j, k). - Mélika Tebni, Feb 02 2023
T(n, k) = binomial(n, k)*(n-1)!/(k-1)! for n, k > 0. - Chai Wah Wu, Nov 30 2023

A157400 A partition product with biggest-part statistic of Stirling_1 type (with parameter k = -2) as well as of Stirling_2 type (with parameter k = -2), (triangle read by rows).

Original entry on oeis.org

1, 1, 2, 1, 6, 6, 1, 24, 24, 24, 1, 80, 180, 120, 120, 1, 330, 1200, 1080, 720, 720, 1, 1302, 7770, 10920, 7560, 5040, 5040, 1, 5936, 57456, 102480, 87360, 60480, 40320, 40320, 1, 26784, 438984, 970704, 1103760, 786240, 544320, 362880, 362880

Offset: 1

Peter Luschny, Mar 09 2009, Mar 14 2009

Partition product of Product_{j=0..n-1} ((k+1)*j - 1) and n! at k = -2, summed over parts with equal biggest part (Stirling_2 type) as well as partition product of Product_{j=0..n-2} (k-n+j+2) and n! at k = -2 (Stirling_1 type).
It shares this property with the signless Lah numbers.
Underlying partition triangle is A130561.
Same partition product with length statistic is A105278.
Diagonal a(A000217) = A000142.
Row sum is A000262.
T(n,k) is the number of nilpotent elements in the symmetric inverse semigroup (partial bijections) on [n] having index k. Equivalently, T(n,k) is the number of directed acyclic graphs on n labeled nodes with every node having indegree and outdegree at most one and the longest path containing exactly k nodes. - Geoffrey Critzer, Nov 21 2021

			Triangle starts:
  1;
  1,   2;
  1,   6,    6;
  1,  24,   24,   24;
  1,  80,  180,  120, 120;
  1, 330, 1200, 1080, 720, 720;
  ...

Peter Luschny, Counting with Partitions.
Peter Luschny, Generalized Stirling_1 Triangles.
Peter Luschny, Generalized Stirling_2 Triangles.

Cf. A157396, A157397, A157398, A157399, A080510, A157401, A157402, A157403, A157404, A157405, A157386, A157385, A157384, A157383, A126074, A157391, A157392, A157393, A157394, A157395.

egf:= k-> exp((x^(k+1)-x)/(x-1))-exp((x^k-x)/(x-1)):
T:= (n, k)-> n!*coeff(series(egf(k), x, n+1), x, n):
seq(seq(T(n, k), k=1..n), n=1..10);  # Alois P. Heinz, Oct 10 2015

egf[k_] := Exp[(x^(k+1)-x)/(x-1)] - Exp[(x^k-x)/(x-1)]; T[n_, k_] := n! * SeriesCoefficient[egf[k], {x, 0, n}]; Table[Table[T[n, k], {k, 1, n}], {n, 1, 10}] // Flatten (* Jean-François Alcover, Oct 11 2015, after Alois P. Heinz *)

T(n,0) = [n = 0] (Iverson notation) and for n > 0 and 1 <= m <= n.
T(n,m) = Sum_{a} M(a)|f^a| where a = a_1,...,a_n such that
1*a_1 + 2*a_2 + ... + n*a_n = n and max{a_i} = m, M(a) = n!/(a_1!*..*a_n!),
f^a = (f_1/1!)^a_1*..*(f_n/n!)^a_n and f_n = Product_{j=0..n-1} (-j-1)
OR f_n = Product_{j=0..n-2} (j-n) since both have the same absolute value n!.
E.g.f. of column k: exp((x^(k+1)-x)/(x-1))-exp((x^k-x)/(x-1)). - Alois P. Heinz, Oct 10 2015

A002868 Largest number in n-th row of triangle of Lah numbers (A008297 and A271703).

Original entry on oeis.org

1, 1, 2, 6, 36, 240, 1800, 15120, 141120, 1693440, 21772800, 299376000, 4390848000, 68497228800, 1133317785600, 19833061248000, 396661224960000, 8299373322240000, 181400588328960000, 4135933413900288000, 98228418580131840000, 2426819753156198400000

Offset: 0

N. J. A. Sloane

N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 0..100
Victor Meally, Comparison of several sequences given in Motzkin's paper "Sorting numbers for cylinders...", letter to N. J. A. Sloane, N. D.
T. S. Motzkin, Sorting numbers for cylinders and other classification numbers, in Combinatorics, Proc. Symp. Pure Math. 19, AMS, 1971, pp. 167-176. [Annotated, scanned copy]
OEIS Wiki, Sorting numbers

Essentially the same as A001286.

Cf. A000262, A008297, A105278, A271703.

a002868 n = if n == 0 then 1 else maximum $ map abs $ a008297_row n
-- Reinhard Zumkeller, Sep 30 2014

with(combinat): for n from 0 to 35 do big := 1: for m from 1 to n do if big < n!*binomial(n-1,m-1)/m! then big := n!*binomial(n-1,m-1)/m! fi: od: printf(`%d,`,big): od:

a[n_] := ( big = 1; For[ m = 1 , m <= n, m++, b = n!*Binomial[n - 1, m - 1]/m!; If[ big < b , big = b ]]; big); Table[a[n], {n, 0, 19}] (* Jean-François Alcover, Sep 21 2012, after Maple *)

For 2 <= n <= 7, equals (n+1)!*n/2. - Alexander R. Povolotsky, Oct 16 2006

More terms from James Sellers, Jan 03 2001

A062137 Coefficient triangle of generalized Laguerre polynomials n!*L(n,3,x) (rising powers of x).

Original entry on oeis.org

1, 4, -1, 20, -10, 1, 120, -90, 18, -1, 840, -840, 252, -28, 1, 6720, -8400, 3360, -560, 40, -1, 60480, -90720, 45360, -10080, 1080, -54, 1, 604800, -1058400, 635040, -176400, 25200, -1890, 70, -1, 6652800, -13305600, 9313920

Offset: 0

Wolfdieter Lang, Jun 19 2001

The row polynomials s(n,x) := n!*L(n,3,x) = Sum_{m=0..n} a(n,m)*x^m have e.g.f. exp(-z*x/(1-z))/(1-z)^4. They are Sheffer polynomials satisfying the binomial convolution identity s(n,x+y) = Sum_{k=0..n} binomial(n,k)*s(k,x)*p(n-k,y), with polynomials p(n,x) = Sum_{m=1..n} |A008297(n,m)|*(-x)^m, n >= 1 and p(0,x)=1 (for Sheffer polynomials see A048854 for S. Roman reference).
These polynomials appear in the radial part of the l=1 (p-wave) eigen functions for the discrete energy levels of the H-atom. See Messiah reference.
The unsigned version of this triangle is the triangle of unsigned 2-Lah numbers A143497. - Peter Bala, Aug 25 2008
This matrix (unsigned) is embedded in that for n!*L(n,-3,-x). Introduce 0,0,0 to each unsigned row and then add 1,-2,1,4,2,1 to the beginning of the array as the first three rows to generate n!*L(n,-3,-x). - Tom Copeland, Apr 21 2014
From Wolfdieter Lang, Jul 07 2014: (Start)
The standard Rodrigues formula for the monic generalized Laguerre polynomials (also called Laguerre-Sonin polynomials) is Lm(n,alpha,x) := (-1)^n*n!*L(n,3,x) is x^(-alpha)*exp(x)*(d/dx)^n(exp(-x)*x^(n+alpha)).
Another Rodrigues type formula is Lm(n,alpha,x) = exp(x)*x^(-alpha+n+1)*(-x^2*d/dx)^n*(exp(-x)*x^(alpha+1)). This is derivable from the differential - difference relation of Lm(n,alpha,x) which yields first the creation operator formula -(x*d/dx + (-x + alpha + n + 1))*Lm(n,alpha,x) = Lm(n+1,alpha,x) or in the variable q = log(x) the operator -(d/dq + alpha + n + 1 - exp(q)).
The identity (due to Christoph Mayer) (d/dq - (d/dq)W(q))*f(q) = exp(W(q)*d/dq(exp(-W(q)*f(q)) is then iterated with W(q) = W(alpha,n,q) = exp(q) - (alpha + n + 1)*q and a further change of variables leads then to the given result. (End)

			The triangle a(n,m) begins:
n\m       0        1       2     3    4   5 ...
0:        1
1:        4       -1
2:       20      -10      1
3:      120      -90     18     -1
4:      840     -840    252    -28    1
5:     6720    -8400   3360   -560   40  -1
... Formatted by _Wolfdieter Lang_, Jul 07 2014
For more rows see the link.
n = 2: 2!*L(2,3,x) = 20 - 10*x + x^2.

A. Messiah, Quantum mechanics, vol. 1, p. 419, eq.(XI.18a), North Holland, 1969.

Indranil Ghosh, Rows 0..125, flattened
Wolfdieter Lang, First eleven rows of the triangle.
Index entries for sequences related to Laguerre polynomials

For m=0..5 the (unsigned) columns give A001715, A061206, A062141-A062144. The row sums (signed) give A062146, the row sums (unsigned) give A062147.
Cf. A143497. - Peter Bala, Aug 25 2008
Cf. A062139, A105278. - Wolfdieter Lang, Jul 07 2014

Flatten[Table[((-1)^m)*n!*Binomial[n+3,n-m]/m!,{n,0,9},{m,0,n}]] (* Indranil Ghosh, Feb 23 2017 *)

row(n) = Vecrev(n!*pollaguerre(n, 3)); \\ Michel Marcus, Feb 06 2021

a(n, m) = ((-1)^m)*n!*binomial(n+3, n-m)/m!.

E.g.f. for m-th column sequence: ((-x/(1-x))^m)/(m!*(1-x)^4), m >= 0.

Previous Showing 11-20 of 43 results. Next

cp's OEIS Frontend

A133314 Coefficients of list partition transform: reciprocal of an exponential generating function (e.g.f.).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Sage

Formula

Extensions

A035342 The convolution matrix of the double factorial of odd numbers (A001147).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

Sage

Formula

Extensions

A021009 Triangle of coefficients of Laguerre polynomials n!*L_n(x) (rising powers of x).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

PARI

Sage

Formula

Extensions

A049029 Triangle read by rows, the Bell transform of the quartic factorial numbers A007696(n+1) without column 0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A111596 The matrix inverse of the unsigned Lah numbers A271703.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Sage

Formula

Extensions

A035469 Triangle read by rows, the Bell transform of the triple factorial numbers A007559(n+1) without column 0.

Views