poj 1177 Picture（线段树周长并）

题目链接：http://poj.org/problem?id=1177

题意：给你n个矩形问你重叠后外边缘总共多长。

周长并与面积并很像只不过是处理的时候是   增加的周长=abs（上一次的线段的长度-更新后线段的长度）

然后分别处理一下竖着的边和横着的边就好了即建两次树就好。

就是一道典型的周长并问题，可以拿来练练周长并的写法。

 

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdio>
using namespace std;
typedef long long ll;
const int M = 1e5 + 10;
struct ss {
    int l , r , h , flag;
}s1[M << 1] , s2[M << 1];
struct TnT {
    int l , r , add , len;
}T[M << 4];
bool cmp(ss a , ss b) {
    if(a.h == b.h)
        return a.flag > b.flag;
    return a.h < b.h;
}
void build(int l , int r , int p) {
    int mid = (l + r) >> 1;
    T[p].l = l , T[p].r = r , T[p].add = T[p].len = 0;
    if(l == r)
        return ;
    build(l , mid , p << 1);
    build(mid + 1 , r , (p << 1) | 1);
}
void pushup(int p) {
    if(T[p].add) {
        T[p].len = T[p].r - T[p].l + 1;
    }
    else if(T[p].l == T[p].r) {
        T[p].len = 0;
    }
    else {
        T[p].len = T[p << 1].len + T[(p << 1) | 1].len;
    }
}
void updata(int l , int r , int p , int ad) {
    int mid = (T[p].l + T[p].r) >> 1;
    if(T[p].l == l && T[p].r == r) {
        T[p].add += ad;
        pushup(p);
        return ;
    }
    if(mid >= r) {
        updata(l , r , p << 1 , ad);
    }
    else if(mid < l) {
        updata(l , r , (p << 1) | 1 , ad);
    }
    else {
        updata(l , mid , p << 1 , ad);
        updata(mid + 1 , r , (p << 1) | 1 , ad);
    }
    pushup(p);
}
int main() {
    int n;
    scanf("%d" , &n);
    for(int i = 1 ; i <= n ; i++) {
        int x1 , y1 , x2 , y2;
        scanf("%d%d%d%d" , &x1 , &y1 , &x2 , &y2);
        x1 += M , x2 += M , y1 += M , y2 += M;
        s1[i].flag = 1;
        s1[i].l = x1;
        s1[i].r = x2;
        s1[i].h = y1;
        s1[i + n].flag = -1;
        s1[i + n].l = x1;
        s1[i + n].r = x2;
        s1[i + n].h = y2;
        s2[i].flag = 1;
        s2[i].l = y1;
        s2[i].r = y2;
        s2[i].h = x1;
        s2[i + n].flag = -1;
        s2[i + n].l = y1;
        s2[i + n].r = y2;
        s2[i + n].h = x2;
    }
    sort(s1 + 1 , s1 + 1 + 2 * n , cmp);
    sort(s2 + 1 , s2 + 1 + 2 * n , cmp);
    int l , r;
    build(1 , 2 * M , 1);
    ll ans = 0;
    for(int i = 1 ; i <= 2 * n ; i++) {
        int last = T[1].len;
        l = s1[i].l;
        r = s1[i].r - 1;
        updata(l , r , 1 , s1[i].flag);
        ans += abs(last - T[1].len);
    }
    build(1 , 2 * M , 1);
    for(int i = 1 ; i <= 2 * n ; i++) {
        int last = T[1].len;
        l = s2[i].l;
        r = s2[i].r - 1;
        updata(l , r , 1 , s2[i].flag);
        ans += abs(last - T[1].len);
    }
    printf("%lld\n" , ans);
    return 0;
}