[HDOJ 4938] Seeing People [排序+二分查找]

有两组人，一组从x轴上以相同速率沿y轴正方向移动，另一组从y轴上以另一相同速率沿x轴正方向运动。

已知每个人的起始位置，开始移动的时刻，视野范围，视野范围为和移动方向垂直的一条线段，问每个人在移动过程中能碰见多少个人？

每个人的开始移动时间和起始位置均不同。

数据范围为10^5。输入的所有的数均为正整数。

首先我们可以知道，同组人是无法在移动时互相看到的。然后我们考虑x组上的一个人看见y组上的一个人的充要条件为Px<=(Py/Vx+Tx-Ty)*Vy<=Px+Wx。

稍作转化即可得到Px-Tx*Vy<=(Py/Vx-Ty)*Vy<=Px+Wx-Tx*Vy。

可以知道Vx和Vy对于任意两对人是恒定不变的，即可以当成常数，中间部分仅与y有关，两边仅与x有关。

所以我们对y进行排序，对于任意一个x，用二分查找判断符合要求的y有多少个，即可知道x与多少人碰见。

y会与多少x碰见也是一样的，两边是对称的。

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

int n,v1,v2;
const long double eps=1e-9;

struct People {
	int i,p,t,w;
	long double v;
	People() {}
	People(long double vv) {
		v=vv;
	}
	People(int ii,int pp,int tt,int ww,int v1,int v2) {
		i=ii; p=pp; t=tt; w=ww;
		v=((long double)p/v2-t)*v1;
	}
	friend bool operator < (const People &a,const People &b) {
		return a.v<b.v;
	}
};

int ans[100000];
People a[100000];
People b[100000];
int an,bn;

int main() {
	int t,tt,i;
	scanf("%d",&t);
	for (tt=1;tt<=t;tt++) {
		scanf("%d%d%d",&n,&v1,&v2);
		an=bn=0;
		memset(ans,0,sizeof(ans));
		for (i=0;i<n;i++) {
			int k,t,p,w;
			scanf("%d%d%d%d",&k,&t,&p,&w);
			if (k==1) a[an++]=People(i,p,t,w,v1,v2);
			else b[bn++]=People(i,p,t,w,v2,v1);
		}
		sort(a,a+an);
		sort(b,b+bn);
		for (i=0;i<an;i++) {
			People upper=People(a[i].p+a[i].w-(long double)a[i].t*v2+eps);
			People lower=People(a[i].p-(long double)a[i].t*v2-eps);
			ans[a[i].i]=upper_bound(b,b+bn,upper)-lower_bound(b,b+bn,lower);
		}
		for (i=0;i<bn;i++) {
			People upper=People(b[i].p+b[i].w-(long double)b[i].t*v1+eps);
			People lower=People(b[i].p-(long double)b[i].t*v1-eps);
			ans[b[i].i]=upper_bound(a,a+an,upper)-lower_bound(a,a+an,lower);
		}
		printf("Case #%d:\n",tt);
		for (i=0;i<n;i++) {
			printf("%d\n",ans[i]);
		}
	}
	return 0;
}