default(realprecision,38);
default(echo,1);
gettime;
test(D) =
{
  for (d = D, D+10^3,
    if (!isfundamental(d),next);
    print(d, " ", quadclassunit(d).cyc)
  );
}
setrand(1); test(10^15); test(-10^15)
