poj1113(凸包卷包裹法)

最新推荐文章于 2021-10-13 14:40:39 发布

原创最新推荐文章于 2021-10-13 14:40:39 发布 · 443 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#凸包

计算几何专栏收录该内容

5 篇文章

订阅专栏

本文介绍了一个基于Jarvis步进法的凸包算法实现过程。该算法首先选择最下端左侧的点作为起始点，然后按逆时针方向确定其余边界点。文中详细解释了如何通过叉积判断点的位置关系，并提供了完整的C++代码实现。

#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;
}