OFFSET
0,5
LINKS
Eric Weisstein's World of Mathematics, Lambert W-Function.
FORMULA
a(n) = n! * Sum_{k=1..floor(n/4)} k^(k-1) * Stirling2(n-3*k,k)/(6^k * (n-3*k)!).
a(n) ~ sqrt(exp(r)*(3+r)/6 - 1/2) * n^(n-1) / (exp(n - 1/2) * r^(n - 3/2)), where r = 1.055946740386774512329072722491942575326676434456... is the root of the equation exp(r+1) - exp(1) = 6/r^3. - Vaclav Kotesovec, Jan 31 2026
PROG
(PARI) my(N=30, x='x+O('x^N)); concat([0, 0, 0, 0], Vec(serlaplace(-lambertw(x^3/6*(1-exp(x))))))
(PARI) a(n) = n!*sum(k=1, n\4, k^(k-1)*stirling(n-3*k, k, 2)/(6^k*(n-3*k)!));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Sep 24 2022
STATUS
approved
