A392994 Expansion of e.g.f. (1/x) * Series_Reversion( x/(1 + x*log(1-x^3)^2) ).

1, 0, 0, 0, 0, 0, 0, 5040, 0, 0, 3628800, 0, 0, 5708102400, 610248038400, 0, 17435658240000, 6046686277632000, 0, 92585437533388800, 68932223565004800000, 3576365952019660800000, 786800509444325376000, 1040543673740120309760000, 171243758878374085263360000

Offset

0,8

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)*log(1-(x*A(x))^3)^2.

a(n) = (n!)^2 * Sum_{k=0..floor(n/3)} (2*(n-3*k))!/((n-3*k)! * (3*k+1)!) * |Stirling1(k,2*(n-3*k))|/k!.

Mathematica

Table[(n!)^2* Sum[(2*(n-3*k))!/((n-3*k)!*(3*k+1)!) *Abs[StirlingS1[k, 2*(n-3*k)]]/k!, {k, 0, Floor[n/3]}], {n, 0, 24}] (* Vincenzo Librandi, Feb 01 2026 *)

Prog

(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(serreverse(x/(1+x*log(1-x^3)^2))/x))
(Magma) [Factorial(n)^2* &+[Factorial(2*(n-3*k))/(Factorial(n-3*k) * Factorial(3*k+1)) * Abs (StirlingFirst(k, 2*(n-3*k)))/Factorial(k): k in [0..Floor(n/3)] ] : n in [0..24] ]; // Vincenzo Librandi, Feb 01 2026

Crossrefs

Cf. A392888, A392938.

Sequence in context: A237690 A269125 A392992 * A111030 A068378 A036458

Adjacent sequences: A392991 A392992 A392993 * A392995 A392996 A392997

Keyword

nonn

Author

Seiichi Manyama, Jan 29 2026

Status

approved