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