第一次做线段成段更新的题,这题需要离散化。
题意:给出几个矩形,求其覆盖面积。
思路:先离散化x坐标,然后加入矩形的上下边,对矩形上下边根据高度进行排序后。建树,从最高的边开始更新线段树,然后就是通过给每条边赋值上边-1下边1,通过节点中
iosum来判断每次线段树中的有效边,每次ans+=有效边*下一层边的高度,最后得出答案。
代码:
#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <memory.h>
#define lson l,m,rt<<1
#define rson m,r,rt<<1|1
using namespace std;
double posx[2001];//对x轴进行离散化
int mx,me;
struct edge
{
double x1,x2,y;
int io;//矩形下边的度为1,矩形上边的度为-1(暂且叫作度)
} ee[2001]; //对变进行高度排序
struct node
{
int iosum;//边的度和
double with;//有效高度
} num[8888];
int cmp1(const void *a,const void *b)
{
return *(double *)a>*(double *)b?1:-1;
}
int cmp2(const void *a,const void *b)
{
return (*(edge *)a).y>(*(edge *)b).y?1:-1;
}
int bin(double x)
{
int low=0,high=mx-1,mid;
while(low<=high)
{
mid=(low+high)>>1;
if(posx[mid]==x) return mid;
else if(posx[mid]>x) high=mid-1;
else low=mid+1;
}
return -1;
}
void updata(int a,int b,int io,int l,int r,int rt)
{
if(a<=l&&r<=b)
{
num[rt].iosum+=io;
if(num[rt].iosum)num[rt].with=posx[r]-posx[l]; //
else if(l+1==r)num[rt].with=0; //
else num[rt].with=num[rt<<1].with+num[rt<<1|1].with;//
return;
}
int m=(l+r)>>1;
if(a<m)updata(a,b,io,lson);
if(b>m)updata(a,b,io,rson);
if(num[rt].iosum)num[rt].with=posx[r]-posx[l]; //
else if(l+1==r)num[rt].with=0; //
else num[rt].with=num[rt<<1].with+num[rt<<1|1].with;//
}
int main()
{
int n,i,ca=1,a,b;
double x1,x2,y1,y2,res;
while(1)
{
scanf("%d",&n);
if(n==0)break;
me=mx=0;
for(i=0; i<n; i++)
{
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
posx[mx++]=x1;
posx[mx++]=x2;
ee[me].x1=x1;
ee[me].x2=x2;
ee[me].y=y1;
ee[me++].io=1;
ee[me].x1=x1;
ee[me].x2=x2;
ee[me].y=y2;
ee[me++].io=-1;
}
qsort(posx,mx,sizeof(posx[0]),cmp1);
for(i=1,mx=1; i<2*n; i++)
if(posx[i]!=posx[i-1])posx[mx++]=posx[i];
qsort(ee,me,sizeof(ee[0]),cmp2);
memset(num,0,sizeof(node)*8888);
res=0;
for(i=0; i<me-1; i++)
{
a=bin(ee[i].x1);
b=bin(ee[i].x2);
updata(a,b,ee[i].io,0,mx-1,1);
res+=num[1].with*(ee[i+1].y-ee[i].y);
}
printf("Test case #%d\n",ca++);
printf("Total explored area: %.2lf\n\n",res);
}
return 0;
}