POJ 1177 Picture（扫描线 + 线段树 矩形覆盖区域周长的并）

题意：

求矩形所覆盖的区域周长的并。

思路：

1. 首先是对于 x 轴线的扫描，这点和上题POJ 1151类似，不过求周长和面积不同的一点是前面一条线和后面一条线的差。

2. 对于高，要另外设置几组参数和标记：ynum[rt] 表示在 rt 表示的区域范围内，y 坐标方向竖线的个数

3. lbd[rt], rbd[rt] 分别表示左右边界是否被扫描线覆盖，这样的话求线段的并的时候变能够判断边界部分是否有重合。

4. 以上操作是在每两条扫描线之间进行的，所以这里需要一个抽象的过程，不过也不算多难理解。

5. 周长和面积的区别是：周长要一直循环 m(扫描线的个数) 次，因为最上层的一条线毕竟需要加入周长范围内。

#include <iostream>
#include <algorithm>
using namespace std;

#define lhs l, m, rt << 1
#define rhs m + 1, r, rt << 1 | 1

const int maxn = 20010;
int xsum[maxn << 2], ynum[maxn << 2];
int cnt[maxn << 2];
bool lbd[maxn << 2], rbd[maxn << 2];

structSegment {
int l, r, h, v;
Segment() { }
Segment(int _l, int _r, int _h, int _v) : l(_l), r(_r), h(_h), v(_v) { }
bool operator < (const Segment& other) { return h < other.h; }
} seg[maxn] ;

void PushUp(int l, int r, int rt)
{
if (cnt[rt])
{
lbd[rt] = rbd[rt] = true;
xsum[rt] = r - l + 1;
ynum[rt] = 2;
}
else if (l == r)
{
lbd[rt] = rbd[rt] = false;
xsum[rt] = ynum[rt] = 0;
}
else
{
lbd[rt] = lbd[rt << 1];
rbd[rt] = rbd[rt << 1 | 1];
xsum[rt] = xsum[rt << 1] + xsum[rt << 1 | 1];
ynum[rt] = ynum[rt << 1] + ynum[rt << 1 | 1];

if (lbd[rt << 1 | 1] && rbd[rt << 1])
ynum[rt] -= 2;
}
}

void Update(int beg, int end, int value, int l, int r, int rt)
{
if (beg <= l && r <= end)
{
cnt[rt] += value;
PushUp(l, r, rt);
return ;
}
int m = (l + r) >> 1;
if (beg <= m)
Update(beg, end, value, lhs);
if (end > m)
Update(beg, end, value, rhs);
PushUp(l, r, rt);
}

int main()
{
int n;
while (~scanf("%d", &n))
{
int m = 0;
int beg = 10000, end = -10000;
for (int i = 0; i < n; ++i)
{
int a, b, c, d;
scanf("%d %d %d %d", &a, &b, &c, &d);
beg = min(a, beg);
end = max(c, end);
seg[m++] = Segment(a, c, b, 1);
seg[m++] = Segment(a, c, d, -1);
}
sort(seg, seg + m);

int ret = 0, pre = 0;
for (int i = 0; i < m; ++i)
{
if (seg[i].l < seg[i].r)
Update(seg[i].l, seg[i].r - 1, seg[i].v, beg, end - 1, 1);
ret += abs(xsum[1] - pre);
pre = xsum[1];
ret += ynum[1] * (seg[i+1].h - seg[i].h);
}
printf("%d\n", ret);
}
return 0;
}