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 的时候,分几种情况:
- cvr[rt] >= 2:即整个区间都被覆盖至少两次,那么有:
one[rt] = two[rt] = x[r] - x[l]
,直接就是区间的长度; - cvr[rt] == 1:整个区间被完全覆盖 1 次,故
one[rt] = x[r] - x[l]
,而对于two[]
,可能存在某些子区间覆盖有不止 1 次,所以还要再分情况:
- l + 1 >= r:已是最小的区间(没有更小的子区间),那么
two[rt] = 0
; - 否则,
two[rt] = one[rt<<1] + one[rt<<1|1]
,因为本区间已完全覆盖 1 次,子区间只需再覆盖一次,就能满足至少两次
- l + 1 >= r:已是最小的区间(没有更小的子区间),那么
- cvr[rt] == 0:区间没有被完全覆盖,还是要看子区间:
- l + 1 >= r:有子区间,那么:
one[rt] = one[rt<<1] + one[rt<<1|1]
,two[rt] = two[rt<<1] + two[rt<<1|1]
- 否则,
one[rt] = two[rt] = 0
- l + 1 >= r:有子区间,那么:
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
- 题目说的是“左上、右下”,但其实它给的是“左下、右上”…
- 这份代码是从上往下扫