HDU 1828 Picture (线段树+扫描线)(周长并)

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

给你n个矩形，让你求出总的周长。

类似面积并，面积并是扫描一次，周长并是扫描了两次，x轴一次，y轴一次。每次加起来的无非都是新加的边(flag为1)或者是新减的边(flag为-1)，即加起来的是此时的总长度(T[1].val)减去上一次扫到的总长度(last)的绝对值(T[1].val - last)。(注意用c++提交，g++会wa)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int MAXN = 2e4 + 10;
struct data {
    int l , r , h , flag;
    bool operator <(const data& cmp) const {
        return h < cmp.h;
    }
}x[MAXN] , y[MAXN];
struct segtree {
    int l , r , val , add;
}T[MAXN << 2];

int f(int a) {
    return (a > 0 ? a : -a);
}

void build(int p , int l , int r) {
    int mid = (l + r) >> 1;
    T[p].l = l , T[p].r = r , T[p].add = T[p].val = 0;
    if(r - l == 1) {
        return ;
    }
    build(p << 1 , l , mid);
    build((p << 1)|1 , mid , r);
}

void pushup(int p) {
    if(T[p].add) {
        T[p].val = T[p].r - T[p].l;
    }
    else if(T[p].r - T[p].l == 1) {
        T[p].val = 0;
    }
    else {
        T[p].val = T[p << 1].val + T[(p << 1)|1].val;
    }
}

void updata(int p , int l , int r , int add) {
    int mid = (T[p].l + T[p].r) >> 1;
    if(l == T[p].l && T[p].r == r) {
        T[p].add += add;
        pushup(p);
        return ;
    }
    if(r <= mid) {
        updata(p << 1 , l , r , add);
    }
    else if(l >= mid) {
        updata((p << 1)|1 , l , r , add);
    }
    else {
        updata(p << 1 , l , mid , add);
        updata((p << 1)|1 , mid , r , add);
    }
    pushup(p);
}

int main()
{
    int n , x1 , x2 , y1 , y2;
    while(~scanf("%d" , &n)) {
        for(int i = 0 ; i < n ; i++) {
            scanf("%d %d %d %d" , &x1 , &y1 , &x2 , &y2);
            x1 += 1e4 , x2 += 1e4 , y1 += 1e4 , y2 += 1e4;
            if(x1 > x2)
                swap(x1 , x2);
            if(y1 > y2)
                swap(y1 , y2);
            int ls = i<<1 , rs = (i<<1)|1;
            x[ls].l = x1 , x[ls].r = x2 , x[ls].h = y1 , x[ls].flag = 1;
            x[rs].l = x1 , x[rs].r = x2 , x[rs].h = y2 , x[rs].flag = -1;
            y[ls].l = y1 , y[ls].r = y2 , y[ls].h = x1 , y[ls].flag = 1;
            y[rs].l = y1 , y[rs].r = y2 , y[rs].h = x2 , y[rs].flag = -1;
        }
        sort(x , x + n * 2);
        sort(y , y + n * 2);
        build(1 , 0 , 2e4);
        int res = 0 , last = 0;
        for(int i = 0 ; i < 2 * n ; i++) {
            updata(1 , x[i].l , x[i].r , x[i].flag);
            res += f(last - T[1].val);
            last = T[1].val;
        }
        for(int i = 0 ; i < 2 * n ; i++) {
            updata(1 , y[i].l , y[i].r , y[i].flag);
            res += f(last - T[1].val);
            last = T[1].val;
        }
        printf("%d
" , res);
    }
}