OFFSET
0,2
FORMULA
E.g.f. A(x) satisfies A(x) = (1 + x*A(x))^3 * exp(3 * x * A(x)).
E.g.f.: B(x)^3, where B(x) is the e.g.f. of A377893.
a(n) = 3 * n! * Sum_{k=0..n} (3*n+3)^(k-1) * binomial(3*n+3,n-k)/k!.
MATHEMATICA
nmax=17; CoefficientList[(1/x)InverseSeries[Series[x*Exp[-3*x]/(1 + x)^3 , {x, 0, nmax}]], x]Range[0, nmax-1]! (* Stefano Spezia, Feb 06 2025 *)
PROG
(PARI) a(n) = 3*n!*sum(k=0, n, (3*n+3)^(k-1)*binomial(3*n+3, n-k)/k!);
CROSSREFS
KEYWORD
nonn,new
AUTHOR
Seiichi Manyama, Feb 06 2025
STATUS
approved