终于做到了这部分,以前看到这种类型的题只能直接放弃的,现在终于AC了第一个这种题
因为坐标都是浮点数,所以需要对坐标进行离散化(需要插入线段树的坐标进行离散化),为了方便,我直接用了map
将每个矩形的上下两条水平边存到数组中(得记录这条边是下边还是上边,为了计算覆盖次数,下边记为1,上边记为-1),按y的大小排序;从y最小的边开始向上扫描,首先将一条边插入线段树,然后得到当前当前扫描线所在位置的覆盖到的边的总长,
第一次是A1B1边插入到书中,扫描线覆盖到的边的总长就是A1B1,再拿他乘以当前高度和下一条边的高度差,得到底下那个矩形的面积
第二次是第二条黑边插入线段树,这时扫描到的总长度为A2B2,再拿他乘以当前高度和下一条边的高度差,得到中间那个长矩形的面积
第三次是A1B1的对边插入线段树,因为这条边是上边,所以覆盖次数会被-1,挺次,原先线段树中被A1B1覆盖的部分会被-1,所以,此时线段树中能够计算得到扫面线扫到的总长为A3B3的长度,再拿他乘以当前高度和下一条边的高度差,得到最上面那个矩形的面积
三者相加就是总面积了
插入的时候也得注意,不是根据坐标计算出其离散值就当做区间直接插入,而是要将右坐标的离散值减1,这样做的原因是:
这样的两条线段其实没有相互覆盖,但是如果r不减1的话,他们就会变成彼此覆盖的了……在这颗线段树中,一个点其实表示一个区间
这么做了以后,在计算时只要再把r加上1就好了
以下是代码:
#include <cstdio>
#include <cstring>
#include <map>
#include <algorithm>
using namespace std;
const int MAX = 210;
struct line
{
double l,r,h;
int f;
line(){}
line(double a,double b,double c,int d):l(a),r(b),h(c),f(d){}
}a[MAX<<1];
double sum[MAX<<2],x0[MAX];
int cov[MAX<<2];
map<double,int>m;
int cmp(const line &a,const line &b)
{
return a.h<b.h;
}
void push_up(int l,int r,int rt)
{
if(cov[rt]) sum[rt] = x0[r+1] - x0[l];
else if(l==r) sum[rt]=0;
else sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
void update(int c,int L,int R,int l,int r,int rt)
{
if(L<=l&&R>=r)
{
cov[rt]+=c;
push_up(l,r,rt);
return;
}
int mid = r+l>>1;
if(L<=mid) update(c,L,R,l,mid,rt<<1);
if(R>mid) update(c,L,R,mid+1,r,rt<<1|1);
push_up(l,r,rt);
}
int main()
{
int n;
double x1,x2,y1,y2;
int cas = 0;
while(scanf("%d",&n)&&n)
{
int cnt = 0;
int flag = 0;
for(int i=0; i<n; i++)
{
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
a[cnt++] = line(x1,x2,y1,1);
a[cnt++] = line(x1,x2,y2,-1);
if(m[x1]==0) m[x1] = ++flag;
if(m[x2]==0) m[x2] = ++flag;
}
map<double,int>::iterator iter;
int k = 0;
for(iter=m.begin(); iter!=m.end(); iter++)
{
x0[++k] = iter->first;
iter->second = k;
}
sort(a,a+cnt,cmp);
double ans = 0;
memset(sum,0,sizeof(sum));
memset(cov,0,sizeof(cov));
for(int i=0; i<cnt-1; i++)
{
int l = m[a[i].l];
int r = m[a[i].r]-1;
update(a[i].f,l,r,1,flag,1);
ans+=(a[i+1].h-a[i].h)*sum[1];
}
printf("Test case #%d\n",++cas);
printf("Total explored area: %.2lf\n\n",ans);
m.clear();
}
return 0;
}