A331793 Expansion of ((1 - 5*x)/sqrt(1 - 10*x + 9*x^2) - 1)/(8*x^2).

1, 10, 87, 740, 6285, 53550, 458115, 3934600, 33913881, 293244050, 2542684463, 22101612780, 192530903461, 1680415209270, 14692052109915, 128653303453200, 1128147127156785, 9905115333850650, 87066787614156807, 766127762539955700, 6747880819438628541

Offset

0,2

Links

Seiichi Manyama, Table of n, a(n) for n = 0..1000

Formula

a(n) = (2/(n+2)) * A331516(n) = Sum_{k=0..n} 4^k * binomial(n+1,k) * binomial(n+1,k+1).

n * (n+2) * a(n) = (n+1) * (5 * (2*n+1) * a(n-1) - 9 * n * a(n-2)) for n>1.

a(n) ~ 3^(2*n + 3) / (2^(5/2) * sqrt(Pi*n)). - Vaclav Kotesovec, Jan 26 2020

From Seiichi Manyama, Aug 23 2025: (Start)

a(n) = Sum_{k=0..floor(n/2)} 4^k * 5^(n-2*k) * binomial(n+1,n-2*k) * binomial(2*k+1,k).

a(n) = Sum_{k=0..n} 2^k * binomial(n+1,k+1) * binomial(2*k+2,k+2). (End)

From Seiichi Manyama, Aug 25 2025: (Start)

a(n) = [x^n] (1+5*x+4*x^2)^(n+1).

E.g.f.: exp(5*x) * BesselI(1, 4*x) / 2, with offset 1. (End)

Mathematica

a[n_] := Sum[4^k * Binomial[n + 1, k] * Binomial[n + 1, k + 1], {k, 0, n}]; Array[a, 21, 0] (* Amiram Eldar, May 05 2021 *)

Prog

(PARI) my(N=30, x='x+O('x^N)); Vec(((1-5*x)/sqrt(1-10*x+9*x^2)-1)/(8*x^2))
(PARI) a(n) = sum(k=0, n, 4^k*binomial(n+1, k)*binomial(n+1, k+1));

Crossrefs

Column 5 of A331791.

Cf. A331516.

Sequence in context: A251193 A198858 A383944 * A121115 A292998 A114648

Adjacent sequences: A331790 A331791 A331792 * A331794 A331795 A331796

Keyword

nonn

Author

Seiichi Manyama, Jan 26 2020

Status

approved