【解题报告】 NYOJ 78 -- 圈水池 Andrew 凸包算法

题目连接: NYOJ 78凸包入门水题

不会凸包的移步这里

// NYOJ 78 -- 圈水池 Andrew 凸包算法 O(nlogn)
// 将凸包的顶点按 从小到大 先x后y的 顺序输出
//
/*test data
9 100
200 400
300 400
300 300
400 300
400 400
500 400
500 200
350 200
200 200
	== 1628

*/

#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;

const int MAXP = 1010;
const double Pzero = 0.00001; 

struct POINT{
	int x,y;
}p[MAXP],CHp[MAXP*2];

double XJ(POINT a, POINT b, POINT c){
	// vector a->b (b.x-a.x, b.y-a.y)
	// vector a->c (c.x-a.x, c.y-a.y)
	return (b.x-a.x)*(c.y-a.y) - (c.x-a.x)*(b.y-a.y);
}

bool comp(POINT a,POINT b){
	return (a.x < b.x) || ((a.x==b.x) && (a.y<b.y));
}

int AndrewConvexHull(POINT *CHp, int n){
	// 计算凸包 CHp[] 逆时针存点 
	sort(p, p + n, comp);
	int nCHp = 0;
	for (int i = 0 ; i < n ; CHp[nCHp++] = p[i++])
		while(nCHp >= 2 && XJ(CHp[nCHp-2], CHp[nCHp-1], p[i]) <= 0)
			nCHp--;
	for (int i = n - 1,m = nCHp + 1 ; i >= 0 ; CHp[nCHp++] = p[i--])
		while(nCHp >= m && XJ(CHp[nCHp-2], CHp[nCHp-1], p[i]) <= 0)
			nCHp--;
	return --nCHp;
}

void init(int n){
	// 输入
	for (int i = 0 ; i < n; i++)
		scanf("%d%d", &p[i].x, &p[i].y);
}

void Output(int nCHp){
	sort(CHp, CHp + nCHp, comp);
	for (int i = 0; i < nCHp ; i++){
		printf("%d %d\n", CHp[i].x, CHp[i].y);
	}
}

int main()
{
//	freopen("in.txt","r",stdin);
	int n;scanf("%d",&n);
	while(~scanf("%d",&n) ){

		init(n);
		int nCHp = AndrewConvexHull(CHp, n);
		Output(nCHp);
	}
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值