Picture 线段树扫描线求轮廓线

http://acm.hdu.edu.cn/showproblem.php?pid=1828

思路来自http://www.cnblogs.com/ka200812/archive/2011/11/13/2247259.html

向大牛膜拜

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<map>
#include<cmath>
using namespace std;
#define MAXN 5005
#define lch p<<1
#define rch p<<1|1
struct node
{
    int l,r;
    int cover,len,num;
    int lcover,rcover;
}t[MAXN*8];
struct edge
{
    int x,yd,yu;
    int lr;
    bool operator < (const edge a) const
    {
        if(x==a.x) return lr>a.lr;
        return x<a.x;
    }
}e[MAXN*2];
int yy[MAXN*2];
void construct(int l,int r,int p)
{
    memset(&t[p],0,sizeof(t[p]));
    t[p].l=l,t[p].r=r;
    if(r-l<=1) return ;
    int m=(l+r)>>1;
    construct(l,m,lch);
    construct(m,r,rch);
}
void pushup(int p)
{
    if(t[p].cover)
    {
        t[p].len=yy[t[p].r]-yy[t[p].l];
        t[p].num=t[p].lcover=t[p].rcover=1;
    }
    else if(t[p].r-t[p].l==1)
    {
        t[p].len=t[p].lcover=t[p].rcover=t[p].num=0;
    }
    else
    {
        t[p].len=t[lch].len+t[rch].len;
        t[p].lcover=t[lch].lcover;
        t[p].rcover=t[rch].rcover;
        t[p].num=t[lch].num+t[rch].num;
        if(t[lch].rcover&&t[rch].lcover) t[p].num--;
    }
}
void modify(int l,int r,int val,int p)
{
    if(t[p].l>r||t[p].r<l) return ;
    if(t[p].l>=l&&t[p].r<=r)
    {
        t[p].cover+=val;
        pushup(p);
        return ;
    }
    modify(l,r,val,lch);
    modify(l,r,val,rch);
    pushup(p);
}
int main()
{
    map<int ,int> hash;
    int n,i;
    int x1,x2,yd,yu;
    while(scanf("%d",&n)!=EOF)
    {
        hash.clear();
        for(i=0;i<2*n;i+=2)
        {
            scanf("%d%d%d%d",&x1,&yd,&x2,&yu);
            e[i].x=x1,e[i].yd=yd,e[i].yu=yu,e[i].lr=1;
            e[i+1].x=x2,e[i+1].yd=yd,e[i+1].yu=yu,e[i+1].lr=-1;
            yy[i]=yd;
            yy[i+1]=yu;
        }
        sort(yy,yy+i);
        sort(e,e+i);
        int up=int(unique(yy,yy+i)-yy);
        for(i=0;i<up;i++)
            hash[yy[i]]=i;
        int pre=0,ans=0;
        construct(0,up-1,1);
        for(i=0;i<2*n;i++)
        {
            modify(hash[e[i].yd],hash[e[i].yu],e[i].lr,1);
            
            ans+=abs(pre-t[1].len);
            
            pre=t[1].len;
            
            if(i!=2*n-1)
                ans+=t[1].num*2*(e[i+1].x-e[i].x);
        }
        printf("%d\n",ans);
    }
}