OFFSET
0,5
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..300
FORMULA
E.g.f. A(x) satisfies A(x) = 1 + (exp((x*A(x))^3) - 1)^2/(x*A(x))^2.
a(n) = (n!)^2 * Sum_{k=0..floor(n/3)} (2*(n-3*k))!/((n-3*k)! * (3*k+1)!) * Stirling2(n-2*k,2*(n-3*k))/(n-2*k)!.
MATHEMATICA
Table[(n!)^2*Sum[(2*(n-3*k))!/((n-3*k)!*(3*k+1)!)*Abs[StirlingS2[n-2*k, 2*(n-3*k)]/(n-2*k)!], {k, 0, Floor[n/3]}], {n, 0, 23}] (* Vincenzo Librandi, Jan 26 2026 *)
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(serreverse(x/(1+(exp(x^3)-1)^2/x^2))/x))
(Magma) [ Factorial(n)^2 * &+[Factorial(2*(n - 3*k)) / (Factorial(n - 3*k) * Factorial(3*k + 1)) * Abs(StirlingSecond(n - 2*k, 2*(n - 3*k)) / Factorial(n - 2*k)): k in [0..Floor(n/3)]] : n in [0..23] ]; // Vincenzo Librandi, Jan 26 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 25 2026
STATUS
approved
