a[1]=1;a[n_]:=a[n]=模[{d,i,j,t,m,m},d=2*n;m=d^2-4;m=数组[0&,{m,m}];
对于[i=1,i<=d-3,i++,M[[i+1,i]]=1];
对于[i=1,i<=d-2,i++,M[[i,i+d-1]]=1];
对于[i=m-d+3,i<=m-1,i++,m[[i,i+1]]=1];
对于[i=m-d+3,i<=m,i++,m[[i-d+1,i]]=1];
对于[i=d-1,i<=m-2*d+2,i++,m[[i,i+d]]=1];
对于[i=1,i<=d-2,i++,
对于[j=1,j<=d-1,j++,t=d*i+j-2;M[[t,t+1]]=If[Mod[i,2]==1,1,-1]]];M=M-转座[M];Sqrt[检测[M]]];
|