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);
}
}