OFFSET
0,6
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..300
FORMULA
E.g.f. A(x) satisfies A(x) = 1 + (x*A(x))^2 * (exp(x*A(x)) - 1)^3.
a(n) = (n!/(n+1)) * Sum_{k=0..floor(n/5)} (3*k)! * binomial(n+1,k) * Stirling2(n-2*k,3*k)/(n-2*k)!.
MATHEMATICA
Table[(n!/(n+1))* Sum[(3*k)!*Binomial[n+1, k]*StirlingS2[n-2*k, 3*k]/(n-2*k)!, {k, 0, Floor[n/5]}], {n, 0, 21}] (* Vincenzo Librandi, Feb 06 2026 *)
PROG
(PARI) a(n) = n!*sum(k=0, n\5, (3*k)!*binomial(n+1, k)*stirling(n-2*k, 3*k, 2)/(n-2*k)!)/(n+1);
(Magma) [(Factorial(n)/(n+1)) * &+[Factorial(3*k)* Binomial(n+1, k) * StirlingSecond(n - 2*k, 3*k) / Factorial(n - 2*k): k in [0..Floor(n/5)]]: n in [0..25] ]; // Vincenzo Librandi, Feb 06 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 24 2026
STATUS
approved
