凸包

算法思想:

1.选则p0作为y坐标最小的点,如果有多个这样的点,则选择x最小的。
2.剩余的点根据他们相对于p0的极角的大小从小到大排序,设排序后的点依次为P[0....n]。
3.  设置一个栈,P0,P1,P2先入栈。
4.对于 P[3..n]的每个点,若它与栈顶的两个点不构成"向左转"的关系,则将栈顶的点出栈,直至没有点需要出栈以后将当前点进栈;
所有点处理完之后栈中保存的点就是凸包了。

模板

#include <iostream>
#include <algorithm>
#include <vector>
#include <cstdio>
using namespace std;
struct Point{
    int x,y;
}p[50005];
const double eps = 1e-8;
int n;
vector<Point >s;
int dist(Point p1,Point p2){
    return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
int xmul(Point p0,Point p1,Point p2){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
//Graham_scan求凸包
bool cmp(Point p1,Point p2){
    if(xmul(p[0],p1,p2)>eps)  return true;
    else if(xmul(p[0],p2,p1)<eps&&dist(p[0],p1)<dist(p[0],p2))  return true;
    return false;
}
void Graham_scan(){
    for(int i=1;i<n;i++)
        if(p[i].y<p[0].y||(p[i].y==p[0].y&&p[i].x<p[0].x))
            swap(p[i],p[0]);
    sort(p+1,p+n,cmp);
    s.clear();
    s.push_back(p[0]);s.push_back(p[1]);
    for(int i=2;i<n;i++){
        while(s.size()>=2&&xmul(s[s.size()-2],s[s.size()-1],p[i])<eps) s.pop_back();
        s.push_back(p[i]);
    }
}


你可能感兴趣的:(凸包)