with(组合):M:=多项式:
b: =proc(n,k)局部l,g;l、 克:=[1,5],
proc(k,m,i,t)选项记忆;局部d,j;d: =l[i];
`如果`(i=1,n^m,加上(m(k,k-(d-t)*j,(d-t,$j)/j*
(d-1)^j*M(M,M-t*j,t$j)*g(k-(d-t)*j,M-t**j,
`如果`(d-t=1,[i-1,0],[i,t+1])[]),j=0..分钟(k/(d-t),
`如果`(t=0,[][],m/t)))
结束;g(k,n-k,nops(l),0)
结束:
a: =n->加上(b(n,j)*stirling2(n,j*二项式(n,j)*j!,j=0..n):
seq(a(n),n=0..20);
|