codevs 3044 矩形面积求并 || hdu 1542

这个线段树的作用其实是维护一组（1维 平面(?) 上的）线段覆盖的区域的总长度，支持加入/删除一条线段。

线段树只能维护整数下标，因此要离散化。

也可以理解为将每一条处理的线段分解为一些小线段，要求每一条要处理的线段都能这么分解

注意端点，线段树维护的是线段，而查询是端点，可能需要稍微变一下

具体的query方法好像类似标记永久化(?)

感觉不优美啊。。。想了一会并没有想到更好的方法。。。。

codevs:

错误记录：

1.打成unordered_map<int,int>

2.70、75行两个tlen误为len

3.要求保留两位小数

 #include<cstdio>
 #include<algorithm>
 #include<tr1/unordered_map>
 using namespace std;
 using namespace tr1;
 struct Q
 {
     int x1,x2;double y1,y2,tx1,tx2;
 }q[];
 struct QQ
 {
     int l,r;double h;int fl;
     friend bool operator<(const QQ &a,const QQ &b)    {return a.h<b.h;}
 }qq[];
 double tt[];
 unordered_map<double,int> ma;
 int len,tlen;double ans;
 int n;
 namespace SegT
 {
 #define lc (num<<1)
 #define rc (num<<1|1)
 #define mid (l+((r-l)>>1))
 //[l,r]表示点l到点r+1之间的线段
     int cover[];double dat[];
     int L,R,x;
     void build(int l,int r,int num)
     {
         cover[num]=;dat[num]=;
         if(l!=r)    build(l,mid,lc),build(mid+,r,rc);
     }
     void addx(int l,int r,int num)
     {
         if(L<=l&&r<=R)
         {
             cover[num]+=x;
             if(cover[num])    dat[num]=tt[r+]-tt[l];
             else    dat[num]=dat[lc]+dat[rc];
             return;
         }
         if(L<=mid)    addx(l,mid,lc);
         if(mid<R)    addx(mid+,r,rc);
         if(cover[num])    dat[num]=tt[r+]-tt[l];
         else    dat[num]=dat[lc]+dat[rc];
     }
 #undef lc
 #undef rc
 #undef mid
 }
 int main()
 {
     int i;
     while()
     {
         scanf("%d",&n);len=tlen=;ans=;ma.clear();
         if(n==)    break;
         for(i=;i<=n;i++)
         {
             scanf("%lf%lf%lf%lf",&q[i].tx1,&q[i].y1,&q[i].tx2,&q[i].y2);
             tt[++tlen]=q[i].tx1;tt[++tlen]=q[i].tx2;
         }
         sort(tt+,tt+tlen+);tlen=unique(tt+,tt+tlen+)-tt-;
         for(i=;i<=tlen;i++)    ma[tt[i]]=i;
         for(i=;i<=n;i++)    q[i].x1=ma[q[i].tx1],q[i].x2=ma[q[i].tx2]-;
         for(i=;i<=n;i++)
         {
             qq[++len]=(QQ){q[i].x1,q[i].x2,q[i].y1,};
             qq[++len]=(QQ){q[i].x1,q[i].x2,q[i].y2,-};
         }
         sort(qq+,qq+len+);SegT::build(,tlen-,);
         for(i=;i<=len;i++)
         {
             ans+=SegT::dat[]*(qq[i].h-qq[i-].h);
             SegT::L=qq[i].l;SegT::R=qq[i].r;SegT::x=qq[i].fl;
             SegT::addx(,tlen-,);
         }
         printf("%.2lf\n",ans);
     }
     return ;
 }

hdu

 #include<cstdio>
 #include<algorithm>
 #include<tr1/unordered_map>
 using namespace std;
 using namespace tr1;
 struct Q
 {
     int x1,x2;double y1,y2,tx1,tx2;
 }q[];
 struct QQ
 {
     int l,r;double h;int fl;
     friend bool operator<(const QQ &a,const QQ &b)    {return a.h<b.h;}
 }qq[];
 double tt[];
 unordered_map<double,int> ma;
 int len,tlen;double ans;
 int n;
 namespace SegT
 {
 #define lc (num<<1)
 #define rc (num<<1|1)
 #define mid (l+((r-l)>>1))
 //[l,r]表示点l到点r+1之间的线段
     int cover[];double dat[];
     int L,R,x;
     void build(int l,int r,int num)
     {
         cover[num]=;dat[num]=;
         if(l!=r)    build(l,mid,lc),build(mid+,r,rc);
     }
     void addx(int l,int r,int num)
     {
         if(L<=l&&r<=R)
         {
             cover[num]+=x;
             if(cover[num])    dat[num]=tt[r+]-tt[l];
             else    dat[num]=dat[lc]+dat[rc];
             return;
         }
         if(L<=mid)    addx(l,mid,lc);
         if(mid<R)    addx(mid+,r,rc);
         if(cover[num])    dat[num]=tt[r+]-tt[l];
         else    dat[num]=dat[lc]+dat[rc];
     }
 #undef lc
 #undef rc
 #undef mid
 }
 int main()
 {
     int i,T=;
     while()
     {
         scanf("%d",&n);len=tlen=;ans=;ma.clear();
         if(n==)    break;
         for(i=;i<=n;i++)
         {
             scanf("%lf%lf%lf%lf",&q[i].tx1,&q[i].y1,&q[i].tx2,&q[i].y2);
             tt[++tlen]=q[i].tx1;tt[++tlen]=q[i].tx2;
         }
         sort(tt+,tt+tlen+);tlen=unique(tt+,tt+tlen+)-tt-;
         for(i=;i<=tlen;i++)    ma[tt[i]]=i;
         for(i=;i<=n;i++)    q[i].x1=ma[q[i].tx1],q[i].x2=ma[q[i].tx2]-;
         for(i=;i<=n;i++)
         {
             qq[++len]=(QQ){q[i].x1,q[i].x2,q[i].y1,};
             qq[++len]=(QQ){q[i].x1,q[i].x2,q[i].y2,-};
         }
         sort(qq+,qq+len+);SegT::build(,tlen-,);
         for(i=;i<=len;i++)
         {
             ans+=SegT::dat[]*(qq[i].h-qq[i-].h);
             SegT::L=qq[i].l;SegT::R=qq[i].r;SegT::x=qq[i].fl;
             SegT::addx(,tlen-,);
         }
         printf("Test case #%d\nTotal explored area: %.2lf\n\n",++T,ans);
     }
     return ;
 }