OFFSET
0,3
COMMENTS
In general, for m >= 0, Sum_{k=0..n} k^m * (n-k)! * binomial(n,k)^2 ~ exp(2*sqrt(n) - n - 1/2) * n^(n + (2*m+1)/4) / sqrt(2). - Vaclav Kotesovec, Feb 20 2021
FORMULA
E.g.f.: x*(1+4*x-5*x^2+x^3)*exp(x/(1-x))/(1-x)^5.
a(n) = n*n!*F([2, 2, 1 - n], [1, 1, 1], -1), where F is the generalized hypergeometric function. - Stefano Spezia, Feb 06 2021
a(n) ~ exp(2*sqrt(n) - n - 1/2) * n^(n + 9/4) / sqrt(2). - Vaclav Kotesovec, Feb 20 2021
MATHEMATICA
a[n_] := Sum[k^4 * (n-k)! * Binomial[n, k]^2, {k, 0, n}]; Array[a, 19, 0] (* Amiram Eldar, Feb 06 2021 *)
PROG
(PARI) a(n) = sum(k=0, n, k^4*(n-k)!*binomial(n, k)^2);
(PARI) my(N=20, x='x+O('x^N)); concat(0, Vec(serlaplace(x*(1+4*x-5*x^2+x^3)*exp(x/(1-x))/(1-x)^5)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Feb 06 2021
STATUS
approved