【HDU1542】Atlantis

题意

给出n个矩形的左下角和右上角的坐标，计算总的面积（相交部分只算一次）。

分析

线段树扫描线的模板题。

将每个矩形都拆成上下两条线段，然后从下网上扫，当遇到底边时就加上这个区间，遇到顶边时，就减去这个区间。这些都很好理解，但是有一个点我感觉很难受！对于普通线段树，先将区间[1,2]+1，再更新区间[2,3]+1的话，2这个点的值应该是2。但是扫描线来说，2应该是1。因为[1,3]是一条线。我们的解决办法是，在线段树[L,R]的区间内储存[L,R+1]的值。

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

 using namespace std;
 const int maxn=;
 struct Node{
     double l,r,h;
     int state;
     bool operator <(const Node & rhs)const{
         return h<rhs.h;
     }
     Node (){}
     Node (double l,double r,double h,int state):l(l),r(r),h(h),state(state){}
 }node[maxn];
 double v[maxn];
 int n,kase;
 double x1,y1,x2,y2,ans;
 double sumv[*maxn];
 int cover[*maxn];

 void maintain(int o,int L,int R){
     if(cover[o])
         sumv[o]=v[R+]-v[L];
     else if(L==R)
         sumv[o]=;
     else
         sumv[o]=sumv[*o]+sumv[*o+];
 }
 int ql,qr,vv;
 void update(int o,int L,int R){
     if(ql<=L&&qr>=R){
         cover[o]+=vv;
         maintain(o,L,R);
         return ;
     }
     int M=L+(R-L)/;
     if(ql<=M)update(*o,L,M);
     if(qr>M)update(*o+,M+,R);
     maintain(o,L,R);
 }
 int main(){
     kase=;
     while(scanf("%d",&n)!=EOF&&n){
         int sz=;
         ++kase;
         for(int i=;i<=n;i++){
             scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
             node[++sz].l=x1,node[sz].r=x2,node[sz].h=y1,node[sz].state=,v[sz]=x1;
             node[++sz].l=x1,node[sz].r=x2,node[sz].h=y2,node[sz].state=-,v[sz]=x2;
         }
         sort(node+,node++sz);
         sort(v+,v++sz);
         int N=unique(v+,v++sz)-v-;
         memset(sumv,,sizeof(sumv));
         memset(cover,,sizeof(cover));
         ans=;
         for(int i=;i<sz;i++){
             ql=lower_bound(v+,v++N,node[i].l)-v,qr=lower_bound(v+,v++N,node[i].r)-v-,vv=node[i].state;
             update(,,N);
             ans+=sumv[]*(node[i+].h-node[i].h);
         }
         printf("Test case #%d\n",kase);
         printf("Total explored area: %.2f\n",ans);
         printf("\n");
     }
 return ;
 }