static double ptLengthSquared(Pt p) {
  ret p.x*p.x+p.y*p.y;
}