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 A367818 #10 Dec 02 2023 07:49:02 %S A367818 1,-4,15,-53,178,-575,1809,-5598,17141,-52113,157724,-475997,1433429, %T A367818 -4311364,12958627,-38909601,116831426,-350844883,1051414421, %U A367818 -3160120038,9491592177,-28218244109,86403627444,-255153772169,722619907385,-2772952748516,4627276967623,-17420488524253 %N A367818 Expansion of e.g.f. exp(1 - 3*x - exp(x)). %F A367818 G.f. A(x) satisfies: A(x) = 1 - 3*x*A(x) - x * A(x/(1 - x)) / (1 - x). %F A367818 a(n) = exp(1) * Sum_{k>=0} (-1)^k * (k-3)^n / k!. %F A367818 a(0) = 1; a(n) = -3*a(n-1) - Sum_{k=1..n} binomial(n-1,k-1) * a(n-k). %F A367818 a(n) = Sum_{k=0..n} binomial(n,k) * (-3)^(n-k) * A000587(k). %t A367818 nmax = 27; CoefficientList[Series[Exp[1 - 3 x - Exp[x]], {x, 0, nmax}], x] Range[0, nmax]! %t A367818 a[0] = 1; a[n_] := a[n] = -3 a[n - 1] - Sum[Binomial[n - 1, k - 1] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 27}] %t A367818 Table[Sum[Binomial[n, k] (-3)^(n - k) BellB[k, -1], {k, 0, n}], {n, 0, 27}] %o A367818 (PARI) my(x='x+O('x^30)); Vec(serlaplace(exp(1 - 3*x - exp(x)))) \\ _Michel Marcus_, Dec 02 2023 %Y A367818 Cf. A000587, A005494, A109747, A153732, A193683, A346738, A367819. %K A367818 sign %O A367818 0,2 %A A367818 _Ilya Gutkovskiy_, Dec 01 2023