问题描述
N(N<5000) 张矩形的海报,照片和其他同样形状的图片贴在墙上。它们的边都是垂直的或水平的。每个矩形可以部分或者全部覆盖其他矩形。所有的矩形组成的集合的轮廓称为周长。写一个程序计算周长。
图 1 是一个有 7 个矩形的例子:
对应的轮廓为图 2 所示的所有线段的集合:
所有矩形的顶点坐标均为整数。所有的坐标都在 [-10000,10000] 的范围内,并且任何一个矩形面积都为整数。结果的值可能需要 32 位有符号整数表示。
输入格式
第1行: N,张贴在墙上的矩形的数目。 第 2..N+1行 接下来的N行中,每行都有两个点的坐标,分别是矩形的左下角坐标和右上角坐标。每一个坐标由 X 坐标和 Y 坐标组成。
输出格式
只有一行,为一个非负整数,表示输入数据中所有矩形集合的轮廓长度。
样例输入
7
-15 0 5 10
-5 8 20 25
15 -4 24 14
0 -6 16 4
2 15 10 22
30 10 36 20
34 0 40 16
样例输出
228
题解
线段树 —— 扫描线
显然整个图形的轮廓线条一定是竖直或水平的。如果我们求出了竖直的,水平的也可用类似方法。
我们先记录下每条边是矩形的左边还是右边,再把边按横坐标排序。
维护一颗线段树,记录下每个点被覆盖的情况。初始时没有点被覆盖。
接下来顺序扫描每一条边。
情况一:
若讨论到一条左边,将左边表示的纵坐标区间覆盖一层,区间被覆盖总长的增加量即为这条边对答案的贡献。
情况二:
若讨论到一条右边,将右边表示的纵坐标区间去掉一层,区间被覆盖总长的减少量即为这条边对答案的贡献。
横边同理。
似乎有些网站上的题加强了数据?离散化即可。
这道题就解决了!解决了?没有!
此题还有一坑点。输入这组数据试试:
2 0 0 4 4 0 4 4 8
答案是32?显然是24 啊!
需要注意,对边进行排序时,若横(纵)坐标相同,左边(讨论到时覆盖一层的边)应排在前面,否则就会尴尬。
具体见代码。
代码
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
struct idata{
int x,y,q;
bool p;
}h[2][10005];
bool cmp(idata a,idata b)
{return a.q==b.q?a.p>b.p:a.q<b.q;}
const int Q=800000;
int n,lazy[Q],ls[Q],rs[Q],co[Q],maxr=10000,maxl=-10000,t,tot;
void pd(int now)
{
if(!ls[now])ls[now]=++tot;
if(!rs[now])rs[now]=++tot;
int v=lazy[now];
co[ls[now]]+=v;
lazy[ls[now]]+=v;
co[rs[now]]+=v;
lazy[rs[now]]+=v;
lazy[now]=0;
}
void xiu(int now,int l,int r,int x,int y,int v)
{
if(x<=l&&y>=r)
{
co[now]+=v;
lazy[now]+=v;
return;
}
if(lazy[now]!=0)pd(now);
int mid=(l+r)>>1;
if(x<mid){
if(!ls[now])ls[now]=++tot;
xiu(ls[now],l,mid,x,y,v);
}
if(y>mid){
if(!rs[now])rs[now]=++tot;
xiu(rs[now],mid,r,x,y,v);
}
}
int gs(int now,int l,int r)
{
if(r-l==1||((!ls[now])&&(!rs[now])))return (co[now]>0)*(r-l);
if(lazy[now]!=0)pd(now);
int mid=(l+r)>>1,temp=0;
if(ls[now])temp+=gs(ls[now],l,mid);
if(rs[now])temp+=gs(rs[now],mid,r);
return temp;
}
int main()
{
int i,ans=0,a,b,c,d,last,haha;
scanf("%d",&n);
for(i=1;i<=n;i++)
{
scanf("%d%d%d%d",&a,&b,&c,&d);
h[0][2*i-1].p=h[1][2*i-1].p=false;
h[0][2*i].p=h[1][2*i].p=true;
h[0][2*i-1].x=h[0][2*i].x=a;
h[0][2*i-1].y=h[0][2*i].y=c;
h[0][2*i].q=b,h[0][2*i-1].q=d;
h[1][2*i].x=h[1][2*i-1].x=b;
h[1][2*i].y=h[1][2*i-1].y=d;
h[1][2*i].q=a,h[1][2*i-1].q=c;
}
n*=2;
sort(h[0]+1,h[0]+n+1,cmp);
sort(h[1]+1,h[1]+n+1,cmp);
for(t=0;t<2;t++)
{
memset(ls,0,sizeof(ls));
memset(rs,0,sizeof(rs));
memset(co,0,sizeof(co));
memset(lazy,0,sizeof(lazy));
tot=1;
last=0;
for(i=1;i<=n;i++)
{
if(h[t][i].p)xiu(1,maxl,maxr,h[t][i].x,h[t][i].y,1);
else xiu(1,maxl,maxr,h[t][i].x,h[t][i].y,-1);
haha=gs(1,maxl,maxr);
ans+=abs(haha-last),last=haha;
}
}
printf("%d",ans);
return 0;
}
最后,感谢机智的何老板提供的反例数据。