凸包模板

#include<bits/stdc++.h>
using namespace std;
const int N=2005;
int n,tot;
struct point
{
    int x,y;
    point(int _x=0,int _y=0){x=_x; y=_y;}
    point operator -(const point &rhs)const
    {
        return point(x-rhs.x,y-rhs.y);
    }
}p[N],cp[N];
typedef point vec;
int cross(const vec &a,const vec &b)
{
    return (a.x*b.y)-(a.y*b.x);
}
int dis(point a,point b)
{
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
bool cmp(const point &a,const point &b)
{
    int z=cross(a-p[0],b-p[0]);
    if(z>0||(z==0 && dis(p[0],a)<dis(p[0],b)))
        return 1;
    else return 0;
}
void get_cp()
{
    int k=0;
    for(int i=0;i<n;i++)
        if(p[i].y<p[k].y || p[i].y==p[k].y && p[i].x<p[k].x) k=i;
    swap(p[0],p[k]);
    sort(p+1,p+n,cmp);
    tot=2,cp[0]=p[0];cp[1]=p[1];
    for(int i=2;i<n;i++)
    {
        while(tot>1 && cross(cp[tot-2]-cp[tot-1],p[i]-cp[tot-1])>0) tot--;
        cp[tot++]=p[i];
    }
}
int main()
{
    scanf("%d",&n);
    for(int i=0;i<n;i++) scanf("%d%d",&p[i].x,&p[i].y);
    get_cp();
    for(int i=0;i<tot;i++) printf("%d %d
",cp[i].x,cp[i].y);
    return 0;
}