平面凸包--Graham

题目描述：

覆盖平面上N个点的最小凸多边形
形象来讲就是平面上有n个柱子，把一根封闭的弹性绳套上这些柱子，绳子绷紧以后形成的多边形就是我们要求的凸包.

算法实现：

上一篇文章介绍了Jarvis算法
复杂度为 O（n*h）
但如果Jarvis的时间效率不够高，我们可以选择另一种更高效的算法–Graham算法.
Gramham算法其实是基于对Jarvis算法的一个优化，让我们可以更高效的找出凸包。
在Jarvis里，我们每次找寻下个点都需要 n 的去扫描
如果我们可以提前预处理出这个序列，那么就可以O（n）的找寻平面凸包。
具体流程如下：
首先找出最靠左的点（相同最靠下），然后将每个点按照到第一个点的逆时针极角排序。
然后搞一个栈，每次判断栈顶的两个元素和新加入的元素形成的角，若向右转则需弹出栈顶元素（因为是凸多边形），到最后留在栈里的元素即为凸包元素。

代码实现：

#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
const int maxm=10005;
const double eps=1e-10;
struct Point{
    double x,y;
    Point operator -(Point &s)
    {return (Point){x-s.x,y-s.y};}
};
Point d[maxm];
double operator * (Point p1,Point p2)
{
    return (p1.x*p2.y-p2.x*p1.y);
} 
int stk[maxm],top,n;
double distance(Point p1,Point p2)
{
    return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
bool judgeleft(Point p0,Point p1,Point p2)
{
    double s=(p1-p0)*(p2-p0);
    return s<0||(s==0&&distance(p1,p0)>=distance(p2,p0));
}
bool comp(Point p1,Point p2)
{
    double s=(p1-d[1])*(p2-d[1]);
    return s<0||(s==0&&distance(p1,d[1])>=distance(p2,d[1]));
}
int main()
{
    scanf("%d",&n);
    int Fis=1;
    for(int i=1;i<=n;i++)
    {
        scanf("%lf%lf",&d[i].x,&d[i].y);
        if(d[i].x<d[Fis].x||(d[i].x==d[Fis].x&&d[i].y<d[Fis].y)) Fis=i;
    }
    std::swap(d[Fis],d[1]);
    std::sort(&d[2],&d[n+1],comp);
    d[n+1]=d[1];
    stk[1]=1,stk[2]=2,top=2;
    for(int i=3;i<=n+1;i++)
    {
        while(top>1&&judgeleft(d[stk[top-1]],d[i],d[stk[top]])) top--; 
        stk[++top]=i;
    }
    top--;
    printf("%d\n",top);
    for(int i=1;i<=top;i++)
     printf("%.2lf %.2lf\n",d[stk[i]].x,d[stk[i]].y);
    //for(int i=1;i<=n;i++) printf("%.2lf %.2lf\n",d[i].x,d[i].y);
    return 0; 
}