%I #30 Apr 20 2023 09:21:51
%S 1,1,2,7,34,216,1696,15898,173468,2161036,30282076,471599316,
%T 8082816160,151218316120,3066890630168,67031194526416,
%U 1570793031033616,39290173530686544,1044871388684004304,29440090627527552976
%N E.g.f. satisfies A(x) = exp(x + x^2/2 * A(x)).
%H Seiichi Manyama, <a href="/A143740/b143740.txt">Table of n, a(n) for n = 0..419</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/LambertW-Function.html">Lambert W-Function</a>.
%F E.g.f.: A(x) = -2*LambertW( -x^2*exp(x)/2 )/x^2.
%F E.g.f.: A(x) = Sum_{n>=0} (n+1)^(n-1)*(x^2/2)^n*exp((n+1)*x)/n!.
%F a(n) ~ sqrt(1+LambertW(1/sqrt(2*exp(1)))) * n^(n-1) /(2^(n+1/2) * exp(n) * (LambertW(1/sqrt(2*exp(1))))^(n+2)). - _Vaclav Kotesovec_, Jul 09 2013
%F Recurrence: a(0)=1, a(1)=1, for n > 1, a(n) = a(n-1) + Sum_{k=0..n-2} (n-1)*(n-k)*binomial(n-2,k)*a(k)*a(n-2-k)/2. - _Tani Akinari_, Nov 01 2017
%F From _Seiichi Manyama_, Apr 20 2023: (Start)
%F E.g.f.: exp(x - LambertW(-x^2/2 * exp(x))).
%F a(n) = n! * Sum_{k=0..floor(n/2)} (1/2)^k * (k+1)^(n-k-1) / (k! * (n-2*k)!). (End)
%e E.g.f.: A(x) = 1 + x + 2*x^2/2! + 7*x^3/3! + 34*x^4/4! + 216*x^5/5! + ...
%t CoefficientList[Series[-2*LambertW[-x^2*E^x/2]/x^2, {x, 0, 20}], x]* Range[0, 20]! (* _Vaclav Kotesovec_, Jul 09 2013 *)
%o (PARI) {a(n)=local(A=1+x*O(x^n));for(i=0,n,A=exp(x+x^2*A/2));(n+0)!*polcoeff(A,n)}
%o (PARI) {a(n)=local(A=sum(m=0,n,(m+1)^(m-1)*(x^2/2)^m*exp((m+1)*x+x*O(x^n))/m!)); n!*polcoeff(A,n)}
%o (PARI) {a(n)=local(Ex=exp(x+x*O(x^n)),W=Ex);for(k=0,n,W=exp(x*W)); n!*polcoeff(subst(W,x,x^2*Ex/2)*Ex,n)}
%o (Maxima) a[n]:=(if n<2 then 1 else a[n-1]+sum((n-1)*(n-k)*binomial(n-2,k)*a[k]*a[n-2-k],k,0,n-2)/2);
%o makelist(a[n],n,0,100); /* _Tani Akinari_, Nov 01 2017 */
%Y Column k=1 of A362377.
%Y Cf. A125500.
%K nonn
%O 0,3
%A _Paul D. Hanna_, Aug 30 2008