1.用处:用来处理矩形的问题(所形成的多边形的周长、面积并等)
2.方法:已知每个矩形的左下端点(x1,y1)与右上端点(x2,y2),对每个端点记录四元组(x,y1,y2,w),每次选择x坐标相邻的两个点进行运算
3.代码(拿面积并举例子)
#include<bits/stdc++.h>
using namespace std;
struct node{
double x,l,r;
int w;
};
node p[1010];
int cmp(node x,node y){
return x.x<y.x;
}
struct segment_tree{
double x,l,r,ans;
int cs;
};
segment_tree tr[1010];//记录每次更新后矩形的高
int n,gs;
double s[1010];
void build(int id,int l,int r);
void insert(int id,double x,double y,int w);
int main(){
scanf("%d",&n);
while(n){
printf("Test case #%d\n",++gs);
for(int i=1;i<=n*2;i++) s[i]=p[i].x=p[i].l=p[i].r=0.00;
for(int i=1;i<=n;i++){
double x,y,xx,yy;
scanf("%lf%lf%lf%lf",&x,&y,&xx,&yy);
p[i].x=x;p[i].l=y;p[i].r=yy;p[i].w=1;//左下的端点记录
p[i+n].x=xx;p[i+n].l=y;p[i+n].r=yy;p[i+n].w=-1;//右上的端点记录
s[i]=y;
s[i+n]=yy;
}
sort(s+1,s+(2*n+1));
sort(p+1,p+(2*n+1),cmp);
int num=unique(s+1,s+(2*n+1))-s-1;//离散化缩小x坐标的范围
build(1,1,num);
double ans=0.00;
for(int i=1;i<=2*n;i++){
double h=p[i].x-p[i-1].x;
ans+=tr[1].ans*h;
insert(1,p[i].l,p[i].r,p[i].w);//每次更新矩形的高
}
cout<<"Total explored area: "<<fixed<<setprecision(2)<<ans<<endl;
scanf("%d",&n);
}
}
void build(int id,int l,int r){
tr[id].l=s[l];
tr[id].r=s[r];
tr[id].ans=0.00;
tr[id].cs=0;
if(l+1==r){
return ;
}
int mid=(l+r)>>1;
build(id<<1,l,mid);
build(id<<1|1,mid,r);
}
void insert(int id,double x,double y,int w){
if(tr[id].l>=x&&tr[id].r<=y){
tr[id].cs+=w;
if(tr[id].cs) tr[id].ans=tr[id].r-tr[id].l;
else tr[id].ans=(tr[id].l+1==tr[id].r?0:tr[id<<1].ans+tr[id<<1|1].ans);
return ;
}
if(tr[id<<1].r>x) insert(id<<1,x,y,w);
if(tr[id<<1|1].l<y) insert(id<<1|1,x,y,w);
tr[id].ans=(tr[id].cs==0?tr[id<<1].ans+tr[id<<1|1].ans:tr[id].ans);
}