P5490 【模板】扫描线

传送门
扫描线模板

分析

逃不掉，总该要学的。学完后，发现思想不难，早干嘛去了

大概说一下
核心要义 扫描&&线段

定义一根线，从下往上扫
因为矩形的话，肯定在竖直方向有一段是连续的（也就是有个高，已知）
如果要知道某一小段的面积，我们需要知道长度为多少？

那么我们就可以使用线段树维护这个长度
长度的定义为，只要有一个矩形贡献了一段长度，这一段的长度就已经有了，不会因为矩形个数增加
把所有长度都统计起来，最后乘以高度即可

接下来如何维护
如果一个矩形的两个x相同，说明没有面积，面积为 $0$ ，没必要统计

我们统计的是线段！！！，这里很重要，也就是说线段树的每个节点也是线段
所以，某个叶子节点即使是最短都会表示一个数学量，长度至少为 $1$
我们举例，此时有一个线段，$[l,r],(l!=r)$
那么我们需要有节点维护 $l\sim l+1$，$l+1\sim l+2$，$l+2\sim l+3$，$......$ ，$r-1\sim r$，总共$r-l$个线段
因为我们线段树维护的是线段，所以，假设我们维护的线段为$x\sim x+1$，那么就让它在线段树的$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;
}