default(realprecision,38);
gettime;

\\ 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)
