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日
