A060821 -id:A060821 - cp's OEIS frontend

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

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

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)

A066325 Coefficients of unitary Hermite polynomials He_n(x).

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

Christian G. Bower, Dec 14 2001

Also number of involutions on n labeled elements with k fixed points times (-1)^(number of 2-cycles).
Also called normalized Hermite polynomials.
He_n(x) := H_n(x/sqrt(2)) / sqrt(2)^n, with the coefficients of H_n(x) given in A060821. See the Maple program. - Wolfdieter Lang, Jan 13 2020

			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;
    0,  -105,   0,  105,   0, -21,  0, 1;
  ...

F. Bergeron, G. Labelle and P. Leroux, Combinatorial Species and Tree-Like Structures, Cambridge, 1998, pp. 89,94 (2.3.41,54).

Robert Israel, Rows n=0..140 of triangle, flattened
P. Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, chapter 8.
P. Diaconis and A. Gamburd, Random matrices, magic squares and matching polynomials, The Electronic Journal of Combinatorics, Volume 11, Issue 2 (2004-6), Research Paper #R2.
E. Elizalde, Cosmology: techniques and observations, arXiv:gr-qc/0409076, 2004.
D. Foata, Une méthode combinatoire pour l'étude des fonctions spéciales, Journées sur les méthodes en mathématiques, Institut Henri Poincaré, Paris 2-3 april 2003.
R. Sazdanovic, A categorification of the polynomial ring, slide presentation, 2011. [_Tom Copeland_, Dec 27 2015]
Index entries for sequences related to Hermite polynomials

Row sums: A001464 (with different signs).
Row sums of absolute values: A000085.
Absolute values are given in A099174.
Cf. A159834, A001147, A060821 (Hermite H_n(x)).

Q:= [seq(orthopoly[H](n,x/sqrt(2))/2^(n/2), n=0..20)]:
seq(seq(coeff(Q[n+1],x,k),k=0..n),n=0..20); # Robert Israel, Jan 01 2016
# Alternative:
T := proc(n,k) option remember; if k > n then 0 elif n = k then 1 else
(T(n, k+2)*(k+2)*(k+1))/(k-n) fi end:
seq(print(seq(T(n, k), k = 0..n)), n = 0..10); # Peter Luschny, Jan 08 2023

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] // Expand; Table[CoefficientList[H[n, x], x], {n, 0, 11}] // Flatten (* Jean-François Alcover, May 11 2015 *)

for(n=0, 12, for(k=0,n, print1(if(Mod(n-k,2)==0, (-2)^((k-n)/2)*n!/(k!*((n-k)/2)!), 0), ", "))) \\ G. C. Greubel, Nov 23 2018

from sympy import Poly
from sympy.abc import x
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 Poly(H(n, x), x).all_coeffs()[::-1]
for n in range(21): print(a(n)) # Indranil Ghosh, May 26 2017

def A066325_row(n):
    T = [0]*(n+1)
    if n==1: return [1]
    for m in (1..n-1):
        a,b,c = 1,0,0
        for k in range(m,-1,-1):
            r = a - (k+1)*c
            if k < m : T[k+2] = u;
            a,b,c = T[k-1],a,b
            u = r
        T[1] = u;
    return T[1:]
for n in (1..11): A066325_row(n)  # Peter Luschny, Nov 01 2012

# uses[riordan_array from A256893]
riordan_array(exp(-x^2/2), x, 8, True) # Peter Luschny, Nov 23 2018

T(n, k) = (-2)^((k-n)/2)*n!/(k!*((n-k)/2)!) for n-k even, 0 otherwise.
E.g.f. of row polynomials {He_n(y)}:  A(x, y) = exp(x*y - x^2/2).
The umbral compositional inverses (cf. A001147) of the polynomials He(n,x) are given by the same polynomials unsigned, A099174. - Tom Copeland, Nov 15 2014
Exp(-D^2/2) x^n = He_n(x) = p_n(x+1) with D = d/dx and p_n(x), the row polynomials of A159834. These are an Appell sequence of polynomials with lowering and raising operators L = D and R = x - D. - Tom Copeland, Jun 26 2018

A059343 Triangle of nonzero coefficients of Hermite polynomials H_n(x) in increasing powers of x.

Original entry on oeis.org

1, 2, -2, 4, -12, 8, 12, -48, 16, 120, -160, 32, -120, 720, -480, 64, -1680, 3360, -1344, 128, 1680, -13440, 13440, -3584, 256, 30240, -80640, 48384, -9216, 512, -30240, 302400, -403200, 161280, -23040, 1024, -665280, 2217600, -1774080, 506880, -56320, 2048, 665280, -7983360, 13305600

Offset: 0

N. J. A. Sloane, Jan 27 2001

			1; 2*x; -2+4*x^2; -12*x+8*x^3; ...

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 801.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 50.

