Poj 1177 Picture

一个经典的扫描线入门题

矩形周长并
题意就不描述了

主要记一下写这个题的时候的一些想法

矩形周长并的扫描线做法中强调的一点是每次要求一个0和1的交换（也就是0变成1或者1变成0）
或者说，每次更新之后我们要统计的是发生0和1交换的区间的长度
但是，为什么呢，为什么一定是0和1的交换呢？

不妨在脑海里建立一个直角坐标系，所有区域的初始值都是0，对于每一个矩形，我们把矩形所占的区域的值都+1

如果你真的按照这个方法做了，你就会发现，这个时候我们需要求的图形的边界，就是0和非0的数的临界点（想一想，为什么

这个时候我们的扫描线算法就呼之欲出了

对于所有的和y轴平行的边，我们扫一遍，然后和x轴平行的边可以类似的扫一遍

具体的过程可以参照代码

值得注意的一点这个线段树的维护，所有的lazy标记是不用pushdown的，因为对于每一个+1的事件，我们总有一个-1的事件和它对应（覆盖的区间相等）。

---

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

const int maxn = 5123*2;

int msp[maxn]; // [L,R)

int arr[maxn*4];
int lazy[maxn*4];



void initS(){
    memset(arr,0,sizeof(arr));
    memset(lazy,0,sizeof(lazy));
}

#define lson o<<1,l,m
#define rson o<<1|1,m+1,r
#define Now int o,int l,int r
#define Mid int m = (l+r)/2

void upd(Now){
    if(lazy[o] == 0)
        arr[o] = arr[o<<1] + arr[o<<1|1];
    else
        arr[o] = msp[r+1]-msp[l];
}

void update(Now,int L,int R,int v){
    if(L <= l && r <= R){
        lazy[o] += v;
        upd(o,l,r);
        return;
    }
    Mid;
    if(L <= m) update(lson,L,R,v);
    if(m+1 <= R) update(rson,L,R,v);
    upd(o,l,r);
}

struct Rect{
    int x1,x2;
    int y1,y2;
    void init(){
        scanf("%d %d %d %d",&x1,&y1,&x2,&y2);
    }
    void swaper(){
        swap(x1,y1);
        swap(x2,y2);
    }
};

Rect rect[maxn];

struct event{
    int l,r;
    int tim;
    int add;
    void init(int L,int R,int t,int a){
        l = L,r = R;
        tim = t;
        add = a;
    }
};

bool operator == (event a,event b){
    return a.tim == b.tim && a.add == b.add;
}

event eve[maxn*2];

bool cmp(event a,event b){
    if(a.tim != b.tim)
        return a.tim < b.tim;
    return a.add > b.add;
}

int getloca(int x,int n){
    return lower_bound(msp,msp+n,x)-msp;
}

int calanser(int n){
    int m = 0;
    int len=0;
    msp[len++] = -10100;
    for(int i=0;i<n;i++){
        eve[m++].init(rect[i].y1,rect[i].y2,rect[i].x1,1);
        eve[m++].init(rect[i].y1,rect[i].y2,rect[i].x2,-1);
        msp[len++] = rect[i].y1;
        msp[len++] = rect[i].y2;
    }
    sort(eve,eve+m,cmp);
    sort(msp,msp+len);
    len = unique(msp,msp+len) - msp;
    initS();
    int ret = 0;
    int bf = 0;
    for(int i=0;i<m;i++){
        if(i>0 && !(eve[i] == eve[i-1])){
            int qa = arr[1];
            ret += abs(bf-qa);
            bf = qa;
        }
        int l = getloca(eve[i].l,len);
        int r = getloca(eve[i].r,len);
        update(1,1,len,l,r-1,eve[i].add);
    }
    return ret+bf;
}

int main(){
    int n;
    while(~scanf("%d",&n)){
        for(int i=0;i<n;i++){
            rect[i].init();
        }
        int ans = 0;
        ans += calanser(n);
        for(int i=0;i<n;i++){
            rect[i].swaper();
        }
        ans += calanser(n);
        printf("%d\n",ans);
    }
    return 0;
}