HDU 5714 树状数组 + 离散化

题意：

Problem Description

小明在旅游的路上看到了一条美丽的河，河上有许多船只，有的船只向左航行，有的船只向右航行。小明希望拍下这一美丽的风景，并且把尽可能多的船只都完整地拍到一张照片中。

小明位于河的边上，并且可以在河边的任意位置进行拍照，照相机的视野恰好为90度角，只能以垂直于河边的方向进行拍照。河上的船只全都可看作是平行于河边的一条线段，跟河边的距离各不相同，有的正在向左移动，有的正在向右移动，但移动速度恰好都是一样的。小明可以等待恰当的时间让尽量多的船只都走进照相机的视野里，你不需要考虑船只之间会互相遮挡视野的情况。

Input

第一行为$T$，表示输入数据组数。

下面$T$组数据，对于每组数据：

第一行是一个数n(1≤n≤104)n(1\leq n\leq 10^{4})n(1≤n≤10​4​​)，表示船只的数量。

接下来$n$行，每行四个整数x,y,z,d(−106≤x<y≤106,1≤z≤104)x,y,z,d(-10^{6}\leq x < y \leq 10^{6},1\leq z\leq 10^{4})x,y,z,d(−10​6​​≤x<y≤10​6​​,1≤z≤10​4​​)，表示船只的左端点位置、右端点位置、距离河边的距离，以及航行的方向。$d$为$-1$表示向左航行，$1$表示向右航行。

Output

对第$i$组数据，输出

Case #i:

然后输出一行，仅包含一个整数，表示最多可以拍到多少完整的船只。

Sample Input

Copy

321 3 1 12 4 1 -121 3 1 -12 4 1 111 4 1 1

Sample Output

Copy

Case #1:2Case #2:1Case #3:0

思路：

首先假定向左方向的船都是静止的，那么船要么向右移动要么静止。

对于区间是[L,R]的船，我们根据视角计算出在岸上可以观察这艘船的位置区间[l,r]，因为拍摄是可以等待船移动的，因为所有船的速度是一样的，如果我们选定了某个位置，那么这个位置能看到一些静止的船，除此之外，还可以等待该位置左边的且向右行驶的船移动到该位置，所以其实就是在当前位置pos的左边找一个位置pos'，使得pos'位置上能看到的向右行驶的船数+pos位置上静止的船数最大即可。

这样就可以利用树状数组来维护，每次区间左端点+1，右端点后一位-1，然后求得的sum(i)就是当前位置能看到的船数，向左和向右的船分别用一个BIT来维护。

另外需要注意到的一点是x,y的区间范围还是很大，需要离散化。

代码：

#include <bits/stdc++.h>using namespace std;const int MAXN = 22222;struct BIT {int n;int C[MAXN];void init(int n) {this->n = n;memset(C, 0, sizeof(C));}int lowbit(int x) {return x & -x;}void add(int x, int y) {while (x <= n) {C[x] += y;x += lowbit(x);}}int sum(int x) {int res = 0;while (x) {res += C[x];x -= lowbit(x);}return res;}} lbit, rbit;struct Boat {int l, r, dir;};vector <Boat> boat;int discre(int n) {int Max = 0;vector <int> vec;for (int i = 0; i < (int)boat.size(); i++) {vec.push_back(boat[i].l);vec.push_back(boat[i].r);}sort(vec.begin(), vec.end());vec.erase(unique(vec.begin(), vec.end()), vec.end());for (int i = 0; i < (int)boat.size(); i++) {boat[i].l = lower_bound(vec.begin(), vec.end(), boat[i].l) - vec.begin() + 1;boat[i].r = lower_bound(vec.begin(), vec.end(), boat[i].r) - vec.begin() + 1;Max = max(Max, max(boat[i].l, boat[i].r));}return Max;}int main() {// freopen("in.txt", "r", stdin);int T, cs = 0;scanf("%d", &T);while (T--) {int n;scanf("%d", &n);boat.clear();for (int i = 1; i <= n; i++) {int l, r, d, dir;scanf("%d%d%d%d", &l, &r, &d, &dir);double len = (d - (r - l) / 2.0), mid = (l + r) / 2.0;if (len >= 0)boat.push_back((Boat) {mid - len, mid + len, dir});}int m = discre(n);lbit.init(m); rbit.init(m); for (int i = 0; i < boat.size(); i++) {if (boat[i].dir == -1) {lbit.add(boat[i].l, 1);lbit.add(boat[i].r + 1, -1);} else {rbit.add(boat[i].l, 1);rbit.add(boat[i].r + 1, -1);} }int maxright = 0, ans = 0;for (int i = 1; i <= m; i++) {int leftcnt = lbit.sum(i);maxright = max(maxright, rbit.sum(i));ans = max(ans, leftcnt + maxright);}printf("Case #%d:\n%d\n", ++cs, ans);}return 0;}