[ABC211F] Rectilinear Polygons 题解

[ABC211F] Rectilinear Polygons 题解

思路什么的上一篇题解已经写的非常明白了,这里只是提供一个补充 & 另一个实现的方法。

思路解析

先说结论:扫描线。顾名思义,扫描线的本质就是用一条线沿着 xy 轴扫过去,每碰到一条边就记录一下加边后是面积是增加还是减少,然后用树状数组或线段树统计。由于我们并不需要知道整张图所有的位置的值,我们就只需要求出来扫描线所在的那一行 or 列每一个点的值即可。

这里就不过多解释,如果没有学过扫描线建议去做一下模板题,推荐看第一篇题解。

接下来就是本题的重点,该如何把不规则图形存下来。如下图,我们先把图上的竖边找出来,同时给点的顺序标一个号。

可见四条边中左边三条的权值都是增加的,只有右边一条是减少的。我们寻找边上的两点和权值的关系,可以发现,由于是点顺序输入的,所以每一条正边权的边的两点中编号较小的点 y 值也是更小的,而每一条边的两点中编号较小的点的编号一定是奇数。因此我们可以给每一条边都判断两点编号的大小来决定是用正边权还是负边权。

实现

我们根据以上所说的方式存好每一条边,然后将边按 x 的值排序,满足扫描线的性质。对于询问我们将询问离线下来,同样按 x 排序,这样我们在进行加边之前判断当前扫描线是否已经处理了所有 x 小于当前询问的 x 的边,若已经处理完就更新答案。

接下来有几点注意事项。

  • 题目中的输入点的顺序不一定满足我们想要的输入顺序,所以我们需要自己找到一个在 x 最小的情况下 y 也最小的点做编号为 1 的点。
  • 由于我们只需要知道一个点被覆盖的次数,所以可以将模板中的线段树换成树状数组区修单查。
  • 由于我们需要先把当前 x 轴上的边加完再统计询问,所以我们要在所有边的最后加上一个 x=inf,val=0 的边,这样就能防止最后有几个点没有处理到。

code

#include<bits/stdc++.h>
using namespace std;
const int N = 1e5 + 10, M = 1e5;
int n, m, xx[N], yy[N], tv[N], q, c[N], ans[N];
struct side {
	int x, ya, yb, val;
};
struct point {
	int x, ya, yb; 
};
struct query {
	int x, y, p;
};
void add(int x, int y) {
	for(; x <= M; x += (x & -x)) {
		c[x] += y;
	}
}
int ask(int x) {
	int sum = 0;
	for(; x > 0; x -= (x & -x)) {
		sum += c[x];
	}
	return sum;
}
int main() {
	cin >> n;
	vector<side> v;
	for(int i = 1; i <= n; i++) {
		cin >> m;
		for(int i = 0; i < m; i++) {
			cin >> xx[i] >> yy[i];
			xx[i]++; yy[i]++;
		}
		int p = 0;
		for(int i = 0; i < m; i++) {
			if(xx[i] == xx[p] && yy[i] < yy[p]) p = i;	//找到起始点,注意 y 轴也要判
			else if(xx[i] < xx[p]) p = i;
		}
		int val = 1;
		for(int i = 0; i < m; i++) {
			int w = (i + p) % m;
			tv[w] = val;	//给点加权
			val = -val;
		}
		for(int i = 0; i < m; i += 2) {	//判断
			if(yy[i] < yy[i + 1]) v.push_back({xx[i], yy[i], yy[i + 1], tv[i]});
			else v.push_back({xx[i], yy[i + 1], yy[i], tv[i + 1]});
		}
	}
	sort(v.begin(), v.end(), [](side x, side y) {	//按 x 排序
		if(x.x != y.x) return x.x < y.x;
		else return x.ya < y.ya;
	});
	v.push_back({M + 1, 1, 1, 0});	//避免结尾有点没判
	vector<query> que;
	cin >> q;
	for(int i = 1; i <= q; i++) {
		cin >> xx[i] >> yy[i];
		xx[i]++; yy[i]++;
		que.push_back({xx[i], yy[i], i});	//离线
	}
	sort(que.begin(), que.end(), [](query x, query y) {	//排序
		if(x.x != y.x) return x.x < y.x;
		else return x.y < y.y;
	});
	int pos = 0;
	for(auto it : v) {
		int x = it.x, ya = it.ya, yb = it.yb, val = it.val;
		for(; que[pos].x < x && pos < (int)que.size(); pos++) {
			ans[que[pos].p] = ask(que[pos].y);	//若已经操作完则更新
		}
		add(ya, val); add(yb, -val);	//加权
	}
	for(int i = 1; i <= q; i++) cout << ans[i] << '\n';
	return 0;
}
posted @   2020luke  阅读(12)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 无需6万激活码!GitHub神秘组织3小时极速复刻Manus,手把手教你使用OpenManus搭建本
· Manus爆火,是硬核还是营销?
· 终于写完轮子一部分:tcp代理 了,记录一下
· 别再用vector<bool>了!Google高级工程师:这可能是STL最大的设计失误
· 单元测试从入门到精通
点击右上角即可分享
微信分享提示