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 }