|
s[1]=x；t[1]=2；z=8；
s[n]：=s[n]=s[n-1]^2-t[n-1]^2；t[n]：=t[n]=s[n-1]*t[n-1]；
系数[n_]：=应用[Riffle，Map[DeleteCase[#，0]&，系数列表[{s[n]，t[n]}，x]]]；
polys=表[Root[Total[Reverse[Ccoeffs[n]]#^（范围[1+（2^（n-1））]-1）]&，1（*2^（n-1）*）]，{n，z}]；
m=地图[NumberFieldDiscriminant，poly](*彼得·J·C·摩西2015年7月30日*）
表[m[[n+1]]/m[[n]]^2，{n，1，z-1}]（*可除性猜想*）
| 