hdu1542

 1 /*
 2 题意：矩形面积并
 3 
 4 分析：离散+线段树+扫描线；
 5 
 6 细节：首先线段记录的信息，len[]表示区间被覆盖的长度,cov[]表示当前区间是否被完全覆盖
 7 其次，线段树的叶子节点[l,r]{l==r}的长度是？区间[l,r]的长度是LX[r]-LX[l]{LX是离散后保存数据的地方}
 8 那叶子节点[l,l]的长度不就变成0，显然这是有问题的
 9 
10 线段树的每一个节点表示一段区间，[l,r]该区间表示LX[r+1]-LX[l]的长度 
11 1___2___3___4___5离散后的状况
12   1   2   3   4  线段树中的每一个节点
13    
14  
15 
16 */
17 #include<cstdio>
18 #include<cstring>
19 #include<cmath>
20 #include<iostream>
21 #include<algorithm>
22 #include<cstdlib>
23 #include<vector>
24 #include<string>
25 #define Find(i) lower_bound(LX.begin(),LX.begin()+n1,i)-LX.begin()
26 #define lson l,m,rt<<1
27 #define rson m+1,r,rt<<1|1
28 using namespace std;
29 const int N=200+10;
30 struct Edge{
31     double hi,ls,rs;
32     int s;
33     Edge(){}
34     Edge(double a,double b,double c,int e):hi(a),ls(b),rs(c),s(e){}
35     bool operator < (const Edge &p)const{
36         return hi<p.hi;
37     }
38 };
39 vector<Edge> E;
40 vector<double> LX;
41 double len[N<<2];
42 int cov[N<<2];//表示区间是否被完全覆盖 
43 int n,n1;
44 void pushup(int l,int r,int rt){
45     if (cov[rt]>=1){
46         len[rt]=(LX[r+1]-LX[l]);
47     }else if (l==r) len[rt]=0;
48     else len[rt]=len[rt<<1]+len[rt<<1|1]; 
49 }
50 void update(int L,int R,int v,int l,int r,int rt){
51     if (L<=l && r<=R){
52         cov[rt]+=v;
53         if (cov[rt]>=1) len[rt]=(LX[r+1]-LX[l]);
54         else if (l==r) len[rt]=0;
55         else len[rt]=len[rt<<1]+len[rt<<1|1];
56         return ;
57     }
58     int m=(l+r)>>1;
59     if (L<=m) update(L,R,v,lson);
60     if (m< R) update(L,R,v,rson);
61     pushup(l,r,rt);
62 }
63 
64 void work(){
65     double ret=0;
66     memset(len,0,sizeof(len));
67     memset(cov,0,sizeof(cov));
68     for (int i=0;i<E.size();i++){
69         int l=Find(E[i].ls),r=Find(E[i].rs)-1;
70         if (l<=r) update(l,r,E[i].s,0,n1-1,1);
71         if (E[i].hi!=E[i+1].hi){
72             ret+=len[1]*(E[i+1].hi-E[i].hi);
73         }
74     }
75     printf("Total explored area: %.2lf\n",ret);
76 
77 }
78 int main(){
79     int cas=0;
80     while (~scanf("%d",&n),n){
81         E.clear();LX.clear();
82         for (int i=0;i<n;i++){
83             double x,y,x2,y2;
84             scanf("%lf%lf%lf%lf",&x,&y,&x2,&y2);
85             E.push_back(Edge(y,x,x2,1));
86             E.push_back(Edge(y2,x,x2,-1));
87             LX.push_back(x);LX.push_back(x2);
88         }    
89         sort(E.begin(),E.end());
90         sort(LX.begin(),LX.end());
91         n1=unique(LX.begin(),LX.end())-LX.begin();
92         printf("Test case #%d\n",++cas);    
93         work();
94         printf("\n");
95     }
96     return 0;
97 }