A001497 -id:A001497 - cp's OEIS frontend

A100861 Triangle of Bessel numbers read by rows: T(n,k) is the number of k-matchings of the complete graph K(n).

1, 1, 1, 1, 1, 3, 1, 6, 3, 1, 10, 15, 1, 15, 45, 15, 1, 21, 105, 105, 1, 28, 210, 420, 105, 1, 36, 378, 1260, 945, 1, 45, 630, 3150, 4725, 945, 1, 55, 990, 6930, 17325, 10395, 1, 66, 1485, 13860, 51975, 62370, 10395, 1, 78, 2145, 25740, 135135, 270270, 135135, 1, 91, 3003, 45045, 315315, 945945, 945945, 135135

Offset: 0

Emeric Deutsch, Jan 08 2005

Row n contains 1 + floor(n/2) terms. Row sums yield A000085. T(2n,n) = T(2n-1,n-1) = (2n-1)!! (A001147).
Inverse binomial transform is triangle with T(2n,n) = (2n-1)!!, 0 otherwise. - Paul Barry, May 21 2005
Equivalently, number of involutions of n with k pairs. - Franklin T. Adams-Watters, Jun 09 2006
From Gary W. Adamson, Dec 09 2009: (Start)
If considered as an infinite lower triangular matrix (cf. A144299),
lim_{n->} A100861^n = A118930: (1, 1, 2, 4, 13, 41, ...).
(End)
Sum_{k=0..floor(n/2)} T(n,k)m^(n-2k)s^(2k) is the n-th non-central moment of the normal probability distribution with mean m and standard deviation s. - Stanislav Sykora, Jun 19 2014
Row n is the list of coefficients of the independence polynomial of the n-triangular graph. - Eric W. Weisstein, Nov 11 2016
Restating the 2nd part of the Name, row n is the list of coefficients of the matching-generating polynomial of the complete graph K_n. - Eric W. Weisstein, Apr 03 2018

			T(4, 2) = 3 because in the graph with vertex set {A, B, C, D} and edge set {AB, BC, CD, AD, AC, BD} we have the following three 2-matchings: {AB, CD},{AC, BD} and {AD, BC}.
Triangle starts:
[0] 1;
[1] 1;
[2] 1,  1;
[3] 1,  3;
[4] 1,  6,   3;
[5] 1, 10,  15;
[6] 1, 15,  45,   15;
[7] 1, 21, 105,  105;
[8] 1, 28, 210,  420, 105;
[9] 1, 36, 378, 1260, 945.
.
From _Eric W. Weisstein_, Nov 11 2016: (Start)
As polynomials:
1,
1,
1 + x,
1 + 3*x,
1 + 6*x + 3*x^2,
1 + 10*x + 15*x^2,
1 + 15*x + 45*x^2 + 15*x^3. (End)

M. Abramowitz and I. Stegun, Handbook of Mathematical Functions (1983 reprint), 10th edition, 1964, expression 22.3.11 in page 775.
C. D. Godsil, Algebraic Combinatorics, Chapman & Hall, New York, 1993.

T. D. Noe, Rows n = 0..100, flattened
Paul Barry, The Gamma-Vectors of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1804.05027 [math.CO], 2018.
Gérard Le Caër, A new family of solvable Pearson-Dirichlet random walks, Journal of Statistical Physics 144:1 (2011), pp. 23-45.
Ji Young Choi and Jonathan D. H. Smith, On the unimodality and combinatorics of Bessel numbers, Discrete Math., 264 (2003), 45-53.
Tom Copeland, Infinitesimal Generators, the Pascal Pyramid, and the Witt and Virasoro Algebras, 2012.
John Engbers, David Galvin, and Clifford Smyth, Restricted Stirling and Lah numbers and their inverses, arXiv:1610.05803 [math.CO], 2016.
Mikael Fremling, On the modular covariance properties of composite fermions on the torus, arXiv:1810.10391 [cond-mat.str-el], 2018.
Gary R. W. Greaves, Jeven Syatriadi, and Charissa I. Utomo, Chromatic polynomials of signed graphs and dominating-vertex deletion formulae, arXiv:2407.00883 [math.CO], 2024. See p. 11.
A. Hernando, R. Hernando, A. Plastino and A. R. Plastino, The workings of the Maximum Entropy Principle in collective human behavior, arXiv preprint arXiv:1201.0905 [stat.AP], 2012.
Robert S. Maier, Boson Operator Ordering Identities from Generalized Stirling and Eulerian Numbers, arXiv:2308.10332 [math.CO], 2023. See p. 18.
Aaron Pollack, Exceptional theta functions, arXiv:2211.05280 [math.NT], Nov 2022. See Lemma 7.5.1.
Eric Weisstein's World of Mathematics, Complete Graph
Eric Weisstein's World of Mathematics, Independence Polynomial
Eric Weisstein's World of Mathematics, Matching-Generating Polynomial
Eric Weisstein's World of Mathematics, Triangular Graph
Wikipedia, Normal distribution, section 'Moments'
Jian Zhou, Quantum deformation theory of the Airy curve and the mirror symmetry of a point, arXiv preprint arXiv:1405.5296 [math.AG], 2014.

Other versions of this same triangle are given in A144299, A001497, A001498, A111924.
Cf. A000085 (row sums).
Cf. A118930, A008619.

a100861 n k = a100861_tabf !! n !! k
a100861_row n = a100861_tabf !! n
a100861_tabf = zipWith take a008619_list a144299_tabl
-- Reinhard Zumkeller, Jan 02 2014

P[0]:=1: for n from 1 to 14 do P[n]:=sort(expand(P[n-1]+(n-1)*t*P[n-2])) od: for n from 0 to 14 do seq(coeff(t*P[n],t^k),k=1..1+floor(n/2)) od; # yields the sequence in triangular form
# Alternative:
A100861 := proc(n,k)
    n!/k!/(n-2*k)!/2^k ;
end proc:
seq(seq(A100861(n,k),k=0..n/2),n=0..10) ; # R. J. Mathar, Aug 19 2014

Table[Table[n!/(i! 2^i (n - 2 i)!), {i, 0, Floor[n/2]}], {n, 0, 10}] // Flatten  (* Geoffrey Critzer, Mar 27 2011 *)
CoefficientList[Table[2^(n/2) (-(1/x))^(-n/2) HypergeometricU[-n/2, 1/2, -1/(2 x)], {n, 0, 10}], x] // Flatten (* Eric W. Weisstein, Apr 03 2018 *)
CoefficientList[Table[(-I)^n Sqrt[x/2]^n HermiteH[n, I/Sqrt[2 x]], {n, 0, 10}], x] // Flatten (* Eric W. Weisstein, Apr 03 2018 *)

T(n,k)=if(k<0 || 2*k>n, 0, n!/k!/(n-2*k)!/2^k) /* Michael Somos, Jun 04 2005 */

