Jackson2:=进程(n)局部m,i,j,m,p,b,s,x;
如果0=(n mod 2),则;
m:=无/2;
M:=矩阵(M,M,0);
对于i从1到m do,对于j从1到mdo;
如果1=igcd(2*i,2*j-1),则M[i,j]:=1;fi;od;od;
s:=线性代数[永久](M);
返回s^2;
其他;
m:=(n+1)/2;
M:=矩阵(M,M,0);
对于i从1到m-1 do,对于j从1到m do;
如果1=igcd(2*i,2*j-1),则M[i,j]:=1;fi;od;od;
对于j到m do
M[M,j]:=x[j];
结束do;
p:=线性代数[永久](M);
b:=[];
对于j到m do
b:=[op(b),系数(p,x[j])];
结束do;
s:=0;
对于i从1到m do,对于j从1到mdo;
如果1=igcd(2*i-1,2*j-1),则s:=s+b[i]*b[j];fi;od;od;fi;
返回s;
结束;
seq(杰克森2(n),n=1..25); #斯蒂芬·洛克2022年2月24日