[ACM_几何] Wall

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=28417#problem/E

题目大意：依次给n个点围成的一个城堡，在周围建围墙，要求围墙离城墙的距离大于一定的值，求围墙最短长度（结果四舍五入

解题思路：求围住所有点的凸包周长+一个圆的周长

#include<iostream>

#include<cmath>

#include<string.h>

#include<string>

#include<stdio.h>

#include<algorithm>

#include<iomanip>



using namespace std;

#define eps 0.0000000001

#define PI acos(-1.0)



//*******************************************************************************

//点和向量

struct Point{

    double x,y;

    Point(double x=0,double y=0):x(x),y(y){}

};

typedef Point Vector;

Vector operator-(Vector a,Vector b){return Vector(a.x-b.x,a.y-b.y);}

bool operator<(const Vector& a,const Vector& b){return a.x<b.x||(a.x==b.x && a.y<b.y);}

int dcmp(double x){

    if(fabs(x)<eps)return 0;

    else return x<0 ? -1:1;

}

double Dot(Vector A,Vector B){return A.x*B.x+A.y*B.y;}//向量点积

double Length(Vector A){return sqrt(Dot(A,A));}//向量模长

double Cross(Vector A,Vector B){return A.x*B.y-A.y*B.x;}

//********************************************************************************

//计算凸包输入点数组p，个数n，输出点数组ch，返回凸包定点数

//输入不能有重复，完成后输入点顺序被破坏

//如果不希望凸包的边上有输入点，把两个<=改成<

//精度要求高时，建议用dcmp比较

//基于水平的Andrew算法-->1、点排序2、删除重复的然后把前两个放进凸包

//3、从第三个往后当新点在凸包前进左边时继续，否则一次删除最近加入的点，直到新点在左边

int ConVexHull(Point* p,int n,Point*ch){

    sort(p,p+n);

    int m=0;

    for(int i=0;i<n;i++){//下凸包

        while(m>1 && Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<=0)m--;

        ch[m++]=p[i];

    }

    int k=m;

    for(int i=n-2;i>=0;i--){//上凸包

        while(m>k && Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<=0)m--;

        ch[m++]=p[i];

    }

    if(n>1)m--;

    return m;

}

//*******************************************************************

int main(){

    int n,L;

    while(cin>>n>>L){

        

        Point p[1005],ch[1005];

        

        for(int i=0;i<n;i++)  //输入点

            cin>>p[i].x>>p[i].y;



        int m=ConVexHull(p,n,ch);   //求解凸包并计算凸包周长

        double sum=Length(ch[0]-ch[m-1]);

        for(int i=1;i<m;i++)

            sum+=Length(ch[i]-ch[i-1]);



        printf("%.0lf\n", sum+2*L*PI);//输出周长+圆周并四舍五入

    }

    return 0;    

}

View Code