T. D. Noe, Rows n=0..100 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].
P. Diaconis and A. Gamburd, Random matrices, magic squares and matching polynomials
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3.
Dana G. Korssjoen, Biyao Li, Stefan Steinerberger, Raghavendra Tripathi, and Ruimin Zhang, Finding structure in sequences of real numbers via graph theory: a problem list, arXiv:2012.04625, Dec 08, 2020.
Eric Weisstein's World of Mathematics, Hermite Polynomial

Cf. A059344.

If initial zeros are included, same as A060821.

with(orthopoly): h:=proc(n) if n mod 2=0 then expand(x^2*H(n,x)) else expand(x*H(n,x)) fi end: seq(seq(coeff(h(n),x^(2*k)),k=1..1+floor(n/2)),n=0..14); # this gives the signed sequence

Flatten[ Table[ Coefficient[ HermiteH[n, x], x, k], {n, 0, 12}, {k, Mod[n, 2], n, 2}]] (* Jean-François Alcover, Jan 23 2012 *)

from sympy import hermite, Poly, Symbol
x = Symbol('x')
def a(n):
    return Poly(hermite(n, x), x).coeffs()[::-1]
for n in range(21): print(a(n)) # Indranil Ghosh, May 26 2017

Edited by Emeric Deutsch, Jun 05 2004

A059344 Triangle read by rows: row n consists of the nonzero coefficients of the expansion of 2^n x^n in terms of Hermite polynomials with decreasing subscripts.

Original entry on oeis.org

1, 1, 1, 2, 1, 6, 1, 12, 12, 1, 20, 60, 1, 30, 180, 120, 1, 42, 420, 840, 1, 56, 840, 3360, 1680, 1, 72, 1512, 10080, 15120, 1, 90, 2520, 25200, 75600, 30240, 1, 110, 3960, 55440, 277200, 332640, 1, 132, 5940, 110880, 831600, 1995840, 665280, 1, 156

Offset: 0

N. J. A. Sloane, Jan 27 2001

			Triangle begins
  1;
  1;
  1,     2;
  1,     6;
  1,    12,    12;
  1,    20,    60;
  1,    30,   180,   120;
  1,    42,   420,   840;
  1,    56,   840,  3360,  1680;
  1,    72,  1512, 10080, 15120;
x^2 = 1/2^2*(Hermite(2,x)+2*Hermite(0,x)); x^3 = 1/2^3*(Hermite(3,x)+6*Hermite(1,x)); x^4 = 1/2^4*(Hermite(4,x)+12*Hermite(2,x)+12*Hermite(0,x)); x^5 = 1/2^5*(Hermite(5,x)+20*Hermite(3,x)+60*Hermite(1,x)); x^6 = 1/2^6*(Hermite(6,x)+30*Hermite(4,x)+180*Hermite(2,x)+120*Hermite(0,x)). - _Vladeta Jovovic_, Feb 21 2003
1 = H(0); 2x = H(1); 4x^2 = H(2)+2H(0); 8x^3 = H(3)+6H(1); etc. where H(k)=Hermite(k,x).

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 801.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 50.

G. C. Greubel, Table of n, a(n) for the first 100 rows, 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].
Paul Barry, The Gamma-Vectors of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1804.05027 [math.CO], 2018.

Cf. A059343, A060821.

Cf. A119275 (signed row reverse).

Flatten[Table[n!/(k! * (n-2k)!), {n, 0, 13}, {k, 0, Floor[n/2]}]]
(* Second program: *)
row[n_] := Table[h[k], {k, n, Mod[n, 2], -2}] /. SolveAlways[2^n*x^n == Sum[h[k]*HermiteH[k, x], {k, Mod[n, 2], n, 2}], x] // First; Table[ row[n], {n, 0, 13}] // Flatten (* Jean-François Alcover, Jan 05 2016 *)

for(n=0,25, for(k=0,floor(n/2), print1(n!/(k!*(n-2*k)!), ", "))) \\ G. C. Greubel, Jan 07 2017

E.g.f.: exp(x^2+y*x). - Vladeta Jovovic, Feb 21 2003

a(n, k) = n!/(k! (n-2k)!). - Dean Hickerson, Feb 24 2003

More terms from Vladeta Jovovic, Feb 21 2003

Edited by Emeric Deutsch, Jun 05 2004

A067147 Triangle of coefficients for expressing x^n in terms of Hermite polynomials.

Original entry on oeis.org

1, 0, 1, 2, 0, 1, 0, 6, 0, 1, 12, 0, 12, 0, 1, 0, 60, 0, 20, 0, 1, 120, 0, 180, 0, 30, 0, 1, 0, 840, 0, 420, 0, 42, 0, 1, 1680, 0, 3360, 0, 840, 0, 56, 0, 1, 0, 15120, 0, 10080, 0, 1512, 0, 72, 0, 1, 30240, 0, 75600, 0, 25200, 0, 2520, 0, 90, 0, 1

