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!.
a(n) ~ sqrt(58/sqrt(65)-6) * 2^(9*n+6) * exp(((sqrt(65)-9)*n + sqrt(65) - 7)/2) * n^(n-1) / (101*sqrt(65)-773)^(n+1). - Vaclav Kotesovec, Jan 31 2026
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
AUTHOR
Seiichi Manyama, Feb 06 2025
STATUS
approved
