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 A143990 #14 Sep 28 2023 04:24:25 %S A143990 1,1,4,42,888,31920,1750320,136115280,14254007040,1934091250560, %T A143990 330078373228800,69199130042380800,17481751115946163200, %U A143990 5237838647954514201600,1836425205487182172262400,744852154338379227748608000,346052141662324885396697088000,182572078442025253754006986752000 %N A143990 a(n) = n!*A001515(n-1) with a(0) = 1. %H A143990 N. J. A. Sloane, <a href="/A143990/b143990.txt">Table of n, a(n) for n = 0..101</a> %F A143990 a(n) = n*(2*n-3)*a(n-1) + n*(n-1)*a(n-2), a(0)=1, a(1)=1. - _Sergei N. Gladkovskii_, May 17 2013 %F A143990 From _G. C. Greubel_, Sep 27 2023: (Start) %F A143990 a(n) = n * A105749(n-1) + [n=0]. %F A143990 a(n) = n! * A144301(n). %F A143990 E.g.f.: 1 + sqrt(Pi*x/2) * exp(-(1-x)^2/(2*x)) * erfi((1-x)/sqrt(2*x)). %F A143990 Sum_{n >= 0} a(n)*x^n/(n!)^2 = exp(1 - sqrt(1-2*x)). (End) %t A143990 With[{m=30}, CoefficientList[Series[Exp[1-Sqrt[1-2*x]], {x,0,m}], x]*(Range[0, m]!)^2] (* _G. C. Greubel_, Sep 27 2023 *) %o A143990 (Magma) [n le 2 select 1 else (n-1)*(2*n-5)*Self(n-1) + (n-1)*(n-2)*Self(n-2): n in [1..30]]; // _G. C. Greubel_, Sep 27 2023 %o A143990 (SageMath) %o A143990 m=30 %o A143990 P.<x> = PowerSeriesRing(QQ, m+2) %o A143990 def A143990(n): return (factorial(n))^2*P( exp(1-sqrt(1-2*x)) ).list()[n] %o A143990 [A143990(n) for n in range(m+1)] # _G. C. Greubel_, Sep 27 2023 %Y A143990 Cf. A001515, A105749, A144301. %K A143990 nonn %O A143990 0,3 %A A143990 _N. J. A. Sloane_, Dec 02 2008