Graham扫描法求凸包

今天刚看了Graham扫描法求凸包，就赶紧来记录一下。以下都是我的个人理解，不保证绝对正确，若是有错误希望大佬能指出，共同进步！

什么是凸包：给了有n个点的点集，让你求出这点集中的m个点，这m个点可以将其余点都包围起来。

求法：先找出左下角的点，就是纵坐标最小的点，若纵坐标相等则取横坐标最小的点。之后将这个点(x0,y0)分别与其他点相连。之后按射线与y=x0所行成的夹角大小排序（升序）。这里要是有图应该会更好理解..不过实在是懒得弄了- -。之后我们维护一个栈stack[maxn]，首先将前两个点压入栈中，之后遍历其余的点。遍历要做的事就是判断栈顶的点是否满足顶点的要求，若满足则将遍历到的点i压入栈中，若不满足则将栈顶弹出。这么说可能有点不清楚，加个额外的描述:

for(i : 2->n){

do : while(stack[top] is not available) top--; 

then : stack[++top] = point[i];

}

这么做的原因是什么呢？要是你画了个图可能就会明白，一开始我们按夹角的大小给点排了个序，那么若点i位于凸包上，那么点j(j>i)必在点i的一侧。我突然想到加上这句可能会更好点，凸包上的点都在点集的最外侧。所以这个循环中做的事就是维护这个栈中的点始终是在"扩张"的。这也是为什么我们为每次都要将point[i]压入栈中，因为对于当前来说，它就是最后一个点，那它现在肯定满足这个条件。不过下一个点加进来后可能就不满足了，而栈中可能有多个点不满足，所以我们要循环判断。

说了那么多，也不知道讲没讲清楚。上个例题代码：（hdu1348）

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<stack>
#include<algorithm>
using namespace std;
#define PI 3.1415926
struct point
{
    double x;
    double y;
}data[1005],fun;
double len_fun;
double getLen(point a,point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(b.y-a.y)*(b.y-a.y));
}
double cha(point a,point b,point c)
{
    return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}
bool cmp(point a,point b)
{
    double ax=atan2(a.y-fun.y,a.x-fun.x);
    double bx=atan2(b.y-fun.y,b.x-fun.x);
    if(ax!=bx) return ax<bx;
    return a.x<b.x;
}
int main()
{
    int cases; scanf("%d",&cases);
    while(cases--)
    {
        int n,l; scanf("%d%d",&n,&l);
        for(int i=0;i<n;i++)
        scanf("%lf%lf",&data[i].x,&data[i].y);
        fun = data[0]; int pos=0;
        for(int i=0;i<n;i++)
        {
            if(data[i].y<fun.y || (data[i].y==fun.y&&data[i].x<fun.x))
            {
                fun=data[i];
                pos=i;
            }
        }
        data[pos]=data[0]; data[0]=fun;
        sort(data+1,data+n,cmp);
        point st[1005]; int top=1;
        st[0]=data[0]; st[1]=data[1];
        for(int i=2;i<n;i++)
        {
            while(cha(st[top-1],data[i],st[top])>0)
            top--;
            st[++top]=data[i];

        }
        double ans=0.0;
        for(int i=1;i<=top;i++)
        {
            ans+=getLen(st[i],st[i-1]);
        }
        ans+=getLen(st[0],st[top]);
        ans+=2*PI*(double)l;
        printf("%.0lf\n",ans);
        if(cases) printf("\n");
    }
    return 0;
}