A385563 Expansion of 1/((1-x) * (1-5*x))^(3/2).
1, 9, 60, 360, 2055, 11403, 62132, 334260, 1781415, 9425295, 49581576, 259601004, 1353939405, 7038232425, 36484340400, 188665670880, 973545780195, 5014258620075, 25783103206100, 132378800689800, 678768332410245, 3476164133573505, 17782899991147500
Offset: 0
Links
- Paolo Xausa, Table of n, a(n) for n = 0..1000
Crossrefs
Programs
-
Mathematica
Module[{a, n}, RecurrenceTable[{a[n] == ((6*n+3)*a[n-1] - 5*(n+1)*a[n-2])/n, a[0] == 1, a[1] == 9}, a, {n, 0, 25}]] (* Paolo Xausa, Aug 21 2025 *)
-
PARI
my(N=30, x='x+O('x^N)); Vec(1/((1-x)*(1-5*x))^(3/2))
Formula
n*a(n) = (6*n+3)*a(n-1) - 5*(n+1)*a(n-2) for n > 1.
a(n) = (1/4)^n * Sum_{k=0..n} 5^k * (2*k+1) * (2*(n-k)+1) * binomial(2*k,k) * binomial(2*(n-k),n-k).
a(n) = Sum_{k=0..n} (2*k+1) * binomial(2*k,k) * binomial(n+2,n-k).
a(n) = Sum_{k=0..n} (-1)^k * 5^(n-k) * (2*k+1) * binomial(2*k,k) * binomial(n+2,n-k).
a(n) = binomial(n+2,2) * A002212(n+1).
a(n) = ((n+2)/2) * Sum_{k=0..floor(n/2)} 3^(n-2*k) * binomial(n+1,n-2*k) * binomial(2*k+1,k).
a(n) = Sum_{k=0..n} (3/2)^k * (-5/6)^(n-k) * (2*k+1) * binomial(2*k,k) * binomial(k,n-k).
a(n) ~ sqrt(n) * 5^(n + 3/2) / (4*sqrt(Pi)). - Vaclav Kotesovec, Aug 21 2025