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]*200e=[0]*200d[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