hdu 1542 Atlantis

求矩形的面积之和。

线段树+离散话+扫描线

 #include<iostream>
 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<string>
 #include<queue>
 #include<algorithm>
 #include<map>
 #include<iomanip>
 #define INF 99999999
 using namespace std;

 const int MAX=+;
 int mark[MAX<<];//记录某个区间的下底边个数
 double sum[MAX<<];//记录某个区间的下底边总长度
 double lisan[MAX];//对x进行离散化,否则x为浮点数且很大无法进行线段树

 //以横坐标作为线段(区间),对横坐标线段进行扫描
 //扫描的作用是每次更新下底边总长度和下底边个数,增加新面积
 struct seg{//线段
     double l,r,h;
     int d;
     seg(){}
     seg(double x1,double x2,double H,int c):l(x1),r(x2),h(H),d(c){}
     bool operator<(const seg &a)const{
         return h<a.h;
     }
 }s[MAX];

 void Upfather(int n,int left,int right){
     if(mark[n])sum[n]=lisan[right+]-lisan[left];//表示该区间整个线段长度可以作为底边
     else if(left == right)sum[n]=;//叶子结点则底边长度为0(区间内线段长度为0)
     else sum[n]=sum[n<<]+sum[n<<|];
 }

 void Update(int L,int R,int d,int n,int left,int right){
     if(L<=left && right<=R){//该区间是当前扫描线段的一部分,则该区间下底边总长以及上下底边个数差更新
         mark[n]+=d;//更新底边相差差个数
         Upfather(n,left,right);//更新底边长
         return;
     }
     int mid=left+right>>;
     if(L<=mid)Update(L,R,d,n<<,left,mid);
     if(R>mid)Update(L,R,d,n<<|,mid+,right);
     Upfather(n,left,right);
 }

 int search(double key,double* x,int n){
     int left=,right=n-;
     while(left<=right){
         int mid=left+right>>;
         if(x[mid] == key)return mid;
         if(x[mid]>key)right=mid-;
         else left=mid+;
     }
     return -;
 }

 int main(){
     int n,num=;
     double x1,x2,y1,y2;
     while(cin>>n,n){
         int k=;
         for(int i=;i<n;++i){
             cin>>x1>>y1>>x2>>y2;
             lisan[k]=x1;
             s[k++]=seg(x1,x2,y1,);
             lisan[k]=x2;
             s[k++]=seg(x1,x2,y2,-);
         }
         sort(lisan,lisan+k);
         sort(s,s+k);
         int m=;
         for(int i=;i<k;++i)//去重复端点
             if(lisan[i] != lisan[i-])lisan[m++]=lisan[i];
         double ans=;
         //memset(mark,0,sizeof mark);
         //memset(sum,0,sizeof sum);//如果下面是i<k-1,则要初始化,因为如果对第k-1条,线段扫描时会使得mark,sum为0才不用初始化的
         for(int i=;i<k;++i){//扫描线段
             int L=search(s[i].l,lisan,m);
             int R=search(s[i].r,lisan,m)-;
             Update(L,R,s[i].d,,,m-);//扫描线段时更新底边长度和底边相差个数
             ans+=sum[]*(s[i+].h-s[i].h);//新增加面积
         }
         printf("Test case #%d\nTotal explored area: %.2lf\n\n",++num,ans);
     }
     return ;
 }