OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..200
MAPLE
C:= binomial:
b:= proc(n, i, k, p) option remember; `if`(n=0, p!, `if`(i<1, 0, add(
b(n-i*j, min(n-i*j, i-1), k, p+j)/j!*C(C(k+i-1, i), j), j=0..n/i)))
end:
a:= n-> add(add(b(n$2, i, 0)*(-1)^(k-i)*C(k, i), i=0..k), k=0..n):
seq(a(n), n=0..23);
MATHEMATICA
c = Binomial;
b[n_, i_, k_, p_] := b[n, i, k, p] = If[n == 0, p!, If[i<1, 0, Sum[b[n - i*j, Min[n - i*j, i-1], k, p+j]/j!*c[c[k+i-1, i], j], {j, 0, n/i}]]];
a[n_] := Sum[Sum[b[n, n, i, 0]*(-1)^(k-i)*c[k, i], {i, 0, k}], {k, 0, n}];
Table[a[n], {n, 0, 23}] (* Jean-François Alcover, Mar 05 2022, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Sep 18 2019
STATUS
approved