旋转卡壳算法

• 点集的直径：指的是点集中最远点对的距离
• 对踵点对：过该点对的两条平行线能把所有凸包夹在之间
  
  D,G就是点对
• 旋转卡壳：
  射线j,k同时逆时针旋转θ角度，绕凸多边形一圈，他们两线之间的最大距离就是点集的直径
  刘汝佳的算法如下：

// 返回点集直径的平方
int diameter2(vector<Point>& points) {
  vector<Point> p = ConvexHull(points);
  int n = p.size();
  if(n == 1) return 0;
  if(n == 2) return Dist2(p[0], p[1]);
  p.push_back(p[0]); // 免得取模
  int ans = 0;
  for(int u = 0, v = 1; u < n; u++) {
    // 一条直线贴住边p[u]-p[u+1]
    for(;;) {
      // 当Area(p[u], p[u+1], p[v+1]) <= Area(p[u], p[u+1], p[v])时停止旋转
      // 即Cross(p[u+1]-p[u], p[v+1]-p[u]) - Cross(p[u+1]-p[u], p[v]-p[u]) <= 0
      // 根据Cross(A,B) - Cross(A,C) = Cross(A,B-C)
      // 化简得Cross(p[u+1]-p[u], p[v+1]-p[v]) <= 0
      int diff = Cross(p[u+1]-p[u], p[v+1]-p[v]);
      if(diff <= 0) {
        ans = max(ans, Dist2(p[u], p[v])); // u和v是对踵点
        if(diff == 0) ans = max(ans, Dist2(p[u], p[v+1])); // diff == 0时u和v+1也是对踵点
        break;
      }
      v = (v + 1) % n;
    }
  }
  return ans;
}

也即是，以凸包的一条边为基边，如上面的L1=p[u]-p[u+1]这条边，然后一直旋转L2=p[v]-p[v+1]这条直线到下一凸包的边

旋转的条件就是当p[v+1]到直线L1的距离比p[v]到直线L1的大时，如果小的话,就更新点p[v]到点p[u]的值,如果相等的话，说明L1与L2平行，答案也有可能是p[v+1]到p[u]的值

判断点到线的距离可以通过判断面积大小，所以可以转化为叉积比较，就是上面代码写得那样