(PARI)permRWNb(a)=n=材料尺寸(a)[1];如果(n==1,则返回(a[1,1]);sg=1;nc=0;in=矢量v(n);x=英寸;x=a[,n]-总和(j=1,n,a[,j])/2;p=触头(i=1,n,x[i]);对于(k=1,2^(n-1)-1,sg=-sg;j=估值(k,2)+1;z=1-2*英寸[j];单位[j]+=z;nc+=z;x+=z*a[,j];p+=触头(i=1,n,x[i],sg);返回(2*(2*(n%2)-1)*p)
对于(n=1,24,a=矩阵(n,n,i,j,isprime((2*i)^2+(2*j-1)^2));print1(permRWNb(a)“,”)\\Herman Jamke(hermanjamke(AT)fastmail.fm),2007年5月13日
|