default(realprecision,38);
\\ all quad imaginary whose class group has exponent 2
v=[-15,-35,-51,-91,-115,-123,-187,-195,-235,-267,-403,-427,-435,-483,-555,-595,-627,-715,-795,-1155,-1435,-1995,-3003,-3315,-20,-24,-40,-52,-84,-88,-120,-132,-148,-168,-228,-232,-280,-312,-340,-372,-408,-420,-520,-532,-660,-708,-760,-840,-1012,-1092,-1320,-1380,-1428,-1540,-1848,-5460];
for(i=1, #v, print(v[i],": ",quadhilbert(v[i])))

quadhilbert(-4036)
quadhilbert(-300003)
quadhilbert(-3628843)

Q(D,f) = lift(quadray(D,f));
Q(-4,31)
Q(-11,2)
Q(-15,3)
Q(-179,2)
Q(-2276,2)
Q(-251,2)
Q(-35,2)
Q(-4,31)
Q(-51,3)
Q(8-48*3,2)
Q(1-48*3,3)
Q(1-48*3,4)
Q(40-48*3,6)
Q(-7,7)

K = bnfinit(y^2+5);
P5=idealprimedec(K,5)[1];
Q(K,P5)

K = bnfinit(y^2+5*12);
P2=idealprimedec(K,2)[1];
P5=idealprimedec(K,5)[1];
Q(K,idealmul(K,P2,P5))

\\#1633
quadray(-11,3)
