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 + log(1-(x*A(x))^3)^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)!) * |Stirling1(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[StirlingS1[n-2*k, 2*(n-3*k)]/(n -2*k)!], {k, 0, Floor[n/3]}], {n, 0, 23}] (* Vincenzo Librandi, Jan 27 2026 *)
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(serreverse(x/(1+log(1-x^3)^2/x^2))/x))
(Magma) seq := [(Factorial(n))^2 * &+[Factorial(2*(n-3*k))/(Factorial(n-3*k) * Factorial(3*k+1)) * Abs(StirlingFirst(n - 2*k, 2*(n - 3*k)) / Factorial(n - 2*k)): k in [0..Floor(n/3)]] : n in [0..25] ]; seq; // Vincenzo Librandi, Jan 27 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 25 2026
STATUS
approved
