cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

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

Views

Author

Keywords

Comments

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

Examples

			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
		

References

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

Crossrefs

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

Programs

  • Haskell
    a001498 n k = a001498_tabl !! n !! k
    a001498_row n = a001498_tabl !! n
    a001498_tabl = map reverse a001497_tabl
    -- Reinhard Zumkeller, Jul 11 2014
    
  • Magma
    /* 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
  • 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;
    # 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
  • Mathematica
    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 *)
  • PARI
    {T(n,k)=if(k<0||k>n, 0, binomial(n, k)*(n+k)!/2^k/n!)} /* Michael Somos, Oct 03 2006 */
    
  • PARI
    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
    

Formula

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

A002119 Bessel polynomial y_n(-2).

Original entry on oeis.org

1, -1, 7, -71, 1001, -18089, 398959, -10391023, 312129649, -10622799089, 403978495031, -16977719590391, 781379079653017, -39085931702241241, 2111421691000680031, -122501544009741683039, 7597207150294985028449, -501538173463478753560673
Offset: 0

Views

Author

Keywords

Comments

Absolute values give denominators of successive convergents to e using continued fraction 1+2/(1+1/(6+1/(10+1/(14+1/(18+1/(22+1/26...)))))).
Absolute values give number of different arrangements of nonnegative integers on a set of n 6-sided dice such that the dice can add to any integer from 0 to 6^n-1. For example when n=2, there are 7 arrangements that can result in any total from 0 to 35. Cf. A273013. The number of sides on the dice only needs to be the product of two distinct primes, of which 6 is the first example. - Elliott Line, Jun 10 2016
Absolute values give number of Krasner factorizations of (x^(6^n)-1)/(x-1) into n polynomials p_i(x), i=1,2,...,n, satisfying p_i(1)=6. In these expressions 6 can be replaced with any product of two distinct primes (Krasner and Ranulac, 1937). - William P. Orrick, Jan 18 2023
Absolute values give number of pairs (s, b) where s is a covering of the 1 X 2n grid with 1 X 2 dimers and equal numbers of red and blue 1 X 1 monomers and b is a bijection between the red monomers and the blue monomers that does not map adjacent monomers to each other. Ilya Gutkovskiy's formula counts such pairs by an inclusion-exclusion argument. The correspondence with Elliott Line's dice problem is that a dimer corresponds to a die containing an arithmetic progression of length 6 and a pair (r, b(r)), where r is a red monomer and b(r) its image under b, corresponds to a die containing the sum of an arithmetic progression of length 2 and an arithmetic progression of length 3. - William P. Orrick, Jan 19 2023

Examples

			Example from _William P. Orrick_, Jan 19 2023: (Start)
For n=2 the Bessel polynomial is y_2(x) = 1 + 3x + 3x^2 which satisfies y_2(-2) = -7.
The |a(2)|=7 dice pairs are
  {{0,1,2,3,4,5}, {0,6,12,18,24,30}},
  {{0,1,2,18,19,20}, {0,3,6,9,12,15}},
  {{0,1,2,9,10,11}, {0,3,6,18,21,24}},
  {{0,1,6,7,12,13}, {0,2,4,18,20,22}},
  {{0,1,12,13,24,25}, {0,2,4,6,8,10}},
  {{0,1,2,6,7,8}, {0,3,12,15,24,27}},
  {{0,1,4,5,8,9}, {0,2,12,14,24,26}}.
The corresponding Krasner factorizations of (x^36-1)/(x-1) are
  {(x^6-1)/(x-1), (x^36-1)/(x^6-1)},
  {((x^36-1)/(x^18-1))*((x^3-1)/(x-1)), (x^18-1)/(x^3-1)},
  {((x^18-1)/(x^9-1))*((x^3-1)/(x-1)), ((x^36-1)/(x^18-1))*((x^9-1)/(x^3-1))},
  {((x^18-1)/(x^6-1))*((x^2-1)/(x-1)), ((x^36-1)/(x^18-1))*((x^6-1)/(x^2-1))},
  {((x^36-1)/(x^12-1))*((x^2-1)/(x-1)), (x^12-1)/(x^2-1)},
  {((x^12-1)/(x^6-1))*((x^3-1)/(x-1)), ((x^36-1)/(x^12-1))*((x^6-1)/(x^3-1))},
  {((x^12-1)/(x^4-1))*((x^2-1)/(x-1)), ((x^36-1)/(x^12-1))*((x^4-1)/(x^2-1))}.
