OFFSET
0,3
COMMENTS
a(206) is negative. - Vaclav Kotesovec, Jan 23 2026
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..300
FORMULA
E.g.f. A(x) satisfies A(x) = 1/(1 - log(1+(x*A(x))^3)/(x*A(x))^2).
a(n) = (1/(n+1)) * Sum_{k=0..floor(n/3)} (2*n-3*k)! * Stirling1(n-2*k,n-3*k)/(n-2*k)!.
MATHEMATICA
Table[1/(n+1) * Sum[(2*n-3*k)! * StirlingS1[n-2*k, n-3*k] / (n-2*k)!, {k, 0, n/3}], {n, 0, 20}]
(* or *)
nmax = 20; CoefficientList[1/x*InverseSeries[Series[x - Log[1 + x^3]/x, {x, 0, nmax+1}], x], x] * Range[0, nmax]! (* Vaclav Kotesovec, Jan 23 2026 *)
PROG
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(serreverse(x-log(1+x^3)/x)/x))
CROSSREFS
KEYWORD
changed,sign
AUTHOR
Seiichi Manyama, Jan 22 2026
STATUS
approved
