【POJ 3304】Segments（计算几何）

Segments

题目链接：POJ 3304

题目大意

给你若干条线段，然后问你是否存在一条直线使得所有线段投影到这条直线之后都有一个公共点。

思路

考虑转化一下题意：
找是否一条直线，使得它穿过了所有的线段。

那我们考虑这样：如果有这么一条直线，那我们可以通过平移，使它会碰到任何两个端点。
那我们就可以枚举任意两个端点，然后判断这条直线是否与所有的线段有交。（复杂度是 $O(n^3)$ 可以通过）

那至于怎么判断直线 AB 和线段 CD 是否有交。
可以用这个方法：看 CA,CD,CB 这三条向量之前的关系。
如果 CD 在 CA 和 CB 的中间或者跟其中一个重合，就是有交的。

然后这个可以用叉积来判断。

代码

#include<cmath>
#include<cstdio>

using namespace std;

int T, n;
struct dot {
	double x, y;
};
struct line {
	dot x, y;
}a[101];
const double eps = 1e-8;

double Abs(double x) {
	return x < 0 ? -x : x;
}

dot operator +(dot x, dot y) {
	return (dot){x.x + y.x, x.y + y.y};
}

dot operator -(dot x, dot y) {
	return (dot){x.x - y.x, x.y - y.y};
}

double operator *(dot x, dot y) {
	return x.x * y.x + x.y * y.y;
}

double operator ^(dot x, dot y) {
	return x.x * y.y - x.y * y.x;
}

double get_dis(dot x, dot y) {
	return sqrt((x - y) * (x - y));
}

double cha(dot z, dot x, dot y) {
	return (x - z) ^ (y - z);
}

int sgn(double x) {
	if (Abs(x) < eps) return 0;
	if (x > eps) return 1;
	return -1;
}

bool line_meet(line x, line y) {//x为直线，y为线段
	if (sgn(cha(y.x, x.x, x.y)) * sgn(cha(y.y, x.x, x.y)) > 0) return 0;//两个叉积符号要不同，也就是岔开来，才能穿过
	return 1;
} 

bool check(dot l, dot r) {
	if (!sgn(get_dis(l, r))) return 0;
	for (int i = 1; i <= n; i++) {
		if (!line_meet((line){l, r}, a[i])) return 0;
	}
	return 1;
}

int main() {
	scanf("%d", &T);
	while (T--) {
		scanf("%d", &n);
		for (int i = 1; i <= n; i++) {
			scanf("%lf %lf %lf %lf", &a[i].x.x, &a[i].x.y, &a[i].y.x, &a[i].y.y);
		}
		
		bool yes = 0;
		for (int i = 1; i <= n; i++) {
			for (int j = 1; j <= n; j++) {
				if (check(a[i].x, a[j].x) || check(a[i].x, a[j].y) || check(a[i].y, a[j].x) || check(a[i].y, a[j].y)) {//四种搭配
					yes = 1;
					break;
				}
			}
			if (yes) break;
		}
		
		if (yes) printf("Yes!\n");
			else printf("No!\n");
	}
	
	return 0;
}