\\ 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,-1435,-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];
foreach(v, p, print(p,": ",quadhilbert(p)))

\\ other quad imaginary whose class group has exponent 2;result depend on arch
{
  v=[[-1155,x^8 + 15*x^6 + 32*x^4 + 15*x^2 + 1],
     [-1995,x^8 + 15*x^6 + 32*x^4 + 15*x^2 + 1],
     [-3003,x^8 - 9*x^6 + 80*x^4 - 9*x^2 + 1],
     [-3315,x^8 + 9*x^6 + 77*x^4 + 36*x^2 + 16]];
  foreach(v, x,
    my([D,q] = x, p = quadhilbert(D));
    if(nfisisom(p,q), print(D,": ",q)));
}

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)
