A107102 -id:A107102 - cp's OEIS frontend

A001497 Triangle of coefficients of Bessel polynomials (exponents in decreasing order).

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

Offset: 0

N. J. A. Sloane

The (reverse) Bessel polynomials P(n,x):=Sum_{m=0..n} a(n,m)*x^m, the row polynomials, called Theta_n(x) in the Grosswald reference, solve x*(d^2/dx^2)P(n,x) - 2*(x+n)*(d/dx)P(n,x) + 2*n*P(n,x) = 0.
With the related Sheffer associated polynomials defined by Carlitz as
B(0,x) = 1
B(1,x) = x
B(2,x) = x + x^2
B(3,x) = 3 x + 3 x^2 + x^3
B(4,x) = 15 x + 15 x^2 + 6 x^3 + x^4
... (see Mathworld reference), then P(n,x) = 2^n * B(n,x/2) are the Sheffer polynomials described in A119274. - Tom Copeland, Feb 10 2008
Exponential Riordan array [1/sqrt(1-2x), 1-sqrt(1-2x)]. - Paul Barry, Jul 27 2010
From Vladimir Kruchinin, Mar 18 2011: (Start)
For B(n,k){...} the Bell polynomial of the second kind we have
B(n,k){f', f'', f''', ...} = T(n-1,k-1)*(1-2*x)^(k/2-n), where f(x) = 1-sqrt(1-2*x).
The expansions of the first few rows are:
1/sqrt(1-2*x);
1/(1-2*x)^(3/2), 1/(1-2*x);
3/(1-2*x)^(5/2), 3/(1-2*x)^2, 1/(1-2*x)^(3/2);
15/(1-2*x)^(7/2), 15/(1-2*x)^3, 6/(1-2*x)^(5/2), 1/(1-2*x)^2. (End)
Also the Bell transform of A001147 (whithout column 0 which is 1,0,0,...). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 19 2016
Antidiagonals of A099174 are rows of this entry. Dividing each diagonal by its first element generates A054142. - Tom Copeland, Oct 04 2016
The row polynomials p_n(x) of A107102 are (-1)^n B_n(1-x), where B_n(x) are the modified Carlitz-Bessel polynomials above, e.g., (-1)^2 B_2(1-x) = (1-x) + (1-x)^2 = 2 - 3 x + x^2 = p_2(x). - Tom Copeland, Oct 10 2016
a(n-1,m-1) counts rooted unordered binary forests with n labeled leaves and m roots. - David desJardins, Feb 23 2019
From Jianing Song, Nov 29 2021: (Start)
The polynomials P_n(x) = Sum_{k=0..n} T(n,k)*x^k satisfy: P_n(x) - (d/dx)P_n(x) = x*P_{n-1}(x) for n >= 1.
{P(n,x)} are related to the Fourier transform of 1/(1+x^2)^(n+1) and x/(1+x^2)^(n+2):
(i) For n >= 0, real number t, we have Integral_{x=-oo..oo} exp(-i*t*x)/(1+x^2)^(n+1) dx = Pi/(2^n*n!) * P_n(|t|) * exp(-|t|);
(ii) For n >= 0, real number t, we have Integral_{x=-oo..oo} x*exp(-i*t*x)/(1+x^2)^(n+2) dx = Pi/(2^(n+1)*(n+1)!) * ((-t)*P_n(-|t|)) * exp(-|t|). (End)
Suppose that f(x) is an n-times differentiable function defined on (a,b) for 0 <= a < b <= +oo, then for n >= 1, the n-th derivative of f(sqrt(x)) on (a^2,b^2) is Sum_{k=1..n} ((-1)^(n-k)*T(n-1,k-1)*f^(k)(sqrt(x))) / (2^n*x^(n-(k/2))), where f^(k) is the k-th derivative of f. - Jianing Song, Nov 30 2023

			Triangle begins
        1,
        1,       1,
        3,       3,      1,
       15,      15,      6,      1,
      105,     105,     45,     10,     1,
      945,     945,    420,    105,    15,    1,
    10395,   10395,   4725,   1260,   210,   21,   1,
   135135,  135135,  62370,  17325,  3150,  378,  28,  1,
  2027025, 2027025, 945945, 270270, 51975, 6930, 630, 36, 1
Production matrix begins
       1,      1,
       2,      2,      1,
       6,      6,      3,     1,
      24,     24,     12,     4,     1,
     120,    120,     60,    20,     5,    1,
     720,    720,    360,   120,    30,    6,   1,
    5040,   5040,   2520,   840,   210,   42,   7,  1,
   40320,  40320,  20160,  6720,  1680,  336,  56,  8, 1,
  362880, 362880, 181440, 60480, 15120, 3024, 504, 72, 9, 1
