平面凸包--Jarvis

问题描述：

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

算法实现：

1.斜率逼近法，不好写而且很慢，不予介绍。
2.jarvis算法
3.Graham算法
第三个是基于对第二个的优化。
本文主要介绍Jarvis算法，Graham算法下篇文章再议。

求平面凸包的数学构造法：

找一条直线L过其中一点A，且其他点都在L的同侧（这样的直线显然存在），则A为所求凸包的一个点。
让L以A为轴点向同一个方向不断选择（默认顺时针），直至L碰到除A以外的第一个点，记为B。若同时碰到多个点B，则求距离最远的。再以B为轴点旋转，重复上述过程，直至找回到A点即求出了凸包，期间找到的点即为凸包上的点。
然而我们在计算机上并不可以实现直线的选择，因此不可以直接套用数学构造法。但是我们可以通过扫描点来达到相同的目的。

代码实现：

#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
const double eps=1e-10;
int n;
struct Point{
    double x,y;
    Point operator -(Point &s)
    {return (Point){x-s.x,y-s.y};}
};
Point d[100051];
int stk[100051];
double operator * (Point p1,Point p2)
{
    return (p1.x*p2.y-p2.x*p1.y);
} 
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));//原本的不能使这个点同侧或者不是最优
}
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;
    }
    int top=0; 
    stk[++top]=Fis; 
    do{
        int ep=1;
        for(int i=2;i<=n;i++)
         if(judgeleft(d[stk[top]],d[i],d[ep])) ep=i;
        stk[++top]=ep;
    }while(stk[top]!=Fis);
    top--;
    printf("%d\n",top);
    for(int i=1;i<=top;i++)
     printf("%.2lf %.2lf\n",d[stk[i]].x,d[stk[i]].y);
    return 0;
}