Offset: 0

Christian G. Bower, Jan 03 2002

x^n = (1/2^n) * Sum_{k=0..n} a(n,k)*H_k(x).

These polynomials, H_n(x), are an Appell sequence, whose umbral compositional inverse sequence HI_n(x) consists of the same polynomials signed with the e.g.f. e^{-t^2} e^{xt}. Consequently, under umbral composition H_n(HI.(x)) = x^n = HI_n(H.(x)). Other differently scaled families of Hermite polynomials are A066325, A099174, and A060821. See Griffin et al. for a relation to the Catalan numbers and matrix integration. - Tom Copeland, Dec 27 2020

			Triangle begins with:
    1;
    0,   1;
    2,   0,   1;
    0,   6,   0,   1;
   12,   0,  12,   0,   1;
    0,  60,   0,  20,   0,   1;
  120,   0, 180,   0,  30,   0,   1;

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 801. (Table 22.12)

G. C. Greubel, Rows n = 0..100 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].
M. Griffin, K. Ono, L. Rolen, and D. Zagier, Jensen polynomials for the Riemann zeta function and other sequences, arXiv:1902.07321 [math.NT], 2019.
Index entries for sequences related to Hermite polynomials

Row sums give A047974. Columns 0-2: A001813, A000407, A001814. Cf. A048854, A060821.

Cf. A060821, A066325, and A099174.

[[Round(Factorial(n)*(1+(-1)^(n+k))/(2*Factorial(k)*Gamma((n-k+2)/2))): k in [0..n]]: n in [0..10]]; // G. C. Greubel, Jun 09 2018

T := proc(n, k) (n - k)/2; `if`(%::integer, (n!/k!)/%!, 0) end:
for n from 0 to 11 do seq(T(n, k), k=0..n) od; # Peter Luschny, Jan 05 2021

Table[n!*(1+(-1)^(n+k))/(2*k!*Gamma[(n-k+2)/2]), {n,0,20}, {k,0,n}]// Flatten (* G. C. Greubel, Jun 09 2018 *)

T(n, k) = round(n!*(1+(-1)^(n+k))/(2*k! *gamma((n-k+2)/2)))
for(n=0,20, for(k=0,n, print1(T(n, k), ", "))) \\ G. C. Greubel, Jun 09 2018

{T(n,k) = if(k<0 || nMichael Somos, Jan 15 2020 */

E.g.f. (rel to x): A(x, y) = exp(x*y + x^2).
Sum_{ k>=0 } 2^k*k!*T(m, k)*T(n, k) = T(m+n, 0) = |A067994(m+n)|. - Philippe Deléham, Jul 02 2005
T(n, k) = 0 if n-k is odd; T(n, k) = n!/(k!*((n-k)/2)!) if n-k is even. - Philippe Deléham, Jul 02 2005
T(n, k) = n!/(k!*2^((n-k)/2)*((n-k)/2)!)*2^((n+k)/2)*(1+(-1)^(n+k))/2^(k+1).
T(n, k) = A001498((n+k)/2, (n-k)/2)2^((n+k)/2)(1+(-1)^(n+k))/2^(k+1). - Paul Barry, Aug 28 2005
Exponential Riordan array (e^(x^2),x). - Paul Barry, Sep 12 2006
G.f.: 1/(1-x*y-2*x^2/(1-x*y-4*x^2/(1-x*y-6*x^2/(1-x*y-8*x^2/(1-... (continued fraction). - Paul Barry, Apr 10 2009
The n-th row entries may be obtained from D^n(exp(x*t)) evaluated at x = 0, where D is the operator sqrt(1+4*x)*d/dx. - Peter Bala, Dec 07 2011
As noted in the comments this is an Appell sequence of polynomials, so the lowering and raising operators defined by L H_n(x) = n H_{n-1}(x) and R H_{n}(x) = H_{n+1}(x) are L = D_x, the derivative, and R = D_t log[e^{t^2} e^{xt}] |{t = D_x} = x + 2 D_x, and the polynomials may also be generated by e^{-D^2} x^n = H_n(x). - _Tom Copeland, Dec 27 2020

cp's OEIS Frontend

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

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Python

Sage

Formula

A066325 Coefficients of unitary Hermite polynomials He_n(x).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Sage

Sage

Formula

A059343 Triangle of nonzero coefficients of Hermite polynomials H_n(x) in increasing powers of x.

Views

Author

Keywords

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Python

Extensions

A059344 Triangle read by rows: row n consists of the nonzero coefficients of the expansion of 2^n x^n in terms of Hermite polynomials with decreasing subscripts.

Views

Author

Keywords

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A067147 Triangle of coefficients for expressing x^n in terms of Hermite polynomials.

Views

Author

A163322 The 3rd Hermite Polynomial evaluated at n: H_3(n) = 8n^3 - 12n.