首页 > 其他 > 详细

poj 1113 Wall (凸包)

时间:2014-08-04 17:58:47      阅读:370      评论:0      收藏:0      [点我收藏+]

链接:poj 1113

题意:给定多边形城堡的n个顶点,绕城堡外面建一个围墙,围住所有点,

并且墙与所有点的距离至少为L,求这个墙最小的长度

思路:最小长度=城堡顶点构成的凸包的总边长+半径为L的圆的周长

先用Graham算法求出凸包,再枚举其顶点求两两之间的边长,记得加上第一个顶点和最后一个顶点的边长

最后要输出四舍五入的整数结果,可以用double存,最后用%.0lf输出

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
struct stu
{
    int x,y;
}p[1010],s[1010];
int m,top;
int chaji(struct stu p1,struct stu p2,struct stu p3)
{
    return (p1.x-p2.x)*(p3.y-p2.y)-(p3.x-p2.x)*(p1.y-p2.y);
}
double dis(struct stu p1,struct stu p2)
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)*1.0+(p1.y-p2.y)*(p1.y-p2.y));
}
int cmp(struct stu p1,struct stu p2)
{
    int k;
    k=chaji(p1,p[0],p2);
    if(k>0||(k==0&&dis(p1,p[0])<dis(p2,p[0])))
        return 1;
    return 0;
}
void graham()
{
    struct stu t;
    int k=0,i;
    for(i=1;i<m;i++)
        if(p[i].y<p[k].y||(p[i].y==p[k].y&&p[i].x<p[k].x))
            k=i;
    t=p[k];
    p[k]=p[0];
    p[0]=t;
    sort(p+1,p+m,cmp);
    s[0]=p[0];
    s[1]=p[1];
    top=1;
    for(i=2;i<m;i++){
        while(top>=1&&chaji(s[top-1],s[top],p[i])>=0)
            top--;
        top++;
        s[top]=p[i];
    }
}
int main()
{
    int n,i;
    double sum;
    while(scanf("%d%d",&m,&n)!=EOF){
        for(i=0;i<m;i++)
            scanf("%d%d",&p[i].x,&p[i].y);
        graham();
        sum=dis(s[0],s[top])+2*3.1415926*n;
        for(i=1;i<=top;i++)
            sum+=dis(s[i-1],s[i]);
        printf("%.0lf\n",sum);
    }
    return 0;
}


poj 1113 Wall (凸包),布布扣,bubuko.com

poj 1113 Wall (凸包)

原文:http://blog.csdn.net/acm_code/article/details/38370049

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!