OFFSET
0,3
FORMULA
a(n) = n! * Sum_{k=0..floor(n/2)} (Product_{j=0..k-1} (3*j+1)) * |Stirling1(n-k,k)|/(n-k)!.
a(n) ~ sqrt(2*Pi) * 3^n * r^(n + 1/3) * n^(n - 1/6) / (Gamma(1/3) * (1/(3*r-1) + log(3*r/(3*r-1)))^(1/3) * exp(n)), where r = 0.67741090606491714108599100770859788... is the root of the equation 3*r*(1 - exp(-r)) = 1. - Vaclav Kotesovec, Jan 23 2026
MATHEMATICA
With[{nn=20}, CoefficientList[Series[1/CubeRoot[1+3x Log[1-x]], {x, 0, nn}], x] Range[0, nn]!] (* Harvey P. Dale, Oct 21 2025 *)
PROG
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(1/(1+3*x*log(1-x))^(1/3)))
(PARI) a(n) = n!*sum(k=0, n\2, prod(j=0, k-1, 3*j+1)*abs(stirling(n-k, k, 1))/(n-k)!);
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Aug 24 2024
STATUS
approved
