Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A001498
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A001498 Triangle of coefficients of Bessel polynomials (exponents in increasing order). +0
27
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 (list; table; graph; listen)
OFFSET

0,5

COMMENT

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 T(n,k) has factorization [C(n,k)][C(k,n-k)]Diag((2n-1)!!) The triangle T(n-k,k) is A100861, which gives coefficients of scaled Hermite polynomials. - Paul Barry (pbarry(AT)wit.ie), May 21 2005

Related to k-matchings of the complete graph K_n by T(n,k)=A100861(n+k,k). Related to the Morgan-Voyce polynomials by T(n,k)=(2k-1)!!*A085478(n,k). - Paul Barry (pbarry(AT)wit.ie), Aug 17 2005

Related to Hermite polynomials by T(n,k)=(-1)^k*A060821(n+k,n-k)/2^n; - Paul Barry (pbarry(AT)wit.ie), Aug 28 2005

The row polynomials, the Bessel polynomials y(n,x):=sum(a(n,m)*x^m,m=0..n) (called y_{n}(x) in the Grosswald reference) solve (x^2)*diff(y(n,x),x$2)+2*(x+1)*diff(y(n,x),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. W. Lang, Sep 14 2007.

REFERENCES

E. Grosswald, Bessel Polynomials, Lecture Notes Math. vol. 698 1978 p. 18.

B. Leclerc, Powers of staircase Schur functions and symmetric analogues of Bessel polynomials, Discrete Math., 153 (1996), 213-227.

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

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

LINKS

T. D. Noe, Rows n=0..50 of triangle, flattened

A. Burstein and T. Mansour, Words restricted by patterns with at most 2 distinct letters.

W. Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.

L. A. Sz\'ekely, P. L. Erd\"os and M. A. Steel, The combinatorics of evolutionary trees

Eric Weisstein's World of Mathematics, Modified Spherical Bessel Function of the Second Kind

Index entries for sequences related to Bessel functions or polynomials

W. Lang, First ten rows.

FORMULA

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) [ Leonard Smiley (smiley(AT)math.uaa.alaska.edu) ]

a(n, m)= A001497(n, n-m) = A001147(m)*binomial(n+m, 2*m) for n >= m >= 0 else 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).

B(r, s) = (r+s)!/[2^s*(r-s)!*s! ]. - R. Stephan, Apr 20 2004

EXAMPLE

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 combinatorics: 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). W. Lang, Sep 14 2007.

MAPLE

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;

PROGRAM

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

CROSSREFS

Columns from left edge include A000217, A050534.

Columns 1-6 from right edge are A001147, A001879, A000457, A001880, A001881, A038121.

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. A001497.

Adjacent sequences: A001495 A001496 A001497 this_sequence A001499 A001500 A001501

Sequence in context: A110640 A094040 A039798 this_sequence A117279 A049323 A084144

KEYWORD

nonn,tabl,nice

AUTHOR

njas

page 1

Search completed in 0.003 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified May 16 23:01 EDT 2008. Contains 139884 sequences.


AT&T Labs Research