【凸包】

例题:cogs圈奶牛

叶子最可爱了

算法流程:

找左下点,以其为原点极角排序,队列维护凸包点。

算法实现:

一、左下点扫一遍

二、极角排序重载运算符,通过向量叉积判断极角大小,若极角相同则与参照点距离小者优先。

三、维护队列,加入点便在队列末端进行叉积判断极角。极角只能越来越大。

叉积

对于两个向量A,B,AxB=Xa*Yb-Xb*Ya;计算结果绝对值为两向量构成的平行四边形面积,不加绝对值称为有向面积。

如果结果大于0,则说明A的极角比B小,反之比B大。我称为“正下负上”。

明显对于一个三角形面积,任取两条边叉积除以2取绝对值即可。

代码

#include<bits/stdc++.h>
using namespace std;
#define in read()
int in{
	int cnt=0,f=1;char ch=0;
	while(!isdigit(ch)){
		ch=getchar();if(ch=='-')f=-1;
	}
	while(isdigit(ch)){
		cnt=cnt*10+ch-48;
		ch=getchar();
	}return cnt*f;
}
struct node{
	double x,y;
	node friend operator -(node a,node b){
		return {a.x-b.x,a.y-b.y};
	}
}a[10003],p[10003];
int n,tail;
double operator ^(node a,node b){
	return a.x*b.y-a.y*b.x;
}
double dis(node a,node b){
	node c=a-b;
	return sqrt(c.x*c.x+c.y*c.y);
}
double get(node a){
	return atan2(a.y,a.x);
}
bool cmpp(node A,node B){
	double c=get(A-a[1]),d=get(B-a[1]);
	if(c==d)return dis(a[1],A)<dis(a[1],B);
	return c<d;
}
int main(){
	freopen("fc.in","r",stdin);
	freopen("fc.out","w",stdout);
	n=in;
	for(int i=1;i<=n;i++){
		scanf("%lf%lf",&a[i].x,&a[i].y);
	}
	if(n==1){
		printf("0.00");return 0;
	}
	if(n==2){
		printf("%.2lf",dis(a[1],a[2]));return 0;
	}
	int k=1;
	for(int i=2;i<=n;i++){
		if(a[i].y<a[k].y||(a[i].y==a[k].y&&a[i].x<a[k].x))k=i;
	}
	swap(a[1],a[k]);
	sort(a+2,a+n+1,cmpp);
	p[++tail]=a[1],p[++tail]=a[2];
	for(int i=3;i<=n;i++){
		while(tail>=2&&((p[tail]-p[tail-1])^(a[i]-p[tail-1]))<=0)tail--;
		p[++tail]=a[i];
	}
	double sum=0;
	for(int i=1;i<tail;i++)sum+=dis(p[i],p[i+1]);
	printf("%.2lf",sum+dis(p[1],p[tail]));
	return 0;
}

 

 

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值