The corresponding monomer-dimer configurations, with dimers, red monomers, and blue monomers represented by the symbols '=', 'R', and 'B', and bijections between red and blue monomers given as sets of ordered pairs, are
  (==, {}),
  (B=R, {(3,1)}),
  (BBRR, {(3,1),(4,2)}),
  (RBBR, {(1,3),(4,2)}),
  (R=B, {(1,3)}),
  (BRRB, {(2,4),(3,1)}),
  (RRBB, {(1,3),(2,4)}).
(End)
		

References

  • L. Euler, 1737.
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 77.
  • 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).

Crossrefs

See also A033815.
Numerators of the convergents of e are A001517, which has a similar interpretation to a(n) in terms of monomer-dimer configurations, but omitting the restriction that adjacent monomers not be mapped to each other by the bijection.
Polynomial coefficients are in A001498.

Programs

  • Maple
    f:=proc(n) option remember; if n <= 1 then 1 else f(n-2)+(4*n-2)*f(n-1); fi; end;
    [seq(f(n), n=0..20)]; # This is for the unsigned version. - N. J. A. Sloane, May 09 2016
    seq(simplify((-1)^n*KummerU(-n, -2*n, -1)), n = 0..17); # Peter Luschny, May 10 2022
  • Mathematica
    Table[(-1)^k (2k)! Hypergeometric1F1[-k, -2k, -1]/k!, {k, 0, 10}] (* Vladimir Reshetnikov, Feb 16 2011 *)
    nxt[{n_,a_,b_}]:={n+1,b,a-2b(2n+1)}; NestList[nxt,{1,1,-1},20][[All,2]] (* Harvey P. Dale, Aug 18 2017 *)
  • PARI
    {a(n)= if(n<0, n=-n-1); sum(k=0, n, (2*n-k)!/ (k!*(n-k)!)* (-1)^(n-k) )} /* Michael Somos, Apr 02 2007 */
    
  • PARI
    {a(n)= local(A); if(n<0, n= -n-1); A= sqrt(1 +4*x +x*O(x^n)); n!*polcoeff( exp((A-1)/2)/A, n)} /* Michael Somos, Apr 02 2007 */
    
  • PARI
    {a(n)= local(A); if(n<0, n= -n-1); n+=2 ; for(k= 1, n, A+= x*O(x^k); A= truncate( (1+x)* exp(A) -1-A) ); A+= x*O(x^n); A-= A^2; -(-1)^n*n!* polcoeff( serreverse(A), n)} /* Michael Somos, Apr 02 2007 */
    
  • Sage
    A002119 = lambda n: hypergeometric([-n, n+1], [], 1)
    [simplify(A002119(n)) for n in (0..17)] # Peter Luschny, Oct 17 2014

Formula

