This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A321965 #14 Jan 01 2019 20:16:01 %S A321965 1,2,8,46,338,2996,30952,364148,4797116,69854968,1113018176, %T A321965 19244304872,358608737368,7160626365296,152458303437728, %U A321965 3446434090192816,82412163484132112,2077739630757428768,55068742629150564736,1530394053934299827168,44490672191650220419616 %N A321965 a(n) = n! [x^n] exp((1/(x - 1)^2 - 1)/2)/(1 - x). %F A321965 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 A321965 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 A321965 egf := exp((1/(x - 1)^2 - 1)/2)/(1 - x): ser := series(egf, x, 22): %p A321965 seq(n!*coeff(ser, x, n), n=0..20); %t A321965 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 A321965 Row sums of A321966. %K A321965 nonn %O A321965 0,2 %A A321965 _Peter Luschny_, Dec 20 2018