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