（PARI）permRWNb（a）=n=材料尺寸（a）[1]；如果（n==1，返回（a[1，1]）；sg=1；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；x+=z*a[，j]；p+=触头（i=1，n，x[i]，sg）；对于（n=1，23，a=矩阵（n，n，i，j，ndiv（abs（i-j））），返回（2*（2*（n%2）-1）*p）ndiv（n）=如果（n>0，numdiv（n），0）；print1（permRWNb（a）“，”）-Herman Jamke（hermanjamke（AT）fastmail.fm），2007年5月14日
