POJ 3335 Rotating Scoreboard 半平面的交

题目描述：http://poj.org/problem?id=3335

解题思路：

把多边形的每条边看做一条直线，只有直线一侧的区域是合法的。有解的情况当且仅当每条直线划分出的区域的交集非空。即半平面的交非空。

求解半平面的交根据需求有不同的算法，求解出区域的凸包有O(n^2)和O(nlogn)的算法，另外用线性规划求解区域中的一点的算法利用随机化期望复杂度可以达到O(n)。

因为题目数据量很小，我写的是O(n^2)算法，简单粗暴。

算法的核心思想就是不断用边来切割已有的解的凸多边形点集。（半平面的交区域必定为凸多边形）

首先将区域初始化为一个很大的矩形区域。比如(-1e6,-1e6),(-1e6,1e6),),(1e6,1e6),(1e6,-1e6)。然后O(n)枚举多边形的边，对于每条边判断点集中的点是否在这条边正确的一侧，如果在则留着；如果不在就删掉，并看与它相连的两条边（线段）是否和这条直线有交点，如果有就添加到解集中。最后得到的点集就是合法区域的顶点。如果顶点数为0了就代表区域为空。

具体实现时，直线由两个点p1,p2确定，并且这两个点是有向的，即直线确定的半平面是p1p2向量的右边（如果定为左边也行，p1p2换个顺序就行了）。这样就可以用叉积来判断点是否在半平面内了。

/*
 *	14.11.16
 *	Problem: 3335		
 *	Memory: 160K		Time: 0MS
 *	Language: C++		Result: Accepted
 *
 */
#include "iostream"
#include "cmath"
#define EPS 1e-8
#define MAXN 207
#define INF 1e6

struct Point{  
    double x,y;  
    Point(){}  
    Point(double _x,double _y):x(_x),y(_y){}  
}p[MAXN], c[MAXN], t[MAXN];;  

int n, cn;

void init() {
	scanf("%d", &n);
	for (int i=1; i<=n; i++) scanf("%lf %lf", &p[i].x, &p[i].y);
	p[n+1] = p[1];
	c[1] = Point(-INF, -INF);
	c[2] = Point(-INF, INF);
	c[3] = Point(INF, INF);
	c[4] = Point(INF, -INF);
	c[0] = c[cn=4]; // 枚举时避免特判1和n的情况
	c[cn+1] = c[1];
}

double cross(Point a, Point b, Point c) {
	double x1 = b.x-a.x, y1 = b.y-a.y;
	double x2 = c.x-a.x, y2 = c.y-a.y;
	return x1*y2-x2*y1;
}

void intersect(Point &t, Point p1, Point p2, Point p3, Point p4) {
	double a = p2.y-p1.y, b = -(p2.x-p1.x), c = p1.x*p2.y-p2.x*p1.y;
	double d = p4.y-p3.y, e = -(p4.x-p3.x), f = p3.x*p4.y-p4.x*p3.y;
	t.x = (c*e-b*f)/(a*e-b*d);
	t.y = (c*d-a*f)/(b*d-a*e);
}

bool cut(Point a, Point b) {
	int tn = 0;
	for (int i=1; i<=cn; i++) {
		if (cross(a, b, c[i])<=EPS) { // 如果点在半平面内
			t[++tn] = c[i];	
		} else {	// 如果点不在半平面内则查看它左右连接的两条边与直线的交点
			if (cross(a, b, c[i-1])<=-EPS) intersect(t[++tn], c[i-1], c[i], a, b);
			if (cross(a, b, c[i+1])<=-EPS) intersect(t[++tn], c[i], c[i+1], a, b);
		}
	}
	if (!tn) return false;
	for (int i=1; i<=tn; i++) c[i] = t[i];
	c[0] = c[cn=tn];
	c[cn+1] = c[1];
	return true;
}

int main() {
	int T;
	scanf("%d", &T);
	while (T--) {
		init();
		bool ok = true;
		for (int i=1; i<=n; i++) {
			if (!cut(p[i], p[i+1])) {
				ok = false;
				break;
			} 
		}
		if (ok) {
			printf("YES\n");
		} else {
			printf("NO\n");
		}
	}
	return 0;
}