OFFSET
0,2
FORMULA
a(n) = Sum_{k=0..floor(n/3)} binomial(n+k-1,k) * binomial(3*n-2*k-1,n-3*k).
The g.f. exp( Sum_{k>=1} a(k) * x^k/k ) has integer coefficients and equals (1/x) * Series_Reversion( x * (1-x) * (1-x-x^3) ). See A368931.
MATHEMATICA
a[n_]:=SeriesCoefficient[1 / ( (1-x) * (1-x-x^3) )^n, {x, 0, n}]; Array[a, 23, 0] (* Stefano Spezia, Jun 02 2026*)
PROG
(PARI) a(n, s=3, t=1, u=1) = sum(k=0, n\s, binomial(t*n+k-1, k)*binomial((t+u+1)*n-(s-1)*k-1, n-s*k));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, May 02 2024
EXTENSIONS
a(23) from Stefano Spezia, Jun 02 2026
STATUS
approved
