login
E.g.f.: 2*exp(x)/(exp(2*x)+1+2*x).
1

%I #26 May 24 2016 23:57:33

%S 1,-1,3,-15,93,-725,6815,-74627,933849,-13148361,205690779,

%T -3539545559,66446203637,-1351309774685,29595401433975,

%U -694475294514315,17382734374217265,-462283425487469585,13017336622169166515,-386916316537712637215,12105656546432789499405,-397693919494074869853285

%N E.g.f.: 2*exp(x)/(exp(2*x)+1+2*x).

%C Empirically, for odd n, n|a(n) and for even n, (n-1)|a(n).

%H Robert Israel, <a href="/A272230/b272230.txt">Table of n, a(n) for n = 0..418</a>

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

%F a(n) = Sum_{k=0..n} Sum_{j=0..k} Sum_{i=0..n} binomial(n,k)*binomial(k+i,i)*binomial(n,i)(i!*(-1)^(i+j)*(2j+1)^(n-i))/(2^k)

%F a(n) ~ n! * (-1)^n*sqrt(LambertW(exp(-1)))*2^(n+1) / (1+LambertW(exp(-1)))^(n+2). - _Vaclav Kotesovec_, May 03 2016

%p S:= series(2*exp(x)/(exp(2*x)+1+2*x),x,31):

%p seq(coeff(S,x,j)*j!,j=0..30); # _Robert Israel_, May 24 2016

%t a[n_]:=Sum[Sum[Sum[Binomial[k, j]*Binomial[k + i, i]*Binomial[n,i]*((i!*(-1)^(i + j)*(2 j + 1)^(n - i))/(2^k)), {i, 0, n}], {j, 0, k}], {k, 0, n}]

%Y Cf. A000364, A122045.

%K sign

%O 0,3

%A _Christopher Ernst_, Apr 22 2016