【PR #5 C】和平共处（整体二分）

和平共处

题目链接：PR #5 C

题目大意

有 n 个黑点 m 个白点，黑点一开始都在，白点按一定顺序加入。
问每次加入之后，你要选一些点删去（只是假设删去，并没有真正删去），使得不存在一个黑点在白点的左下方。
问你每次加入之后，最少要删多少个点。

思路

首先考虑白点都加入了要怎么搞，就是要删哪些，或者留那些。
那比如选了一个黑点，它右上方就不会有白点，那它右上方的黑点就可以全部保留。
那我们不难想象出最后的样子，一个阶梯状作为分解，而且是左上右下的那种。
上面包括阶梯线上的黑色+下面不包括阶梯的白色就是我们能保留的。

那一个显然的事情是随着白色点的加入，这个阶梯只会往右上移动。
那就是单调的，我们可以用整体二分来弄（一个分界线确定之后，它两边的点只会分别在 \(<mid\) 和 \(>mid\) 有是否贡献的不确定性）

那我们就继续考虑要怎么求一次的分界线。
考虑从左往右看，每个黑色的地方看是否要从当前的高度降到他这个高度。
那降不降就看那个缩小的区间中，黑白点的数量。
考虑黑白点进行一个匹配，黑色跟右上方的点匹配， 理解为作用抵消，然后没有匹配的话肯定是黑色多，折。
那如果两个都在先下面（注意上面是白色所以上面不能压线），那我们肯定优先不要在这里往下（因为我们是尽量不往下的）。
那如果黑色的在下面，白色的在上面，那我们就要黑色往下了，不然就不优了。

不过会发现匹配不可以随便匹配。
思考会发现应该从后往前贪心的匹配，黑点每次选还没有匹配的右边的第一个比它高的白点跟它匹配。

代码

#include<set>
#include<cstdio>
#include<vector>
#include<algorithm>
#define INF 0x3f3f3f3f3f3f3f3f

using namespace std;

const int N = 2e5 + 100;
struct dian {
	int x, y, op;
};
int n, m, ans[N], nxt[N];
vector <dian> p;

bool cmp(dian x, dian y) {
	if (x.x != y.x) return x.x < y.x;
	if (x.y != y.y) return x.y < y.y;
	return x.op < y.op;
}

void slove(int l, int r, vector <dian> &p, int an) {
	if (l > r) return ;
	int mid = (l + r) >> 1;
	set <pair<int, int> > s;
	for (int i = 0; i < p.size(); i++) nxt[i] = -1;
	for (int i = p.size() - 1; i >= 0; i--) {
		if (p[i].op) {
			if (p[i].op <= mid) s.insert(make_pair(p[i].y, i));
		}
		else {
			set <pair<int, int> > ::iterator it = s.lower_bound(make_pair(p[i].y, -1));
			if (it == s.end()) continue;
			nxt[(*it).second] = i; nxt[i] = (*it).second;
			s.erase(it);
		}
	}
	vector <dian> pl, pr; int Up = INF, lcnt = 0, rcnt = 0;
	for (int i = 0; i < p.size(); i++) {
		if (p[i].op) {
			if (p[i].op <= mid) {
				if (p[i].y < Up) pl.push_back(p[i]), rcnt++;
					else pr.push_back(p[i]);
			}
			else pr.push_back(p[i]);
		}
		else {
			if (nxt[i] == -1 || p[nxt[i]].y >= Up) Up = min(Up, p[i].y);
			if (p[i].y >= Up) pr.push_back(p[i]), lcnt++;
				else pl.push_back(p[i]);
		}
	}
	ans[mid] = n + mid - (an + lcnt + rcnt);
	slove(l, mid - 1, pl, an + lcnt);
	slove(mid + 1, r, pr, an + rcnt);
}

int main() {
	scanf("%d", &n);
	for (int i = 1; i <= n; i++) {
		dian x; scanf("%d %d", &x.x, &x.y); x.op = 0;
		p.push_back(x);
	}
	scanf("%d", &m);
	for (int i = 1; i <= m; i++) {
		dian x; scanf("%d %d", &x.x, &x.y); x.op = i;
		p.push_back(x);
	}
	sort(p.begin(), p.end(), cmp);
	
	slove(1, m, p, 0);
	for (int i = 1; i <= m; i++) printf("%d\n", ans[i]);
	
	return 0;
}