凸包模板

#include<iostream>
#include<cmath>
#include<stack>
using namespace std;
#define Max 1000
int multi(int p1,int p2,int p3);
int Point[Max][2] , a[Max][2] , n;
int multi(int p1,int p2,int p3)
{
    int v;
    int v1x=Point[p2][0]-Point[p1][0],
    v1y=Point[p2][1]-Point[p1][1],
    v2x=Point[p3][0]-Point[p1][0],
    v2y=Point[p3][1]-Point[p1][1];
    if(p1==p2) return p3;
    if(p1==p3) return p2;
    else{
        v=v1x*v2y-v1y*v2x;
        if(v>0) return p2;
        if(v<0) return p3;
        if(v==0){
        if((v1x*v1x+v1y*v1y)>(v2x*v2x+v2y*v2y))
        return p2;
        else return p3;
   }

}
}
int qiutubao(int point[Max][2],int &n,int a[Max][2])//point是输入的点集 求出的凸包存到a中
{                                                    //返回 的值表示凸包上顶点个数
int i;
    int m=Point[0][1],t=0;
    for(int k=1;k<n;k++){
        if (Point[k][1]<m)
        {
            m=Point[k][1];t=k;
        }
        else if(Point[k][1]==m)
        {
            if (Point[k][0]<Point[t][0])
            {m=Point[k][1];t=k;}
        }
}
    Point[n][0] = Point[t][0];
    Point[n][1] = Point[t][1];
    int q=n;
    stack<int> p;
    p.push(t);
    int kk=t,g=0;
    while(true)
    {
        for(int i=1;i<=n;i++)
            g= multi(kk,g,i);
        if(g==q)   break;
        p.push(g); kk=g; g=0;
}
    int get; i=0;
    while(!p.empty())     //出栈并存到a中
    {        
        get = p.top();
   a[i][0]=Point[get][0];a[i][1]=Point[get][1];
   i++;   p.pop();
    }
    return i;
}
int main()
{
    int i=0,l;
    scanf("%d",&n);
    for(i=0;i<n;i++)
        scanf("%d%d",&Point[i][0],&Point[i][1]);
l=qiutubao(Point,n,a);
cout<<endl;
for(i=0;i<l;i++)
   cout<<a[i][0]<<' '<<a[i][1]<<endl;
return 0;
}