hdoj 1828 （poj 1177） 线段树求矩形并的周长

做这个题还是费了很大心思的，因为之前没有怎么接触过扫描线的方法，另外这个题的线段树节点的设计也是蛮经典的。

具体讲解我就不说了，我也唠叨不清楚。推荐两个博客：

http://blog.csdn.net/wuxinliulei/article/details/9052725

http://blog.csdn.net/k1246195917/article/details/5966523

如果把这两个博客搞懂了，这题就能水之~

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAX 5555
int n;
int ans=0;
struct Node
{
    int l,r;
    int c,lc,rc;
    int line;
    int sum;
};
struct Line
{
    int x,y1,y2;
    int f;
};
Node node[MAX*4];
Line line[MAX*3];
int Y[MAX*3];
int cmp(Line a,Line b)
{
    if(a.x==b.x) return a.f>b.f;
    return a.x<b.x;
}
void build_tree(int l,int r,int k)
{
    node[k].l=l;node[k].r=r;
    node[k].c=node[k].lc=node[k].rc=0;
    node[k].sum=node[k].line=0;
    if(l+1==r) return ;
    int mid=(l+r)/2;
    build_tree(l,mid,k*2);
    build_tree(mid,r,k*2+1);
}
void callen(int k)
{
    if(node[k].c>0)
        node[k].sum=Y[node[k].r]-Y[node[k].l];
    else if(node[k].l+1==node[k].r) node[k].sum=0;
    else node[k].sum=node[k*2].sum+node[k*2+1].sum;
}
void calline(int k)
{
    if(node[k].c>0)
    {
        node[k].line=node[k].rc=node[k].lc=1;
    }
    else if(node[k].l+1==node[k].r)
    {
        node[k].line=node[k].rc=node[k].lc=0;
    }
    else  {
        node[k].rc=node[k*2+1].rc;
        node[k].lc=node[k*2].lc;
        node[k].line=node[k*2].line+node[k*2+1].line-node[k*2].rc*node[k*2+1].lc;
    }
}
void update(Line L,int k)
{
    if(L.y1>=Y[node[k].r]||L.y2<=Y[node[k].l])
        return ;
    if(L.y1<=Y[node[k].l]&&Y[node[k].r]<=L.y2)
    {
        node[k].c+=L.f;
        callen(k);
        calline(k);
        return ;
    }
    update(L,k*2);
    update(L,k*2+1);
    callen(k);
    calline(k);
}
int main()
{
    int i,j;
    int n;
    while(scanf("%d",&n)!=EOF){
        int t=1;
        int x1,y1,x2,y2;
        for(i=1;i<=n;i++)
        {
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            line[t].x=x1;
            line[t].y1=y1;
            line[t].y2=y2;
            line[t].f=1;
            Y[t++]=y1;
            line[t].x=x2;
            line[t].y1=y1;
            line[t].y2=y2;
            line[t].f=-1;
            Y[t++]=y2;
        }
        sort(line+1,line+t,cmp);
        sort(Y+1,Y+t);
        int cnt=1;
        for(i=1;i<t-1;i++)
        {
            if(Y[i]!=Y[i+1])
                Y[cnt++]=Y[i];
        }
        Y[cnt++]=Y[t-1];
        build_tree(1,cnt-1,1);
        ans=0;
        int nowsum=0,nowline=0;
        for(i=1;i<t;i++)
        {
            update(line[i],1);
            if(i>1){
                ans+=2*nowline*(line[i].x-line[i-1].x);
            }
            ans+=abs(node[1].sum-nowsum);
            nowsum=node[1].sum;
            nowline=node[1].line;
        }
        printf("%d\n",ans);
    }
    return 0;
}