OFFSET
0,2
FORMULA
a(n) = Sum_{k=0..n} binomial(n+3*k+1,k) * binomial(3*n+3*k+1,n-k)/(n+3*k+1).
MATHEMATICA
terms = 21; A[_] = 0; Do[A[x_] = 1/((1-x*A[x]^3)*(1 -x*A[x])^2) + O[x]^terms // Normal, terms]; CoefficientList[A[x], x] (* Stefano Spezia, Jun 14 2025 *)
PROG
(PARI) a(n) = sum(k=0, n, binomial(n+3*k+1, k)*binomial(3*n+3*k+1, n-k)/(n+3*k+1));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Dec 17 2024
STATUS
approved