T(n, k) = n!/(k!(n-2k)!*2^k).
E.g.f.: exp(z+tz^2/2).
G.f.: g(t, z) satisfies the differential equation g = 1 + zg + tz^2*(d/dz)(zg).
Row generating polynomial = P[n] = [-i*sqrt(t/2)]^n*H(n, i/sqrt(2t)), where H(n, x) is a Hermite polynomial and i=sqrt(-1). Row generating polynomials P[n] satisfy P[0]=1, P[n] = P[n-1] + (n-1)tP[n-2].
T(n, k) = binomial(n, 2k)(2k-1)!!. - Paul Barry, May 21 2002 [Corrected by Roland Hildebrand, Mar 06 2009]
T(n,k) = (n-2k+1)*T(n-1,k-1) + T(n-1,k). - Franklin T. Adams-Watters, Jun 09 2006
E.g.f.: 1 + (x+y*x^2/2)/(E(0)-(x+y*x^2/2)), where E(k) = 1 + (x+y*x^2/2)/(1 + (k+1)/E(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Nov 08 2013
T(n,k) = A144299(n,k), k=0..n/2. - Reinhard Zumkeller, Jan 02 2014

A001498 Triangle a(n,k) (n >= 0, 0 <= k <= n) of coefficients of Bessel polynomials y_n(x) (exponents in increasing order).

Original entry on oeis.org

1, 1, 1, 1, 3, 3, 1, 6, 15, 15, 1, 10, 45, 105, 105, 1, 15, 105, 420, 945, 945, 1, 21, 210, 1260, 4725, 10395, 10395, 1, 28, 378, 3150, 17325, 62370, 135135, 135135, 1, 36, 630, 6930, 51975, 270270, 945945, 2027025, 2027025, 1, 45, 990, 13860, 135135, 945945, 4729725, 16216200, 34459425, 34459425

Offset: 0

N. J. A. Sloane

The row polynomials with exponents in increasing order (e.g., third row: 1+3x+3x^2) are Grosswald's y_{n}(x) polynomials, p. 18, Eq. (7).
Also called Bessel numbers of first kind.
The triangle a(n,k) has factorization [C(n,k)][C(k,n-k)]Diag((2n-1)!!) The triangle a(n-k,k) is A100861, which gives coefficients of scaled Hermite polynomials. - Paul Barry, May 21 2005
Related to k-matchings of the complete graph K_n by a(n,k)=A100861(n+k,k). Related to the Morgan-Voyce polynomials by a(n,k)=(2k-1)!!*A085478(n,k). - Paul Barry, Aug 17 2005
Related to Hermite polynomials by a(n,k)=(-1)^k*A060821(n+k, n-k)/2^n. - Paul Barry, Aug 28 2005
The row polynomials, the Bessel polynomials y(n,x):=Sum_{m=0..n} (a(n,m)*x^m) (called y_{n}(x) in the Grosswald reference) satisfy (x^2)*(d^2/dx^2)y(n,x) + 2*(x+1)*(d/dx)y(n,x) - n*(n+1)*y(n,x) = 0.
a(n-1, m-1), n >= m >= 1, enumerates unordered n-vertex forests composed of m plane (aka ordered) increasing (rooted) trees. Proof from the e.g.f. of the first column Y(z):=1-sqrt(1-2*z) (offset 1) and the Bergeron et al. eq. (8) Y'(z)= phi(Y(z)), Y(0)=0, with out-degree o.g.f. phi(w)=1/(1-w). See their remark on p. 28 on plane recursive trees. For m=1 see the D. Callan comment on A001147 from Oct 26 2006. - Wolfdieter Lang, Sep 14 2007
The asymptotic expansions of the higher order exponential integrals E(x,m,n), see A163931 for information, lead to the Bessel numbers of the first kind in an intriguing way. For the first four values of m these asymptotic expansions lead to the triangles A130534 (m=1), A028421 (m=2), A163932 (m=3) and A163934 (m=4). The o.g.f.s. of the right hand columns of these triangles in their turn lead to the triangles A163936 (m=1), A163937 (m=2), A163938 (m=3) and A163939 (m=4). The row sums of these four triangles lead to A001147, A001147 (minus a(0)), A001879 and A000457 which are the first four right hand columns of A001498. We checked this phenomenon for a few more values of m and found that this pattern persists: m = 5 leads to A001880, m=6 to A001881, m=7 to A038121 and m=8 to A130563 which are the next four right hand columns of A001498. So one by one all columns of the triangle of coefficients of Bessel polynomials appear. - Johannes W. Meijer, Oct 07 2009
a(n,k) also appear as coefficients of (n+1)st degree of the differential operator D:=1/t d/dt, namely D^{n+1}= Sum_{k=0..n} a(n,k) (-1)^{n-k} t^{1-(n+k)} (d^{n+1-k}/dt^{n+1-k}. - Leonid Bedratyuk, Aug 06 2010
a(n-1,k) are the coefficients when expanding (xI)^n in terms of powers of I. Let I(f)(x) := Integral_{a..x} f(t) dt, and (xI)^n := x Integral_{a..x} [ x_{n-1} Integral_{a..x_{n-1}} [ x_{n-2} Integral_{a..x_{n-2}} ... [ x_1 Integral_{a..x_1} f(t) dt ] dx_1 ] .. dx_{n-2} ] dx_{n-1}. Then: (xI)^n = Sum_{k=0..n-1} (-1)^k * a(n-1,k) * x^(n-k) * I^(n+k)(f)(x) where I^(n) denotes iterated integration. - Abdelhay Benmoussa, Apr 11 2025

			The triangle a(n, k), n >= 0, k = 0..n, begins:
  1
  1  1
  1  3   3
  1  6  15    15
  1 10  45   105    105
  1 15 105   420    945    945
  1 21 210  1260   4725  10395   10395
  1 28 378  3150  17325  62370  135135   135135
  1 36 630  6930  51975 270270  945945  2027025  2027025
  1 45 990 13860 135135 945945 4729725 16216200 34459425 34459425
  ...
And the first few Bessel polynomials are:
  y_0(x) = 1,
  y_1(x) = x + 1,
  y_2(x) = 3*x^2 + 3*x + 1,
  y_3(x) = 15*x^3 + 15*x^2 + 6*x + 1,
  y_4(x) = 105*x^4 + 105*x^3 + 45*x^2 + 10*x + 1,
  y_5(x) = 945*x^5 + 945*x^4 + 420*x^3 + 105*x^2 + 15*x + 1,
  ...
Tree counting: a(2,1)=3 for the unordered forest of m=2 plane increasing trees with n=3 vertices, namely one tree with one vertex (root) and another tree with two vertices (a root and a leaf), labeled increasingly as (1, 23), (2,13) and (3,12). - _Wolfdieter Lang_, Sep 14 2007

J. Riordan, Combinatorial Identities, Wiley, 1968, p. 77.

T. D. Noe, Rows n=0..50 of triangle, flattened
Alexander Alldridge, Joachim Hilgert, and Martin R. Zirnbauer, Chevalley's restriction theorem for reductive symmetric superpairs, arXiv:0812.3530 [math.RT], 2008-2009; J. Alg. 323 (4) (2010) 1159-1185 doi:10.1016/j.jalgebra.2009.11.014, Remark 3.17.
David Applegate and N. J. A. Sloane, The Gift Exchange Problem, arXiv:0907.0513 [math.CO], 2009.
Juan Antonio Barcelo and Anthony Carbery, On the magnitudes of compact sets in Euclidean spaces, arXiv preprint arXiv:1507.02502 [math.MG], 2015.
François Bergeron, Philippe Flajolet, and Bruno Salvy, Varieties of Increasing Trees, in Lecture Notes in Computer Science vol. 581, ed. J.-C. Raoult, Springer 1922, pp. 24-48.
Alexander W. Boldyreff, Decomposition of Rational Fractions into Partial Fractions, Nat. Math. Mag. 17 (6) (1943), 261-267; coefficients (m)N(r).
Alexander Burstein and Toufik Mansour, Words restricted by patterns with at most 2 distinct letters, arXiv:math/0110056 [math.CO], 2001.
Roudy El Haddad, Repeated Integration and Explicit Formula for the n-th Integral of x^m*(ln x)^m', arXiv:2102.11723 [math.GM], 2021.
Andrew Francis and Michael Hendriksen, Counting spinal phylogenetic networks, arXiv:2502.14223 [q-bio.PE], 2025. See p. 9.
Emil Grosswald, Bessel Polynomials: Recurrence Relations, Lecture Notes Math. vol. 698, 1978, p. 18.
Cameron Jakub and Mihai Nica, Depth Degeneracy in Neural Networks: Vanishing Angles in Fully Connected ReLU Networks on Initialization, arXiv:2302.09712 [stat.ML], 2023.
Taekyun Kim, and Dae San Kim, Identities involving Bessel polynomials arising from linear differential equations, arXiv:1602.04106 [math.NT], 2016.
H. L. Krall and Orrin Frink, A new class of orthogonal polynomials, Trans. Amer. Math. Soc. 65, 100-115, 1949.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
Wolfdieter Lang, First ten rows.
B. Leclerc, Powers of staircase Schur functions and symmetric analogues of Bessel polynomials, Discrete Math., 153 (1996), 213-227.
Shi-Mei Ma, Toufik Mansour, and Matthias Schork. Normal ordering problem and the extensions of the Stirling grammar, arXiv preprint arXiv:1308.0169 [math.CO], 2013.
Shi-Mei Ma, Toufik Mansour, Jean Yeh, and Yeong-Nan Yeh, Normal ordered grammars, arXiv:2404.15119 [math.CO], 2024. See p. 11.
Guillermo Navas-Palencia, On the computation of the cumulative distribution function of the Normal Inverse Gaussian distribution, arXiv:2502.16015 [math.NA], 2025. See p. 25.
Andrew Elvey Price and Alan D. Sokal, Phylogenetic trees, augmented perfect matchings, and a Thron-type continued fraction (T-fraction) for the Ward polynomials, arXiv:2001.01468 [math.CO], 2020.
John Riordan, Notes to N. J. A. Sloane, Jul. 1968
Florian Stober, Average case considerations for MergeInsertion, Master's Thesis, University of Stuttgart, Institute of Formal Methods in Computer Science, 2018.
Florian Stober and Armin Weiß, On the Average Case of MergeInsertion, arXiv:1905.09656 [cs.DS], 2019.
Laszlo A. Székely, Pál L. Erdős, and M. A. Steel, The combinatorics of evolutionary trees, Séminaire Lotharingien de Combinatoire, B28e (1992), 15 pp.
Juan G. Triana, Bessel polynomials by context-free grammars (Polinomios de Bessel mediante gramáticas independientes del contexto), Bistua, Univ. de Pamplona (Colombia, 2024) Vol 22, No. 2. See p. 3.
Jonas Wahl, Traces on diagram algebras II: Centralizer algebras of easy groups and new variations of the Young graph, arXiv:2009.08181 [math.RT], 2020.
Eric Weisstein's World of Mathematics, Modified Spherical Bessel Function of the Second Kind
Index entries for sequences related to Bessel functions or polynomials

Cf. A001497 (same triangle but rows read in reverse order). Other versions of this same triangle are given in A144331, A144299, A111924 and A100861.
Columns from left edge include A000217, A050534.
Columns 1-6 from right edge are A001147, A001879, A000457, A001880, A001881, A038121.
Bessel polynomials evaluated at certain x are A001515 (x=1, row sums), A000806 (x=-1), A001517 (x=2), A002119 (x=-2), A001518 (x=3), A065923 (x=-3), A065919 (x=4). Cf. A043301, A003215.
Cf. A245066 (central terms). A113025 (y_n(2*x)).

a001498 n k = a001498_tabl !! n !! k
a001498_row n = a001498_tabl !! n
a001498_tabl = map reverse a001497_tabl
-- Reinhard Zumkeller, Jul 11 2014

/* As triangle: */ [[Factorial(n+k)/(2^k*Factorial(n-k)*Factorial(k)): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Feb 15 2016

Bessel := proc(n,x) add(binomial(n+k,2*k)*(2*k)!*x^k/(k!*2^k),k=0..n); end; # explicit Bessel polynomials
Bessel := proc(n) option remember; if n <=1 then (1+x)^n else (2*n-1)*x*Bessel(n-1)+Bessel(n-2); fi; end; # recurrence for Bessel polynomials
bessel := proc(n,x) add(binomial(n+k,2*k)*(2*k)!*x^k/(k!*2^k),k=0..n); end;
f := proc(n) option remember; if n <=1 then (1+x)^n else (2*n-1)*x*f(n-1)+f(n-2); fi; end;
# Alternative:
T := (n,k) -> pochhammer(n+1,k)*binomial(n,k)/2^k:
for n from 0 to 9 do seq(T(n,k), k=0..n) od; # Peter Luschny, May 11 2018
T := proc(n, k) option remember; if k = 0 then 1 else if k = n then T(n, k-1)
else (n - k + 1)* T(n, k - 1) + T(n - 1, k) fi fi end:
for n from 0 to 9 do seq(T(n, k), k = 0..n) od;  # Peter Luschny, Oct 02 2023

max=50; Flatten[Table[(n+k)!/(2^k*(n-k)!*k!), {n, 0, Sqrt[2 max]//Ceiling}, {k, 0, n}]][[1 ;; max]] (* Jean-François Alcover, Mar 20 2011 *)

{T(n,k)=if(k<0||k>n, 0, binomial(n, k)*(n+k)!/2^k/n!)} /* Michael Somos, Oct 03 2006 */

A001497_ser(N,t='t) = {
  my(x='x+O('x^(N+2)));
  serlaplace(deriv(exp((1-sqrt(1-2*t*x))/t),'x));
};
concat(apply(Vecrev, Vec(A001497_ser(9)))) \\ Gheorghe Coserea, Dec 27 2017

a(n, k) = (n+k)!/(2^k*(n-k)!*k!) (see Grosswald and Riordan). - Ralf Stephan, Apr 20 2004
a(n, 0)=1; a(0, k)=0, k > 0; a(n, k) = a(n-1, k) + (n-k+1) * a(n, k-1) = a(n-1, k) + (n+k-1) * a(n-1, k-1). - Len Smiley
a(n, m) = A001497(n, n-m) = A001147(m)*binomial(n+m, 2*m) for n >= m >= 0, otherwise 0.
G.f. for m-th column: (A001147(m)*x^m)/(1-x)^(2*m+1), m >= 0, where A001147(m) = double factorials (from explicit a(n, m) form).
Row polynomials y_n(x) are given by D^(n+1)(exp(t)) evaluated at t = 0, where D is the operator 1/(1-t*x)*d/dt. - Peter Bala, Nov 25 2011
G.f.: conjecture: T(0)/(1-x), where T(k) = 1 - x*y*(k+1)/(x*y*(k+1) - (1-x)^2/T(k+1)); (continued fraction). - Sergei N. Gladkovskii, Nov 13 2013
Recurrence from Grosswald, p. 18, eq. (5), for the row polynomials: y_n(x) = (2*n-1)*x*y_{n-1} + y_{n-2}(x), y_{-1}(x) = 1 = y_{0} = 1, n >= 1. This becomes, for n >= 0, k = 0..n: a(n, k) = 0 for n < k (zeros not shown in the triangle), a(n, -1) = 0, a(0, 0) = 1 = a(1, 0) and otherwise a(n, k) = (2*n-1)*a(n-1, k-1) + a(n-2, k). Compare with the above given recurrences. - Wolfdieter Lang, May 11 2018
T(n, k) = Pochhammer(n+1,k)*binomial(n,k)/2^k = A113025(n,k)/2^k. - Peter Luschny, May 11 2018
a(n, k) = Sum_{i=0..min(n-1, k)} (n-i)(k-i) * a(n-1, i) where x(n) = x*(x-1)*...*(x-n+1) is the falling factorial, this equality follows directly from the operational formula we wrote in Apr 11 2025.- Abdelhay Benmoussa, May 18 2025

A105278 Triangle read by rows: T(n,k) = binomial(n,k)*(n-1)!/(k-1)!.

Original entry on oeis.org

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

Offset: 1

Miklos Kristof, Apr 25 2005

T(n,k) is the number of partially ordered sets (posets) on n elements that consist entirely of k chains. For example, T(4, 3)=12 since there are exactly 12 posets on {a,b,c,d} that consist entirely of 3 chains. Letting ab denote a<=b and using a slash "/" to separate chains, the 12 posets can be given by a/b/cd, a/b/dc, a/c/bd, a/c/db, a/d/bc, a/d/cb, b/c/ad, b/c/da, b/d/ac, b/d/ca, c/d/ab and c/d/ba, where the listing of the chains is arbitrary (e.g., a/b/cd = a/cd/b =...cd/b/a). - Dennis P. Walsh, Feb 22 2007
Also the matrix product |S1|.S2 of Stirling numbers of both kinds.
This Lah triangle is a lower triangular matrix of the Jabotinsky type. See the column e.g.f. and the D. E. Knuth reference given in A008297. - Wolfdieter Lang, Jun 29 2007
The infinitesimal matrix generator of this matrix is given in A132710. See A111596 for an interpretation in terms of circular binary words and generalized factorials. - Tom Copeland, Nov 22 2007
Three combinatorial interpretations: T(n,k) is (1) the number of ways to split [n] = {1,...,n} into a collection of k nonempty lists ("partitions into sets of lists"), (2) the number of ways to split [n] into an ordered collection of n+1-k nonempty sets that are noncrossing ("partitions into lists of noncrossing sets"), (3) the number of Dyck n-paths with n+1-k peaks labeled 1,2,...,n+1-k in some order. - David Callan, Jul 25 2008
Given 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 = D where D(n,k) = T(n,k)*[a(.)+b(.)]^(n-k), umbrally. - Tom Copeland, Aug 21 2008
An e.g.f. for the row polynomials of A(n,k) = T(n,k)*a(n-k) is exp[a(.)* D_x * x^2] exp(x*t) = exp(x*t) exp[(.)!*Lag(.,-x*t,1)*a(.)*x], umbrally, where [(.)! Lag(.,x,1)]^n = n! Lag(n,x,1) is a normalized Laguerre polynomial of order 1. - Tom Copeland, Aug 29 2008
Triangle of coefficients from the Bell polynomial of the second kind for f = 1/(1-x). B(n,k){x1,x2,x3,...} = B(n,k){1/(1-x)^2,...,(j-1)!/(1-x)^j,...} = T(n,k)/(1-x)^(n+k). - Vladimir Kruchinin, Mar 04 2011
The triangle, with the row and column offset taken as 0, is the generalized Riordan array (exp(x), x) with respect to the sequence n!*(n+1)! 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!^2 is A021009 unsigned). - Peter Bala, Aug 15 2013
For a relation to loop integrals in QCD, see p. 33 of Gopakumar and Gross and Blaizot and Nowak. - Tom Copeland, Jan 18 2016
Also the Bell transform of (n+1)!. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 27 2016
Also the number of k-dimensional flats of the n-dimensional Shi arrangement. - Shuhei Tsujie, Apr 26 2019
The numbers T(n,k) appear as coefficients when expanding the rising factorials (x)^k = x(x+1)...(x+k-1) in the basis of falling factorials (x)k = x(x-1)...(x-k+1). Specifically, (x)^n = Sum{k=1..n} T(n,k) (x)k. - _Jeremy L. Martin, Apr 21 2021

			T(1,1) = C(1,1)*0!/0! = 1,
T(2,1) = C(2,1)*1!/0! = 2,
T(2,2) = C(2,2)*1!/1! = 1,
T(3,1) = C(3,1)*2!/0! = 6,
T(3,2) = C(3,2)*2!/1! = 6,
T(3,3) = C(3,3)*2!/2! = 1,
Sheffer a-sequence recurrence: T(6,2)= 1800 = (6/3)*120 + 6*240.
B(n,k) =
   1/(1-x)^2;
   2/(1-x)^3,  1/(1-x)^4;
   6/(1-x)^4,  6/(1-x)^5,  1/(1-x)^6;
  24/(1-x)^5, 36/(1-x)^6, 12/(1-x)^7, 1/(1-x)^8;
The triangle T(n,k) begins:
  n\k      1       2       3      4      5     6    7  8  9 ...
  1:       1
  2:       2       1
  3:       6       6       1
  4:      24      36      12      1
  5:     120     240     120     20      1
  6:     720    1800    1200    300     30     1
  7:    5040   15120   12600   4200    630    42    1
  8:   40320  141120  141120  58800  11760  1176   56  1
  9:  362880 1451520 1693440 846720 211680 28224 2016 72  1
  ...
Row n=10: [3628800, 16329600, 21772800, 12700800, 3810240, 635040, 60480, 3240, 90, 1]. - _Wolfdieter Lang_, Feb 01 2013
From _Peter Bala_, Feb 24 2025: (Start)
The array factorizes as an infinite product (read from right to left):
  /  1                \        /1             \^m /1           \^m /1           \^m
  |  2    1            |      | 0   1          |  |0  1         |  |1  1         |
  |  6    6   1        | = ...| 0   0   1      |  |0  1  1      |  |0  2  1      |
  | 24   36  12   1    |      | 0   0   1  1   |  |0  0  2  1   |  |0  0  3  1   |
  |120  240 120  20   1|      | 0   0   0  2  1|  |0  0  0  3  1|  |0  0  0  4  1|
  |...                 |      |...             |  |...          |  |...          |
where m = 2. Cf. A008277 (m = 1), A035342 (m = 3), A035469 (m = 4), A049029 (m = 5) A049385 (m = 6), A092082 (m = 7), A132056 (m = 8), A223511 - A223522 (m = 9 through 20), A001497 (m = -1), A004747 (m = -2), A000369 (m = -3), A011801 (m = -4), A013988 (m = -5). (End)

Reinhard Zumkeller, Rows n = 1..100 of triangle, flattened
Peter Bala, Factorising (r,b)-Stirling arrays
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.
Paul Barry, Eulerian polynomials as moments, via exponential Riordan arrays, arXiv preprint arXiv:1105.3043 [math.CO], 2011, J. Int. Seq. 14 (2011) # 11.9.5
Jean-Paul Blaizot and Maciej A. Nowak, Large N_c confinement and turbulence, arXiv:0801.1859 [hep-th], 2008.
David Callan, Sets, Lists and Noncrossing Partitions, arXiv:0711.4841 [math.CO], 2007-2008.
Pietro Codara, Ottavio M. D'Antona, and Pavol Hell, A simple combinatorial interpretation of certain generalized Bell and Stirling numbers, arXiv preprint arXiv:1308.1700 [cs.DM], 2013.
Tom Copeland, Mathemagical Forests, Addendum to Mathemagical Forests, The Inverse Mellin Transform, Bell Polynomials, a Generalized Dobinski Relation, and the Confluent Hypergeometric Functions, A Class of Differential Operators and the Stirling Numbers
Siad Daboul, Jan Mangaldan, Michael Z. Spivey and Peter Taylor, The Lah Numbers and the n-th Derivative of exp(1/x), Math. Mag., 86 (2013), 39-47.
Bérénice Delcroix-Oger and Clément Dupont, Lie-operads and operadic modules from poset cohomology, arXiv:2505.06094 [math.CO], 2025. See p. 33.
G. H. E. Duchamp et al., Feynman graphs and related Hopf algebras, J. Phys. (Conf Ser) 30 (2006) 107-118.
Rajesh Gopakumar and David J. Gross, Mastering the master field, arXiv:hep-th/9411021, 1994.
Gábor Hetyei, Meixner polynomials of the second kind and quantum algebras representing su(1,1), arXiv preprint arXiv:0909.4352 [math.QA], 2009, p. 4. - From _Tom Copeland_, Oct 01 2015
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3.
Donald E. Knuth, Convolution polynomials, The Mathematica J., 2 (1992), 67-78.
Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104v2 [math.CO], 2012. - From _N. J. A. Sloane_, Aug 21 2012
MacTutor History of Mathematics archive: Biography of Ivo Lah.
Robert S. Maier, Boson Operator Ordering Identities from Generalized Stirling and Eulerian Numbers, arXiv:2308.10332 [math.CO], 2023. See p. 19.
Norihiro Nakashima and Shuhei Tsujie, Enumeration of Flats of the Extended Catalan and Shi Arrangements with Species, arXiv:1904.09748 [math.CO], 2019. See p. 18.
Michael Penn, Lah Numbers and an appearance of exponential generating functions, YouTube video, 2025.
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. See p. 4.
Tilman Piesk, Illustration of the first four rows
Kornelia Ufniarz and Grzegorz Siudem, Combinatorial origins of the canonical ensemble, arXiv:2008.00244 [math-ph], 2020. See p. 5.
Weiping Wang and Tianming Wang, Generalized Riordan arrays, Discrete Mathematics, Vol. 308, No. 24, 6466-6500.
Wikipedia, Lah number

Cf. A000262, A103194, A105220.
Triangle of Lah numbers (A008297) unsigned.
Cf. A111596 (signed triangle with extra n=0 row and m=0 column).
Cf. A130561 (for a natural refinement).
Cf. A094638 (for differential operator representation).
Cf. A248045 (central terms), A002868 (row maxima).
Cf, A059110.
Cf. A001263, A008277.
Cf. A089231 (triangle with mirrored rows).
Cf. A271703 (triangle with extra n=0 row and m=0 column).

Flat(List([1..10],n->List([1..n],k->Binomial(n,k)*Factorial(n-1)/Factorial(k-1)))); # Muniru A Asiru, Jul 25 2018

a105278 n k = a105278_tabl !! (n-1) !! (k-1)
a105278_row n = a105278_tabl !! (n-1)
a105278_tabl = [1] : f [1] 2 where
   f xs i = ys : f ys (i + 1) where
     ys = zipWith (+) ([0] ++ xs) (zipWith (*) [i, i + 1 ..] (xs ++ [0]))
-- Reinhard Zumkeller, Sep 30 2014, Mar 18 2013

/* As triangle */ [[Binomial(n,k)*Factorial(n-1)/Factorial(k-1): k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Oct 31 2014

The triangle: for n from 1 to 13 do seq(binomial(n,k)*(n-1)!/(k-1)!,k=1..n) od;
the sequence: seq(seq(binomial(n,k)*(n-1)!/(k-1)!,k=1..n),n=1..13);
# The function BellMatrix is defined in A264428.
# Adds (1, 0, 0, 0, ...) as column 0.
BellMatrix(n -> (n+1)!, 9); # Peter Luschny, Jan 27 2016

nn = 9; a = x/(1 - x); f[list_] := Select[list, # > 0 &]; Flatten[Map[f, Drop[Range[0, nn]! CoefficientList[Series[Exp[y a], {x, 0, nn}], {x, y}], 1]]] (* Geoffrey Critzer, Dec 11 2011 *)
nn = 9; Flatten[Table[(j - k)! Binomial[j, k] Binomial[j - 1, k - 1], {j, nn}, {k, j}]] (* Jan Mangaldan, Mar 15 2013 *)
rows = 10;
t = Range[rows]!;
T[n_, k_] := BellY[n, k, t];
Table[T[n, k], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 23 2018, after Peter Luschny *)
T[n_, n_] := 1; T[n_, k_] /;0Oliver Seipel, Dec 06 2024 *)

use ntheory ":all"; say join ", ", map { my $n=$; map { stirling($n,$,3) } 1..$n; } 1..9; # Dana Jacobsen, Mar 16 2017

T(n,k) = Sum_{m=n..k} |S1(n,m)|*S2(m,k), k>=n>=1, with Stirling triangles S2(n,m):=A048993 and S1(n,m):=A048994.
T(n,k) = C(n,k)*(n-1)!/(k-1)!.
Sum_{k=1..n} T(n,k) = A000262(n).
n*Sum_{k=1..n} T(n,k) = A103194(n) = Sum_{k=1..n} T(n,k)*k^2.
E.g.f. column k: (x^(k-1)/(1-x)^(k+1))/(k-1)!, k>=1.
Recurrence from Sheffer (here Jabotinsky) a-sequence [1,1,0,...] (see the W. Lang link under A006232): T(n,k)=(n/k)*T(n-1,m-1) + n*T(n-1,m). - Wolfdieter Lang, Jun 29 2007
The e.g.f. is, umbrally, exp[(.)!* L(.,-t,1)*x] = exp[t*x/(1-x)]/(1-x)^2 where L(n,t,1) = Sum_{k=0..n} T(n+1,k+1)*(-t)^k = Sum_{k=0..n} binomial(n+1,k+1)* (-t)^k / k! is the associated Laguerre polynomial of order 1. - Tom Copeland, Nov 17 2007
For this Lah triangle, the n-th row polynomial is given umbrally by
  n! C(B.(x)+1+n,n) = (-1)^n C(-B.(x)-2,n), where C(x,n)=x!/(n!(x-n)!),
  the binomial coefficient, and B_n(x)= exp(-x)(xd/dx)^n exp(x), the n-th Bell / Touchard / exponential polynomial (cf. A008277). E.g.,
  2! C(-B.(-x)-2,2) = (-B.(x)-2)(-B.(x)-3) = B_2(x) + 5*B_1(x) + 6 = 6 + 6x + x^2.
  n! C(B.(x)+1+n,n) = n! e^(-x) Sum_{j>=0} C(j+1+n,n)x^j/j! is a corresponding Dobinski relation. See the Copeland link for the relation to inverse Mellin transform. - Tom Copeland, Nov 21 2011
The row polynomials are given by D^n(exp(x*t)) evaluated at x = 0, where D is the operator (1+x)^2*d/dx. Cf. A008277 (D = (1+x)*d/dx), A035342 (D = (1+x)^3*d/dx), A035469 (D = (1+x)^4*d/dx) and A049029 (D = (1+x)^5*d/dx). - Peter Bala, Nov 25 2011
T(n,k) = Sum_{i=k..n} A130534(n-1,i-1)*A008277(i,k). - Reinhard Zumkeller, Mar 18 2013
Let E(x) = Sum_{n >= 0} x^n/(n!*(n+1)!). Then a generating function is exp(t)*E(x*t) = 1 + (2 + x)*t + (6 + 6*x + x^2)*t^2/(2!*3!) + (24 + 36*x + 12*x^2 + x^3)*t^3/(3!*4!) + ... . - Peter Bala, Aug 15 2013
P_n(x) = L_n(1+x) = n!*Lag_n(-(1+x);1), where P_n(x) are the row polynomials of A059110; L_n(x), the Lah polynomials of A105278; and Lag_n(x;1), the Laguerre polynomials of order 1. These relations follow from the relation between the iterated operator (x^2 D)^n and ((1+x)^2 D)^n with D = d/dx. - Tom Copeland, Jul 23 2018
Dividing each n-th diagonal by n!, where the main diagonal is n=1, generates the Narayana matrix A001263. - Tom Copeland, Sep 23 2020
T(n,k) = A089231(n,n-k). - Ron L.J. van den Burg, Dec 12 2021
T(n,k) = T(n-1,k-1) + (n+k-1)*T(n-1,k). - Bérénice Delcroix-Oger, Jun 25 2025

Stirling comments and e.g.f.s from Wolfdieter Lang, Apr 11 2007

A054142 Triangular array binomial(2*n-k, k), k=0..n, n >= 0.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 6, 1, 1, 7, 15, 10, 1, 1, 9, 28, 35, 15, 1, 1, 11, 45, 84, 70, 21, 1, 1, 13, 66, 165, 210, 126, 28, 1, 1, 15, 91, 286, 495, 462, 210, 36, 1, 1, 17, 120, 455, 1001, 1287, 924, 330, 45, 1, 1, 19, 153, 680, 1820, 3003, 3003, 1716, 495, 55, 1

Offset: 0

Clark Kimberling

Row sums are odd-indexed Fibonacci numbers.
T(n,k) is the number of nondecreasing Dyck paths of semilength n+1, having k double rises. Mirror image of A085478. - Emeric Deutsch, May 31 2004
Diagonal sums are A052535. - Paul Barry, Jan 21 2005
Matrix inverse is the triangle of Salie numbers A098435. - Paul Barry, Jan 21 2005
Coefficients of Morgan-Voyce polynomial b(n,x); e.g., b(3,x)=x^3+5x^2+6x+1. See A172431 for coefficients of Morgan-Voyce polynomial B(n,x). - Clark Kimberling, Feb 13 2010
T(n,k) is the number of stack polyominoes of perimeter 2n+4 with k+1 columns. - Emanuele Munarini, Apr 07 2011
Roots of signed n-th polynomials are chaotic with respect to the operation (-2, x^2), with cycle lengths A003558(n). Example: starting with a root to x^3 - 5x^2 + 6x - 1 = 0; (2 + 2*cos(2*Pi/N) = 3.24697... = A116415; we obtain the trajectory (3.24697...-> 1.55495...-> 0.198062...; the 3 roots to the polynomial with cycle length 3 matching A003558(3) = 3. The operation (-2, x^2) is the reversal of the well known chaotic operation (x^2 - 2) [Kappraff, Adamson, 2004] starting with seed 2*cos(2*Pi/N). Check: given 2*cos(2*Pi/7) = 1.24697..., we obtain the 3-cycle using (x^2 - 2): (1.24697...-> -0.445041...-> 1.801937...; where the terms in either set are intermediate terms in the other, irrespective of sign. - Gary W. Adamson, Sep 22 2011
A054142 is jointly generated with A172431 as an array of coefficients of polynomials u(n,x): initially, u(1,x)=v(1,x)=1; for n>1, u(n,x)=x*u(n-1,x)+v(n-1,x) and v(n,x)=x*u(n-1,x)+(x+1)*v(n-1,x). See the Mathematica section of A172431. - Clark Kimberling, Mar 09 2012
Subtriangle of the triangle given by (0, 1, 0, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, 1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Apr 01 2012
The o.g.f. for row n of the array A(n, k) = binomial(2*n-k,k), k >= 0, n >= 0 is G(n,x) = Sum_{k=0..n} T(n, k)*x^k + (-x)^(2*n+1) * c(-x)^(2*n+1) / sqrt(1-4*(-x)), for n >= 0. Here c(x) is the o.g.f. of A000108 (Catalan). For powers of c(x) see the W. Lang link in A115139. For the alternating sign case replace x by -x. - Wolfdieter Lang, Sep 12 2016
Multiplying the n-th diagonal by A001147(n) generates A001497. - Tom Copeland, Oct 04 2016

			Triangle begins:
  1;
  1,  1;
  1,  3,  1;
  1,  5,  6,   1;
  1,  7, 15,  10,   1;
  1,  9, 28,  35,  15,   1;
  1, 11, 45,  84,  70,  21,   1;
  1, 13, 66, 165, 210, 126,  28,  1;
  1, 15, 91, 286, 495, 462, 210, 36, 1; ...
...
(0, 1, 0, 0, 0, 0, ...) DELTA (1, 0, 1, 0, 0, 0, ...) begins:
  1;
  0, 1;
  0, 1, 1;
  0, 1, 3,  1;
  0, 1, 5,  6,  1;
  0, 1, 7, 15, 10,  1;
  0, 1, 9, 28, 35, 15, 1. _Philippe Deléham_, Apr 01 2012

Vincenzo Librandi, Table of n, a(n) for n = 0..495
E. Barcucci, A. Del Lungo, S. Fezzi and R. Pinzani, Nondecreasing Dyck paths and q-Fibonacci numbers, Discrete Math., 170, 1997, 211-217.
D. Dumont and J. Zeng, Polynomes d'Euler et les fractions continues de Stieltjes-Rogers, Ramanujan J. 2 (1998) 3, 387-410.
Molly Fenn and Eric Sommers, A transitivity result for ad-nilpotent ideals in type A, arXiv:2101.04091 [math.RT], 2021.
J. L. Jacobsen, and J. Salas, Transfer Matrices and Partition-Function Zeros for Antiferromagnetic Potts Models IV. Chromatic polynomial with cyclic boundary conditions, J. Stat. Phys. 122 (2006) 705-760, arXiv:cond-mat/0407444 See Eq. 2.27. Mentions this sequence. - _N. J. A. Sloane_, Mar 14 2014
Jay Kappraff and Gary W. Adamson, Polygons and Chaos, 5th Interdispl Symm. Congress and Exh. Jul 8-14, Sydney, 2001 - [with commercial pop-ups].
Jay Kappraff and Gary W. Adamson, Polygons and Chaos, Journal of Dynamical Systems and Geometric Theories, Vol. 2 pp. 79-94, (Nov 2004).

These are the even-indexed rows of A011973, the odd-indexed rows form A053123.
Cf. A000108, A003558, A027989, A052535, A054142, A076756, A084938, A085478, A098435, A115139, A172431, A172991, A188648.
Cf. A001147, A001497.

Flat(List([0..12], n-> List([0..n], k-> Binomial(2*n-k,k) ))); # G. C. Greubel, Aug 01 2019

[Binomial(2*n-k,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Aug 01 2019

T:=(n,k)->binomial(2*n-k,k): seq(seq(T(n,k), k=0..n), n=0..11);

Flatten[Table[Binomial[2n - k, k], {n, 0, 11}, {k, 0, n}]] (* Emanuele Munarini, Apr 07 2011 *)

create_list(binomial(2*n-k,k),n,0,10,k,0,n); /* Emanuele Munarini, Apr 07 2011 */

T(n,k)=if(n<0,0,polcoeff(charpoly(matrix(n,n,i,j,-min(i,j))),k))

[[binomial(2*n-k,k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Aug 01 2019

G.f.: (1-t*z)/((1-t*z)^2-z). - Emeric Deutsch, May 31 2004
Column k has g.f.: (Sum_{j=0..k+1} binomial(k+1, 2j)*x^j)*x^k/(1-x)^(k+1). - Paul Barry, Jun 22 2005
Recurrence: T(n+2,k+2) = T(n+1,k+2) + 2*T(n+1,k+1) - T(n,k). - Emanuele Munarini, Apr 07 2011
T(n, k) = binomial(2*n-k, k) = A085478(n, n-k), for n >= 0, k = 0..n. - Wolfdieter Lang, Mar 25 2020

A099174 Triangle read by rows: coefficients of modified Hermite polynomials.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 0, 3, 0, 1, 3, 0, 6, 0, 1, 0, 15, 0, 10, 0, 1, 15, 0, 45, 0, 15, 0, 1, 0, 105, 0, 105, 0, 21, 0, 1, 105, 0, 420, 0, 210, 0, 28, 0, 1, 0, 945, 0, 1260, 0, 378, 0, 36, 0, 1, 945, 0, 4725, 0, 3150, 0, 630, 0, 45, 0, 1, 0, 10395, 0, 17325, 0, 6930, 0, 990, 0, 55, 0, 1

Offset: 0

Ralf Stephan, on a suggestion of Karol A. Penson, Oct 13 2004

Absolute values of A066325.
T(n,k) is the number of involutions of {1,2,...,n}, having k fixed points (0 <= k <= n). Example: T(4,2)=6 because we have 1243,1432,1324,4231,3214 and 2134. - Emeric Deutsch, Oct 14 2006
Riordan array [exp(x^2/2),x]. - Paul Barry, Nov 06 2008
Same as triangle of Bessel numbers of second kind, B(n,k) (see Cheon et al., 2013). - N. J. A. Sloane, Sep 03 2013
The modified Hermite polynomial h(n,x) (as in the Formula section) is the numerator of the rational function given by f(n,x) = x + (n-2)/f(n-1,x), where f(x,0) = 1. - Clark Kimberling, Oct 20 2014
Second lower diagonal T(n,n-2) equals positive triangular numbers A000217 \ {0}. - M. F. Hasler, Oct 23 2014
From James East, Aug 17 2015: (Start)
T(n,k) is the number of R-classes (equivalently, L-classes) in the D-class consisting of all rank k elements of the Brauer monoid of degree n.
For n < k with n == k (mod 2), T(n,k) is the rank (minimal size of a generating set) and idempotent rank (minimal size of an idempotent generating set) of the ideal consisting of all rank <= k elements of the Brauer monoid. (End)
This array provides the coefficients of a Laplace-dual sequence H(n,x) of the Dirac delta function, delta(x), and its derivatives, formed by taking the inverse Laplace transform of these modified Hermite polynomials. H(n,x) = h(n,D) delta(x) with h(n,x) as in the examples and the lowering and raising operators L = -x and R = -x + D = -x + d/dx such that L H(n,x) = n * H(n-1,x) and R H(n,x) = H(n+1,x). The e.g.f. is exp[t H(.,x)] = e^(t^2/2) e^(t D) delta(x) = e^(t^2/2) delta(x+t). - Tom Copeland, Oct 02 2016
Antidiagonals of this entry are rows of A001497. - Tom Copeland, Oct 04 2016
This triangle is the reverse of that in Table 2 on p. 7 of the Artioli et al. paper and Table 6.2 on p. 234 of Licciardi's thesis, with associations to the telephone numbers. - Tom Copeland, Jun 18 2018 and Jul 08 2018
See A344678 for connections to a Heisenberg-Weyl algebra of differential operators, matching and independent edge sets of the regular n-simplices with partially labeled vertices, and telephone switchboard scenarios. - Tom Copeland, Jun 02 2021

			h(0,x) = 1
h(1,x) = x
h(2,x) = x^2 + 1
h(3,x) = x^3 + 3*x
h(4,x) = x^4 + 6*x^2 + 3
h(5,x) = x^5 + 10*x^3 + 15*x
h(6,x) = x^6 + 15*x^4 + 45*x^2 + 15
From _Paul Barry_, Nov 06 2008: (Start)
Triangle begins
   1,
   0,  1,
   1,  0,  1,
   0,  3,  0,  1,
   3,  0,  6,  0,  1,
   0, 15,  0, 10,  0,  1,
  15,  0, 45,  0, 15,  0,  1
Production array starts
  0, 1,
  1, 0, 1,
  0, 2, 0, 1,
  0, 0, 3, 0, 1,
  0, 0, 0, 4, 0, 1,
  0, 0, 0, 0, 5, 0, 1 (End)

Alois P. Heinz, Rows n = 0..150, flattened
M. Artioli, G. Dattoli, S. Licciardi, and S. Pagnutti, Motzkin Numbers: an Operational Point of View, arXiv:1703.07262 [math.CO], 2017.
Paul Barry, Riordan array, orthogonal polynomials as moments, and Hankel transforms, arXiv:1102.0921 [math.CO], 2011.
G.-S. Cheon, J.-H. Jung and L. W. Shapiro, Generalized Bessel numbers and some combinatorial settings, Discrete Math., 313 (2013), 2127-2138.
T. Copeland, Juggling Zeros in the Matrix (Example II), 2020.
James East and Robert D. Gray, Diagram monoids and Graham-Houghton graphs: idempotents and generating sets of ideals, arXiv:1404.2359 [math.GR], 2014. See Theorem 8.4 and Table 7. - _James East_, Aug 17 2015
A. Horzela, P. Blasiak, G. E. H. Duchamp, K. A. Penson and A. I. Solomon, A product formula and combinatorial field theory, arXiv:quant-ph/0409152, 2004.
Alexander Kreinin, Combinatorial Properties of Mills' Ratio, arXiv:1405.5852 [math.CO], 2014. See Table 2. - _N. J. A. Sloane_, May 29 2014
S. Licciardi, Umbral Calculus, a Different Mathematical Language, arXiv:1803.03108 [math.CA], 2018.
R. Paris, A uniform asymptotic expansion for the incomplete gamma function, Journal of Computational and Applied Mathematics, 148 (2002), p. 223-239 (See p. 329 and A137286. From _Tom Copeland_, Jan 03 2016).
R. Sazdanovic, A categorification of the polynomial ring, slide presentation, 2011
S. Yang and Z. Qiao, The Bessel numbers and Bessel matrices, Jrn. Math. Rsch. and Exposition, July 2011, Vol. 31, No. 4, pp.627-636. DOI:10.3770/j.issn:1000-341X.2011.04.006.

Row sums (polynomial values at x=1) are A000085.
Polynomial values: A005425 (x=2), A202834 (x=3), A202879(x=4).
Cf. A000217, A001147, A059343, A060821, A066325.
Cf. A000384, A014105, A034839, A049403, A096713, A100861, A104556, A122848, A130757, A176230, A176231.
Cf. A137286.
Cf. A001497.
Cf. A111062, A344678.

T:=proc(n,k) if n-k mod 2 = 0 then n!/2^((n-k)/2)/((n-k)/2)!/k! else 0 fi end: for n from 0 to 12 do seq(T(n,k),k=0..n) od; # yields sequence in triangular form; Emeric Deutsch, Oct 14 2006

nn=10;a=y x+x^2/2!;Range[0,nn]!CoefficientList[Series[Exp[a],{x,0,nn}],{x,y}]//Grid  (* Geoffrey Critzer, May 08 2012 *)
H[0, x_] = 1; H[1, x_] := x; H[n_, x_] := H[n, x] = x*H[n-1, x]-(n-1)* H[n-2, x]; Table[CoefficientList[H[n, x], x], {n, 0, 11}] // Flatten // Abs (* Jean-François Alcover, May 23 2016 *)
T[ n_, k_] := If[ n < 0, 0, Coefficient[HermiteH[n, x I/Sqrt[2]] (Sqrt[1/2]/I)^n, x, k]]; (* Michael Somos, May 10 2019 *)

T(n,k)=if(k<=n && k==Mod(n,2), n!/k!/(k=(n-k)/2)!>>k) \\ M. F. Hasler, Oct 23 2014

import sympy
from sympy import Poly
from sympy.abc import x, y
def H(n, x): return 1 if n==0 else x if n==1 else x*H(n - 1, x) - (n - 1)*H(n - 2, x)
def a(n): return [abs(cf) for cf in Poly(H(n, x), x).all_coeffs()[::-1]]
for n in range(21): print(a(n)) # Indranil Ghosh, May 26 2017

def Trow(n: int) -> list[int]:
    row: list[int] = [0] * (n + 1); row[n] = 1
    for k in range(n - 2, -1, -2):
        row[k] = (row[k + 2] * (k + 2) * (k + 1)) // (n - k)
    return row  # Peter Luschny, Jan 08 2023

def A099174_triangle(dim):
    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]+(k+1)*M[n-1,k+1]
    return M
A099174_triangle(9)  # Peter Luschny, Oct 06 2012

h(k, x) = (-I/sqrt(2))^k * H(k, I*x/sqrt(2)), H(n, x) the Hermite polynomials (A060821, A059343).
T(n,k) = n!/(2^((n-k)/2)*((n-k)/2)!k!) if n-k >= 0 is even; 0 otherwise. - Emeric Deutsch, Oct 14 2006
G.f.: 1/(1-x*y-x^2/(1-x*y-2*x^2/(1-x*y-3*x^2/(1-x*y-4*x^2/(1-... (continued fraction). - Paul Barry, Apr 10 2009
E.g.f.: exp(y*x + x^2/2). - Geoffrey Critzer, May 08 2012
Recurrence: T(0,0)=1, T(0,k)=0 for k>0 and for n >= 1 T(n,k) = T(n-1,k-1) + (k+1)*T(n-1,k+1). - Peter Luschny, Oct 06 2012
T(n+2,n) = A000217(n+1), n >= 0. - M. F. Hasler, Oct 23 2014
The row polynomials P(n,x) = (a. + x)^n, umbrally evaluated with (a.)^n = a_n = aerated A001147, are an Appell sequence with dP(n,x)/dx = n * P(n-1,x). The umbral compositional inverses (cf. A001147) of these polynomials are given by the same polynomials signed, A066325. - Tom Copeland, Nov 15 2014
From Tom Copeland, Dec 13 2015: (Start)
The odd rows are (2x^2)^n x n! L(n,-1/(2x^2),1/2), and the even, (2x^2)^n n! L(n,-1/(2x^2),-1/2) in sequence with n= 0,1,2,... and L(n,x,a) = Sum_{k=0..n} binomial(n+a,k+a) (-x)^k/k!, the associated Laguerre polynomial of order a. The odd rows are related to A130757, and the even to A176230 and A176231. Other versions of this entry are A122848, A049403, A096713 and A104556, and reversed A100861, A144299, A111924. With each non-vanishing diagonal divided by its initial element A001147(n), this array becomes reversed, aerated A034839.
Create four shift and stretch matrices S1,S2,S3, and S4 with all elements zero except S1(2n,n) = 1 for n >= 1, S2(n,2n) = 1 for n >= 0, S3(2n+1,n) = 1 for n >= 1, and S4(n,2n+1) = 1 for n >= 0. Then this entry's lower triangular matrix is T = Id + S1 * (A176230-Id) * S2 + S3 * (unsigned A130757-Id) * S4 with Id the identity matrix. The sandwiched matrices have infinitesimal generators with the nonvanishing subdiagonals A000384(n>0) and A014105(n>0).
As an Appell sequence, the lowering and raising operators are L = D and R = x + dlog(exp(D^2/2))/dD = x + D, where D = d/dx, L h(n,x) = n h(n-1,x), and R h(n,x) = h(n+1,x), so R^n 1 = h(n,x). The fundamental moment sequence has the e.g.f. e^(t^2/2) with coefficients a(n) = aerated A001147, i.e., h(n,x) = (a. + x)^n, as noted above. The raising operator R as a matrix acting on o.g.f.s (formal power series) is the transpose of the production matrix P below, i.e., (1,x,x^2,...)(P^T)^n (1,0,0,...)^T = h(n,x).
For characterization as a Riordan array and associations to combinatorial structures, see the Barry link and the Yang and Qiao reference. For relations to projective modules, see the Sazdanovic link.
(End)
From the Appell formalism, e^(D^2/2) x^n = h_n(x), the n-th row polynomial listed below, and e^(-D^2/2) x^n = u_n(x), the n-th row polynomial of A066325. Then R = e^(D^2/2) * x * e^(-D^2/2) is another representation of the raising operator, implied by the umbral compositional inverse relation h_n(u.(x)) = x^n. - Tom Copeland, Oct 02 2016
h_n(x) = p_n(x-1), where p_n(x) are the polynomials of A111062, related to the telephone numbers A000085. - Tom Copeland, Jun 26 2018
From Tom Copeland, Jun 06 2021: (Start)
In the power basis x^n, the matrix infinitesimal generator M = A132440^2/2, when acting on a row vector for an o.g.f., is the matrix representation for the differential operator D^2/2.
e^{M} gives the coefficients of the Hermite polynomials of this entry.
The only nonvanishing subdiagonal of M, the second subdiagonal (1,3,6,10,...), gives, aside from the initial 0, the triangular numbers A000217, the number of edges of the n-dimensional simplices with (n+1) vertices. The perfect matchings of these simplices are the aerated odd double factorials A001147 noted above, the moments for the Hermite polynomials.
The polynomials are also generated from A036040 with x[1] = x, x[2] = 1, and the other indeterminates equal to zero. (End)

A122848 Exponential Riordan array (1, x(1+x/2)).

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 3, 1, 0, 0, 3, 6, 1, 0, 0, 0, 15, 10, 1, 0, 0, 0, 15, 45, 15, 1, 0, 0, 0, 0, 105, 105, 21, 1, 0, 0, 0, 0, 105, 420, 210, 28, 1, 0, 0, 0, 0, 0, 945, 1260, 378, 36, 1, 0, 0, 0, 0, 0, 945, 4725, 3150, 630, 45, 1, 0, 0, 0, 0, 0, 0, 10395, 17325, 6930, 990, 55, 1, 0, 0

Offset: 0

Paul Barry, Sep 14 2006

Entries are Bessel polynomial coefficients. Row sums are A000085. Diagonal sums are A122849. Inverse is A122850. Product of A007318 and A122848 gives A100862.
T(n,k) is the number of self-inverse permutations of {1,2,...,n} having exactly k cycles. - Geoffrey Critzer, May 08 2012
Bessel numbers of the second kind. For relations to the Hermite polynomials and the Catalan (A033184 and A009766) and Fibonacci (A011973, A098925, and A092865) matrices, see Yang and Qiao. - Tom Copeland, Dec 18 2013.
Also the inverse Bell transform of the double factorial of odd numbers Product_{k= 0..n-1} (2*k+1) (A001147). For the definition of the Bell transform see A264428 and for cross-references A265604. - Peter Luschny, Dec 31 2015

			Triangle begins:
    1
    0    1
    0    1    1
    0    0    3    1
    0    0    3    6    1
    0    0    0   15   10    1
    0    0    0   15   45   15    1
    0    0    0    0  105  105   21    1
    0    0    0    0  105  420  210   28    1
    0    0    0    0    0  945 1260  378   36    1
From _Gus Wiseman_, Jan 12 2021: (Start)
As noted above, a(n) is the number of set partitions of {1..n} into k singletons or pairs. This is also the number of set partitions of subsets of {1..n} into n - k pairs. In the first case, row n = 5 counts the following set partitions:
  {{1},{2,3},{4,5}}  {{1},{2},{3},{4,5}}  {{1},{2},{3},{4},{5}}
  {{1,2},{3},{4,5}}  {{1},{2},{3,4},{5}}
  {{1,2},{3,4},{5}}  {{1},{2,3},{4},{5}}
  {{1,2},{3,5},{4}}  {{1,2},{3},{4},{5}}
  {{1},{2,4},{3,5}}  {{1},{2},{3,5},{4}}
  {{1},{2,5},{3,4}}  {{1},{2,4},{3},{5}}
  {{1,3},{2},{4,5}}  {{1},{2,5},{3},{4}}
  {{1,3},{2,4},{5}}  {{1,3},{2},{4},{5}}
  {{1,3},{2,5},{4}}  {{1,4},{2},{3},{5}}
  {{1,4},{2},{3,5}}  {{1,5},{2},{3},{4}}
  {{1,4},{2,3},{5}}
  {{1,4},{2,5},{3}}
  {{1,5},{2},{3,4}}
  {{1,5},{2,3},{4}}
  {{1,5},{2,4},{3}}
In the second case, we have:
  {{1,2},{3,4}}  {{1,2}}  {}
  {{1,2},{3,5}}  {{1,3}}
  {{1,2},{4,5}}  {{1,4}}
  {{1,3},{2,4}}  {{1,5}}
  {{1,3},{2,5}}  {{2,3}}
  {{1,3},{4,5}}  {{2,4}}
  {{1,4},{2,3}}  {{2,5}}
  {{1,4},{2,5}}  {{3,4}}
  {{1,4},{3,5}}  {{3,5}}
  {{1,5},{2,3}}  {{4,5}}
  {{1,5},{2,4}}
  {{1,5},{3,4}}
  {{2,3},{4,5}}
  {{2,4},{3,5}}
  {{2,5},{3,4}}
(End)

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
Peter Bala, Generalized Dobinski formulas
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, Infinitesimal Generators, the Pascal Pyramid, and the Witt and Virasoro Algebras
H. Han and S. Seo, Combinatorial proofs of inverse relations and log-concavity for Bessel numbers, Eur. J. Combinat. 29 (7) (2008) 1544-1554. [From _R. J. Mathar_, Mar 20 2009]
Robert S. Maier, Boson Operator Ordering Identities from Generalized Stirling and Eulerian Numbers, arXiv:2308.10332 [math.CO], 2023. See p. 18.
S. Yang and Z. Qiao, The Bessel Numbers and Bessel Matrices, Journal of Mathematical Research & Exposition, July, 2011, Vol. 31, No. 4, pp. 627-636. [From _Tom Copeland_, Dec 18 2013]

Cf. A001497, A008275, A039755, A096713, A104556, A113278, A130757.
Row sums are A000085.
Column sums are A001515.
Same as A049403 but with a first column k = 0.
The same set partitions counted by number of pairs are A100861.
Reversing rows gives A111924 (without column k = 0).
A047884 counts standard Young tableaux by size and greatest row length.
A238123 counts standard Young tableaux by size and least row length.
A320663/A339888 count unlabeled multiset partitions into singletons/pairs.
A322661 counts labeled covering half-loop-graphs.
A339742 counts factorizations into distinct primes or squarefree semiprimes.
Cf. A000110, A000258, A320732, A321729, A339741, A339887.

# The function BellMatrix is defined in A264428.
BellMatrix(n -> `if`(n<2,1,0), 9); # Peter Luschny, Jan 27 2016

t[n_, k_] := k!*Binomial[n, k]/((2 k - n)!*2^(n - k)); Table[ t[n, k], {n, 0, 11}, {k, 0, n}] // Flatten
(* Second program: *)
rows = 12;
t = Join[{1, 1}, Table[0, rows]];
T[n_, k_] := BellY[n, k, t];
Table[T[n, k], {n, 0, rows}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 23 2018,after Peter Luschny *)
sbs[{}]:={{}};sbs[set:{i_,_}]:=Join@@Function[s,(Prepend[#1,s]&)/@sbs[Complement[set,s]]]/@Cases[Subsets[set],{i}|{i,_}];
Table[Length[Select[sbs[Range[n]],Length[#]==k&]],{n,0,6},{k,0,n}] (* Gus Wiseman, Jan 12 2021 *)

{T(n,k)=if(2*kn, 0, n!/(2*k-n)!/(n-k)!*2^(k-n))} /* Michael Somos, Oct 03 2006 */

# uses[inverse_bell_transform from A265605]
multifact_2_1 = lambda n: prod(2*k + 1 for k in (0..n-1))
inverse_bell_matrix(multifact_2_1, 9) # Peter Luschny, Dec 31 2015

Number triangle T(n,k) = k!*C(n,k)/((2k-n)!*2^(n-k)).
T(n,k) = A001498(k,n-k). - Michael Somos, Oct 03 2006
E.g.f.: exp(y(x+x^2/2)). - Geoffrey Critzer, May 08 2012
Triangle equals the matrix product A008275*A039755. Equivalently, the n-th row polynomial R(n,x) is given by the Type B Dobinski formula R(n,x) = exp(-x/2)*Sum_{k>=0} P(n,2*k+1)*(x/2)^k/k!, where P(n,x) = x*(x-1)*...*(x-n+1) denotes the falling factorial polynomial. Cf. A113278. - Peter Bala, Jun 23 2014
From Daniel Checa, Aug 28 2022: (Start)
E.g.f. for the m-th column: (x^2/2+x)^m/m!.
T(n,k) = T(n-1,k-1) + (n-1)*T(n-2,k-1) for n>1 and k=1..n, T(0,0) = 1. (End)

A001879 a(n) = (2n+2)!/(n!*2^(n+1)).

Original entry on oeis.org

1, 6, 45, 420, 4725, 62370, 945945, 16216200, 310134825, 6547290750, 151242416325, 3794809718700, 102776096548125, 2988412653476250, 92854250304440625, 3070380543400170000, 107655217802968460625, 3989575718580595893750, 155815096120119939628125

Offset: 0

N. J. A. Sloane

From Wolfdieter Lang, Oct 06 2008: (Start)
a(n) is the denominator of the n-th approximant to the continued fraction 1^2/(6+3^2/(6+5^2/(6+... for Pi-3. W. Lang, Oct 06 2008, after an e-mail from R. Rosenthal. Cf. A142970 for the corresponding numerators.
The e.g.f. g(x)=(1+x)/(1-2*x)^(5/2) satisfies (1-4*x^2)*g''(x) - 2*(8*x+3)*g'(x) -9*g(x) = 0 (from the three term recurrence given below). Also g(x)=hypergeom([2,3/2],[1],2*x). (End)
Number of descents in all fixed-point-free involutions of {1,2,...,2(n+1)}. A descent of a permutation p is a position i such that p(i) > p(i+1). Example: a(1)=6 because the fixed-point-free involutions 2143, 3412, and 4321 have 2, 1, and 3 descents, respectively. - Emeric Deutsch, Jun 05 2009
First differences of A193651. - Vladimir Reshetnikov, Apr 25 2016
a(n-2) is the number of maximal elements in the absolute order of the Coxeter group of type D_n. - Jose Bastidas, Nov 01 2021

J. Riordan, Combinatorial Identities, Wiley, 1968, p. 77 (Problem 10, values of Bessel polynomials).
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..200
Reinis Cirpons, James East, and James D. Mitchell, Transformation representations of diagram monoids, arXiv:2411.14693 [math.RA], 2024. See pp. 3, 33.
Selden Crary, Richard Diehl Martinez and Michael Saunders, The Nu Class of Low-Degree-Truncated Rational Multifunctions. Ib. Integrals of Matern-correlation functions for all odd-half-integer class parameters, arXiv:1707.00705 [stat.ME], 2017, Table 1.
Alexander Kreinin, Integer Sequences and Laplace Continued Fraction, Preprint 2016.
J. Riordan, Notes to N. J. A. Sloane, Jul. 1968

Cf. A002544, A001814, A001876, A001877, A001878.
Second column of triangle A001497. Equals (A001147(n+1)-A001147(n))/2.
Equals row sums of A163938.

[Factorial(2*n+2)/(Factorial(n)*2^(n+1)): n in [0..20]]; // Vincenzo Librandi, Nov 22 2011

restart: G(x):=(1-x)/(1-2*x)^(1/2): f[0]:=G(x): for n from 1 to 29 do f[n]:=diff(f[n-1],x) od:x:=0:seq(f[n],n=2..20); # Zerinvary Lajos, Apr 04 2009

Table[(2n+2)!/(n!2^(n+1)),{n,0,20}] (* Vincenzo Librandi, Nov 22 2011 *)

PARI
```
a(n)=if(n<0,0,(2*n+2)!/n!/2^(n+1))
```

E.g.f.: (1+x)/(1-2*x)^(5/2).
a(n)*n = a(n-1)*(2n+1)*(n+1); a(n) = a(n-1)*(2n+4)-a(n-2)*(2n-1), if n>0. - Michael Somos, Feb 25 2004
From Wolfdieter Lang, Oct 06 2008: (Start)
a(n) = (n+1)*(2*n+1)!! with the double factorials (2*n+1)!!=A001147(n+1).
D-finite with recurrence a(n) = 6*a(n-1) + ((2*n-1)^2)*a(n-2), a(-1)=0, a(0)=1. (End)
With interpolated 0's, e.g.f.: B(A(x)) where B(x)= x exp(x) and A(x)=x^2/2.
E.g.f.: -G(0)/2 where G(k) = 1 - (2*k+3)/(1 - x/(x - (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Dec 06 2012
G.f.: (1-x)/(2*x^2*Q(0)) - 1/(2*x^2), where Q(k) = 1 - x*(k+1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 20 2013
From Karol A. Penson, Jul 12 2013: (Start)
Integral representation as n-th moment of a signed function w(x) of bounded variation on (0,infinity),
  w(x) = -(1/4)*sqrt(2)*sqrt(x)*(1-x)*exp(-x/2)/sqrt(Pi):
  a(n) = Integral_{x>=0} x^n*w(x), n>=0.
  For x>1, w(x)>0. w(0)=w(1)=limit(w(x),x=infinity)=0. For x<1, w(x)<0.
Asymptotics: a(n)->(1/576)*2^(1/2+n)*(1152*n^2+1680*n+505)*exp(-n)*(n)^(n), for n->infinity. (End)
G.f.: 2F0(3/2,2;;2x). - R. J. Mathar, Aug 08 2015

Entry revised Aug 31 2004 (thanks to Ralf Stephan and Michael Somos)

E.g.f. in comment line corrected by Wolfdieter Lang, Nov 21 2011

A000457 Exponential generating function: (1+3x)/(1-2x)^(7/2).

Original entry on oeis.org

1, 10, 105, 1260, 17325, 270270, 4729725, 91891800, 1964187225, 45831035250, 1159525191825, 31623414322500, 924984868933125, 28887988983603750, 959493919812553125, 33774185977401870000, 1255977541034632040625

Offset: 0

N. J. A. Sloane

			G.f. = 1 + 10*x + 105*x^2 + 1260*x^3 + 17325*x^4 + 270270*x^5 + ... - _Michael Somos_, Dec 15 2023

L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 256.
F. N. David and D. E. Barton, Combinatorial Chance. Hafner, NY, 1962, p. 296.
C. Jordan, Calculus of Finite Differences. Eggenberger, Budapest and Röttig-Romwalter, Sopron 1939; Chelsea, NY, 1965, p. 172.
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).

G. C. Greubel, Table of n, a(n) for n = 0..200
Selden Crary, Richard Diehl Martinez, and Michael Saunders, The Nu Class of Low-Degree-Truncated Rational Multifunctions. Ib. Integrals of Matern-correlation functions for all odd-half-integer class parameters, arXiv:1707.00705 [stat.ME], 2017, Table 1.
H. W. Gould, Harris Kwong, and Jocelyn Quaintance, On Certain Sums of Stirling Numbers with Binomial Coefficients, J. Integer Sequences, 18 (2015), #15.9.6.
C. Jordan, On Stirling's Numbers, Tohoku Math. J., 37 (1933), 254-278.
Alexander Kreinin, Integer Sequences Connected to the Laplace Continued Fraction and Ramanujan's Identity, Journal of Integer Sequences, 19 (2016), #16.6.2.
J. Riordan, Notes to N. J. A. Sloane, Jul. 1968
Eric Weisstein's World of Mathematics, Stirling Number of the First Kind.

Equals (1/2)*A000906.
Third column of triangle A001497.
Second column (m=1) of unsigned Laguerre-Sonin a=1/2 triangle |A130757|.
Diagonal k=n-1 of triangle A134991.
Cf. A160473, A163939.

[Factorial(2*n+3)/(6*Factorial(n)*2^n): n in [0..30]]; // G. C. Greubel, May 15 2018

Table[(2n+3)!/(3!*n!*2^n), {n,0,30}] (* G. C. Greubel, May 15 2018 *)

for(n=0, 30, print1((2*n+3)!/(3!*n!*2^n), ", ")) \\ G. C. Greubel, May 15 2018

a(n) = (2n+3)!/( 3!*n!*2^n ).
a(n) = (n+1)*(2*n+3)!!/3, n>=0, with (2*n+3)!! = A001147(n+2).
a(n) = Sum_{j=0..n} (j + 1) * Eulerian2(n + 2, n - j). - Peter Luschny, Feb 13 2023

More terms from Sascha Kurz, Aug 15 2002

A004747 Triangle read by rows: the Bell transform of the triple factorial numbers A008544 without column 0.

Original entry on oeis.org

1, 2, 1, 10, 6, 1, 80, 52, 12, 1, 880, 600, 160, 20, 1, 12320, 8680, 2520, 380, 30, 1, 209440, 151200, 46480, 7840, 770, 42, 1, 4188800, 3082240, 987840, 179760, 20160, 1400, 56, 1, 96342400, 71998080, 23826880, 4583040, 562800, 45360, 2352, 72, 1

Offset: 1

Wolfdieter Lang

Previous name was: Triangle of numbers related to triangle A048966; generalization of Stirling numbers of second kind A008277, Bessel triangle A001497.
T(n,m) = S2p(-2; n,m), a member of a sequence of triangles including S2p(-1; n,m) = A001497(n-1,m-1) (Bessel triangle) and ((-1)^(n-m))*S2p(1; n,m) = A008277(n, m) (Stirling 2nd kind). T(n,1)= A008544(n-1).
T(n,m), n>=m>=1, enumerates unordered n-vertex m-forests composed of m plane (aka ordered) increasing (rooted) trees where vertices of out-degree r>=0 come in r+1 different types (like an (r+1)-ary vertex). Proof from the e.g.f. of the first column Y(z) = 1 - (1-3*x)^(1/3) and the F. Bergeron et al. eq. (8) Y'(z)= phi(Y(z)), Y(0) = 0, with out-degree o.g.f. phi(w)=1/(1-w)^2. - Wolfdieter Lang, Oct 12 2007
Also the Bell transform of the triple factorial numbers A008544 which adds a first column (1,0,0 ...) on the left side of the triangle. For the definition of the Bell transform see A264428. See A051141 for the triple factorial numbers A032031 and A203412 for the triple factorial numbers A007559 as well as A039683 and A132062 for the case of double factorial numbers. - Peter Luschny, Dec 21 2015

			Triangle begins:
       1;
       2,      1;
      10,      6,     1;
      80,     52,    12,    1;
     880,    600,   160,   20,   1;
   12320,   8680,  2520,  380,  30,  1;
  209440, 151200, 46480, 7840, 770, 42, 1;
Tree combinatorics for T(3,2)=6: Consider first the unordered forest of m=2 plane trees with n=3 vertices, namely one vertex with out-degree r=0 (root) and two different trees with two vertices (one root with out-degree r=1 and a leaf with r=0). The 6 increasing labelings come then from the forest with rooted (x) trees x, o-x (1,(3,2)), (2,(3,1)) and (3,(2,1)) and similarly from the second forest x, x-o (1,(2,3)), (2,(1,3)) and (3,(1,2)).

G. C. Greubel, Rows n = 1..50 of the triangle, flattened
F. Bergeron, Ph. Flajolet and B. Salvy, Varieties of increasing trees, 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.
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
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, Combinatorial Interpretation of Generalized Stirling Numbers, J. Int. Seqs. Vol. 12 (2009) #09.3.3.
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.
Index entries for sequences related to Bessel functions or polynomials

Cf. A015735 (row sums).
Cf. A007559, A008277, A008544, A032031, A039683.
Cf. A048966, A051141, A132062, A203412, A264428.
Triangles with the recurrence T(n,k) = (m*(n-1)-k)*T(n-1,k) + T(n-1,k-1): A010054 (m=1), A001497 (m=2), this sequence (m=3), A000369 (m=4), A011801 (m=5), A013988 (m=6).

function T(n,k) // T = A004747
  if k eq 0 then return 0;
  elif k eq n then return 1;
  else return (3*(n-1)-k)*T(n-1,k) + T(n-1,k-1);
  end if;
end function;
[T(n,k): k in [1..n], n in [1..12]]; // G. C. Greubel, Oct 03 2023

T := (n, m) -> 3^n/m!*(1/3*m*GAMMA(n-1/3)*hypergeom([1-1/3*m, 2/3-1/3*m, 1/3-1/3*m], [2/3, 4/3-n], 1)/GAMMA(2/3)-1/6*m*(m-1)*GAMMA(n-2/3)*hypergeom( [1-1/3*m, 2/3-1/3*m, 4/3-1/3*m], [4/3, 5/3-n], 1)/Pi*3^(1/2)*GAMMA(2/3)):
for n from 1 to 6 do seq(simplify(T(n,k)),k=1..n) od;
# Karol A. Penson, Feb 06 2004
# The function BellMatrix is defined in A264428.
# Adds (1,0,0,0, ..) as column 0.
BellMatrix(n -> mul(3*k+2, k=(0..n-1)), 9); # Peter Luschny, Jan 29 2016

(* First program *)
T[1,1]= 1; T[, 0]= 0; T[0, ]= 0; T[n_, m_]:= (3*(n-1)-m)*T[n-1, m]+T[n-1, m-1];
Flatten[Table[T[n, m], {n,12}, {m,n}] ][[1 ;; 45]] (* Jean-François Alcover, Jun 16 2011, after recurrence *)
(* Second program *)
f[n_, m_]:= m/n Sum[Binomial[k, n-m-k] 3^k (-1)^(n-m-k) Binomial[n+k-1, n-1], {k, 0, n-m}]; Table[n! f[n, m]/(m! 3^(n-m)), {n,12}, {m,n}]//Flatten (* Michael De Vlieger, Dec 23 2015 *)
(* Third program *)
rows = 12;
T[n_, m_]:= BellY[n, m, Table[Product[3k+2, {k, 0, j-1}], {j, 0, rows}]];
Table[T[n, m], {n,rows}, {m,n}]//Flatten (* Jean-François Alcover, Jun 22 2018 *)

# uses [bell_transform from A264428]
triplefactorial = lambda n: prod(3*k+2 for k in (0..n-1))
def A004747_row(n):
    trifact = [triplefactorial(k) for k in (0..n)]
    return bell_transform(n, trifact)
[A004747_row(n) for n in (0..10)] # Peter Luschny, Dec 21 2015

T(n, m) = n!*A048966(n, m)/(m!*3^(n-m));

T(n+1, m) = (3*n-m)*T(n, m)+ T(n, m-1), for n >= m >= 1, with T(n, m) = 0, for n

E.g.f. of m-th column: ( 1 - (1-3*x)^(1/3) )^m/m!.

Sum_{k=1..n} T(n, k) = A015735(n).

For a formula expressed as special values of hypergeometric functions 3F2 see the Maple program below. - Karol A. Penson, Feb 06 2004

T(n,1) = A008544(n-1). - Peter Luschny, Dec 23 2015

New name from Peter Luschny, Dec 21 2015

A163936 Triangle related to the o.g.f.s. of the right-hand columns of A130534 (E(x,m=1,n)).

Original entry on oeis.org

1, 1, 0, 2, 1, 0, 6, 8, 1, 0, 24, 58, 22, 1, 0, 120, 444, 328, 52, 1, 0, 720, 3708, 4400, 1452, 114, 1, 0, 5040, 33984, 58140, 32120, 5610, 240, 1, 0, 40320, 341136, 785304, 644020, 195800, 19950, 494, 1, 0, 362880, 3733920, 11026296, 12440064, 5765500, 1062500

Offset: 1

Johannes W. Meijer, Aug 13 2009

The asymptotic expansions of the higher-order exponential integral E(x,m=1,n) lead to triangle A130524, see A163931 for information on E(x,m,n). The o.g.f.s. of the right-hand columns of triangle A130534 have a nice structure: gf(p) = W1(z,p)/(1-z)^(2*p-1) with p = 1 for the first right-hand column, p = 2 for the second right-hand column, etc. The coefficients of the W1(z,p) polynomials lead to the triangle given above, n >= 1 and 1 <= m <= n. Our triangle is the same as A112007 with an extra right-hand column, see also the second Eulerian triangle A008517. The row sums of our triangle lead to A001147.
We observe that the row sums of the triangles A163936 (m=1), A163937 (m=2), A163938 (m=3) and A163939 (m=4) for z=1 lead to A001147, A001147 (minus a(0)), A001879 and A000457 which are the first four left-hand columns of the triangle of the Bessel coefficients A001497 or, if one wishes, the right-hand columns of A001498. We checked this phenomenon for a few more values of m and found that this pattern persists: m = 5 leads to A001880, m=6 to A001881, m=7 to A038121 and m=8 to A130563 which are the next left- (right-) hand columns of A001497 (A001498). An interesting phenomenon.
If one assumes the triangle not (1,1) based but (0,0) based, one has T(n, k) = E2(n, n-k), where E2(n, k) are the second-order Eulerian numbers A340556. - Peter Luschny, Feb 12 2021

			Triangle starts:
[ 1]      1;
[ 2]      1,       0;
[ 3]      2,       1,      0;
[ 4]      6,       8,      1,      0;
[ 5]     24,      58,     22,      1,      0;
[ 6]    120,     444,    328,     52,      1,     0;
[ 7]    720,    3708,   4400,   1452,    114,     1,   0;
[ 8]   5040,   33984,  58140,  32120,   5610,   240,   1,  0;
[ 9]  40320,  341136, 785304, 644020, 195800, 19950, 494,  1, 0;
The first few W1(z,p) polynomials are
W1(z,p=1) = 1/(1-z);
W1(z,p=2) = (1 + 0*z)/(1-z)^3;
W1(z,p=3) = (2 + 1*z + 0*z^2)/(1-z)^5;
W1(z,p=4) = (6 + 8*z + 1*z^2 + 0*z^3)/(1-z)^7.

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened

Row sums equal A001147.
A000142, A002538, A002539, A112008, A112485 are the first few left hand columns.
A000007, A000012, A005803(n+2), A004301, A006260 are the first few right hand columns.
Cf. A163931 (E(x,m,n)), A048994 (Stirling1) and A008517 (Euler).
Cf. A112007, A163937 (E(x,m=2,n)), A163938 (E(x,m=3,n)) and A163939 (E(x,m=4,n)).
Cf. A001497 (Bessel), A001498 (Bessel), A001147 (m=1), A001147 (m=2), A001879 (m=3) and A000457 (m=4), A001880 (m=5), A001881 (m=6) and A038121 (m=7).
Cf. A340556.

with(combinat): a := proc(n, m): add((-1)^(n+k+1)*binomial(2*n-1, k)*stirling1(m+n-k-1, m-k), k=0..m-1) end: seq(seq(a(n, m), m=1..n), n=1..9);  # Johannes W. Meijer, revised Nov 27 2012

Table[Sum[(-1)^(n + k + 1)*Binomial[2*n - 1, k]*StirlingS1[m + n - k - 1, m - k], {k, 0, m - 1}], {n, 1, 10}, {m, 1, n}] // Flatten (* G. C. Greubel, Aug 13 2017 *)

for(n=1,10, for(m=1,n, print1(sum(k=0,m-1,(-1)^(n+k+1)* binomial(2*n-1,k)*stirling(m+n-k-1,m-k, 1)), ", "))) \\ G. C. Greubel, Aug 13 2017

\\ assuming offset = 0:
E2poly(n,x) = if(n == 0, 1, x*(x-1)^(2*n)*deriv((1-x)^(1-2*n)*E2poly(n-1,x)));
{ for(n = 0, 9, print(Vec(E2poly(n,x)))) } \\ Peter Luschny, Feb 12 2021

a(n, m) = Sum_{k=0..(m-1)} (-1)^(n+k+1)*binomial(2*n-1,k)*Stirling1(m+n-k-1,m-k), for 1 <= m <= n.

Assuming offset = 0 the T(n, k) are the coefficients of recursively defined polynomials. T(n, k) = [x^k] x^n*E2poly(n, 1/x), where E2poly(n, x) = x*(x - 1)^(2*n)*d_{x}((1 - x)^(1 - 2*n)*E2poly(n - 1, x))) for n >= 1 and E2poly(0, x) = 1. - Peter Luschny, Feb 12 2021

Previous Showing 11-20 of 43 results. Next

cp's OEIS Frontend

A100861 Triangle of Bessel numbers read by rows: T(n,k) is the number of k-matchings of the complete graph K(n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A001498 Triangle a(n,k) (n >= 0, 0 <= k <= n) of coefficients of Bessel polynomials y_n(x) (exponents in increasing order).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Formula

A105278 Triangle read by rows: T(n,k) = binomial(n,k)*(n-1)!/(k-1)!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

Perl

Formula

Extensions

A054142 Triangular array binomial(2*n-k, k), k=0..n, n >= 0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

Maxima

PARI

Sage

Formula

A099174 Triangle read by rows: coefficients of modified Hermite polynomials.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

A000457 Exponential generating function: (1+3x)/(1-2x)^(7/2).