(PARI)g(m)={my(f);如果(m==0,返回(1));f=向量(m+1);f[1]=1;
对于(j=1,m,f[j+1]=1/j*和(k=1,j,sumdiv(k,d,d*f[d])*f[j-k+1]);f[m+1]};
全局(最大值=130);A000081号=向量(max_n,n,g(n-1));
F(n,t)={my(s=0,D,c,PŤ1);forpart(P_1=n,D=Set(P_1);c=向量(#D);
对于(k=1,#D,c[k]=#选择(x->x==D[k],Vec(P#1));
s+=生产(k=1,#D,二项式(A000081号[D[k]-1]+c[k]-1,c[k]))
,[2,n],[t,t]);s};
seq(n)=和(t=1,n\2,t*F(n,t));\华盛顿博菲姆2020年7月8日
|