旋转卡壳算法

double rotating_calipers(POINT *ch, int n){  
    // 旋转卡壳 求点对最长距离     效率 O(n)  
    // ch[] 凸包 为逆时针给出的  
    // 计算最长距离的时候，枚举了每一条边 |p p+1| 及 离该边最远的点q 的情况， 每一次计算点到边两端的距离，并更新一下ans（取最大值）  
    // 如果遇到下一个点与该线段围成的三角形面积(即叉积) 比这一次变大了，则比较下一个点 （即q++）  
    // 如果遇到下一个点与该线段围成的三角形面积(即叉积) 比这一次变小了，则计算两个距离（因为会出现平行的情况）并更新ans  
    // POJ 2187 暴力可过 肯定凸包的顶点数不多  
  
    int q = 1;double ans = 0;  
    ch[n] = ch[0]; // ch[q+1] 可能会访问到  
    for(int p = 0; p < n; p++){  
        // 此处遍历的线段为 |p->p+1| 点为 q   
        // 所以叉积的意义为 线段|p->p+1| 和 点q 围成的三角形的有向面积，逆时针为正，因此在这里均为正  
        while(XJ(ch[p+1], ch[q], ch[p]) > XJ(ch[p+1], ch[q+1], ch[p]))  
            q = (q + 1) % n; // 移到下一个点，因为q可能会转几圈，因此对n取余  
  
        // 这里要比较q点到线段 |p->p+1| 两个端点的距离，因为有可能在旋转的时候两个边平行  
        ans = max( ans, max(unsqrtDistance(ch[p], ch[q]), unsqrtDistance(ch[p+1], ch[q+1])) );              
    }  
    return ans;   
}  
  

  直接按照这个描述可以实现旋转卡壳算法，但是代码肯定相当冗长。逆向思考，如果qa，qb是凸包上最远两点，必然可以分别过qa，qb画出一对平行线。通过旋转这对平行线，我们可以让它和凸包上的一条边重合,如图中蓝色直线，可以注意到，qa是凸包上离p和qb所在直线最远的点。于是我们的思路就是枚举凸包上的所有边，对每一条边找出凸包上离该边最远的顶点，计算这个顶点到该边两个端点的距离，并记录最大的值。直观上这是一个O(n2)的算法，和直接枚举任意两个顶点一样了。但是注意到当我们逆时针枚举边的时候，最远点的变化也是逆时针的，这样就可以不用从头计算最远点，而可以紧接着上一次的最远点继续计算(详细的证明可以参见上面链接中的论文)。于是我们得到了O(n)的算法。

int rotating_calipers(Point *ch,int n)
{
    int q=1,ans=0;
    ch[n]=ch[0];
    for(int p=0;p<n;p++)
    {
        while(cross(ch[p+1],ch[q+1],ch[p])>cross(ch[p+1],ch[q],ch[p]))
            q=(q+1)%n;
        ans=max(ans,max(dist2(ch[p],ch[q]),dist2(ch[p+1],ch[q+1])));            
    }
    return ans; 
}

完全是搬运的：

http://www.cppblog.com/staryjy/archive/2009/11/19/101412.html

http://blog.csdn.net/x314542916/article/details/7949096