HDU 1255 覆盖的面积 （线段树+扫描线+离散化）

题意：给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积。

题解：线段树+扫描线+离散化
有一道题跟这差不多，不过比这简单，是求总的面积，只要记录覆盖1次及以上的长度len即可。

这题因为要求多次覆盖的面积，所以我们得新增一个变量len2来维护被覆盖两次及以上的区域长度，考虑4种情况：

①cnt>1，被覆盖两次及以上长度就是区间长度。

②叶子节点，为0。

③cnt=1，区间被完全覆盖了1次，而有没有被完全覆盖两次或以上则无法确定，只要加上只要加上左右子结点的len，看看左右孩子区间被覆盖了一次或以上的长度，那么叠加在双亲上就是双亲被覆盖两次或以上的长度。

④其他情况，直接用左右孩子的len2更新根结点。

len正常维护即可。

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<fstream>
#include<set>
#include<map>
#include<sstream>
#include<iomanip>
#define ll long long
using namespace std;

//扫描线 从下至上
const int MAXN = 2111;
struct Seg {
    double l, r, h;  //入边、出边的左右横坐标，以及高度
    int flag;  //标记入边1 出边-1
    Seg() {}
    Seg(double a, double b, double c, int d) :l(a), r(b), h(c), flag(d) {}
    bool operator < (const Seg& cmp) const {
        return h < cmp.h;
    }
}e[MAXN];
struct node {
    int cnt;
    double len, len2;
} tree[MAXN << 2];
double X[MAXN];
void pushUp(int l, int r, int rt) {
    if (tree[rt].cnt)//当前的边被标记，就把当前的长度加上
        tree[rt].len = X[r + 1] - X[l];
    else if (l == r)//当为一个点的时候长度为0
        tree[rt].len = 0;
    else//其他情况把左右两个区间的值加上
        tree[rt].len = tree[rt << 1].len + tree[rt << 1 | 1].len;

    if (tree[rt].cnt > 1)
        tree[rt].len2 = X[r + 1] - X[l];
    else if (l == r)
        tree[rt].len2 = 0;
    else if (tree[rt].cnt == 1)
        tree[rt].len2 = tree[rt << 1].len + tree[rt << 1 | 1].len;
    else
        tree[rt].len2 = tree[rt << 1].len2 + tree[rt << 1 | 1].len2;
}
void update(int L, int R, int l, int r, int rt, int val) {
    if (L <= l && r <= R) {
        tree[rt].cnt += val;  //加上标记的值
        pushUp(l, r, rt);
        return;
    }
    int m = (l + r) >> 1;
    if (L <= m) update(L, R, l, m, rt << 1, val);
    if (R > m) update(L, R, m + 1, r, rt << 1 | 1, val);
    pushUp(l, r, rt);
}
int n;
double xa, ya, xb, yb;
int main() {
    int t;
    scanf("%d", &t);
    while (t--) {
        scanf("%d", &n);
        memset(tree, 0, sizeof(tree));
        int num = 0;
        for (int i = 1; i <= n; i++) {
            scanf("%lf%lf%lf%lf", &xa, &ya, &xb, &yb);
            X[num] = xa;
            e[num++] = Seg(xa, xb, ya, 1); //ya 入边
            X[num] = xb;
            e[num++] = Seg(xa, xb, yb, -1); //yb 出边
        }
        sort(X, X + num);
        sort(e, e + num);
        int len = unique(X, X + num) - X;
        double ans = 0;
        for (int i = 0; i < num; i++) {
            int l = lower_bound(X, X + len, e[i].l) - X;
            int r = lower_bound(X, X + len, e[i].r) - X;
            update(l, r - 1, 0, len, 1, e[i].flag);
            ans += tree[1].len2  * (e[i + 1].h - e[i].h);
        }
        printf("%.2f\n", ans);
    }
    return 0;
}