(PARI)getT(n,k,T)=如果(!T[n+1,k+1],T[n+1,k+1]=总和(j=0,n-k-1,getT(j+k,k,T)*getT(n-j,k+1,T));T[n+1,k+1];
tabl(nn)={my(T=矩阵(nn+1,nn+1));对于(i=1,nn+1,T[i,i]=i);针对(i=0,nn,对于(j=0,i,T[i+1,j+1]=getT(i,j,T););T;}/*A127059号*/
列表(nn)={my(T=tabl(nn”));向量(nn,k,vecsum(T[k,]));}
(鼠尾草)
@缓存函数
定义T(n,k):
如果(k==n):返回n+1
else:(0..n-k-1)中j的返回和(T(j+k,k)*T(n-j,k+1))
def a(n):返回和(T(n,j)for j in(0..n))