D = -5291
T = [BinaryQF([1,1,1323])]
a = [BinaryQF([3,1,441]),BinaryQF([5,3,265]),BinaryQF([7,1,189])]
r = [0,0,0]
s = [0,0,0]
for i in range(len(a)):
r[i],s[i]=extend_table(T,a[i])
assert len(T) == QuadraticField(D).class_number()
print "a =",a
print "r =",r
print "s =",s
a = [3*x^2 + x*y + 441*y^2, 5*x^2 + 3*x*y + 265*y^2, 7*x^2 + x*y + 189*y^2]
r = [9, 2, 2]
s = [0, 4, 6]