1015 - 计算几何之Graham扫描法求凸包 - Cows（POJ 3348）

传送门

 

题意

给你一堆点，求这些点的凸包，并求出面积

 

分析

很久之前就做过的一道题了，还记得那是凯爷（凯爷好厉害好厉害的）讲的，是Jarris步进法：按照横纵坐标对所有的点进行排序（横坐标优先）

然后就是和Graham类似的方法了，边扫描边弹出栈，只是Jarris算法要进行两遍扫描，一遍解决了下面一半的凸包，一遍解决上面一半的凸包，拼在一起就是一个完整的凸包咯

但这不是今天要讲的重点，重点是Graham算法

其实也就是在开始预处理的时候排序的方式与Jarris不一样而已，

先找到横坐标最小的那一个点，这个点肯定是在凸包上的

然后依次对每个点进行极角排序，

然后再边扫描边弹

但实际上，在操作的时候我们并没有使用极角排序，我们直接按照每个点的偏离程度来排，越往下的越靠前

要是遇到奇葩情况：三点共线，我们就把离第一个点（横坐标最小的）最近的放在前面。因为这样在后面扫描的时候就可以把这个点弹掉了（在凸包上为了方便，我们都保留同一条直线上最远的那一个点，这也就是为什么在判弹出条件的时候是>=0了）

 

代码

 

Jarris算法

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int eps=1e-8;
//Your task:构造一个凸包并且计算它的面积然后除以50 
struct point{
	double x,y;
}p[10009],a[10009],init;
inline double cross(const point &a,const point &b,const point &c){
	return ((c.x-a.x)*(b.y-a.y)-(c.y-a.y)*(b.x-a.x))/2;
}
int n;
int cmp(const point &a,const point &b){
	if(a.x<b.x) return 1;
	return 0;
}
int main(){
	scanf("%d",&n);
	int i,j;
	for(i=1;i<=n;++i)
		scanf("%lf%lf",&a[i].x,&a[i].y);
	sort(a+1,a+n+1,cmp);
	int t=0;
	p[++t]=a[1];
	p[++t]=a[2];
	for(i=3;i<=n;++i){
		p[++t]=a[i];
		while(cross(p[t-1],p[t-2],p[t])<=0&&t>2){
			p[t-1]=p[t];
			t--;
		}
	}
	for(i=n-1;i>=1;--i){
		p[++t]=a[i];
		while(cross(p[t-1],p[t-2],p[t])<=0&&i!=n-1){
			p[t-1]=p[t];
			t--;
		}
	}
	init.x=0,init.y=0;
	double s=0;
	t--;
	for(i=1;i<=t;++i){
		s+=cross(p[i],p[i+1],init);
	}
	s=fabs(s);
	
	s=s/50;
	printf("%d",(int)s);
	return 0;
}

 

 

Graham算法

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#define N 10009
using namespace std;
struct point{
	int x,y;
	point(){}
	point(int _x,int _y):x(_x),y(_y){}
	friend inline point operator +(const point &a,const point &b)
	{return point(a.x+b.x,a.y+b.y);}
	friend inline point operator -(const point &a,const point &b)
	{return point(a.x-b.x,a.y-b.y);}
	friend inline point operator *(double k,const point &a)
	{return point(a.x*k,a.y*k);}
	friend inline int dot(const point &a,const point &b){
		return a.x*b.x+a.y*b.y;
	}
	friend inline int cross(const point &a,const point &b){
		return a.x*b.y-b.x*a.y;
	}
}a[N],b[N],S;
int n,m;
bool cmp1(const point &a,const point &b){
	if(a.x!=b.x) return a.y<b.y;
        return a.x<b.x;
    //上面是错误的，下面才是正确的
    //上面那个能A纯粹是数据水，不存在a.x==b.x的情况，这样就能找到最低点
    //if(a.x==b.x) return a.y<b.y;
    //return a.x<b.x;
	
}
bool cmp2(const point &a,const point &b){
	int hh=cross(a-S,b-S);
	if(hh==0) return dot(a-S,a-S)<dot(b-S,b-S);/ 如果直接写dot(a,a)求的就是到原点的距离了 
	return hh>0;
}
void Graham(){
	sort(a+1,a+n+1,cmp1);
	S=a[1];
	sort(a+2,a+n+1,cmp2);
	m=1;b[1]=a[1];
	for(int i=2;i<=n;++i)
	{
		while(m>=2&&cross(a[i]-b[m-1],b[m]-b[m-1])>=0) --m;///是a[i]啊！！ 
		b[++m]=a[i];/就是a[i]
	}
}
int main(){
	scanf("%d",&n);
	int i,j,k;
	for(i=1;i<=n;++i) scanf("%d%d",&a[i].x,&a[i].y);
	Graham();
	int res=0;b[m+1]=b[1];
	for(i=1;i<=m;++i)
		res+=cross(b[i],b[i+1]);
	printf("%d",abs(res/100));
	return 0;
}