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

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

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

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

;
}
void down(int w)
{
    if(s[w]!=-1)
    {
        s[w<<1] = s[w<<1|1] = s[w];
        s[w] = -1;
    }
}
void build(int l,int r,int w)
{
    s[w] = 0;
    if(l==r)
    {
        return ;
    }
    int m = (l+r)>>1;
    build(l,m,w<<1);
    build(m+1,r,w<<1|1);
    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)>>1;
    int k = 0;
    if(a<=m)
    k = update(a,b,d,l,m,w<<1);
    if(b>m)
    k = update(a,b,d,m+1,r,w<<1|1);
    up(w);
    return k;
}
int main()
{
    int n,w,h,i,j;
    scanf("%d%d%d",&n,&w,&h);
    int g = 0;
    for(i = 1;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+1,p+g+1);
    sort(a+1,a+g+1);
    int o = 0;
    f[a[1]] = ++o;
    for(i = 2; i <= g;i++)
    {
        if(a[i]!=a[i-1])
        {
            f[a[i]] = ++o;
            val[o] =a[i];
        }
    }
    build(1,o,1);
    for(i = 1; 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>0)
            {
                int k = update(l,r,p[i].d,1,o,1);
                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],1,o,1);
            }
            //fresh(l,r,1,o,1);
        }
    }
    sum[0] =(LL)w*h;
    for(i = 0 ;i <= n; i++)
    ans[i] = sum[i];
    for(i = 0; i <= n; i++)
    {
        for(j = 0 ;j < ed[i].size() ; j++)
        {
            int v = ed[i][j];
            ans[i]-=sum[v];
        }
    }
    sort(ans,ans+n+1);
    for(i = 0 ;i < n; i++)
    printf("%I64d ",ans[i]);
    printf("%I64d\n",ans[n]);
    return 0;
}