hdu 1542 Atlantis 矩形面积并

//	hdu 1542 Atlantis 矩形面积并
//
//	解题思路:
//
//		将矩形沿着x轴或者y轴切割.然后求和
//		具体做法:
//			将每个矩形看成是两条线段.将x方向
//		坐标离散化.给每个线段一个值,下面线段
//		为1,上面线段为-1.将这些线段映射到线段
//		树上.具体做法是,只要保存线段的左端点就
//		可以了.然后是线段树的基本操作了.更新成
//		段更新,查询是整个区间的值就是sum[1].仔细
//		理解倒也不是很难.强调一下,线段树保存的是
//		线段,并不是单个节点.计算sum的时候要多加
//		注意.

#include <cstring>
#include <algorithm>
#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <vector>
#include <queue>
#define For(x,a,b,c) for (int x = a; x <= b; x += c)
#define Ffor(x,a,b,c) for (int x = a; x >= b; x -= c)
#define cls(x,a) memset(x,a,sizeof(x))
using namespace std;
typedef long long ll;

const double PI = acos(-1.0);
const double EPS = 1e-10;
const int MAX_N = 500 + 8;
const double INF = 1e9;
int N,M;

double pos[MAX_N];

struct Seg{
double l,r,h;
int cov;
Seg(){

}
Seg(double a,double b,double c,int cov):l(a),r(b),h(c),cov(cov){

}

bool operator < (const Seg& other) const{
return h < other.h;
}

}s[MAX_N];

struct InterverTree{

#define ls(x) (x << 1)
#define rs(x) (x << 1 | 1)
int cov[MAX_N<<2];
double sum[MAX_N<<2];

void push_up(int rt,int L,int R){
if (cov[rt])	sum[rt] = pos[R + 1] - pos[L]; // 这条线段已经插入
else if (L == R)	sum[rt] = 0; // 到了叶子节点,并且没有被插入到线段树中
else sum[rt] = sum[ls(rt)] + sum[rs(rt)];
}

void update(int rt,int L,int R,int ql,int qr,int v){
if (ql <= L && R <= qr){
cov[rt] += v;
push_up(rt,L,R); // 如果更新的是单个节点不加这里是会出错的
return ;
}

int M = (L + R) >> 1;
if (ql <= M)	update(ls(rt),L,M,ql,qr,v);
if (M < qr)		update(rs(rt),M+1,R,ql,qr,v);
push_up(rt,L,R);
}
}it;

void print(int n){
for (int i = 1 ;i < n;i ++){
printf("%lf ",pos[i]);
}
puts("");
}

void input(){
int k = 1;
for (int i = 1 ;i <= N;i ++){
double a,b,c,d;
scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
pos[k] = a;
s[k++] = Seg(a,c,b,1);
pos[k] = c;
s[k++] = Seg(a,c,d,-1);
}

sort(pos+1,pos+k);
sort(s+1,s+k);

int m = unique(pos + 1, pos + k) - pos; // 离散化

//print(m);
//it.build(1,1,k-1);
cls(it.sum,0);
cls(it.cov,0);
double res = 0;
for (int i = 1 ;i < k - 1; i++){
int l = lower_bound(pos+1,pos + m,s[i].l) - pos; // 离散化操作
int r = lower_bound(pos+1,pos + m,s[i].r) - pos - 1; // 为什么-1,这里是插入线段的左端点
//cout << l << " " << r << " " << s[i].l << " " << s[i].r << endl;

if (l <= r)	it.update(1,1,k-1,l,r,s[i].cov);
//printf("sum[1] = %lf\n",it.sum[1]);
res += it.sum[1] *(s[i + 1].h - s[i].h);
}

printf("Total explored area: %.2lf\n",res);
puts("");
}

int main(){
int t;
//freopen("1.in","r",stdin);

int kase = 1;
while(scanf("%d",&N)!=EOF && N){
printf("Test case #%d\n",kase++);
input();
}
}