static double distance(Vec2D a, Vec2D b) { double x = b[0] - a[0]; double y = b[1] - a[1]; return sqrt(x * x + y * y); }