HDU 1392 Surround the Trees（求凸包周长）

题意：有n棵树，每棵树给出它的坐标，现在以一些树为端点构造围墙将所有树都围起来，求出围墙的总长度。

思路：凸包问题，将树看成点求出凸包，然后求出凸包的周长即可。

注意：当n=1和n=2时要特判，n=1时为0，n=2时为两点间的距离。(吐槽一下为毛我的C++不能编译G++就能AC)

#include<stdio.h>
#include<math.h>
#include<algorithm>
#include<iostream>
using namespace std;
const int MAXN=1000;

struct point
{
    int x,y;
};
point list[MAXN];
int stack[MAXN],top;

int cross(point p0,point p1,point p2) //计算叉积  p0p1 X p0p2
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
}
double dis(point p1,point p2)  //计算 p1p2的 距离
{
    return sqrt((double)(p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y));
}
bool cmp(point p1,point p2) //极角排序函数 ， 角度相同则距离小的在前面
{
    int tmp=cross(list[0],p1,p2);
    if(tmp>0) return true;
    else if(tmp==0&&dis(list[0],p1)<dis(list[0],p2)) return true;
    else return false;
}
void init(int n) //输入，并把  最左下方的点放在 list[0]  。并且进行极角排序
{
    int i,k;
    point p0;
    scanf("%d%d",&list[0].x,&list[0].y);
    p0.x=list[0].x;
    p0.y=list[0].y;
    k=0;
    for(i=1;i<n;i++)
    {
        scanf("%d%d",&list[i].x,&list[i].y);
        if( (p0.y>list[i].y) || ((p0.y==list[i].y)&&(p0.x>list[i].x)) )
        {
            p0.x=list[i].x;
            p0.y=list[i].y;
            k=i;
        }
    }
    list[k]=list[0];
    list[0]=p0;

    sort(list+1,list+n,cmp);
}

void graham(int n)
{
    int i;
    if(n==1)
    {
        top=0;
        stack[0]=0;
    }
    if(n==2)
    {
        top=1;
        stack[0]=0;
        stack[1]=1;
    }
    if(n>2)
    {
        for(i=0;i<=1;i++) stack[i]=i;
        top=1;

        for(i=2;i<n;i++)
        {
            while(top>0&&cross(list[stack[top-1]],list[stack[top]],list[i])<=0)
                top--;
            top++;
            stack[top]=i;
        }
    }
}

int main()
{
    int n;
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    while(scanf("%d",&n)&&n!=0)
    {
        if(n==1)
            printf("0.00\n");
        else if(n==2)
        {
            scanf("%d%d%d%d",&list[0].x,&list[0].y,&list[1].x,&list[1].y);
            printf("%.2lf\n",dis(list[0],list[1]));
        }
        else
        {
            double ans=0;
            init(n);
            graham(n);
            //求凸包的周长
            for(int i=0;i<=top-1;i++)//注意凸包中的点是从0保存到top的
            {
                int xx1=list[stack[i]].x,yy1=list[stack[i]].y;
                int xx2=list[stack[i+1]].x,yy2=list[stack[i+1]].y;
                ans+=sqrt((yy1-yy2)*(yy1-yy2)+(xx1-xx2)*(xx1-xx2));
            }
            int xx1=list[0].x,yy1=list[0].y;
            int xx2=list[stack[top]].x,yy2=list[stack[top]].y;
            ans+=sqrt((yy1-yy2)*(yy1-yy2)+(xx1-xx2)*(xx1-xx2));
            printf("%.2lf\n",ans);
        }
    }
    return 0;
}