OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..200
FORMULA
a(n) = n! * [x^n] (exp(x)-x^3/3!)^n.
a(n) ~ c * d^n * n^n / exp(n), where d = 2.52566039645910026750819504865..., c = 1.031458655073968039932844239... . - Vaclav Kotesovec, Jul 24 2014
MAPLE
b:= proc(n, i) option remember; `if`(n=0 and i=0, 1, `if`(i<1, 0,
add(`if`(j=3, 0, b(n-j, i-1) *binomial(n, j)), j=0..n)))
end:
a:= n-> b(n$2):
seq(a(n), n=0..25);
MATHEMATICA
Table[n!*SeriesCoefficient[(E^x-x^3/6)^n, {x, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Jul 23 2014 *)
With[{k=3}, Flatten[{1, Table[Sum[Binomial[n, j]*Binomial[n, k*j]*(-1)^j*(n-j)^(n-k*j)*(k*j)!/(k!)^j, {j, 0, n/k}], {n, 1, 20}]}]] (* Vaclav Kotesovec, Jul 24 2014 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jul 21 2014
STATUS
approved