OFFSET
0,4
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..500
FORMULA
a(n) = (1/(n+1)) * Sum_{k=0..floor(n/2)} binomial(n+1,k) * binomial(n+k-1,n-2*k).
a(n) = (1/(n+1)) * [x^n] (1 + x^2 / (1 - x)^3)^(n+1).
D-finite with recurrence -3*(3*n+2)*(3*n+1) *(5329638178135*n -6740091898537)*(n+1)*a(n) +(957297182797754*n^4 -2683976978115556*n^3 +2789459719303369*n^2 -911914488473519*n +40440551391222) *a(n-1) +(-2136238900546091*n^4 +15699572589884723*n^3 -39252401421363832*n^2 +39901025026117972*n -14020651307189502) *a(n-2) +6*(574225490711666*n^4 -4669351681023436*n^3 +14992495899086774*n^2 -21603111733355204*n +10995922036219955) *a(n-3) -(n-4) *(4457503378649747*n^3 -27859307256243307*n^2 +51102685535867403*n -25763046371721693) *a(n-4) +(n-4) *(n-5) *(4482834765302906*n^2 -23916955442793690*n +21414530310542049) *a(n-5) -31*(n-4) *(n-5) *(n-6) *(15703570527019*n -18045386290574) *a(n-6)=0. - R. J. Mathar, Sep 29 2025
MATHEMATICA
Table[Coefficient[(1+x)^(n-1) (1+x^2+x^3)^(n+1), x, n]/(n+1), {n, 0, 25}] (* Vincenzo Librandi, Oct 01 2025 *)
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serreverse(x/(1+x^2/(1-x)^3))/x)
(Magma) a := [ n eq 0 select 1 else (n eq 1 select 0 else &+[Binomial(n+1, k)*Binomial(n+k-1, 3*k-1) : k in [1..Floor(n/2)]]/(n+1)) : n in [0..35] ]; a; // Vincenzo Librandi, Oct 01 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Sep 26 2025
STATUS
approved