D-finite with recurrence a(n) = -2(2n-1)*a(n-1) + a(n-2). - T. D. Noe, Oct 26 2006
If y = x + Sum_{k>=2} A005363(k)*x^k/k!, then y = x + Sum_{k>=2} a(k-2)(-y)^k/k!. - Michael Somos, Apr 02 2007
a(-n-1) = a(n). - Michael Somos, Apr 02 2007
a(n) = (1/n!)*Integral_{x>=-1} (-x*(1+x))^n*exp(-(1+x)). - Paul Barry, Apr 19 2010
G.f.: 1/Q(0), where Q(k) = 1 - x + 2*x*(k+1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 17 2013
Expansion of exp(x) in powers of y = x*(1 + x): exp(x) = 1 + y - y^2/2! + 7*y^3/3! - 71*y^4/4! + 1001*y^5/5! - .... E.g.f.: (1/sqrt(4*x + 1))*exp(sqrt(4*x + 1)/2 - 1/2) = 1 - x + 7*x^2/2! - 71*x^3/3! + .... - Peter Bala, Dec 15 2013
a(n) = hypergeom([-n, n+1], [], 1). - Peter Luschny, Oct 17 2014
a(n) = sqrt(Pi/exp(1)) * BesselI(1/2+n, 1/2) + (-1)^n * BesselK(1/2+n, 1/2) / sqrt(exp(1)*Pi). - Vaclav Kotesovec, Jul 22 2015
a(n) ~ (-1)^n * 2^(2*n+1/2) * n^n / exp(n+1/2). - Vaclav Kotesovec, Jul 22 2015
From G. C. Greubel, Aug 16 2017: (Start)
G.f.: (1/(1-t))*hypergeometric2f0(1, 1/2; -; -4*t/(1-t)^2).
E.g.f.: (1+4*t)^(-1/2) * exp((sqrt(1+4*t) - 1)/2). (End)
a(n) = Sum_{k=0..n} (-1)^k*binomial(n,k)*binomial(n+k,k)*k!. - Ilya Gutkovskiy, Nov 24 2017
a(n) = (-1)^n*KummerU(-n, -2*n, -1). - Peter Luschny, May 10 2022

Extensions

More terms from Vladeta Jovovic, Apr 03 2000

A232692 E.g.f. satisfies: A(x) = exp( 1/A(x)^3 * Integral A(x)^8 dx ).

Original entry on oeis.org

1, 1, 3, 24, 213, 3096, 46071, 967608, 20251809, 555747048, 15004870731, 508165972056, 16810393586733, 677183788645704, 26523956467895103, 1238567261126084856, 56056407696184372281, 2976966230117448265128, 152872356339113679491859, 9098430770913969095416728
Offset: 0

Views

Author

Paul D. Hanna, Dec 06 2013

Keywords

Comments

Compare e.g.f. to: B(x) = exp( 1/B(x)^3 * Integral B(x)^3 dx ) where B(y) = Bessel polynomial y_n(-3) (cf. A065923).
Note that G(x) = exp(1/G(x)^3 * Integral G(x)^7 dx) has negative coefficients.
CONJECTURE:
Given G(x,n,k) = G such that G = exp( 1/G^n * Integral G^k dx ) then G(x,n,k) consists solely of positive coefficients when k >= A047399(n) where A047399 lists numbers that are congruent to {0,3,6} mod 8.

Examples

			E.g.f.: A(x) = 1 + x + 3*x^2/2! + 24*x^3/3! + 213*x^4/4! + 3096*x^5/5! +...
Related expansions:
log(A(x)) = x + 2*x^2/2! + 17*x^3/3! + 120*x^4/4! + 1905*x^5/5! + 23640*x^6/6! +...
Integral A(x)^8 dx = x + 8*x^2/2! + 80*x^3/3! + 1032*x^4/4! + 16320*x^5/5! +...
1/A(x)^3 = 1 - 3*x + 3*x^2/2! - 24*x^3/3! + 117*x^4/4! - 2088*x^5/5! +...
		

Crossrefs

Programs

  • Maple
    seq(n! * coeff(series((3*LambertW(-1, (25*x-8)/3*exp(-8/3))/(25*x-8))^(1/5), x, n+1), x, n), n=0..20) # Vaclav Kotesovec, Jan 05 2014
  • Mathematica
    m = 20; A[] = 1; Do[A[x] = Exp[1/A[x]^3 Integrate[A[x]^8 + O[x]^m, x]] + O[x]^m // Normal, {m}]; CoefficientList[A[x], x] Range[0, m-1]! (* Jean-François Alcover, Nov 03 2019 *)
  • PARI
    {a(n)=local(A=1+x);for(i=1,n,A=exp(1/A^3*intformal(A^8+x*O(x^n))));n!*polcoeff(A,n)}
    for(n=0,30,print1(a(n),", "))

Formula

E.g.f.: (3*LambertW(-1, (25*x-8)/3*exp(-8/3))/(25*x-8))^(1/5). - Vaclav Kotesovec, Jan 05 2014

A078977 Denominator of n-th convergent to e^(2/3).

Original entry on oeis.org

1, 1, 19, 134, 153, 287, 3023, 163529, 2619487, 2783016, 5402503, 105430573, 9494154073, 237459282398, 246953436471, 484412718869, 13810509564803, 1740608617884047, 59194503517622401, 60935112135506448
Offset: 1

Views

Author

Benoit Cloitre, Dec 19 2002

Keywords

Crossrefs

Programs

  • Mathematica
    Denominator[Convergents[E^(2/3),20]] (* Harvey P. Dale, Dec 01 2013 *)
  • PARI
    a(n)=component(component(contfracpnqn(contfrac(exp(2/3),n)),1),2) \\ (Warning: this will give only a limited number of correct terms, depending on the precision used. - The Editors, Oct 13 2009. See A078976 for better code.)

Formula

Special cases : a(5k+1)=abs(A065923(3k)); a(5k+3)=abs(A065923(3k+2)) where A065923(n)=y(n, -3) where y(n, x)=sum (k=0, n, (n+k)!*(x/2)^k/((n-k)!*k!))
Showing 1-4 of 4 results.