nmax:=15:与(组合):cfn1:=proc(n,k):和(-1)^j*stirling1(n+1,n+1-k+j)*stirling 1(n+1,n+1-k-j),j=-k.k)结束过程:Omega(0):=1:对于从1到nmax的n,做Omega 1),k1=1..n))/(2*n-1)!end do:对于n从1到nmax do d(n):=2^(2*n-1)*Omega(n)end do:对于n从2到nmmax do Zc(n-1):=d(n-1 1):=c(n)end do:seq(p(n),n=2..nmax);