扫描线，矩形覆盖两次及以上的面积hdu1255

给一些矩形的顶点，求它们重叠了两次及以上的部分的面积。

这个题需要在扫描线求矩形面积的基础上修改一些，不是简单的求出所有矩形的面积，再减去实际的面积。因为由的部分是覆盖了三次及以上的直接减去是错误的。

需要额外增加一个数组用来保存覆盖两次及以上的面积的长，维护这个数组的话只需要每次更新的时候计算一下。

其他的地方和简单的扫描线求矩形面积并没有什么区别。

#pragma warning(disable:4996)
#include<iostream>
#include<cstring>
#include<cstdio>
#include<set>
#include<algorithm>
#include<climits>
using namespace std;
typedef long long ll;
const int maxn = 20005;
struct node
{
	double l, r;
	double h;
	ll val;
	bool operator<(const node& x)const
	{
		return h < x.h;
	}
};
node line[maxn];
double sum[maxn], x[maxn], ss[maxn];
//ss代表被覆盖两次及以上的长度，sum代表覆盖一次及以上的长度
ll mark[maxn];
void pushup(ll l, ll r, ll k)
{
	//先计算sum
	if (mark[k])
	{
		sum[k] = x[r + 1] - x[l];
	}
	else if (l == r)sum[k] = 0;
	else sum[k] = sum[2 * k] + sum[2 * k + 1];
	//大于1直接计算出长度
	if (mark[k] > 1)
	{
		ss[k] = x[r + 1] - x[l];
	}
	else if (l == r)ss[k] = 0;
	//等于1的时候，相当于原来被覆盖一次的部分被覆盖了两次
	else if (mark[k] == 1)ss[k] = sum[2 * k] + sum[2 * k + 1];
	else ss[k] = ss[2 * k] + ss[2 * k + 1];
}
void update(ll k, ll L, ll R, ll val, ll l, ll r)
{
	if (l >= L && r <= R)
	{
		mark[k] += val;
		pushup(l, r, k);
		return;
	}
	ll mid = (l + r) / 2;
	if (L <= mid)update(2 * k, L, R, val, l, mid);
	if (R > mid)update(2 * k + 1, L, R, val, mid + 1, r);
	pushup(l, r, k);
}
int main()
{
	ll n, i, k, k1,T;
	double x1, x2, y1, y2, ans,tol;
	scanf("%lld", &T);
	while (T--)
	{
		scanf("%lld", &n);
		memset(sum, 0, sizeof(sum));
		memset(mark, 0, sizeof(mark));
		memset(ss, 0, sizeof(ss));
		k = 0; k1 = 0; ans = 0; tol = 0;
		for (i = 1; i <= n; i++)
		{
			scanf("%lf%lf%lf%lf", &x1, &y1, &x2, &y2);
			tol += (y2 - y1) * (x2 - x1);
			line[++k].l = x1; line[k].r = x2; line[k].val = 1;
			line[k].h = y1; line[++k].l = x1;
			line[k].r = x2; line[k].h = y2; line[k].val = -1;
			x[++k1] = x1; x[++k1] = x2;
		}
		//排序去重
		sort(line + 1, line + 1 + k);
		sort(x + 1, x + 1 + k1);
		k1 = unique(x + 1, x + 1 + k1) - x - 1;
		ll L, R;
		for (i = 1; i <k; i++)
		{
			L = lower_bound(x + 1, x + 1 + k1, line[i].l) - x;
			R = lower_bound(x + 1, x + 1 + k1, line[i].r) - x;
			R--;
			//update中应该是k1-1，因为只有k1-1个区间
			update(1, L, R, line[i].val, 1, k1 - 1);
			ans += (line[i + 1].h - line[i].h) * ss[1];//算出被覆盖两次及以上的面积
		}
		printf("%.2lf\n", ans);
	}
	return 0;
}