poj 1113 Wall（凸包）

题目：http://poj.org/problem?id=1113

题意：给定多边形城堡的n个顶点，绕城堡外面建一个围墙，围住所有点，并且墙与所有点的距离至少为L，求这个墙最小的长度。

公式：

城堡围墙长度最小值 = 城堡顶点坐标构成的散点集的凸包总边长 + 半径为L的圆周长

#include <iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define PI 3.1415926535
using namespace std;
const double eps=1e-6;
typedef struct node
{
    int x,y;
}point;
int n;
point pt[1100];
point res[1100];
bool cmp(point a,point b)
{
    if(a.x==b.x)
    return a.y<b.y;
    else
    return a.x<b.x;
}
int dcml(double x)
{
    if(fabs(x)<eps)
    return 0;
    if(x<0)
    return -1;
    else
    return 1;
}
double dot(int x1,int y1,int x2,int y2)
{
    return x1*y2-x2*y1;
}
int cross(point a,point b,point c)
{
    return dcml(dot(a.x-c.x,a.y-c.y,b.x-c.x,b.y-c.y));
}
int dist(point a,point b)
{
    return ((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int main()
{
    int L;
    scanf("%d%d",&n,&L);
    int i;
    for(i=0;i<n;i++)
    {
        scanf("%d%d",&pt[i].x,&pt[i].y);
    }
    sort(pt,pt+n,cmp);
    int m;
    m=0;
    for(i=0;i<n;i++)
    {
        while(m>1&&cross(res[m-1],pt[i],res[m-2])<=0)
        m--;
        res[m++]=pt[i];
    }
    int k=m;
    for(i=n-2;i>=0;i--)
    {
        while(m>k&&cross(res[m-1],pt[i],res[m-2])<=0)
        m--;
        res[m++]=pt[i];
    }
    if(n>1)
    m--;
    double dis=0;
    for(i=0;i<m-1;i++)
    {
        dis+=sqrt(dist(res[i],res[i+1]));
    }
    dis+=sqrt(dist(res[0],res[m-1]));
    dis+=2*PI*L;
    printf("%d",(int)(dis+0.5));
    return 0;
}