【POJ 2318】TOYS（计算几何）

TOYS

题目链接：POJ 2318

题目大意

给你一个长方形，然后有若干条不相交的线把它分成若干个区域。
然后有一些点，保证不会再线上，然后问你每个区域有多少个点。

思路

首先我们发现它线互不相交，而且点是不会在线上。

那我们可以考虑直接枚举每个区间看它在不在这个区间里面。
然后怎么判断呢？我们可以直接就要求它在左边线的右边，在右边先的左边即可。
（具体来讲我们可以求出某条线段在某个 $y$ 值，即你要判断点的 $y$ 值的位置上它的 $x$ 值是多少，然后跟要判断的点的 $x$ 值来判断即可）

然后这样搞是 $O(nm)$ 的，勉勉强强能过。
然后其实你可以就直接变成找这个点右边的第一个线，然后直接用二分来判断。（$O(m\log n)$，这里就没有写了）

代码

#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

int n, m, x, y, xx, yy;
int ans[5001], X, Y;
double k[5001];
struct node {
	int x, y;
}a[5002];

bool check(int num) {
	double x = 1.0 * a[num].x + k[num] * (Y - y);
	return x > X;
}

int main() {
	scanf("%d", &n);
	while (n) {
		scanf("%d %d %d %d %d", &m, &x, &y, &xx, &yy);
		
		memset(ans, 0, sizeof(ans));
		a[0] = (node){x, x}; a[n + 1] = (node){xx, xx};
		k[0] = k[n + 1] = 0;
		for (int i = 1; i <= n; i++) {
			scanf("%d %d", &a[i].x, &a[i].y);
			k[i] = 1.0 * (a[i].y - a[i].x) / (yy - y);
		}
		
		for (int i = 1; i <= m; i++) {
			scanf("%d %d", &X, &Y);
			for (int j = 0; j <= n; j++) {//其实这里可以二分（即找到它第一个右边的线）
				if (!check(j) && check(j + 1)) {
					ans[j]++;
					break;
				}
			}
		}
		
		for (int i = 0; i <= n; i++) {
			printf("%d: %d\n", i, ans[i]);
		}
		printf("\n");
		
		scanf("%d", &n);
	}
	
	return 0; 
}