OFFSET
0,4
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
FORMULA
a(n) = (1/(n+1)) * Sum_{k=0..floor(n/3)} binomial(n+k,k) * binomial(n+3*k+1,n-3*k).
a(n) = (1/(n+1)) * [x^n] ((1 + x) / (1 - x^3 * (1 + x)^3))^(n+1).
MATHEMATICA
Table[SeriesCoefficient[((1+x)/(1-x^3*(1+x)^3))^(n+1), {x, 0, n}]/(n+1), {n, 0, 30}] (* Vincenzo Librandi, Oct 17 2025 *)
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serreverse(x*(1-x^3*(1+x)^3)/(1+x))/x)
(Magma) [1/(n+1)*&+[Binomial(n+k, k)*Binomial(n+3*k+1, n-3*k): k in [0..Floor(n/3)]]: n in [0..30]]; // Vincenzo Librandi, Oct 17 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Oct 04 2025
STATUS
approved
