HDU 1542 Atlantics 线段树+离散化扫描

将 x 轴上的点进行离散化，扫描线沿着 y 轴向上扫描

每次添加一条边不断找到当前状态有效边的长度 ， 根据这个长度和下一条边形成的高度差得到一块合法的矩形的面积

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>

 using namespace std;

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

 struct Seg{//线段
     double x1,x2,y;
     int d;
     Seg(){}
     Seg(double x1 , double x2 , double y , int d):x1(x1),x2(x2),y(y),d(d){}
     bool operator<(const Seg &a)const{
         return y<a.y;
     }
 }s[MAX];

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

 void update(int o , int l , int r , int s , int t , int d){
     if(s<=l && r<=t){//该区间是当前扫描线段的一部分,则该区间下底边总长以及上下底边个数差更新
         flag[o] += d;//更新底边相差差个数
         updateCur(o , l , r);//更新底边长
         return;
     }
     int mid = (l + r)>> , ls = o<< , rs = o<<|;
     if(mid >= s) update(ls , l , mid , s , t , d);
     if(mid+ <= t) update(rs , mid+ , r , s , t , d);
     updateCur(o , l , r);
 }

 int bin_search(double key , int n){
     int l =  , r = n-;
     while(l<=r){
         int mid = (l+r) >>;
         if(x[mid] == key) return mid;
         if(x[mid] > key) r=mid-;
         else l=mid+;
     }
     return -;
 }

 int main()
 {
    // freopen("a.in" , "r" , stdin);
     int n,cas=;
     double x1,x2,y1,y2;
     while(scanf("%d" , &n) , n){
         int k=;
         for(int i=;i<n;++i){
             cin>>x1>>y1>>x2>>y2;
             x[k]=x1;
             s[k++]=Seg(x1,x2,y1,);
             x[k]=x2;
             s[k++]=Seg(x1,x2,y2,-);
         }
         sort(x,x+k);
         sort(s,s+k);
         int t = unique(x , x+k) - x;
         double ans=;

         for(int i=;i<k;++i){//扫描线段
             int pos1 = bin_search(s[i].x1 , t);
             /*
             这里因为是把它理解成合法线段，所以最右端第 i 个点 ， 对应第 i - 1条边
             同样道理，在updateCur中，计算有效线段长度时 ， 需要x[r+1]
             */
             int pos2 = bin_search(s[i].x2 , t)-;
             update( ,  , t- , pos1 , pos2 , s[i].d);//扫描线段时更新底边长度和底边相差个数
             ans += sum[]*(s[i+].y-s[i].y);//新增加面积
         }
         printf("Test case #%d\nTotal explored area: %.2lf\n\n",++cas,ans);
     }
     return ;
 }