OFFSET
0,5
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..500
FORMULA
a(n) = (1/(n+1)) * Sum_{k=0..floor(n/3)} binomial(n+1,k) * binomial(n-k-1,n-3*k).
a(n) = (1/(n+1)) * [x^n] (1 + x^3 / (1 - x)^2)^(n+1).
MATHEMATICA
a[n_]:=SeriesCoefficient[(1+x^3/(1-x)^2)^(n+1), {x, 0, n}]/(n+1); Table[a[n], {n, 0, 25}] (* Vincenzo Librandi, Oct 01 2025 *)
PROG
(PARI) my(N=40, x='x+O('x^N)); Vec(serreverse(x/(1+x^3/(1-x)^2))/x)
(Magma) a := [ n eq 0 select 1 else &+[Binomial(n+1, k) * Binomial(n-k-1, n-3*k) : k in [0..Floor(n/3)] ] div (n+1): n in [0..30] ]; a; // Vincenzo Librandi, Oct 01 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Sep 26 2025
STATUS
approved
