CDOJ 1058 秋实大哥与家 线段树+扫描线

有一个W*H的矩阵 ，里面有n个家具，相互之间互不重叠。现在有一个1*M的家具，想问有多少种放法（横放，竖放都可以）

对于每个矩形（家具），我们把它横向向x轴正方向延长M-1个单位，然后，所有剩下的没有被覆盖的点就对应一个可以放新家具的位置，有多少个点就有多少种横放的放法。所以问题就变成了，总面积-矩形面积并的问题了。对于竖放，就是纵向延长M-1个单位。

所以重点就是求矩形面积并，就是扫描线的思想。

吭哧吭哧地写了好久，也没有感觉理解的很清晰。网上的扫描线讲解都搞着一堆什么离散化之类的东西，看着也有点头疼。

有一篇比较朴素的代码 点击打开链接 ，里面的最后一个题，贴的代码比较朴素，感觉也方便理解一点，至于思想，网上有很多。

重点就是它是一个区间覆盖问题，所以update和push_down会有些不同，主要不同在push_down上。有时间还是要继续理解。

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define maxn 100003
#define lid (id<<1)
#define rid (id<<1|1)
int W, H, n, M;
struct segment
{
	int l, r;
	long long lazy, sum;
}segtree[maxn * 8];
struct Edgex
{
	int x, l, r, val;
	void make(int x, int l, int r, int val)
	{
		this->x = x;
		this->l = l;
		this->r = r;
		this->val = val;
	}
}edgex[maxn * 2];
bool cmpx(Edgex a, Edgex b)
{
	return a.x < b.x;
}
struct Edgey
{
	int y, l, r, val;
	void make(int y, int l, int r, int val)
	{
		this->y = y;
		this->l = l;
		this->r = r;
		this->val = val;
	}
}edgey[maxn * 2];
bool cmpy(Edgey a, Edgey b)
{
	return a.y < b.y;
}
void bulid(int id, int l, int r)
{
	segtree[id].l = l;
	segtree[id].r = r;
	segtree[id].sum = 0;
	segtree[id].lazy = 0;
	if (l == r)
		return;
	int mid = (l + r) >> 1;
	bulid(lid, l, mid);
	bulid(rid, mid + 1, r);

}
void push_down(int id)
{
	if (segtree[id].lazy > 0)
		segtree[id].sum = segtree[id].r - segtree[id].l + 1;
	else if (segtree[id].l == segtree[id].r)
		segtree[id].sum = 0;
	else
		segtree[id].sum = segtree[lid].sum + segtree[rid].sum;
}
void update(int id, int l, int r, int v)
{
	if (l == segtree[id].l&&segtree[id].r == r)
	{
		segtree[id].lazy += v;
		push_down(id);
		return;
	}
	int mid = (segtree[id].l + segtree[id].r) >> 1;
	if (r <= mid) update(lid, l, r, v);
	else if (l > mid) update(rid, l, r, v);
	else
	{
		update(lid, l, mid, v);
		update(rid, mid + 1, r, v);
	}
	push_down(id);

}
/*long long query(int id, int l, int r)
{
	if (l == segtree[id].l&&r == segtree[id].r)
		return segtree[id].sum;
	push_down(id);
	int mid = (segtree[id].l + segtree[id].r) >> 1;
	if (r <= mid)
		return query(lid, l, r);
	else if (l > mid)
		return query(rid, l, r);
	else
		return query(lid, l, mid) + query(rid, mid + 1, r);
}*/

int main()
{
	//freopen("input.txt", "r", stdin);
	scanf("%d%d%d%d", &W, &H, &n, &M);
	int x1, y1, x2, y2;
	for (int i = 0; i < n; ++i)
	{
		scanf("%d%d%d%d", &x1, &y1, &x2, &y2);
		edgex[i * 2].make(x1, y1, y2, 1);
		edgex[i * 2 + 1].make(x2 + M, y1, y2, -1);
		edgey[i * 2].make(y1, x1, x2, 1);
		edgey[i * 2 + 1].make(y2 + M, x1, x2, -1);
	}
	sort(edgex, edgex + 2 * n, cmpx);
	sort(edgey, edgey + 2 * n, cmpy);
	bulid(1, 1, 2 * maxn);
	long long ans = 0;
	for (int i = 1, t = 0; i <= W; ++i)
	{
		while (edgex[t].x == i&&t < 2 * n)
		{
			update(1, edgex[t].l, edgex[t].r, edgex[t].val);
			++t;
		}
		if (i >= M)
			ans += H - segtree[1].sum;
	}
	//printf("%lld\n", ans);
	if (M != 1)
	{
		bulid(1, 1, 2 * maxn);
		for (int i = 1, t = 0; i <= H; ++i)
		{
			while (edgey[t].y == i&&t < 2 * n)
			{
				update(1, edgey[t].l, edgey[t].r, edgey[t].val);
				++t;
			}
			if (i >= M)
				ans += W - segtree[1].sum;
		}
	}
	printf("%lld\n", ans);
	//system("pause");
	//while (1);
	return 0;
}