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