static double ptDistanceSquared(Pt a, Pt b) {
  ret sqr(a.x-b.x) + sqr(a.y-b.y);
}

static double ptDistanceSquared(double x1, double y1, double x2, double y2) {
  ret sqr(x1-x2) + sqr(y1-y2);
}