POJ 1912 凸包

题意

传送门 POJ 1912

题解

直线在全部点的一侧，考虑凸包即可，凸包算法需要至少 $2$ 个点，需要特别处理点数小于 $2$ 的情况。最简单地，枚举点判断是否在直线同一侧，复杂度 $O(nm)$，需要优化；考虑与当前处理直线斜率相同的直线，与凸包的 $2$ 个切点即可；二分凸包上大于直线斜率正负值对应弧度的最小弧度的边，其端点即为切点，复杂度 $O(nlogm)$。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <queue>
#include <vector>
#include <cmath>
#define min(a,b)    (((a) < (b)) ? (a) : (b))
#define max(a,b)    (((a) > (b)) ? (a) : (b))
#define abs(x)    ((x) < 0 ? -(x) : (x))
#define INF 0x3f3f3f3f
#define delta 0.85
using namespace std;

#define PI 3.14159265358979323846
#define EPS 1e-10

// 考虑误差的加法运算
double add(double a, double b){
	if(abs(a + b) < EPS * (abs(a) + abs(b))) return 0;
	return a + b;
}

// 二维向量结构体
struct P{
	double x, y;
	P(){}
	P(double x, double y) : x(x), y(y) {}
	P operator + (P p){
		return P(add(x, p.x), add(y, p.y));
	}
	P operator - (P p){
		return P(add(x, -p.x), add(y, -p.y));
	}
	P operator * (double d){
		return P(x * d, y * d);
	}
	double dot(P p){	// 内积
		return add(x * p.x, y * p.y);
	}
	double det(P p){	// 外积
		return add(x * p.y, -y * p.x);
	}
};

// 字典序比较
bool cmp_x(const P& p, const P& q){
	if(p.x != q.x) return p.x < q.x;
	return p.y < q.y;
}

// 求凸包
vector<P> convex_hull(P* ps, int n){
	sort(ps, ps + n, cmp_x);
	int k = 0;				// 凸包的顶点数
	vector<P> qs(n * 2);	// 构造中的凸包
	// 构造凸包的下侧
	for(int i = 0; i < n; i++){
		while(k > 1 && (qs[k - 1] - qs[k - 2]).det(ps[i] - qs[k - 1]) <= 0) k--;
		qs[k++] = ps[i];
	}
	// 构造凸包的上侧
	//i=n-2;
	//i = n - 2, 
	for(int i = n - 2, t = k; i >= 0; i--){
		while(k > t && (qs[k - 1] - qs[k - 2]).det(ps[i] - qs[k - 1]) <= 0) k--;
		qs[k++] = ps[i];
	}
	qs.resize(k - 1);
	return qs;
}

// 计算点 q 是否在直线下方
int below_line(P q, P p1, P p2){
	if(!cmp_x(p1, p2)) swap(p1, p2);
	return (p2 - p1).det(q - p1) < 0 ? 1 : 0;
}

// 计算向量 p 的弧度
double atan2(const P& p){
	return atan2(p.y, p.x);
}

#define MAX_N 100000
typedef pair<double, int> PP;
int N;
P house[MAX_N];
PP angle[MAX_N];

// 考虑误差的比较
bool cmp_PP(PP a, PP b){
	return add(a.first, -b.first) < 0;
}

void solve(){
	vector<P> qs;
	int sz;
	// 凸包算法保证至少 2 个点
	if(N > 1){
		qs = convex_hull(house, N);
		sz = qs.size();
		qs.push_back(qs[0]);
		for(int i = 0; i < sz; i++){
			angle[i] = PP(atan2(qs[i + 1] - qs[i]), i);
		}
		sort(angle, angle + sz, cmp_PP);
	}
	
	P q, p;
	while(~scanf("%lf%lf%lf%lf", &q.x, &q.y, &p.x, &p.y)){
		if(N <= 1){
			puts("GOOD");
			continue;
		}
		int i1 = upper_bound(angle, angle + sz, PP(atan2(p - q), 0), cmp_PP) - angle;
		int i2 = upper_bound(angle, angle + sz, PP(atan2(q - p), 0), cmp_PP) - angle;
		// 处理边界情况
		i1 = i1 == sz ? angle[0].second : angle[i1].second;
		i2 = i2 == sz ? angle[0].second : angle[i2].second;
		puts(1 & (below_line(qs[i1], q, p) ^ below_line(qs[i2], q, p)) ? "BAD" : "GOOD");
	}
}

int main(){
	while(~scanf("%d", &N)){
		for(int i = 0; i < N; i++) scanf("%lf%lf", &house[i].x, &house[i].y);
		solve();
	}
	return 0;
}