hdu 1255 覆盖的面积

Problem

acm.hdu.edu.cn/showproblem.php?pid=1255

Reference

hdu 1255 覆盖的面积
矩形面积并、矩形面积交、矩形周长并（线段树、扫描线总结）

Meaning

给出 n 个矩形，求它们相交部分的总面积。

Analysis

扫描线求矩形面积交，与求面积并类似。
求面积并时有个 cvr 数组（CoVeR），表示区间被完全覆盖的次数。相交部分，也即被覆盖至少两次。但不能简单就把判断条件从cvr[rt] > 0改成cvr[rt] > 1就结束，因为整个区间没有完全覆盖多于一次，但某些子区间可能有。
算面积并时只要一个数组统计总长，算面积交时要两个：一个和算面积并的那个相同，统计覆盖至少一次的区间总长（one[]）；另一个统计覆盖至少两次的区间总长（two[]）。
push_up 的时候，分几种情况：

1. cvr[rt] >= 2：即整个区间都被覆盖至少两次，那么有：one[rt] = two[rt] = x[r] - x[l]，直接就是区间的长度；
2. cvr[rt] == 1：整个区间被完全覆盖 1 次，故one[rt] = x[r] - x[l]，而对于two[]，可能存在某些子区间覆盖有不止 1 次，所以还要再分情况：
   
   1. l + 1 >= r：已是最小的区间（没有更小的子区间），那么two[rt] = 0；
   2. 否则，two[rt] = one[rt<<1] + one[rt<<1|1]，因为本区间已完全覆盖 1 次，子区间只需再覆盖一次，就能满足至少两次
3. cvr[rt] == 0：区间没有被完全覆盖，还是要看子区间：
   
   1. l + 1 >= r：有子区间，那么：one[rt] = one[rt<<1] + one[rt<<1|1]，two[rt] = two[rt<<1] + two[rt<<1|1]
   2. 否则，one[rt] = two[rt] = 0

Code

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 1000;

struct segment
{
    double l, r, h;
    int v;

    segment() {}

    segment(double _l, double _r, double _h, int _v):
        l(_l), r(_r), h(_h), v(_v) {}

    bool operator < (const segment &rhs) const
    {
        return h > rhs.h;
    }
} s[N<<1];

double x[N<<1], one[N<<3], two[N<<3];
int cvr[N<<3];

void pushup(int l, int r, int rt)
{
    if(cvr[rt] > 1) // 区间被完全覆盖超过 1 次
        one[rt] = two[rt] = x[r] - x[l];
    else if(cvr[rt] == 1) // 被完全覆盖只有 1 次
    {
        one[rt] = x[r] - x[l];
        if(l + 1 >= r) // 没有子区间
            two[rt] = 0.0;
        else // 只需子区间再覆盖 1 次
            two[rt] = one[rt<<1] + one[rt<<1|1];
    }
    else if(l + 1 >= r) // 没被完全覆盖过，又没有子区间
        one[rt] = two[rt] = 0.0;
    else // 没被完全覆盖，但有子区间
    {
        one[rt] = one[rt<<1] + one[rt<<1|1];
        two[rt] = two[rt<<1] + two[rt<<1|1];
    }
}

void update(int ul, int ur, int v, int l, int r, int rt)
{
    if(ul <= l && r <= ur)
    {
        cvr[rt] += v;
        pushup(l, r, rt);
        return;
    }
    int m = l + r >> 1;
    if(ul < m)
        update(ul, ur, v, l, m, rt<<1);
    if(ur > m)
        update(ul, ur, v, m, r, rt<<1|1);
    pushup(l, r, rt);
}

int main()
{
    int T;
    scanf("%d", &T);
    while(T--)
    {
        int n;
        scanf("%d", &n);
        for(int i = 0; i < n; ++i)
        {
            double l, r, u, d;
            scanf("%lf%lf%lf%lf", &l, &d, &r, &u);
            s[i] = segment(l, r, u, 1);
            s[i+n] = segment(l, r, d, -1);
            x[i] = l;
            x[i+n] = r;
        }
        n <<= 1;
        sort(x, x + n);
        sort(s, s + n);
        int m = unique(x, x + n) - x;
        memset(cvr, 0, sizeof cvr);
        memset(one, 0, sizeof one);
        memset(two, 0, sizeof two);
        double ans = 0.0;
        for(int i = 0, l, r; i < n - 1; ++i)
        {
            l = lower_bound(x, x + m, s[i].l) - x;
            r = lower_bound(x, x + m, s[i].r) - x;
            update(l, r, s[i].v, 0, m-1, 1);
            ans += two[1] * (s[i].h - s[i+1].h);
        }
        printf("%.2f\n", ans);
    }
    return 0;
}

Post Script

• 题目说的是“左上、右下”，但其实它给的是“左下、右上”…
• 这份代码是从上往下扫