|
%I #8 Mar 27 2019 10:02:05
%S 1,-1,2,-9,68,-705,9234,-146209,2717000,-57986433,1397949830,
%T -37576332321,1114326129564,-36141571087297,1272713716466906,
%U -48360394499269665,1972269941821097744,-85929979225787811585,3983422470176606823054,-195765982110500512129057
%N Expansion of e.g.f. 2/(1 + sqrt(1 + 4*x*exp(x))).
%F E.g.f.: 1/(1 + x*exp(x)/(1 + x*exp(x)/(1 + x*exp(x)/(1 + x*exp(x)/(1 + ...))))), a continued fraction.
%F a(n) ~ sqrt(2*(1+LambertW(-1/4))) * n^(n-1) / (exp(n) * (LambertW(-1/4))^n). - _Vaclav Kotesovec_, Nov 18 2017
%p a:=series(2/(1+sqrt(1+4*x*exp(x))),x=0,20): seq(n!*coeff(a,x,n),n=0..19); # _Paolo P. Lava_, Mar 27 2019
%t nmax = 19; CoefficientList[Series[2/(1 + Sqrt[1 + 4 x Exp[x]]), {x, 0, nmax}], x] Range[0, nmax]!
%t nmax = 19; CoefficientList[Series[1/(1 + ContinuedFractionK[x Exp[x], 1, {k, 1, nmax}]), {x, 0, nmax}], x] Range[0, nmax]!
%t Table[Sum[(-1)^(n - k) Binomial[n, k] k! Sum[(-1)^m (m + 1)^(k - m - 1) Binomial[2 m, m]/(k - m)!, {m, 0, k}], {k, 0, n}], {n, 0, 19}]
%Y Cf. A000108, A006531, A052895, A295238.
%K sign
%O 0,3
%A _Ilya Gutkovskiy_, Nov 18 2017
|
