HDOJ 1542 线段树 矩形面积并

链接：

http://acm.hdu.edu.cn/showproblem.php?pid=1542

题意：

给你n个矩形，问你这个矩形覆盖的面积

题解：

就是线段树+扫描线+坐标离散化，详情在代码里面

代码：

#include <map>
#include <set>
#include <cmath>
#include <queue>
#include <stack>
#include <cstdio>
#include <string>
#include <vector>
#include <cstdlib>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
#define rep(i,a,n) for (int i=a;i<n;i++)
#define per(i,a,n) for (int i=n-1;i>=a;i--)
#define all(x) (x).begin(),(x).end()
#define pb push_back
#define mp make_pair
#define lson l,m,rt<<1  
#define rson m+1,r,rt<<1|1 
typedef long long ll;
typedef vector<int> VI;
typedef pair<int, int> PII;
const ll MOD = 1e9 + 7;
const int INF = 0x3f3f3f3f;
const int MAXN = 210;
// head

int cover[MAXN << 2];   //记录某个区间被覆盖次数
double x[MAXN << 1];    //对x进行离散化
double sum[MAXN << 2];  //记录某个区间的下底边总长度 

struct seg {
    double l, r, h;
    int d;
    const bool operator<(const seg &t) const {
        return h < t.h;
    }
}s[MAXN << 1];

void pushup(int l, int r, int rt) {
    if (cover[rt]) sum[rt] = x[r + 1] - x[l];　　//表示该区间至少被覆盖一次
    else if (l == r) sum[rt] = 0;         　　   //叶子结点则底边长度为0(区间内线段长度为0)   
    else sum[rt] = sum[rt << 1] + sum[rt << 1 | 1];
}

void update(int L, int R, int d, int l, int r, int rt) {
    if (L <= l && r <= R) {  //更新该区间下底边总长以及覆盖次数
        cover[rt] += d;      //更新覆盖次数，因为先加后减，所以不会是负数 
        pushup(l, r, rt);    //更新底边长   
        return;
    }
    int m = (l + r) >> 1;
    if (L <= m) update(L, R, d, lson);
    if (R > m) update(L, R, d, rson);
    pushup(l, r, rt);
}

int main() {
    int n, cas = 1;
    while (cin >> n, n) {
        int k = 0;
        rep(i, 0, n) {
            double x1, y1, x2, y2;
            scanf("%lf%lf%lf%lf", &x1, &y1, &x2, &y2);
            x[k] = x1;
            s[k++] = seg{ x1, x2, y1, 1 };
            x[k] = x2;
            s[k++] = seg{ x1, x2, y2, -1 };
        }
        sort(x, x + k);
        sort(s, s + k);
        int m = 1;
        rep(i, 1, k) if (x[i] != x[i - 1]) x[m++] = x[i];    //去重
        double ans = 0;
        rep(i, 0, k) {
            int L = lower_bound(x, x + m, s[i].l) - x;
            int R = lower_bound(x, x + m, s[i].r) - x - 1;
            update(L, R, s[i].d, 0, m - 1, 1);        //更新底边长度和覆盖次数，注意区间为[0,m-1]
            ans += sum[1] * (s[i + 1].h - s[i].h);    //sum[1]存的就是总长度
        }
        printf("Test case #%d
Total explored area: %.2f

", cas++, ans);
    }
    return 0;
}