OFFSET
0,3
FORMULA
E.g.f.: ( (1/x) * Series_Reversion(x/(1 + 3*x*exp(x))) )^(1/3).
a(n) = (n!/(3*n+1)) * Sum_{k=0..n} 3^k * k^(n-k) * binomial(n+1/3,k)/(n-k)!.
PROG
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace((serreverse(x/(1+3*x*exp(x)))/x)^(1/3)))
(PARI) a(n) = n!*sum(k=0, n, 3^k*k^(n-k)*binomial(n+1/3, k)/(n-k)!)/(3*n+1);
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 10 2025
STATUS
approved