poj 1113Walls解题报告

链接：http://poj.org/problem?id=1113

求凸包的模板题，只是这个精度卡死人，而且我的g++死活过不了，换成c++就过了，唉，poj上太多这样的题了，实在让人无语啊

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define N 1005
#define pi 3.141592653
struct point
{
    int x,y;
};
int n;
int dis(point p1,point p2)
{
    return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
point p[N],res[N];
int mul(point p1,point p2,point p3)
{
    return (p2.x-p1.x)*(p3.y-p1.y)-(p2.y-p1.y)*(p3.x-p1.x);
}
int cmp(const void *a,const void *b)//极角排序，要注意快速排序的排序函数
{
    point c=*(point *)a;
    point d=*(point *)b;
    int t=mul(p[0],c,d);
    if(t)
    return -t;
    return dis(p[0],d)-dis(p[0],c);
}
int graham()
{
    int i,j,k,top=2;
    j=0;
    for(i=1;i<n;i++)
    {
        if(p[i].y<p[j].y||(p[i].y==p[j].y&&p[i].x<p[j].x))
        j=i;
    }
    if(j)
    {
        point temp=p[j];
        p[j]=p[0];
        p[0]=temp;
    }
    qsort(p+1,n-1,sizeof(p[0]),cmp);
    //for(i=0;i<n;i++)
    //printf("%d %d\n",p[i].x,p[i].y);
    res[0]=p[0];
    res[1]=p[1];
    res[2]=p[2];
    for(i=3;i<n;i++)
    {
        while(top>1&&mul(res[top-1],res[top],p[i])<=0)
        --top;
        res[++top]=p[i];
    }
    return top;
}
int main()
{
    int i,j,k;
    int l;
    while(scanf("%d%d",&n,&l)!=EOF)
    {
        for(i=0;i<n;i++)
        scanf("%d%d",&p[i].x,&p[i].y);
        int t=graham();
        res[t+1]=res[0];
        //for(i=0;i<=t+1;i++)
        //printf("%.lf %.lf\n",res[i].x,res[i].y);
        double ans=0;
        for(i=0;i<=t;i++)
        ans+=sqrt((double)dis(res[i],res[i+1]));
        ans+=2*pi*l;
        printf("%.lf\n",ans);
    }
    return 0;
}