POJ1177--Picture（扫描线）

题目大意：给出n个矩形，求周长并

分析：线段树+扫描线。这里给出的是整数，而且范围不大，所以，不需要离散化。

和面积并差不多。不过，还是有点区别的。

numseg[i]数组记录有几条横边需要增加

lbd[i],ubd[i]表示当前区间的上下（也可以理解为左右，因为我是按Y轴建树的，所以自然就是上下啦~）两个端点是否被覆盖

cover[i]表示当前区间被几条边覆盖

len[i]表示当前区间的被覆盖的长度

另外，由于这里有负数，所以，计算区间的一半，一定要用二进制来做，即(l+r) >> 1，如果(l+r)/2的话，就会出现错误。这是为什么呢？我们可以举一个例子来看一下。[-6, 1]把它分为两半的话，明显应该是[-6,-3]，[-2, 1]，一边三个数，但是用(l+r)/2，得到的却是[-6,-2]。那为什么(l+r)>>1就能得到-3，这里就需要计算机组成原理的知识了（因为没好好学，这个错误花了3个小时才发现。呜呜呜...），计算机对负数进行计算的时候，使用补码，首先看-5的原码10000101，负数的补码等于反码+1，所以，补码是11111011，然后右移一位，11111101，然后取反加一，就是它的绝对值，00000011，为3，所以最终结果为-3。

代码：

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;

const int maxn = 20010;

struct Seg{
    int x, y1, y2, left;    //left表示是否为矩形的左边，左边则为1，右边则为-1
    Seg(){}
    Seg(int x, int y1, int y2, int left):x(x), y1(y1), y2(y2), left(left){}
    bool operator < (const Seg & cmp) const {
        if(x == cmp.x) return left > cmp.left;  //这里很重要！当两条竖边重合，如果先更新前一个矩形的右边，那么len[0] = 0，
                                                //再更新后一个矩形的左边，则会多加了重合的那一部分，导致最后结果出错
        else  return x < cmp.x;
    }
}seg[maxn];

int lbd[4*maxn], ubd[4*maxn];
int len[4*maxn];
int numseg[4*maxn];
int cover[4*maxn];

void PushUP(int root, int l, int r) {
    if(cover[root]) {
        lbd[root] = ubd[root] = 1;
        len[root] = r-l+1;
        numseg[root] = 2;   //当前区间被一条竖边覆盖，那么接下来要计算周长时，就多了两个横边
    }
    else if(l == r) {
        len[root] = numseg[root] = lbd[root] = ubd[root] = 0;
    }
    else {
        lbd[root] = lbd[2*root+1];
        ubd[root] = ubd[2*root+2];
        len[root] = len[2*root+1]+len[2*root+2];
        numseg[root] = numseg[2*root+1] + numseg[2*root+2];
        if(ubd[2*root+1] && lbd[2*root+2]) numseg[root] -= 2;   //两个矩形相交，则要减少2条横边
    }
}

void Update(int root, int l, int r, int L, int R, int c) {
    if(L <= l && r <= R) {
        cover[root] += c;
        PushUP(root ,l , r);
        return;
    }
    int m = (l+r) >> 1 ;
    if(L <= m) Update(2*root+1, l, m, L, R, c);
    if(m < R) Update(2*root+2, m+1, r, L, R, c);
    PushUP(root, l, r);
    return;
}

int main() {
    int n;
    printf("%d\n", (-6+1) >> 1);
    while(~scanf("%d", &n)) {
        int m = 0;
        int ymin = 10000, ymax = -10000;
        for(int i = 0; i < n; i++) {
            int a, b, c, d;
            scanf("%d%d%d%d", &a, &b, &c, &d);
            ymin = min(ymin, b);
            ymax = max(ymax, d);
            seg[m++] = Seg(a, b, d ,1);
            seg[m++] = Seg(c, b, d, -1);
        }
        sort(seg, seg+m);
        int ans = 0, last = 0;  //last记录上一条边的情况
        for(int i = 0; i < m; i++) {
            Update(0, ymin, ymax-1, seg[i].y1, seg[i].y2-1, seg[i].left);
            ans += numseg[0]*(seg[i+1].x-seg[i].x);
            ans += abs(len[0]-last);
            last = len[0];
        }
        printf("%d\n", ans);
    }
    return 0;
}