旋转卡壳，三种实现方式。

上述中说的三种是方式分别是

1. kuangbin的-容易理解
2. 挑战书上的-容易实现
3. 我的-我的理解

第一种kuangbin板子

class Point
{
public:
    double x,y;
    Point() {}
    Point(double x,double y):x(x),y(y)
    {
    }
    Point operator+ (Point p)
    {
        return Point(add(x,p.x),add(y,p.y));
    }
    Point operator -(Point p)
    {
        return Point(add(x,-p.x),add(y,-p.y));
    }
    Point operator *(double d)
    {
        return Point(x*d,y*d);
    }
    double operator *(Point p)
    {
        return add(x*p.x,y*p.y);//外积
    }
    double operator ^(Point p) //内积
    {
        return add(x*p.y,-y*p.x);
    }
    double det(Point p)
    {
        return add(x*p.y,-y*p.x);
    }
    double len()
    {
        return sqrt(add(x*x,y*y));
    }
};
}
//旋转卡壳，求两点间距离平方的最大值
int rotating_calipers(Point p[],int n)//要求p序列为凸包序列
{
    int ans = 0;
    Point v;
    int cur = 1;
    for(int i = 0; i < n; i++)
    {
        v = p[i]-p[(i+1)%n];
        while((v^(p[(cur+1)%n]-p[cur])) < 0)             cur = (cur+1)%n;
        ans = max(ans,max(dist2(p[i],p[cur]),dist2(p[(i+1)%n],p[(cur+1)%n])));
    }
    return ans;
}

第二种挑战书上的-，不过我现在还是不是非常理解。

double add(double a,double b){//考虑误差的加法运算
    if(abs(a+b) qs = convex_hull(ps,pn);//求凸包
    int n=qs.size();
    if(n==2){
        return  dist(qs[0], qs[1]);
    }
    int i=0,j=0;//某个方向的对踵点对
    //求出x轴方向的对踵点对
    for(int k=0;k

第三种 根据自己理解写出的，旋转卡壳求出每个边对应的对踵点，更新点对最大距离。

double rotaing_calipers(P* ps,int pn)// 保证输入的已经是凸包且按顺序给出点
    vector
 qs(ps,ps+pn);
    int n=qs.size();
    if(n==2){
        return  dist(qs[0], qs[1]);
    }
    int j=1;
    double ans=0.0;
    for(int i=0;i=0 ) j=(j+1)%n;//找到i-i+1边对应的对踵点j
        ans=max(ans, max(dist(qs[j],qs[i]), dist(qs[j],qs[(i+1)%n] )) );
    }
    return ans;
}