OFFSET
0,5
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..300
FORMULA
E.g.f.: (1/x) * Series_Reversion( x/(1 - x*log(1-x^3)) ).
a(n) = (n!)^2 * Sum_{k=0..floor(n/3)} 1/(3*k+1)! * |Stirling1(k,n-3*k)|/k!.
MATHEMATICA
Table[(n!)^2* Sum[1/(3*k+1)!*Abs[StirlingS1[k, n-3*k]]/k!, {k, 0, Floor[n/3]}], {n, 0, 20}] (* Vincenzo Librandi, Feb 04 2026 *)
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(serreverse(x/(1-x*log(1-x^3)))/x))
(Magma) [Factorial(n)^2*&+[1/Factorial(3*k+1)* Abs(StirlingFirst(k, n-3*k))/Factorial(k): k in [0..Floor(n/3)] ] : n in [0..24] ]; // Vincenzo Librandi, Feb 04 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 27 2026
STATUS
approved
