(PARI)见Benyi/Claesson/Dukes参考第18页的公式(5)
N=40;
M=矩阵(N,N,r,c,-1);\\记忆
a(n,k)=
{
如果(n==0&&k==0,返回(1));
如果(k==0,则返回(0));
如果(n==0,则返回(0));
如果(M[n,k]!=-1,返回(M[n,k]);
本人;
s=总和(i=0,n,总和(j=0,k-1,
(-1)^j*二项式(k-j,i)*二项(i,j)*a(n-i,k-j-1));
M[n,k]=s;
申报表;
}
对于(n=0,n,print1(sum(k=1,n,a(n,k)),“,”););
\\对于(n=0,n,对于(k=0,n,打印1(a(n,k),“,”););打印(););