题意:
给定n个矩形,求这些矩形所构成图形的周长。
题解:
跟hdu1542很像,那道题是求构成图形的面积,我一开始也以为是面积,最后求了之后才发现时周长。不过那题的代码还是有用的,这题相当于将那道题代码的二倍化。
我们先离散化矩形各条竖直线的y坐标,然后根据x坐标排序,我们可以发现相邻的x坐标间的图形是一个或者多个矩形,那么图形的上下方向上的周长怎么求呢?我们很容易想到ans=∑(xi-xj)*w*2(其中j=i-1,而w表示被两个横坐标分割区域的矩形数目),那么我们就能转换下,每次碰到左边的时候维护线段树[yi,yj](这条边所覆盖的离散化区间)区域加一,碰到右边的时候减一。那么矩形个数就相当于线段树上非0段的个数了。同理,我们将x坐标离散化,排序y坐标可以求得左右方向上的周长。
线段树结点介绍:
p表示该段区间左边是否有非0数,q表示该段区间右边是否有非0数,w表示该段区间的非0段个数,cov表示覆盖情况。
代码:
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <vector>
#include <map>
#include <queue>
#include <stack>
using namespace std;
const int maxn=1e4+10;
struct node{
int l,r,cov;
int p,q,w;
}e[maxn*4];
struct segment{
int x,up,down;
int cov;
}f[maxn],g[maxn];
int y[maxn],x[maxn];
int cmp(segment a,segment b)
{
return a.x<b.x;
}
void build(int a,int b,int c)
{
e[c].l=a;e[c].r=b;
e[c].w=e[c].p=e[c].q=e[c].cov=0;
if(a==b)return ;
int mid=(a+b)/2;
build(a,mid,2*c);
build(mid+1,b,2*c+1);
}
void push_up(int c)
{
if(e[2*c].q==1&&e[2*c+1].p==1)e[c].w=e[2*c].w+e[2*c+1].w-1;
else e[c].w=e[2*c].w+e[2*c+1].w;
e[c].p=e[2*c].p;e[c].q=e[2*c+1].q;
}
void update(int a,int b,int c,int val)
{
//printf("%d\n",c);
if(e[c].l==a&&e[c].r==b)
{
e[c].cov+=val;
if(e[c].cov!=0)e[c].w=e[c].p=e[c].q=1;
else
{
if(e[c].l==e[c].r)e[c].w=e[c].p=e[c].q=0;
else push_up(c);
}
return ;
}
int mid=(e[c].l+e[c].r)/2;
if(b<=mid)update(a,b,2*c,val);
else if(a>mid)update(a,b,2*c+1,val);
else
{
update(a,mid,2*c,val);
update(mid+1,b,2*c+1,val);
}
if(e[c].cov==0)
push_up(c);
}
int main()
{
int n,tt=0;
while(scanf("%d",&n)!=EOF)
{
if(n==0)break;
int i,j,k,t=0,tt=0,l,r;
int x1,y1,x2,y2;
for(i=0;i<n;i++)
{
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
y[++t]=y1;f[t].x=x1;f[t].up=y2;f[t].down=y1,f[t].cov=1;
y[++t]=y2;f[t].x=x2;f[t].up=y2;f[t].down=y1,f[t].cov=-1;
x[++tt]=x1;g[tt].x=y1;g[tt].up=x2;g[tt].down=x1;g[tt].cov=1;
x[++tt]=x2;g[tt].x=y2;g[tt].up=x2;g[tt].down=x1;g[tt].cov=-1;
}
int ans=0;
sort(y+1,y+t+1);
sort(f+1,f+t+1,cmp);
build(1,t,1);
for(i=1;i<=t;i++)
{
//printf("%d:%d\n",i,e[1].w);
ans+=(f[i].x-f[i-1].x)*e[1].w*2;
l=lower_bound(y+1,y+t+1,f[i].down)-y;
r=lower_bound(y+1,y+t+1,f[i].up)-y-1;
update(l,r,1,f[i].cov);
}
//printf("%d\n",ans);
sort(x+1,x+tt+1);
sort(g+1,g+tt+1,cmp);
build(1,tt,1);
for(i=1;i<=tt;i++)
{
//printf("%d:%d\n",i,e[1].w);
ans+=(g[i].x-g[i-1].x)*e[1].w*2;
l=lower_bound(x+1,x+tt+1,g[i].down)-x;
r=lower_bound(x+1,x+tt+1,g[i].up)-x-1;
update(l,r,1,g[i].cov);
}
printf("%d\n",ans);
}
return 0;
}