login
E.g.f. satisfies A(x) = x * exp( 2*(x + A(x)) ).
3

%I #17 Mar 12 2023 10:52:44

%S 0,1,8,96,1792,46080,1511424,60325888,2837970944,153778913280,

%T 9432255692800,646039266656256,48874810528235520,4047655951598092288,

%U 364221261622538141696,35384754572803304325120,3691411033400626898796544,411569264258973944034361344

%N E.g.f. satisfies A(x) = x * exp( 2*(x + A(x)) ).

%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) = (-1/2) * LambertW(-2*x * exp(2*x)).

%F a(n) = Sum_{k=1..n} (2*k)^(n-1) * binomial(n,k) = 4^(n-1) * A100526(n).

%F a(n) ~ sqrt(1 + LambertW(exp(-1))) * 2^(n-1) * n^(n-1) / (LambertW(exp(-1))^n * exp(n)). - _Vaclav Kotesovec_, Feb 17 2023

%p A360548 := proc(n)

%p add((2*k)^(n-1)*binomial(n,k),k=1..n) ;

%p end proc:

%p seq(A360548(n),n=0..60) ; # _R. J. Mathar_, Mar 12 2023

%o (PARI) my(N=20, x='x+O('x^N)); concat(0, Vec(serlaplace(-lambertw(-2*x*exp(2*x))/2)))

%o (PARI) a(n) = sum(k=1, n, (2*k)^(n-1)*binomial(n, k));

%Y Cf. A100526, A216857, A360545.

%Y Cf. A038049, A201595, A214225, A360547.

%K nonn

%O 0,3

%A _Seiichi Manyama_, Feb 11 2023