HDU 1392 凸包模板题，求凸包周长

１、HDU 1392　Surround the Trees　　

２、题意：就是求凸包周长

３、总结：第一次做计算几何，没办法，还是看了大牛的博客

#include<iostream>
#include<cstring>
#include<cmath>
#include<queue>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#define F(i,a,b) for (int i=a;i<b;i++)
#define FF(i,a,b) for (int i=a;i<=b;i++)
#define mes(a,b) memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
#define LL long long
using namespace std;
const int N=,MAX=;
const double eps=1e-;

struct Point
{
    double x,y;
}p[N],stack[N];

int isZero(double x) {      //x是否为０，还是不明白为什么要引入eps
    return (x>?x:-x)<eps;
}

double crossProd(Point A,Point B,Point C) {     //叉积，这样写即A->B到A->C逆时针为正
    return (B.x-A.x)*(C.y-A.y)-(B.y-A.y)*(C.x-A.x);
}

double dis(Point A,Point B) {
    return sqrt((B.x-A.x)*(B.x-A.x)+(B.y-A.y)*(B.y-A.y));
}

//以最左下的点为基准点，其他各点（逆时针方向）以极角从小到大的排序规则
int cmp(const void *a,const void *b) {    //这里表示*a在*b后
    Point *c=(Point *)a, *d=(Point *)b;
    double k=crossProd(p[],*c,*d);     //极角大小转化为求叉乘
    if(k<eps||isZero(k)&&dis(p[],*d)<dis(p[],*c)) return ;
    return -;
}

int n;
double Graham()
{
    int mi=;
    double mx=p[].x, my=p[].y;
    FF(i,,n) {     //找到最左下点
        if(my>p[i].y||my==p[i].y&&mx>p[i].x) {
            mi=i, my=p[i].y, mx=p[i].x;
        }
    }
    p[]=p[mi], p[mi]=p[], p[]=p[];  //最左下点要换到p开头，不能让它和它自己比
    qsort(p+,n-,sizeof(Point),cmp);
    stack[]=p[], stack[]=p[], stack[]=p[];
    p[n+]=p[];    //在结尾加最左下点为结束点，一开始没加，直接在len加了一个dis(p[1],p[n])，这样忽略了很多
    int top=;
    for(int i=; i<=n+; i++) {
        //加入一个点后,向右偏拐或共线，则上一个点不在凸包内，则--top，该过程直到不向右偏拐或没有三点共线的点
        while(crossProd(stack[top-],stack[top],p[i])<eps&&top>=) top--;
        stack[++top]=p[i];
    }
    double len=;
    F(i,,top) {
        len+=dis(stack[i],stack[i+]);
    }
    return len;
}

int main()
{
    while(scanf("%d",&n)!=EOF,n) {
        FF(i,,n) {
            scanf("%lf%lf",&p[i].x,&p[i].y);
        }
        if(n==) printf("0.00\n");
        else if(n==) printf("%.2lf\n",dis(p[],p[]));
        else {
            printf("%.2lf\n",Graham());       //这里printf用G++11竟然输出0
        }
    }

    return ;
}