poj 1151 Atlantis(矩形面积并)

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=500;
struct node{
double x;
int l,r,t;   //t为上下边标志
}f[maxn];  //水平条
int n;
double q[maxn],x1[maxn],yy1[maxn],x2[maxn],yy2[maxn];
//q存储排序后的x坐标
struct segment{
int mark;  //并区间标志
double len;//并区间长度
}tree[maxn*20];//线段树
int cmp(node a,node b){
return a.x<b.x;
}
int ins(int k,int l,int r,int lc,int rc,int t){
if(lc<=l&&r<=rc){  //将标志为t的水平条[l,r]插入根为k,区间为[lc,rc]的线段树
tree[k].mark+=t;//若覆盖，则计算并区域标志；否则插入左右子树
}
else{
if((l+r)/2>=lc) ins(k*2,l,(l+r)/2,lc,rc,t);
if((l+r)/2<rc) ins(k*2+1,(l+r)/2+1,r,lc,rc,t);
}
//并区间结束，则k节点的区间长度设为左右子区间的长度之和；否则为q中第r个x坐标与第l个x坐标之间的距离
if(tree[k].mark==0)  tree[k].len=tree[k*2].len+tree[k*2+1].len;
else tree[k].len=q[r+1]-q[l];
return 0;
}
int main(){
int t=0;
while(scanf("%d",&n)!=EOF){
if(n==0) break;
double ans=0;
for(int i=1;i<=n;i++){
scanf("%lf%lf%lf%lf",&x1[i],&yy1[i],&x2[i],&yy2[i]);
if(x1[i]>x2[i]) swap(x1[i],x2[i]);
if(yy1[i]>yy2[i]) swap(yy1[i],yy2[i]);
q[i*2-2]=x1[i];
q[i*2-1]=x2[i];//存储x坐标
}
sort(q,q+n*2);
int m=unique(q,q+n*2)-q;
for(int i=1;i<=n;i++){
//将矩阵i的上边左右端点在q表的指针、y坐标、上边标志存入f[i*2-2];将矩阵i的下边左右端点在q表的指针、y坐标、下边标志存入f[i*2-1]
f[i*2-2].l=lower_bound(q,q+m,x1[i])-q;
f[i*2-2].r=lower_bound(q,q+m,x2[i])-q;
f[i*2-2].x=yy1[i];
f[i*2-2].t=1;
f[i*2-1].l=lower_bound(q,q+m,x1[i])-q;
f[i*2-1].r=lower_bound(q,q+m,x2[i])-q;
f[i*2-1].x=yy2[i];
f[i*2-1].t=-1;
}
sort(f,f+n*2,cmp);
for(int i=0;i<n*2;i++){
if(i) ans+=tree[1].len*(f[i].x-f[i-1].x);
ins(1,0,m,f[i].l,f[i].r-1,f[i].t);//将当前水平条插入线段树
}
printf("Test case #%d\n",++t);
printf("Total explored area: %.2f\n\n",ans);
}
return 0;
}