【bzoj1670】[Usaco2006 Oct]Building the Moat护城河的挖掘凸包

最新推荐文章于 2018-12-19 14:56:00 发布

原创最新推荐文章于 2018-12-19 14:56:00 发布 · 517 阅读

0 ·

CC 4.0 BY-SA版权

计算几何专栏收录该内容

11 篇文章

订阅专栏

本文详细介绍了裸的凸包算法的实现过程，包括结构定义、运算符重载、距离计算、排序和Graham扫描算法的实现，最终输出凸包周长。

裸的凸包啦，第一次写。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iostream>
#define maxn 5010 

using namespace std;

struct yts
{
	long long x,y;
}a[maxn],s[maxn];

long long operator*(yts x,yts y)
{
	return x.x*y.y-x.y*y.x;
}

yts operator-(yts x,yts y)
{
	yts ans;
	ans.x=y.x-x.x;ans.y=y.y-x.y;
	return ans;
}

double dis(yts x,yts y)
{
	return sqrt((double)(x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y));
}

int n,m;
double ans;

bool cmp(yts x,yts y)
{
	int t=(x-a[1])*(y-a[1]);
	if (t==0) return dis(x,a[1])<dis(y,a[1]);
	return t<0;
}

void graham()
{
	int t=1;
	for (int i=2;i<=n;i++)
	  if (a[i].x<a[t].x || (a[i].x==a[t].x && a[i].y<a[t].y)) t=i;
	swap(a[t],a[1]);
	sort(a+2,a+n+1,cmp);
	int top=1;
	s[top]=a[1];
	for (int i=2;i<=n;i++)
	{
		while (top>1 && (s[top]-s[top-1])*(a[i]-s[top])>=0) top--;
		s[++top]=a[i];
	}
	s[top+1]=a[1];
	for (int i=1;i<=top;i++) ans+=dis(s[i],s[i+1]);
}

int main()
{
	scanf("%d",&n);
	for (int i=1;i<=n;i++) scanf("%lld%lld",&a[i].x,&a[i].y);
	graham();
	printf("%.2lf\n",ans);
	return 0;
}