OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..430
FORMULA
a(n) = Sum_{k=0..n} binomial(n-1,k-1)*A000217(k)*n!/k!.
a(n) ~ n^(n + 3/4) / (2^(3/2) * exp(n - 2*sqrt(n) + 1/2)). - Vaclav Kotesovec, Jan 27 2019
a(n) = n!*Hypergeometric2F2([1-n, 3], [1, 2], -1). - G. C. Greubel, Mar 05 2021
MAPLE
seq(n!*coeff(series(x*(2-x)*exp(x/(1-x))/(2*(1-x)^2), x=0, 22), x, n), n=0..21); # Paolo P. Lava, Jan 29 2019
MATHEMATICA
nmax = 21; CoefficientList[Series[x (2 - x) Exp[x/(1 - x)]/(2 (1 - x)^2), {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[Binomial[n - 1, k - 1] Binomial[k + 1, 2] n!/k!, {k, 0, n}], {n, 0, 21}]
PROG
(Sage) [0]+[sum(binomial(n-1, k)*binomial(k+2, 2)*factorial(n)/factorial(k+1) for k in (0..n-1)) for n in [1..20]] # G. C. Greubel, Mar 05 2021
(Magma) [n eq 0 select 0 else (&+[Binomial(n-1, k)*Binomial(k+2, 2)* Factorial(n)/Factorial(k+1): k in [0..n-1]]): n in [0..20]]; // G. C. Greubel, Mar 05 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 27 2019
STATUS
approved