OFFSET
0,3
COMMENTS
Generally, A(n,k) ~ n! * n * sum(t>=1, t^k*(t^2+t-1)/(t+2)!) = n! * n * ((Bell(k) - Bell(k+1) + sum(j=0..k, (-1)^j*(2^j*((2*k-j+1)/(j+1))-1) *Bell(k-j)*C(k,j)))*exp(1) - (-1)^k*(2^k-1)), where Bell(k) are Bell numbers A000110. Set k=8 for this sequence. - Vaclav Kotesovec, Sep 12 2013
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..200
FORMULA
a(n) ~ n! * (2914*exp(1)-255)*n. - Vaclav Kotesovec, Sep 12 2013
MAPLE
a:= proc(n) option remember; `if`(n<3, [0, 1, 258][n+1],
((56*n^7-644*n^6+3332*n^5-9590*n^4+16016*n^3-14588*n^2
+5546*n+127)*a(n-1) -(n-1)*(28*n^7-280*n^6+1414*n^5
-4060*n^4+6748*n^3-5992*n^2+2017*n+254)*a(n-2) +(n-1)*(n-2)*
(28*n^6-168*n^5+490*n^4-840*n^3+868*n^2-504*n+127)*a(n-3))/
(28*n^6-336*n^5+1750*n^4-5040*n^3+8428*n^2-7728*n+3025))
end:
seq(a(n), n=0..30);
MATHEMATICA
k=8; Table[n^k+Sum[t^k*n!*(n*(t^2+t-1)-t*(t^2-4)+1)/(t+2)!+Floor[t/n]*(1/(t*(t+3)+2)), {t, 1, n-1}], {n, 0, 20}] (* Vaclav Kotesovec, Sep 12 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Sep 10 2013
STATUS
approved