【HDU】1225 覆盖的面积（线段树 扫描线）

思路：

做到这一题的时候才真正的理解了上一题矩形并面积的思路

col[rt]表示这个区间被覆盖了几次，sum[rt]表示这个区间有效的长度。

这边需要引入第三个变量sum2[rt]表示这个区间内被覆盖两次以上的有效部分大小,sum[rt]就变成了覆盖一次的长度。

这边需要思考的便是uprt函数，即更新父节点的函数。

col[rt]假如不等于0，则表示的就是有一个矩形刚好等于或者包含这段区间。

所以更新父节点的函数应该是

void uprt(int l,int r,int rt)
{
	//col==1的时候表示这个区间正好被覆盖
	//所以这个区间内合法的覆盖区间为，里面被覆盖一次的，还有被覆盖两次即以上的。
	//因为里面小区间被覆盖了一次，再加上外面大区间被覆盖了一次，所以原来的小区间就被覆盖了两次
	//所以sum2[rt]=sum[ls]+sum[rs]+sum2[ls]+sum2[rs];没能够充分理解col就没办法想出来
	if (col[rt] == 1)
	{
		sum2[rt] = sum[ls] + sum[rs] + sum2[ls] + sum2[rs];
		sum[rt] = posx[r + 1] - posx[l] - sum2[rt];
		return;
	}
	if (col[rt] == 2)
	{
		sum2[rt] = posx[r + 1] - posx[l];
		sum[rt] = 0;
		return;
	}
	if (col[rt] == 0)
	{
		sum2[rt] = sum2[ls] + sum2[rs];
		sum[rt] = sum[ls] + sum[rs];
		return;
	}
}

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAX 3005
#define ls rt<<1
#define rs ls|1
#define m (l+r)>>1

double sum[MAX << 2];
double sum2[MAX << 2];
int col[MAX << 2];
double posx[MAX ];
struct pos
{
	double l, r, h;
	int s;
	pos(){}
	pos(double _l, double _r, double _h, int _s)
	{
		l = _l;
		r = _r;
		h = _h;
		s = _s;
	}
	bool operator<(pos b) const
	{
		return h < b.h;
	}
}p[MAX];

void uprt(int l,int r,int rt)
{
	//col==1的时候表示这个区间正好被覆盖
	//所以这个区间内合法的覆盖区间为，里面被覆盖一次的，还有被覆盖两次即以上的。
	//因为里面小区间被覆盖了一次，再加上外面大区间被覆盖了一次，所以原来的小区间就被覆盖了两次
	//所以sum2[rt]=sum[ls]+sum[rs]+sum2[ls]+sum2[rs];没能够充分理解col就没办法想出来
	if (col[rt] == 1)
	{
		sum2[rt] = sum[ls] + sum[rs] + sum2[ls] + sum2[rs];
		sum[rt] = posx[r + 1] - posx[l] - sum2[rt];
		return;
	}
	if (col[rt] == 2)
	{
		sum2[rt] = posx[r + 1] - posx[l];
		sum[rt] = 0;
		return;
	}
	if (col[rt] == 0)
	{
		sum2[rt] = sum2[ls] + sum2[rs];
		sum[rt] = sum[ls] + sum[rs];
		return;
	}
}

void updata(int L, int R, int c, int l, int r, int rt)
{
	if (L <= l&&r <= R)
	{
		col[rt] += c;
		uprt(l, r, rt);
		return;
	}
	int mid = m;
	if (L <= mid)
		updata(L, R, c, l, mid, ls);
	if (mid < R)
		updata(L, R, c, mid + 1, r, rs);
	uprt(l, r, rt);
}

int main()
{
	int t;
	cin >> t;
	while (t--)
	{
		int n;
		int cnt = 0;
		scanf("%d", &n);
		double x1, x2, y1, y2;
		for (int i = 0; i < n; i++)
		{
			scanf("%lf%lf%lf%lf", &x1, &y1, &x2,&y2);
			p[cnt] = pos(x1, x2, y1, 1);
			posx[cnt++] = x1;
			p[cnt] = pos(x1, x2, y2, -1);
			posx[cnt++] = x2;
		}
		sort(posx, posx + cnt);
		sort(p, p + cnt);
		int cnt2 = unique(posx,posx + cnt) - posx;
		memset(col, 0, sizeof(col));//之前是用build函数来重置，但一直wa
		memset(sum, 0, sizeof(sum));
		memset(sum2, 0, sizeof(sum2));
		double ans = 0;
		for (int i = 0; i < cnt-1; i++)
		{
			int curl = lower_bound(posx, posx+cnt2, p[i].l)-posx;
			int curr = lower_bound(posx, posx+cnt2, p[i].r)-posx-1;
			updata(curl, curr, p[i].s, 0, cnt2-1, 1);
			ans += (sum2[1]) * (p[i + 1].h - p[i].h);
		}
		printf("%.2lf\n", ans);
	}
	return 0;
}