A037953 a(n) = binomial(n, floor((n-5)/2)).

0, 0, 0, 0, 0, 1, 1, 7, 8, 36, 45, 165, 220, 715, 1001, 3003, 4368, 12376, 18564, 50388, 77520, 203490, 319770, 817190, 1307504, 3268760, 5311735, 13037895, 21474180, 51895935, 86493225, 206253075

Offset

0,8

Links

G. C. Greubel, Table of n, a(n) for n = 0..1000

Formula

a(n) = binomial(n, A004526(n-5)). - Wesley Ivan Hurt, Nov 28 2013

(n+6)*(n-5)*a(n) = 2*(n)*a(n-1) + 4*(n-1)*n*a(n-2). - R. J. Mathar, Jul 26 2015

From G. C. Greubel, Jun 21 2022: (Start)

G.f.: ((1 +x -6*x^2 -5*x^3 +9*x^4 +5*x^5 -2*x^6) - (1 +x -4*x^2 -3*x^3 +3*x^4 +x^5)*sqrt(1-4*x^2))/(2*x^6*sqrt(1-4*x^2)).

E.g.f.: BesselI(5, 2*x) + BesselI(6, 2*x). (End)

Maple

A037953:=n->binomial(n, floor((n-5)/2)); seq(A037953(n), n=0..50); # Wesley Ivan Hurt, Nov 28 2013

Mathematica

Table[Binomial[n, Floor[(n-5)/2]], {n, 0, 40}] (* Harvey P. Dale, Oct 11 2012 *)

Prog

(Magma) [Binomial(n, Floor((n-5)/2)): n in [0..40]]; // G. C. Greubel, Jun 21 2022
(SageMath) [binomial(n, (n-5)//2) for n in (0..40)] # G. C. Greubel, Jun 21 2022
(PARI) a(n)=binomial(n, (n-5)\2) \\ Charles R Greathouse IV, Oct 23 2023

Crossrefs

Cf. A035951, A035952, A035954, A035955, A035956, A035957.

Cf. A004526, A089940, A101491.

Sequence in context: A154745 A048064 A122605 * A296636 A041106 A303732

Adjacent sequences: A037950 A037951 A037952 * A037954 A037955 A037956

Keyword

nonn

Author

N. J. A. Sloane

Status

approved