题目链接:http://poj.org/problem?id=1177
题意:
题目就是说n个矩形组成的新图形的周长,每个矩形会给出左下角和右上角的坐标。
分析:
线段树+扫描线+离散化,我是从左到右扫描的,需要将纵坐标离散化,然后每扫描到一根线,求出其覆盖y方向的总长度,减去上次求得的长度,即为本条线增加或减少的长度,同时可以求出每两条线之间的距离,即为横坐标方向的周长的一部分,由于是矩形构成的,横坐标方向是为周长的部分一定成对出现,同时需要用到一个很重要的记录(num),因为可能出现两条竖线之间有4,6…甚至更多条外围的线段。
源代码:
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<algorithm>
using namespace std;
const int N=5005;
typedef __int64 ll;
struct seg
{
int y1,y2;
int x;
int flag;
seg(){}
seg(int _y1,int _y2,int _x,int f):y1(_y1),y2(_y2),x(_x),flag(f){}
friend bool operator < (seg a,seg b)
{
if(a.x==b.x)
return a.flag>b.flag;
return a.x<b.x;
}
}Line[N<<1];
struct TNode
{
int l,r;
int num;//子树中有多少条线段
int cover;//是否被覆盖,入边加1,出边-1
int len;//当前节点覆盖y方向的总长度
bool lb,rb;//左右节点是否被覆盖
}tre[N<<3];
vector<int>vec;
void build(int o,int l,int r)
{
tre[o].l=l;
tre[o].r=r;
tre[o].num=0;
if(l==r-1)
return;
int mid=(l+r)/2;
build(o<<1,l,mid);
build(o<<1|1,mid,r);
}
void update_len(int o)
{
if(tre[o].cover>0)
{
tre[o].len=vec[tre[o].r]-vec[tre[o].l];
tre[o].lb=tre[o].rb=true;
tre[o].num=1;
return;
}
if(tre[o].l==tre[o].r-1)
{
tre[o].cover=tre[o].lb=tre[o].rb=0;
tre[o].num=tre[o].len=0;
return;
}
tre[o].lb=tre[o<<1].lb;
tre[o].rb=tre[o<<1|1].rb;
tre[o].len=tre[o<<1].len+tre[o<<1|1].len;
tre[o].num=tre[o<<1].num+tre[o<<1|1].num-(tre[o<<1].rb&tre[o<<1|1].lb);
}
void update(int o,seg p)
{
if(p.y1<=vec[tre[o].l]&&p.y2>=vec[tre[o].r])
{
tre[o].cover+=p.flag;
update_len(o);
return;
}
if(tre[o].l==tre[o].r-1)
return;
int mid=(tre[o].l+tre[o].r)/2;
if(p.y1<=vec[mid])
update(o<<1,p);
if(p.y2>vec[mid]) update(o<<1|1,p);
update_len(o);
}
int main()
{
int n;
while(~scanf("%d",&n))
{
int k=0;
for(int i=0;i<n;i++)
{
int x1,y1,x2,y2;
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
Line[k++]=seg(y1,y2,x1,1);//进
//y[k]=y1;
Line[k++]=seg(y1,y2,x2,-1);//出
//y[k]=y2;
vec.push_back(y1);
vec.push_back(y2);
}
sort(vec.begin(),vec.end());
vec.erase(unique(vec.begin(),vec.end()),vec.end());
sort(Line,Line+k);
build(1,0,vec.size()-1);
//printf("1111\n");
int ans=0;
int old=0;
int lines=0;
for(int i=0;i<k;i++)
{
//printf("%d\n",Line[i].x);
update(1,Line[i]);
if(i>=1)
{
ans+=lines*2*(Line[i].x-Line[i-1].x);
}
ans+=abs(tre[1].len-old);
old=tre[1].len;
lines=tre[1].num;
}
printf("%d\n",ans);
}
return 0;
}