(帕里))克(问)={我的(V(V)=Vec公司(问),D类=设置(V(V)),d日=#D类);如果(d日==4,返回(三*如果[D类[1]]*如果[D类[2]]*如果[D类[三]]*如果[D类[4]]));
克(问) = {我的(V(V)=Vec公司(问),D类=设置(V(V)),d日=#D类);如果(d日==4,返回(三*如果[D类[1]]*如果[D类[2]]*如果[D类[三]] *如果[D类[4]]));如果(d==1,则返回((f[d[1]]^4+2*f[d[1]]^3+3*f[D1]]^2+2*f[1]])/8);
我的(k=1,,米=#=================================================================#选择(x->x====D[k],V),),t) ;而(m==1,,k++;米=#=================================================================#选择(x->x====D[千],],V(V)));));t=D[1];D[1]=D[k];D[k]=t;
序列(maxn)={my(s,a=向量(maxn));),U型);f=矢量(max_n);f[1]=1;
对于(j=1,max_n-1,如果(j个%100==0,打印(j个));f[j+1]=1/j*总和(k=1,j,sumdiv(k,d,d*f[d])*f[j-k+1]);
对于(n=4,max_n,s=0;对于部分(Q=n,秒+=克(如果( (问), [1,n个][4,[4]);一] >问[三]) && (问[n个] =秒);一三]-1>问[4..最大值_n个] }; \\ _华盛顿 邦芬_,七月 19 2012 和 12月 22 20202]),
U=U/(f[Q[4]+1]*f[Q[3]-1])*f[Q[4]*f[Q[3]],U=g(Q));s+=U,
[1,n],[4,4]);a[n]=s;如果(n%100==0,打印(n“:”s));[4..max_n]};
\\华盛顿·邦菲姆2012年7月19日和2020年12月22日
|