凸包(安德鲁算法)

凸包是指包含点集合P中所有点的最小凸多边形。

凸包可以理解为在木板上钉了许多钉子，然后用一根橡皮筋框住所有钉子时所得到的多边形。

安德鲁算法：

1.将给定集合的点按x坐标升序排列，x相同的按y坐标升序排列。

2.按下述流程创建凸包的上部

将排序后的点按照x坐标从小到大的顺序加入凸包U。如果新加入的点使得U不再是凸多边形，则逆序删除之前加入U的点，直到U重新成为凸多边形为止。

3.按下述流程创建凸包的下部

将排序后的点按照x坐标从大到小的顺序加入凸包L。如果新加入的点使得L不再是凸多边形，则逆序删除之前加入L的点，直到L重新成为凸多边形为止。

 

具体步骤：

第一步：将最开始的2个点加入凸包，把加入当前凸包的操作称为push_back，从第二步开始，按照下标顺序将S中的点依次push_back到凸包U中。不过，在将Si加入凸包之前，要先进行下面循环处理

第二步:对于当前的凸包U，只要点Si（即将加入U的点）位于倒数第二个点和倒数第一个点所构成的向量的逆时针位置，就从凸包U中删除倒数第一个点。将删除点的操作称为pop_back。

struct point{
    double x,y;
    point(){}
    point(double x,double y):x(x),y(y){}
    point operator + (point &p){
        return point(x+p.x,y+p.y);
    }
    point operator - (point &p){
        return point(x-p.x,y-p.y);
    }
    point operator * (double k){
        return point(x*k,y*k);
    }
    point operator / (double k){
        return point(x/k,y/k);
    }
};

typedef vector<point>polygon;

int cmp(point a,point b)
{
    if(a.x!=b.x){
        return a.x<b.x;
    }else{
        return a.y<b.y;
    }
}

double cross(point a,point b)
{
    return a.x*b.y-a.y*b.x;
}

#define eps (1e-10)
#define CLOCKWISE 1
#define COUNTER_CLOCKWISE -1
int ccw(point a,point b,point c)
{
    point base=b-a;
    point linshi=c-a;
    if(cross(base,linshi)<eps){
        return CLOCKWISE;
    }else{
        return COUNTER_CLOCKWISE;
    }
}

polygon andrewscanf(polygon s)
{
    int n=s.size();
    if(n<3){
        return s;
    }
    sort(s.begin(),s.end(),cmp);
    polygon u,l;
    u.push_back(s[0]);
    u.push_back(s[1]);

    l.push_back(s[n-1]);
    l.push_back(s[n-2]);
    for(int i=2;i<n;i++){
        for(int j=u.size();j>=2&&(ccw(u[j-2],u[j-1],s[i])!=CLOCKWISE);j--){
            u.pop_back();
        }
        u.push_back(s[i]);
    }
    for(int i=n-3;i>=0;i--){
        for(int j=l.size();j>=2&&(ccw(l[j-2],l[j-1],s[i])!=CLOCKWISE);j--){
           // cout<<s[i].x<<" "<<s[i].y<<" "<<l[j-2].x<<" "<<l[j-2].y<<" "<<l[j-1].x<<" "<<l[j-1].y<<endl;
            l.pop_back();
        }
        l.push_back(s[i]);
    }
    reverse(l.begin(),l.end());
    for(int i=u.size()-2;i>=1;i--){
        l.push_back(u[i]);
    }
    return l;
}