weixin_47616793的博客

def f(p,q) : K = NumberField(x^3 +p*x+q, ‘a’); return K.integral_basis() [print(p,q,f(p,q)) for p in range(1,7,2) for q in range(1,7,2)] 1 1 [1, a, a^2] 1 3 [1, a, a^2] 1 5 [1, a, a^2] 3 1 [1, a, a^2] 3 3 [1, a, a^2] 35[1/3∗a2+1/3∗a+1/3,a,a2]3 5 [1/3*a^2 +

def f®: d=[0]*200 e=[0]*200 d[0]=ceil®; e[0]=1/(d[0]-r); for i in range(1,8): d[i]=ceil(e[i-1]); e[i]=1/(d[i]-e[i-1]); M=MatrixSpace(ZZ,8,1) d=M(d) return transpose(d) f(N(sqrt(163))) 二次域得到类数与上面的连分数有关，周期为n1，。。。，nr时，其和/3-r为二次域的类数，在p为模4余3时成立 def f§: K.=Quadr