sgu316Kalevich Strikes Back（线段树+扫描线）

做法：总体想法是求出一个矩形的面积以及它所包含的矩形，然后用自己的面积减掉所包含的。主要问题是怎样求解它所包含的矩形。

因为是没有相交点的，可以利用扫描线的方法去做，类似染色，当前段如果是x色，也就是第x个矩形，那么再把他染成y颜色时，说明x包含y，而当扫到y的上边时，这一段又恢复到x色。标记一下被包含的矩形，记录所包含的矩形。

因为会有恢复染色操作，up需要时时更新，左儿子和右儿子一样颜色时就可以合并为一段。

 ;
 }
 void down(int w)
 {
     if(s[w]!=-)
     {
         s[w<<] = s[w<<|] = s[w];
         s[w] = -;
     }
 }
 void build(int l,int r,int w)
 {
     s[w] = ;
     if(l==r)
     {
         return ;
     }
     int m = (l+r)>>;
     build(l,m,w<<);
     build(m+,r,w<<|);
     up(w);
 }
 int update(int a,int b,int d,int l,int r,int w)
 {
     if(a<=l&&b>=r)
     {
         int k = s[w];
        // cout<<l<<" "<<r<<" "<<s[w]<<endl;
         s[w] = d;
         return k;
     }
     down(w);
     int m  = (l+r)>>;
     int k = ;
     if(a<=m)
     k = update(a,b,d,l,m,w<<);
     if(b>m)
     k = update(a,b,d,m+,r,w<<|);
     up(w);
     return k;
 }
 int main()
 {
     int n,w,h,i,j;
     scanf("%d%d%d",&n,&w,&h);
     int g = ;
     for(i = ;i <= n; i++)
     {
         int x1,x2,y1,y2;
         scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
         p[++g] = node(x1,x2,min(y1,y2),i);
         a[g] = x1;
         p[++g] = node(x1,x2,max(y1,y2),-i);
         a[g] = x2;
         sum[i] = (LL)abs(x1-x2)*abs(y1-y2);
     }
     sort(p+,p+g+);
     sort(a+,a+g+);
     int o = ;
     f[a[]] = ++o;
     for(i = ; i <= g;i++)
     {
         if(a[i]!=a[i-])
         {
             f[a[i]] = ++o;
             val[o] =a[i];
         }
     }
     build(,o,);
     for(i = ; i <= g;i++)
     {
         int l = min(f[p[i].x1],f[p[i].x2]);
         int r = max(f[p[i].x2],f[p[i].x1]);
         //cout<<l<<" "<<r<<endl;
         if(l<=r)
         {
             if(p[i].d>)
             {
                 int k = update(l,r,p[i].d,,o,);
                 ed[k].push_back(p[i].d);
                 pre[p[i].d] = k;
                 //cout<<k<<" "<<p[i].d<<endl;
             }
             else
             {
                 update(l,r,pre[-p[i].d],,o,);
             }
             //fresh(l,r,1,o,1);
         }
     }
     sum[] =(LL)w*h;
     for(i =  ;i <= n; i++)
     ans[i] = sum[i];
     for(i = ; i <= n; i++)
     {
         for(j =  ;j < ed[i].size() ; j++)
         {
             int v = ed[i][j];
             ans[i]-=sum[v];
         }
     }
     sort(ans,ans+n+);
     for(i =  ;i < n; i++)
     printf("%I64d ",ans[i]);
     printf("%I64d\n",ans[n]);
     return ;
 }