|
%I #14 Jan 01 2019 20:16:01
%S 1,2,8,46,338,2996,30952,364148,4797116,69854968,1113018176,
%T 19244304872,358608737368,7160626365296,152458303437728,
%U 3446434090192816,82412163484132112,2077739630757428768,55068742629150564736,1530394053934299827168,44490672191650220419616
%N a(n) = n! [x^n] exp((1/(x - 1)^2 - 1)/2)/(1 - x).
%F a(n + 3) = (n + 1)^2*(n + 2)*a(n) - (5 + 3*n)*(n + 2)*a(n + 1) + (8 + 3*n)*a(n + 2). - _Robert Israel_, Dec 20 2018
%F a(n) ~ exp(-1/3 + n^(1/3)/2 + 3*n^(2/3)/2 - n) * n^(n + 1/6) / sqrt(3). - _Vaclav Kotesovec_, Dec 20 2018
%p egf := exp((1/(x - 1)^2 - 1)/2)/(1 - x): ser := series(egf, x, 22):
%p seq(n!*coeff(ser, x, n), n=0..20);
%t CoefficientList[Exp[(1/(x - 1)^2 - 1)/2]/(1 - x) + O[x]^21, x] Range[0, 20]! (* _Jean-François Alcover_, Jan 01 2019 *)
%Y Row sums of A321966.
%K nonn
%O 0,2
%A _Peter Luschny_, Dec 20 2018
|
