A387367 Expansion of 1/(1 - 10*x + x^2)^(5/2).

1, 25, 435, 6475, 88270, 1137150, 14081970, 169370850, 1991916135, 23013193775, 262062237437, 2948690451525, 32845189782860, 362721036253100, 3975956599494420, 43300257350934900, 468875116313950845, 5051523021827188725, 54177811767428268535, 578700942412768257775

Offset

0,2

Links

Vincenzo Librandi, Table of n, a(n) for n = 0..1000

Formula

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

a(n) = (binomial(n+4,2)/6) * A387369(n).

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

Mathematica

CoefficientList[Series[1/(1-10*x+x^2)^(5/2), {x, 0, 33}], x] (* Vincenzo Librandi, Aug 29 2025 *)

Prog

(PARI) my(N=20, x='x+O('x^N)); Vec(1/(1-10*x+x^2)^(5/2))
(Magma) R<x> := PowerSeriesRing(Rationals(), 34); f := 1/(1 - 10*x + x^2)^(5/2); coeffs := [ Coefficient(f, n) : n in [0..33] ]; coeffs; // Vincenzo Librandi, Aug 29 2025
(Python)
from sympy import Symbol, series, Rational
x = Symbol('x')
def A387367(n): return int(series(1/(1 - 10*x + x**2)**Rational(5, 2), x, 0, n+1).removeO().coeff(x, n)) # Aitzaz Imtiaz, Oct 03 2025

Crossrefs

Cf. A000332, A387341.

Cf. A387369.

Sequence in context: A019742 A019722 A180800 * A004346 A021324 A092430

Adjacent sequences: A387364 A387365 A387366 * A387368 A387369 A387370

Keyword

nonn

Author

Seiichi Manyama, Aug 27 2025

Status

approved