poj1113(凸包 卷包裹法)

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;

struct Point{
int x,y;
}s[1005];

bool cmp(Point a,Point b) //选择最下左的点作为起点
{
	if(a.y==b.y) return a.x<b.x; 
	else return a.y<b.y;
}

int cross(Point p1,Point p2,Point p3) //叉积
{
	int tmp=(p3.x-p2.x)*(p1.y-p2.y)-(p1.x-p2.x)*(p3.y-p2.y);
	if(tmp>=0) return 1;
	else return -1;
}

int n,l,st[1005],ans[1005],tail,cnt;

void Jarvis()
{
	sort(s,s+n,cmp);
	tail=-1,cnt=0;
	st[++tail]=0,st[++tail]=1;
	for(int i=2;i<n;i++) //因为排过序故到达凸包最顶端的点停止,即包住一半
	{
		while(tail>0&&cross(s[st[tail-1]],s[st[tail]],s[i])<0) //右转不行
			tail--;
		st[++tail]=i;
	}
	for(int i=0;i<=tail;i++)
		ans[cnt++]=st[i];

	tail=-1;
	st[++tail]=n-1,st[++tail]=n-2; //从顶端开始包住另一半，把起点又算了进去
	for(int i=n-3;i>=0;i--)
	{
		while(tail>0&&cross(s[st[tail-1]],s[st[tail]],s[i])<0)
			tail--;
		st[++tail]=i;
	}
	for(int i=0;i<=tail;i++)
		ans[cnt++]=st[i];

}

double dist(Point a,Point b) //计算两点间的距离
{
	double tmp=(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
	return sqrt(tmp);
}
int main()
{
	cin>>n>>l;
	for(int i=0;i<n;i++)
		cin>>s[i].x>>s[i].y;
	Jarvis();
	//cout<<cnt<<endl;
	double re=2*acos(-1.0)*l;
	for(int i=0;i<cnt-1;i++)
	{
		re+=dist(s[ans[i]],s[ans[i+1]]);
	}
	printf("%.0f\n",re);
	return 0;
}