POJ - 1113 Wall 凸包求周长

题目链接：https://cn.vjudge.net/problem/POJ-1113

题解：求下凸包上的点，然后计算围起来的周长，加上2*PI*L即可

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
#define eps 1e-8
#define PI acos(-1)
struct Point{
	int x, y;
	Point (int x_=0,int y_=0)
	{
		x=x_;y=y_;
	}
	Point operator + (const Point &xx)const 
	{
		return Point(x+xx.x,y+xx.y);
	}
	Point operator - (const Point &xx)const 
	{
		return Point(x-xx.x,y-xx.y);
	}
	
	
}a[50010],b[50010];
int cross(Point xx,Point yy)
{
	return (xx.x*yy.y-xx.y*yy.x);
}
double dis(Point xx,Point yy)
{
	return sqrt(1.0*(xx.x-yy.x)*(xx.x-yy.x) + (xx.y-yy.y)*(xx.y-yy.y));
}
int dist(Point xx,Point yy)
{
	return (xx.x-yy.x)*(xx.x-yy.x) + (xx.y-yy.y)*(xx.y-yy.y);
}
int n,m,l;
bool cmp(Point xx,Point yy)
{
	int O=cross(xx-a[1],yy-a[1]);
	if(O>0 || O==0&&dis(a[1],xx)<dis(a[1],yy)) return 1;
	else return 0;
}
void solve()
{
	for(int i=1;i<=min(2,n);i++)
		b[m++]=a[i];
	for(int i=3;i<=n;i++)
	{
		while(m>1 && cross(b[m-1]-b[m-2],a[i]-b[m-2])<=0)m--;
		b[m++]=a[i];
	}
	double ans = 2*PI*l;
	for(int i = 0;i < m;i ++)
		ans+=dis(b[i],b[(i+1)%m]);
	printf("%d\n",(int)(ans+0.5));
}
int main()
{
	while(~scanf("%d%d",&n,&l))
	{
		int k=1;
		m=0;
		for(int i=1;i<=n;i++)
		{
			scanf("%d%d",&a[i].x,&a[i].y);
			if(a[i].y<a[k].y|| (a[i].y==a[k].y&&a[i].x<a[k].x))
			{
				k=i;
			}
		}
		if(n==1)
		{
			printf("0\n");
			continue;
		}
		swap(a[1].x,a[k].x);swap(a[1].y,a[k].y);
		sort(a+2,a+1+n,cmp);
		solve();
	}
	return 0;
}