This is the exponential Riordan array A094587, or [1/(1-x),x], beheaded.
- _Paul Barry_, Mar 18 2011

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

T. D. Noe, Rows n = 0..50 of triangle, flattened
Peter Bala, Generalized Dobinski formulas
Paul Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, chapter 8.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers
E. Deutsch, L. Ferrari and S. Rinaldi, Production Matrices, Advances in Applied Mathematics, 34 (2005) pp. 101-122.
O. Frink and H. L. Krall, A new class of orthogonal polynomials, Trans. Amer. Math. Soc. 65,100-115, 1945. [From _Roger L. Bagula_, Feb 15 2009]
E. Grosswald, Bessel Polynomials, Lecture Notes Math. vol. 698 1978 p. 18.
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.
B. Leclerc, Powers of staircase Schur functions and symmetric analogues of Bessel polynomials, Discrete Math., 153 (1996), 213-227.
Robert S. Maier, Boson Operator Ordering Identities from Generalized Stirling and Eulerian Numbers, arXiv:2308.10332 [math.CO], 2023. See p. 19.
Toufik Mansour, Matthias Schork and Mark Shattuck, The Generalized Stirling and Bell Numbers Revisited, Journal of Integer Sequences, Vol. 15 (2012), #12.8.3.
Toufik Mansour, Matthias Schork and Mark Shattuck, On the Stirling numbers associated with the meromorphic Weyl algebra, Applied Mathematics Letters, Volume 25, Issue 11, November 2012, Pages 1767-1771. - From _N. J. A. Sloane_, Sep 15 2012
W. Mlotkowski and A. Romanowicz, A family of sequences of binomial type, Probability and Mathematical Statistics, Vol. 33, Fasc. 2 (2013), pp. 401-408.
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.
Feng Qi and Bai-Ni Guo, "Some Properties and Generalizations of the Catalan, Fuss, and Fuss-Catalan Numbers", Mathematical Analysis and Applications : Selected Topics (2018), Wiley, Ch. 5, 101-133.
Feng Qi, X.-T. Shi and F.-F. Liu, Several formulas for special values of the Bell polynomials of the second kind and applications, Preprint 2015.
Alexander Stoimenow, On the number of chord diagrams, Discr. Math. 218 (2000), 209-233. Lemma 2.2.
Eric Weisstein's World of Mathematics, Bessel Polynomial
Index entries for sequences related to Bessel functions or polynomials

Reflected version of A001498 which is considered the main entry.
Other versions of this same triangle are given in A144299, A111924 and A100861.
Row sums give A001515. a(n, 0)= A001147(n) (double factorials).
Cf. A104556 (matrix inverse). A039683, A122850.
Cf. A245066 (central terms).
Cf. A054142, A099174, A107102.

a001497 n k = a001497_tabl !! n !! k
a001497_row n = a001497_tabl !! n
a001497_tabl = [1] : f [1] 1 where
   f xs z = ys : f ys (z + 2) where
     ys = zipWith (+) ([0] ++ xs) (zipWith (*) [z, z-1 ..] (xs ++ [0]))
-- Reinhard Zumkeller, Jul 11 2014

