凸包：HDU1392(ZJU1453)

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<iostream>
#include<math.h>
#include<algorithm>
using namespace std;
/* 
PointSet[]：输入的点集 
ch[]：输出的凸包上的点集，按照逆时针方向排列 
n：PointSet中的点的数目 
len：输出的凸包上的点的个数 
*/  
struct Point  
{  
    int  x,y;  
};  
//小于0,说明向量p0p1的极角大于p0p2的极角  
double multiply(Point p1,Point p2,Point p0)  
{  
    return((p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y));  
}  
  
double dis(Point p1,Point p2)  
{  
    return(sqrt((double)(p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y)));  
}  
  
  
void Graham_scan(Point PointSet[],Point ch[],int n,int &len)  
{  
    int i,j,k=0,top=2;  
    Point tmp;  
  
    //找到最下且偏左的那个点  
    for(i=1;i<n;i++)  
        if ((PointSet[i].y<PointSet[k].y)||((PointSet[i].y==PointSet[k].y)&&(PointSet[i].x<PointSet[k].x)))  
            k=i;  
    //将这个点指定为PointSet[0]  
    tmp=PointSet[0];  
    PointSet[0]=PointSet[k];  
    PointSet[k]=tmp;  
  
    //按极角从小到大,距离偏短进行排序  
    for (i=1;i<n-1;i++)  
    {  
        k=i;  
        for (j=i+1;j<n;j++)  
            if( (multiply(PointSet[j],PointSet[k],PointSet[0])>0)  
                ||((multiply(PointSet[j],PointSet[k],PointSet[0])==0)  
                    &&(dis(PointSet[0],PointSet[j])<dis(PointSet[0],PointSet[k]))) )  
                k=j;//k保存极角最小的那个点,或者相同距离原点最近  
        tmp=PointSet[i];  
        PointSet[i]=PointSet[k];  
        PointSet[k]=tmp;  
    }  
    //第三个点先入栈  
    ch[0]=PointSet[0];  
    ch[1]=PointSet[1];  
    ch[2]=PointSet[2];  
    //判断与其余所有点的关系  
    for (i=3;i<n;i++)  
    {  
        //不满足向左转的关系,栈顶元素出栈  
        while(multiply(PointSet[i],ch[top],ch[top-1])>=0) top--;  
        //当前点与栈内所有点满足向左关系,因此入栈.  
        ch[++top]=PointSet[i];  
    }  
    len=top+1;  
}  
  
const int maxN=1010;  
Point PointSet[maxN];  
Point ch[maxN];  
int n;  
int len;  
  
int main()  
{  
    double ans;  
    while(scanf("%d",&n)&&n)
    {
        ans=0;
        for(int i=0;i<n;i++)  
            cin>>PointSet[i].x>>PointSet[i].y;//input the data;  
        if(n==1)
        {
            printf("0.00\n");
            continue;
        }
        if(n==2)
        {
           printf("%.2lf\n",dis(PointSet[0],PointSet[1]));
           continue;
        }
        Graham_scan(PointSet,ch,n,len);  
        ch[len]=ch[0];
        for(int i=0;i<len;i++)  
            ans+=dis(ch[i],ch[(i+1)]);
        printf("%.2lf\n",ans);
    }
    return 0;  
}