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 A356559 #8 Feb 16 2025 08:34:03 %S A356559 1,0,0,1,7,43,281,2056,17004,157809,1622515,18245335,222004597, %T A356559 2898508416,40343356184,595578837205,9287308741827,152459628788599, %U A356559 2627373030049669,47425289731038656,895098852673047772,17644305594671247141,363065584549610882703,7799894520723959486795 %N A356559 a(n) = exp(-1) * n! * Sum_{k>=0} Laguerre(n,k) / k!. %H A356559 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/LaguerrePolynomial.html">Laguerre Polynomial</a> %F A356559 E.g.f.: exp(exp(-x/(1 - x)) - 1) / (1 - x). %F A356559 a(n) = Sum_{k=0..n} (-1)^(n-k) * binomial(n,k)^2 * k! * Bell(n-k). %t A356559 Table[Exp[-1] n! Sum[LaguerreL[n, k]/k!, {k, 0, Infinity}], {n, 0, 23}] %t A356559 nmax = 23; CoefficientList[Series[Exp[Exp[-x/(1 - x)] - 1]/(1 - x), {x, 0, nmax}], x] Range[0, nmax]! %t A356559 Table[Sum[(-1)^(n - k) Binomial[n, k]^2 k! BellB[n - k], {k, 0, n}], {n, 0, 23}] %o A356559 (PARI) my(x='x+O('x^30)); Vec(serlaplace(exp(exp(-x/(1 - x)) - 1) / (1 - x))) \\ _Michel Marcus_, Aug 12 2022 %Y A356559 Cf. A000110, A009940, A101053, A317362, A317366. %K A356559 nonn %O A356559 0,5 %A A356559 _Ilya Gutkovskiy_, Aug 12 2022