[HDOJ1255] 覆盖的面积（线段树，扫描线，矩形面积交）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1255

题意：求矩形面积交

上信安的时候撕烤了很久如何求交，发现只需要在线段树里额外维护一个tlen，代表覆盖两次以上的长度，在pushup的时候分情况更新就行了。

一直调不出来，最后发现是pushup的时候犯傻B了：当前节点代表的线段覆盖过一次的时候，那么覆盖两次的长度应该从左右儿子的覆盖1次的线段更新上来，而不是覆盖2次的（从覆盖2次的更新上来就是覆盖三次的啦喂）。

#include <bits/stdc++.h>
using namespace std;

#define lrt rt << 1
#define rrt rt << 1 | 1
typedef struct Seg {
    double len, tlen;
    int sign;
}Seg;
typedef struct Event {
    double l, r, h;
    int sign;
}Event;
const int maxn = 2100;
int n, hcnt;
vector<Event> event;
double h[maxn];
Seg seg[maxn<<2];

bool cmp(Event a, Event b) {
    if(a.h != b.h) return a.h < b.h;
    return a.sign > b.sign;
}

int id(double x) {
    return lower_bound(h, h+hcnt, x) - h + 1;
}

void build(int l, int r, int rt) {
    seg[rt].len = seg[rt].tlen = .0;
    seg[rt].sign = 0;
    if(l == r) return;
    int mid = (l + r) >> 1;
    build(l, mid, lrt);
    build(mid+1, r, rrt);
}

void pushup(int l, int r, int rt) {
    if(seg[rt].sign) {
        seg[rt].len = h[r] - h[l-1];
        if(seg[rt].sign > 1) seg[rt].tlen = seg[rt].len;
        else seg[rt].tlen = seg[lrt].len + seg[rrt].len;
    }
    else {
        if(l == r) seg[rt].len = seg[rt].tlen = 0;
        else {
            seg[rt].len = seg[lrt].len + seg[rrt].len;
            seg[rt].tlen = seg[lrt].tlen + seg[rrt].tlen;
        }
    }
}

void update(int L, int R, int sign, int l, int r, int rt) {
    if(L <= l && r <= R) {
        seg[rt].sign += sign;
        pushup(l, r, rt);
        return;
    }
    int mid = (l + r) >> 1;
    if(L <= mid) update(L, R, sign, l, mid, lrt);
    if(mid < R) update(L, R, sign, mid+1, r, rrt);
    pushup(l, r, rt);
}

int main() {
    // freopen("in", "r", stdin);
    int T;
    double ax, ay, bx, by;
    scanf("%d", &T);
    while(T--) {
        event.clear(); hcnt = 0;
        scanf("%d", &n);
        for(int i = 0; i < n; i++) {
            scanf("%lf%lf%lf%lf",&ax,&ay,&bx,&by);
            event.push_back(Event{ax,bx,ay,1});
            event.push_back(Event{ax,bx,by,-1});
            h[hcnt++] = ax, h[hcnt++] = bx;
        }
        sort(event.begin(), event.end(), cmp);
        sort(h, h+hcnt); hcnt = unique(h, h+hcnt) - h;
        build(1, hcnt, 1);
        double ret = .0;
        for(int i = 0; i < event.size(); i++) {
            int l = id(event[i].l);
            int r = id(event[i].r) - 1;
            int sign = event[i].sign;
            update(l, r, sign, 1, hcnt, 1);
            ret += (event[i+1].h - event[i].h) * seg[1].tlen;
        }
        printf("%.2f
", ret);
    }
    return 0;
}