传送门
扫描线模板
分析
逃不掉,总该要学的。学完后,发现思想不难,早干嘛去了
大概说一下
核心要义 扫描&&线段
定义一根线,从下往上扫
因为矩形的话,肯定在竖直方向有一段是连续的(也就是有个高,已知)
如果要知道某一小段的面积,我们需要知道长度为多少?
那么我们就可以使用线段树维护这个长度
长度的定义为,只要有一个矩形贡献了一段长度,这一段的长度就已经有了,不会因为矩形个数增加
把所有长度都统计起来,最后乘以高度即可
接下来如何维护
如果一个矩形的两个x相同,说明没有面积,面积为
0
0
0 ,没必要统计
我们统计的是线段!!!,这里很重要,也就是说线段树的每个节点也是线段
所以,某个叶子节点即使是最短都会表示一个数学量,长度至少为
1
1
1
我们举例,此时有一个线段,
[
l
,
r
]
,
(
l
!
=
r
)
[l,r],(l!=r)
[l,r],(l!=r)
那么我们需要有节点维护
l
∼
l
+
1
l \sim l+1
l∼l+1,
l
+
1
∼
l
+
2
l+1 \sim l+2
l+1∼l+2,
l
+
2
∼
l
+
3
l+2 \sim l+3
l+2∼l+3,
.
.
.
.
.
.
... \space ...
... ... ,
r
−
1
∼
r
r-1 \sim r
r−1∼r,总共
r
−
l
r-l
r−l个线段
因为我们线段树维护的是线段,所以,假设我们维护的线段为
x
∼
x
+
1
x \sim x+1
x∼x+1,那么就让它在线段树的
x
x
x位置上维护
那么所有线段都能对应线段树的一段区间,此时线段树维护已经解决了
OK,如果上面没问题了,这里又再说一下他们常用的玄学优化(不用下传标记)
每次我们加减线段,是不是都是一整个区间 加 或者 减 啊
这个意思是说,只要我们碰到要减的时候,肯定也是找到这个区间,不会触碰到这个区间往下的区间
所以加的时候加到包含的区间为止,减的时候也是减到包含的区间为止,任何一个线段的参与,都不会涉及到深层次的叶子节点
所以这里没必要下传标记(线段树下传标记的思想,感性理解一下)
代码
//P5490
/*
@Author: YooQ
*/
#include <bits/stdc++.h>
using namespace std;
#define sc scanf
#define pr printf
#define ll long long
#define int long long
#define FILE_OUT freopen("out", "w", stdout);
#define FILE_IN freopen("in", "r", stdin);
#define debug(x) cout << #x << ": " << x << "\n";
#define max3(a, b, c) max(a, max(b, c))
#define min3(a, b, c) min(a, min(b, c))
#define AC 0
#define WA 1
#define INF 0x3f3f3f3f
const ll MAX_N = 1e6+5;
const ll MOD = 1e9+7;
int N, M, K;
struct Line {
int l, r, h;
int opt;
bool operator < (const Line& B) const {
return h < B.h;
}
}line[MAX_N<<1];
int uniarr[MAX_N];
int unicnt = 0;
struct Tr {
int l, r, sum, len;
}tr[MAX_N<<2];
void push_up(int rt) {
int l = tr[rt].l;
int r = tr[rt].r;
if (tr[rt].sum) {
tr[rt].len = uniarr[r+1] - uniarr[l];
} else {
tr[rt].len = tr[rt<<1].len + tr[rt<<1|1].len;
}
}
void build(int rt, int l, int r) {
tr[rt] = {l, r, 0, 0};
if (l == r) return;
int mid = l + ((r-l)>>1);
build(rt<<1, l, mid);
build(rt<<1|1, mid+1, r);
}
void update(int rt, int l, int r, int x, int y, int k) {
if (x <= l && r <= y) {
tr[rt].sum += k;
push_up(rt);
return;
}
int mid = l + ((r-l)>>1);
if (x <= mid) update(rt<<1, l, mid, x, y, k);
if (y > mid) update(rt<<1|1, mid+1, r, x, y, k);
push_up(rt);
}
void solve() {
sc("%lld", &N);
int x, y, a, b;
for (int i = 1; i <= N; ++i) {
sc("%lld%lld%lld%lld", &a, &b, &x, &y);
uniarr[++unicnt] = a;
uniarr[++unicnt] = x;
line[2*i-1] = {a, x, b, 1};
line[2*i] = {a, x, y, -1};
}
N <<= 1;
sort(line+1, line+1+N);
sort(uniarr+1, uniarr+1+unicnt);
unicnt = unique(uniarr+1, uniarr+1+unicnt) - uniarr - 1;
build(1, 1, unicnt);
int ans = 0;
for (int i = 1; i <= N; ++i) {
line[i].l = lower_bound(uniarr+1, uniarr+1+unicnt, line[i].l) - uniarr;
line[i].r = lower_bound(uniarr+1, uniarr+1+unicnt, line[i].r) - uniarr;
update(1, 1, unicnt, line[i].l, line[i].r - 1, line[i].opt);
ans += tr[1].len * (line[i + 1].h - line[i].h);
}
pr("%lld\n", ans);
}
signed main() {
#ifndef ONLINE_JUDGE
//FILE_IN
FILE_OUT
#endif
int T = 1;//cin >> T;
while (T--) solve();
return AC;
}