OFFSET
1,3
COMMENTS
In general, for k > 0, if e.g.f. A(x) satisfies A'(x) = exp(A(x)*A'(x)^k), then a(n) ~ (k+1)^(2*n-3) * exp(n/k - 1) * n^(n-2) / (sqrt(k) * (k + exp(1 + 1/k))^(n - 3/2)). - Vaclav Kotesovec, Jun 18 2026
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..320
FORMULA
For k >= 0, e.g.f. A(x) satisfies A'(x) = exp( A(x)*A'(x)^k ), with A(0)=0.
Let R(x) = Series_Reversion( (1-(1-k*(k+1)*x) * exp(-(k+1)*x))/(k+1)^2 ).
A(x) = R(x) * exp(-k*R(x)).
For k > 0, with V(x) = LambertW(e^(1/k) * (1-(k+1)^2*x)/k), we have R(x) = (1-k*V(x))/(k*(k+1)) and A(x) = (1-k*V(x))/(k*(k+1)) * exp((k*V(x)-1)/(k+1)).
a(n) ~ 5^(2*n-3) * exp(n/4 - 1) * n^(n-2) / (2 * (4 + exp(5/4))^(n - 3/2)). - Vaclav Kotesovec, Jun 18 2026
PROG
(PARI) my(k=4, N=20, x='x+O('x^N), R=serreverse((1-(1-k*(k+1)*x)*exp(-(k+1)*x))/(k+1)^2)); Vec(serlaplace(R*exp(-k*R)))
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Seiichi Manyama, Jun 16 2026
STATUS
approved
