POJ 1177 Picture

/*
 矩形的周长并
 解决思路和矩形的面积并类似,采用扫描线解决,每一次周长的增量也就是扫描线线段长的变化量
 由于一些边界情况比较难处理,考虑x方向和y方向上恰好有两个矩形相接的情形,为了简化处理,可以从x方向上扫一遍,再从y方向上扫一遍
 */
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;

const int maxn = 11111;

int calc(int x){ return x<0 ? -x:x;}

struct edge
{
    int a,b,h;
    int s;
};

bool cmp(const edge &p,const edge &q)
{
    if(p.h==q.h)  return p.s<q.s;
    return p.h<q.h;
}

edge e[maxn];
edge f[maxn];
int tx1[maxn],ty1[maxn],tx2[maxn],ty2[maxn];
int x[maxn];
int sum[maxn<<2];
int lazy[maxn<<2];
int n,m;
int cnt;

void add_rectangel(int x1,int y1,int x2,int y2)
{
    x[m]=x1;
    e[m].a=x1,e[m].b=x2,e[m].h=y1,e[m++].s=-1;
    x[m]=x2;
    e[m].a=x1,e[m].b=x2,e[m].h=y2,e[m++].s=1;
}


void pushup(int l,int r,int rt)
{
    if(lazy[rt]!=0)  sum[rt]=x[r+1]-x[l];
    else  if(l==r)  sum[rt]=0;
    else sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}

void update(int L,int R,int s,int l,int r,int rt)
{
    if(L<=l&&r<=R)
    {
        lazy[rt]+=s;
        pushup(l,r,rt);
        return;
    }
    int mid = (l+r)>>1;
    if(L<=mid)  update(L,R,s,l,mid,rt<<1);
    if(R>mid)  update(L,R,s,mid+1,r,rt<<1|1);
    pushup(l,r,rt);
}

int bst(int v,int x[],int size)
{
    int mid;
    int l,r;
    l=0,r=size-1;
    while(l<=r)
    {
        mid = (l+r)>>1;
        if(x[mid]==v)  return mid;
        if(x[mid]<v)  l=mid+1;
        else  r=mid-1;
    }
    return -1;
}

int input()
{
    if(scanf("%d",&n)==EOF)  return 0;
    for(int i=0;i<n;i++)
    {
        scanf("%d%d%d%d",&tx1[i],&ty1[i],&tx2[i],&ty2[i]);
    }
    return 1;
}

void solve()
{
    m=0;
    for(int i=0;i<n;i++)
    {
        add_rectangel(tx1[i],ty1[i],tx2[i],ty2[i]);
    }
    sort(e,e+m,cmp);
    sort(x,x+m);
    cnt=1;
    for(int i=1;i<m;i++)
    {
        if(x[i]!=x[i-1])  x[cnt++]=x[i];
    }
    memset(sum,0,sizeof(sum));
    memset(lazy,0,sizeof(lazy));
    int ans = 0;
    int L,R;
    for(int i=0;i<m;i++)
    {
        L=bst(e[i].a,x,cnt);
        R=bst(e[i].b,x,cnt)-1;
        if(L<=R)
        {
            int tem=sum[1];
            update(L,R,e[i].s,0,cnt-1,1);
            ans+=calc(sum[1]-tem);
        }
    }
    
    m=0;
    for(int i=0;i<n;i++)
    {
        add_rectangel(ty1[i],tx1[i],ty2[i],tx2[i]);
    }
    sort(e,e+m,cmp);
    sort(x,x+m);
    cnt=1;
    for(int i=1;i<m;i++)
    {
        if(x[i]!=x[i-1])  x[cnt++]=x[i];
    }
    memset(sum,0,sizeof(sum));
    memset(lazy,0,sizeof(lazy));
    for(int i=0;i<m;i++)
    {
        L=bst(e[i].a,x,cnt);
        R=bst(e[i].b,x,cnt)-1;
        if(L<=R)
        {
            int tem=sum[1];
            update(L,R,e[i].s,0,cnt-1,1);
            ans+=calc(sum[1]-tem);
        }
    }
    printf("%d\n",ans);
}

int main()
{
    while(input()!=0)
    {
        solve();
    }
}