cp's OEIS Frontend

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.

A308946 Expansion of e.g.f. 1/(1 - x*(1 + x/2)*exp(x)).

This page as a plain text file.
%I A308946 #7 Aug 09 2021 03:21:26
%S A308946 1,1,5,30,244,2485,30351,432502,7043660,129050649,2627117875,
%T A308946 58829021416,1437117395946,38032508860177,1083932872119839,
%U A308946 33098858988564090,1078083456543449416,37309607437056658129,1367138649165397662627,52879280631976735387588
%N A308946 Expansion of e.g.f. 1/(1 - x*(1 + x/2)*exp(x)).
%F A308946 E.g.f.: 1 / (1 - Sum_{k>=1} (k*(k + 1)/2)*x^k/k!).
%F A308946 a(0) = 1; a(n) = Sum_{k=1..n} binomial(n,k) * A000217(k) * a(n-k).
%F A308946 a(n) ~ n! * (2 + r) / ((2 + 4*r + r^2) * r^n), where r = 0.49122518354447387971550543450091640839121607... is the root of the equation  exp(r)*r*(2 + r) = 2. - _Vaclav Kotesovec_, Aug 09 2021
%t A308946 nmax = 19; CoefficientList[Series[1/(1 - x (1 + x/2) Exp[x]), {x, 0, nmax}], x] Range[0, nmax]!
%t A308946 a[0] = 1; a[n_] := a[n] = Sum[Binomial[n, k] Binomial[k + 1, 2] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 19}]
%Y A308946 Cf. A000217, A001793, A006153, A052529, A279361, A308861.
%K A308946 nonn
%O A308946 0,3
%A A308946 _Ilya Gutkovskiy_, Jul 02 2019