POJ 1113 && HDU 1348 Wall（凸包）

Description
给定多边形城堡的n个顶点，绕城堡外面建一个围墙，围住所有点，并且墙与所有点的距离至少为L，求这个墙最小的长度
Input
第一行两个整数n，l（n,l<=1000）分别表示定点个数和墙与城堡最短距离，以后n行每行两个整数x，y（-10000<=x,y<=10000）
Output
墙的最短长度（结果四舍五入）
Sample Input
9 100
200 400
300 400
300 300
400 300
400 400
500 400
500 200
350 200
200 200
Sample Output
1628
Solution
城墙最短=凸包的周长 + 以L为半径的圆的周长。
首先，我们得出城墙最短=凸包的周长+相应顶点的弧长（半径都为L）之和。
再证明：相应顶点的弧长（半径都为L）之和=以L为半径的圆的周长。
事实上，设凸包顶点为n,n个顶点组成n个周角，角度为360*n=2*180*n，凸包的内角和为180*(n-2)，作了2*n条垂线，和为2*n*90=180*n,所以总圆周角为2*180*n-180*(n-2)-180*n=360，组成圆
Code

#include<stdio.h>
#include<stdlib.h>
#include<math.h>

typedef struct 
{
    double x;//横坐标 
    double y;//纵坐标 
    double k;//与第一个定点连线的正切值 
}point;
point a[50005];
point b[50005];
int n,res;
double multiply(point x1,point x2,point x3)//叉乘 
{
    return ((x3.y-x2.y)*(x2.x-x1.x)-(x2.y-x1.y)*(x3.x-x2.x));
}
double distance(point x1,point x2)//两点间间距 
{
    double c;
    c=(x1.x-x2.x)*(x1.x-x2.x)+(x1.y-x2.y)*(x1.y-x2.y);
    return sqrt(c);
}
int compare1(const void*x1,const void*x2)//按坐标排序 
{
    point *x3,*x4;
    x3=(point*)x1;
    x4=(point*)x2;
    if(x3->y==x4->y)
        return x3->x>x4->x?1:-1;
    else
        return x3->y>x4->y?1:-1;
}
int compare2(const void*x1,const void*x2)//按斜率排序 
{
    point *x3,*x4;
    x3=(point*)x1;
    x4=(point*)x2;
    if(x3->k==x4->k)
    {
        if(x3->x==x4->x)
            return x4->y>x3->y?1:-1;
        else
            return x4->x>x3->x?1:-1;
    }
    else
        return x3->k>x4->k?1:-1;
}
int main()
{
    int i,j,t,flag1,flag2;
    double dis,len=0;
    float r;
    while(scanf("%d",&n)!=EOF)
    {
        scanf("%f",&r);
        for(i=0;i<n;i++)
            scanf("%lf%lf",&a[i].x,&a[i].y);
        qsort(&a[0],n,sizeof(double)*3,compare1);//对顶点排序 
        for(i=1;i<n;i++)
        {
            a[i].k=atan((double)((a[i].y-a[0].y)/(a[i].x-a[0].x))); 
            if(a[i].k<(double(0)))//若正切值小于0则加PI 
                a[i].k+=acos((double)-1);
        }
        qsort(&a[1],n-1,sizeof(double)*3,compare2);//对斜率排序 
        b[0]=a[0];
        b[1]=a[1];
        b[2]=a[2];
        res=3;
        for(i=3;i<n;i++)
        {
            while(res>2&&multiply(b[res-2],b[res-1],a[i])<=0)//若加入新点后破坏凸性则将凹的部分的点从末尾除去 
                res--;
            b[res++]=a[i];//加入新点 
        }
        for(i=0;i<res-1;i++)//求凸包周长 
            len+=distance(b[i],b[i+1]);
        len+=distance(b[res-1],b[0]);//加上起点与终点间距 
        len+=2*r*acos((double)-1);//以r为半径的圆的周长 
        printf("%.f\n",len);//四舍五入输出 
    }
    return 0;
}