OFFSET
0,4
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..300
FORMULA
E.g.f.: 1/(1+LambertW(-x)^3)^(1/3). - Vaclav Kotesovec, Sep 01 2014
a(n) ~ Gamma(5/6) * (n^(n-1/3) / (12^(1/3) * sqrt(Pi))) * (1 - 2^(7/6) * Gamma(1/3)^3 / (9 * Pi^(3/2) * sqrt(n))). - Vaclav Kotesovec, Sep 01 2014
MAPLE
with(combinat):
b:= proc(n, i) option remember; `if`(n=0, 1,
`if`(i>n, 0, add(b(n-i*j, i+3)*(i-1)!^j*
multinomial(n, n-i*j, i$j)/j!, j=0..n/i)))
end:
a:= a->add(b(j, 3)*n^(n-j)*binomial(n-1, j-1), j=0..n):
seq(a(n), n=0..20);
MATHEMATICA
CoefficientList[Series[1/(1+LambertW[-x]^3)^(1/3), {x, 0, 20}], x] * Range[0, 20]! (* Vaclav Kotesovec, Sep 01 2014 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Aug 31 2014
STATUS
approved