POJ-1177-Picture

这个题要求求图形的周长，再次用到扫描线+离散化+线段树

参考了下这个人的博客：http://www.cppblog.com/abilitytao/archive/2010/07/21/120927.html

并且看了陈宏的论文，理解中~

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cstdlib>
#include<algorithm>
using namespace std;
const int maxn=10100;
struct Point
{
int x;
int y1;
int y2;
int flag;
bool operator <(const Point &a)const
{
return x<a.x;
}
}p[maxn];
struct Tree
{
int l;//记录左边
int r;//记录右边
int len;//区间长度
int segnum;//区间被分成的段数
int cover;//区间覆盖总数
int sum;//区间覆盖总长
bool lcover,rcover;
void Clear()
{
len=segnum=cover=sum=lcover=rcover=0;
}
}t[maxn*3];
int n,m,len,py[maxn],y[maxn];
int Search(int val)
{
int l=1,r=len;
while(l<=r)
{
int mid=(l+r)>>1;
if(y[mid]==val)
return mid;
if(y[mid]<val)
l=mid+1;
else
r=mid-1;
}
}
void Build(int l,int r,int index)
{
t[index].Clear();
t[index].l=l;
t[index].r=r;
t[index].len=y[r]-y[l];
if(l+1==r)
return;
int mid=(l+r)>>1;
Build(l,mid,index<<1);
Build(mid,r,index<<1|1);
}
void UpdateNode(int index)
{
if(t[index].cover>0)
{
t[index].sum=t[index].len;
t[index].lcover=t[index].rcover=true;
t[index].segnum=1;
}
else if(t[index].r-t[index].l>1)
{
t[index].sum=t[index<<1].sum+t[index<<1|1].sum;
t[index].lcover=t[index<<1].lcover;
t[index].rcover=t[index<<1|1].rcover;
t[index].segnum=t[index<<1].segnum+t[index<<1|1].segnum-t[index<<1].rcover*t[index<<1|1].lcover;//最后减去的是中间可能的重复计算边
}
else
{
t[index].lcover=t[index].rcover=0;
t[index].segnum=0;
t[index].sum=0;
}
}
void Update(int l,int r,int val,int index)
{
if(t[index].l==l&&t[index].r==r)
{
t[index].cover+=val;
UpdateNode(index);
return;
}
int mid=(t[index].l+t[index].r)>>1;
if(r<=mid)
Update(l,r,val,index<<1);
else if(l>=mid)
Update(l,r,val,index<<1|1);
else
{
Update(l,mid,val,index<<1);
Update(mid,r,val,index<<1|1);
}
UpdateNode(index);
}
int main()
{
while(scanf("%d",&n)!=EOF)
{
m=1;
for(int i=0;i<n;i++)
{
int x1,y1,x2,y2;
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
p[m].x=x1;
p[m].y1=y1;
p[m].y2=y2;
p[m].flag=1;
py[m++]=y1;
p[m].x=x2;
p[m].y1=y1;
p[m].y2=y2;
p[m].flag=-1;
py[m++]=y2;
}
sort(py+1,py+m);
sort(p+1,p+m);
len=1;
y[1]=py[1];
for(int i=2;i<m;i++)
if(py[i]!=py[i-1])
y[++len]=py[i];
Build(1,len,1);
int ans=0,last=0;
for(int i=1;i<m-1;i++)
{
int a=Search(p[i].y1);
int b=Search(p[i].y2);
if(p[i].flag==1)
Update(a,b,1,1);
else
Update(a,b,-1,1);
ans+=t[1].segnum*(p[i+1].x-p[i].x)*2;
ans+=abs(t[1].sum-last);
last=t[1].sum;
}
int a=Search(p[m-1].y1);
int b=Search(p[m-1].y2);
Update(a,b,-1,1);
ans+=abs(t[1].sum-last);
printf("%d\n",ans);
}
return 0;
}