http://acm.hdu.edu.cn/showproblem.php?pid=1542
1 #include <iostream> 2 #include <cstdio> 3 #include <cstring> 4 #include <algorithm> 5 #define lson l,mid,i<<1 6 #define rson mid+1,r,i<<1|1 7 using namespace std; 8 const int Ni = 210; 9 int n; 10 struct node { 11 double l,r,h; 12 int s; 13 bool operator < (const node & a) const 14 { 15 return h<a.h; 16 } 17 }arr[Ni*2]; 18 double sum[Ni*3]; 19 int col[Ni*3]; 20 double x[Ni*2]; 21 void pushup(int l,int r,int i) 22 { 23 if(col[i]) sum[i]=x[r+1]-x[l];//这一段被覆盖，一个点代表这个点后面的这一段 24 else if(l==r) sum[i]=0;//没有覆盖的叶子结点 25 else sum[i]=sum[i<<1]+sum[i<<1|1];//不是叶子等于左右儿子之和 26 } 27 void update(double ql,double qr,int s,int l=0,int r=n,int i=1) 28 { 29 if(ql<=x[l]&&x[r]<qr)//[ql,qr)这一段 30 { 31 col[i]+=s; 32 pushup(l,r,i); 33 return; 34 } 35 int mid=(l+r)>>1; 36 if(ql<=x[mid]) update(ql,qr,s,lson); 37 if(qr>x[mid+1]) update(ql,qr,s,rson); 38 pushup(l,r,i); 39 } 40 int main() 41 { 42 int cs=1; 43 while(cin>>n,n) 44 { 45 int m=0; 46 for(int i=0;i<n;i++) 47 { 48 double a,b,c,d; 49 scanf("%lf%lf%lf%lf",&a,&b,&c,&d); 50 x[m]=a; 51 arr[m].l=a;arr[m].r=c;arr[m].h=b;arr[m++].s=1; 52 x[m]=c; 53 arr[m].l=a;arr[m].r=c;arr[m].h=d;arr[m++].s=-1; 54 } 55 sort(x,x+m); 56 sort(arr,arr+m); 57 int q=0; 58 for(int i=1;i<m;i++) 59 { 60 if(x[q]!=x[i]) x[++q]=x[i]; 61 } 62 n=q; 63 memset(col,0,sizeof(col)); 64 memset(sum,0,sizeof(sum)); 65 double ret=0; 66 for(int i=0;i<m-1;i++) 67 { 68 update(arr[i].l,arr[i].r,arr[i].s); 69 ret+=sum[1]*(arr[i+1].h-arr[i].h); 70 } 71 printf("Test case #%d\n",cs++); 72 printf("Total explored area: %.2lf\n\n",ret); 73 } 74 return 0; 75 }