HDU1255覆盖的面积(线段树+离散化+扫描线)

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

线段树求面积并升级版

题目大意：

被矩形覆盖过两次的地方的面积

来自大佬的分析：http://www.cnblogs.com/scau20110726/archive/2013/04/14/3020998.html

这题的代码在一般的线段树求面积并的基础上进行了修改，但是所用的思想是一样的，所以不难理解

回忆一下一般的求矩形覆盖面积，线段树节点里面有一个重要的变量，cnt。这个变量表示了该节点表示的区间被完全覆盖，如果cnt=0，说明没有被完全覆盖（但不代表没有被覆盖），要算出该节点所代表的区间被覆盖的长度，需要由它左右孩子节点被覆盖的长度相加所得。如果cnt=1，表示被完全覆盖，覆盖长度就是该区间长度。如果cnt>1说明也是被完全覆盖，不过不止覆盖了一次，在算覆盖长度的时候，和cnt=1的计算方法是一样的。注意一点，节点里还有另一个变量len，就是该区间被覆盖的长度，但是我们注意一下，这个len准确的意义应该是，被覆盖了一次或以上的长度，只是这个意义在一般的求面积问题中，不需要过分强调

 

而在这题中我们要计算被覆盖两次或以上的部分面积，我们在线段树节点中增设了一个变量，ss，其中s表示该该区间内被覆盖了1次或以上的长度，ss表示被覆盖了2次或以上的长度

我们是怎么计算最后的面积的？一样的道理，从下往上扫描矩形，每次添加一条矩形上下边，然后看看t[1].ss是多少，再乘上高度差。因为t[1]表示了总区间，而ss表示被覆盖两次或以上的长度，即计算时我们忽略掉只被覆盖一次的长度

问题的关键变为怎么计算一个节点的ss

分情况讨论

1.cnt>1 : 说明该区间被覆盖两次或以上，那么长度就可以直接计算，就是该区间的长度

剩下的情况就是cnt=1或cnt=0

2.先看叶子节点，因为是叶子没有孩子了，所以被覆盖两次货以上的长度就是0（无论cnt=1或cnt=0都是0，因为是叶子。。。）

3.不是叶子节点 ，且cnt=1.注意这里，cnt=1确切的意义是什么，应该是，可以确定，这个区间被完全覆盖了1次，而有没有被完全覆盖两次或以上则不知道无法确定，那么怎么怎么办了，只要加上t[lch].s + t[rch].s  即，看看左右孩子区间被覆盖了一次或以上的长度，那么叠加在双亲上就是双亲被覆盖两次或以上的长度

3.不是叶子节点，且cnt=0，确切的意义应该是不完全不知道被覆盖的情况（不知道有没有被覆盖，被覆盖了几次，长度是多少都不知道），这种情况，只能由其左右孩子的信息所得

t[lch].ss + t[rch].ss  , 即直接将左右孩子给覆盖了两次或以上的长度加起来，这样才能做到不重不漏

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
const int maxn = 1010;
double pos[maxn*2];
struct segment
{
    double l, r, h;
    int v;
}ss[maxn<<1];
struct Node
{
    int l, r, cnt;
    double s, ss;
}segTree[maxn<<3];
bool cmp(segment p, segment q)
{
    return p.h < q.h;
}
void build(int l, int r, int rt)
{
    segTree[rt].l = l;
    segTree[rt].r = r;
    segTree[rt].cnt = segTree[rt].s = segTree[rt].ss = 0;
    if(l == r)
        return ;
    int mid = (segTree[rt].l + segTree[rt].r) >> 1;
    build(l, mid, rt<<1);
    build(mid+1, r, (rt<<1)|1);
}
int binary(double key, int low, int high)
{
    while(low <= high)
    {
        int mid = (low + high) >> 1;
        if(pos[mid] == key)
            return mid;
        else if(pos[mid] < key)
            low = mid + 1;
        else
            high = mid - 1;
    }
}
void cal(int rt)
{
    if(segTree[rt].cnt)
        segTree[rt].s = pos[segTree[rt].r+1] - pos[segTree[rt].l];
    else if(segTree[rt].l == segTree[rt].r)
        segTree[rt].s = 0;
    else
        segTree[rt].s = segTree[rt<<1].s + segTree[(rt<<1)|1].s;

    if(segTree[rt].cnt > 1)
        segTree[rt].ss = pos[segTree[rt].r+1] - pos[segTree[rt].l];
    else if(segTree[rt].l == segTree[rt].r)
        segTree[rt].ss = 0;
    else if(segTree[rt].cnt == 1)
        segTree[rt].ss = segTree[rt<<1].s + segTree[(rt<<1)|1].s;
    else
        segTree[rt].ss = segTree[rt<<1].ss + segTree[(rt<<1)|1].ss;
}
void update(int l, int r, int val, int rt)
{
    if(segTree[rt].l == l && segTree[rt].r == r)
    {
        segTree[rt].cnt += val;
        cal(rt);
        return ;
    }
    int mid = (segTree[rt].l + segTree[rt].r) >> 1;
    if(r <= mid)
        update(l, r, val, rt<<1);
    else if(l > mid)
        update(l, r, val, (rt<<1)|1);
    else
    {
        update(l, mid, val, rt<<1);
        update(mid+1, r, val, (rt<<1)|1);
    }
    cal(rt);
}
int main()
{
    int T, n;
    scanf("%d", &T);
    while(T--)
    {
        scanf("%d", &n);
        int nums = 0;
        for(int i = 0; i < n; i++, nums += 2)
        {
            double x1, y1, x2, y2;
            scanf("%lf%lf%lf%lf", &x1, &y1, &x2, &y2);
            pos[nums] = x1;pos[nums+1] = x2;
            ss[nums].l = x1;ss[nums].r = x2;ss[nums].h = y1;ss[nums].v = 1;
            ss[nums+1].l = x1;ss[nums+1].r = x2;ss[nums+1].h = y2;ss[nums+1].v = -1;
        }
        sort(pos, pos+nums);
        sort(ss, ss+nums, cmp);
        int m = 1;
        for(int i = 1; i < nums; i++)
            if(pos[i] != pos[i-1])
                pos[m++] = pos[i];
        build(0, m-1, 1);
        //cout << "yes\n";
        double ans = 0;
        for(int i = 0; i < nums-1; i++)
        {
            int l = binary(ss[i].l, 0, m-1);
            int r = binary(ss[i].r, 0, m-1) - 1;
            update(l, r, ss[i].v, 1);
            ans += segTree[1].ss * (ss[i+1].h - ss[i].h);
        }
        printf("%.2f\n", ans);
    }
    return 0;
}