OFFSET
0,2
FORMULA
G.f.: B(x)^3 where B(x) is the g.f. of A006319.
a(n) = 3 * Sum_{k=0..n} binomial(2*k+3,k) * binomial(n+k-1,n-k)/(2*k+3).
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec((1+x*((1-x-sqrt(1-6*x+x^2))/(2*x))^2)^3)
(PARI) a(n, r=3, s=2, t=2, u=0) = r*sum(k=0, n, binomial(t*k+u*(n-k)+r, k)*binomial(n+(s-1)*k-1, n-k)/(t*k+u*(n-k)+r));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Mar 31 2024
STATUS
approved