BZOJ 1670: [Usaco2006 Oct]Building the Moat护城河的挖掘|二维凸包

纠结之后还是写的 水平序的凸包

极点序也不好写 可能有隐含错误

凸包裸题

#include <cmath>
#include <queue>
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define LL long long
#define ULL unsigned long long
const int mo=1e9+7;
const int MAXN=5e3+10;
const int MAXM=1e6+10;
using namespace std;

struct P
{
	int x,y;
}num[MAXN];

bool cmp(P a,P b){if(a.x==b.x) return a.y<b.y;return a.x<b.x;}

int q[MAXN],top;
bool b[MAXN];

int ff(int x1,int y1,int x2,int y2){return x1*y2-y1*x2;}
int f(int z,int y,int x)
{
	return ff(num[y].x-num[x].x,num[y].y-num[x].y,num[z].x-num[y].x,num[z].y-num[y].y);
}
double dis(int x1,int y1,int x2,int y2)
{
	return sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));
}
int main()
{
	int n;
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d %d",&num[i].x,&num[i].y);
	sort(num+1,num+n+1,cmp);
	
	top=0;
	for(int i=1;i<=n;i++){while(top>1&&f(i,q[top],q[top-1])>0) top--;q[++top]=i;}
	for(int i=2;i<=top;i++) b[q[i]]=true;
	for(int i=n;i>=1;i--) if(!b[i]){while(f(i,q[top],q[top-1])>0) top--;q[++top]=i;}
	double ans=0;
	for(int i=2;i<=top;i++) ans+=dis(num[q[i]].x,num[q[i]].y,num[q[i-1]].x,num[q[i-1]].y);
	printf("%.2lf\n",ans);
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值