/* As triangle */ [[Factorial(2*n-k)/(Factorial(k)*Factorial(n-k)*2^(n-k)): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Aug 12 2015

f := proc(n) option remember; if n <=1 then (1+x)^n else expand((2*n-1)*x*f(n-1)+f(n-2)); fi; end;
row := n -> seq(coeff(f(n), x, n - k), k = 0..n): seq(row(n), n = 0..9);

m = 9; Flatten[ Table[(n + k)!/(2^k*k!*(n - k)!), {n, 0, m}, {k, n, 0, -1}]] (* Jean-François Alcover, Sep 20 2011 *)
y[n_, x_] := Sqrt[2/(Pi*x)]*E^(1/x)*BesselK[-n-1/2, 1/x]; t[n_, k_] := Coefficient[y[n, x], x, k]; Table[t[n, k], {n, 0, 9}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Mar 01 2013 *)

T(k, n) = if(n>k||k<0||n<0,0,(2*k-n)!/(n!*(k-n)!*2^(k-n))) /* Ralf Stephan */

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

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

a(n, m) = (2*n-m)!/(m!*(n-m)!*2^(n-m)) if n >= m >= 0 else 0 (from Grosswald, p. 7).
a(n, m)= 0, n= m >= 0 (from Grosswald p. 23, (19)).
E.g.f. for m-th column: ((1-sqrt(1-2*x))^m)/(m!*sqrt(1-2*x)).
G.f.: 1/(1-xy-x/(1-xy-2x/(1-xy-3x/(1-xy-4x/(1-.... (continued fraction). - Paul Barry, Jan 29 2009
T(n,k) = if(k<=n, C(2n-k,2(n-k))*(2(n-k)-1)!!,0) = if(k<=n, C(2n-k,2(n-k))*A001147(n-k),0). - Paul Barry, Mar 18 2011
Row polynomials for n>=1 are given by 1/t*D^n(exp(x*t)) evaluated at x = 0, where D is the operator 1/(1-x)*d/dx. - Peter Bala, Nov 25 2011
The matrix product A039683*A008277 gives a signed version of this triangle. Dobinski-type formula for the row polynomials: R(n,x) = (-1)^n*exp(x)*Sum_{k = 0..inf} k*(k-2)*(k-4)*...*(k-2*(n-1))*(-x)^k/k!. Cf. A122850. - Peter Bala, Jun 23 2014

A100862 Triangle read by rows: T(n,k) is the number of k-matchings of the corona K'(n) of the complete graph K(n) and the complete graph K(1); in other words, K'(n) is the graph constructed from K(n) by adding for each vertex v a new vertex v' and the edge vv'.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 6, 6, 1, 1, 10, 21, 10, 1, 1, 15, 55, 55, 15, 1, 1, 21, 120, 215, 120, 21, 1, 1, 28, 231, 665, 665, 231, 28, 1, 1, 36, 406, 1736, 2835, 1736, 406, 36, 1, 1, 45, 666, 3990, 9891, 9891, 3990, 666, 45, 1, 1, 55, 1035, 8310, 29505, 45297, 29505, 8310, 1035, 55, 1, 1, 66, 1540, 16005, 77715, 173712, 173712, 77715, 16005, 1540, 66, 1

Offset: 0

Emeric Deutsch, Jan 08 2005

Row n has n+1 terms.

Row sums yield A005425.

			T(3,2)=6 because in the graph with vertex set {A,B,C,a,b,c} and edge set {AB,AC,BC,Aa,Bb,Cc} we have the following six 2-matchings: {Aa,BC},{Bb,AC},{Cc,AB},{Aa,Bb},{Aa,Cc} and {Bb,Cc}.
Triangle starts:
  1;
  1,  1;
  1,  3,  1;
  1,  6,  6,  1;
  1, 10, 21, 10,  1;
  1, 15, 55, 55, 15,  1;

Paul Barry, On a Generalization of the Narayana Triangle, J. Int. Seq., Vol. 14 (2011), Article 11.4.5.
Paul Barry, Constructing Exponential Riordan Arrays from Their A and Z Sequences, Journal of Integer Sequences, Vol. 17 (2014), Article 14.2.6.
Paul Barry, The Gamma-Vectors of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1804.05027 [math.CO], 2018.

Cf. A005425.

Cf. A107102 (matrix inverse).

P[0]:=1: for n from 1 to 11 do P[n]:=sort(expand((1+t)*P[n-1]+(n-1)*t*P[n-2])) od: for n from 0 to 11 do seq(coeff(t*P[n],t^k),k=1..n+1) od; # gives the sequence in triangular form

P[0] = 1; P[1] = 1+t; P[n_] := P[n] = (1+t) P[n-1] + (n-1) t P[n-2];
Table[CoefficientList[P[n], t], {n, 0, 11}] // Flatten (* Jean-François Alcover, Jul 23 2018 *)

{T(n,k)=local(X=x+x*O(x^n),Y=y+y*O(y^k)); n!*polcoeff(polcoeff(exp(X+Y*X^2/2+X*Y),n,x),k,y)} \\ Paul D. Hanna, Jul 18 2005

T(n,k) = T(n,n-k).
E.g.f.: exp(z+t*z+t*z^2/2).
Row generating polynomial P[n] = [ -i*sqrt(t/2)]^n*H(n, i(1+t)/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]=(1+t)P[n-1]+(n-1)tP[n-2].
From Fabián Pereyra, Jan 05 2022: (Start)
T(n,k) = T(n-1,k) + (n-1)*T(n-2,k-1) + T(n-1,k-1), n>=0, 0<=k<=n; T(0,0) = 1.
T(n,k) = Sum_{j=k..n} C(n,j)*B(j,j-k), where B are the Bessel numbers A100861. (End)

cp's OEIS Frontend

A107104 Absolute row sums of triangle A107102, which is the matrix inverse of A100862.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

PARI

Formula

A107103 Column 1 of triangle A107102, which is the matrix inverse of A100862.

Views

Author

Keywords

Comments

Crossrefs

Programs

PARI

Formula

A001497 Triangle of coefficients of Bessel polynomials (exponents in decreasing order).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

A100862 Triangle read by rows: T(n,k) is the number of k-matchings of the corona K'(n) of the complete graph K(n) and the complete graph K(1); in other words, K'(n) is the graph constructed from K(n) by adding for each vertex v a new vertex v' and the edge vv'.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula