【USACO5.5.1】Picture矩形周长

问题描述

N(N<5000) 张矩形的海报，照片和其他同样形状的图片贴在墙上。它们的边都是垂直的或水平的。每个矩形可以部分或者全部覆盖其他矩形。所有的矩形组成的集合的轮廓称为周长。写一个程序计算周长。
图 1 是一个有 7 个矩形的例子：

对应的轮廓为图 2 所示的所有线段的集合：

所有矩形的顶点坐标均为整数。所有的坐标都在 [-10000，10000] 的范围内，并且任何一个矩形面积都为整数。结果的值可能需要 32 位有符号整数表示。

输入格式

第1行: N，张贴在墙上的矩形的数目。 第 2..N+1行 接下来的N行中，每行都有两个点的坐标，分别是矩形的左下角坐标和右上角坐标。每一个坐标由 X 坐标和 Y 坐标组成。

输出格式

只有一行，为一个非负整数，表示输入数据中所有矩形集合的轮廓长度。

样例输入

7
-15 0 5 10
-5 8 20 25
15 -4 24 14
0 -6 16 4
2 15 10 22
30 10 36 20
34 0 40 16

样例输出

228

---

题解

线段树 —— 扫描线
显然整个图形的轮廓线条一定是竖直或水平的。如果我们求出了竖直的，水平的也可用类似方法。
我们先记录下每条边是矩形的左边还是右边，再把边按横坐标排序。
维护一颗线段树，记录下每个点被覆盖的情况。初始时没有点被覆盖。
接下来顺序扫描每一条边。

情况一：

若讨论到一条左边，将左边表示的纵坐标区间覆盖一层，区间被覆盖总长的增加量即为这条边对答案的贡献。

情况二：

若讨论到一条右边，将右边表示的纵坐标区间去掉一层，区间被覆盖总长的减少量即为这条边对答案的贡献。

横边同理。
似乎有些网站上的题加强了数据？离散化即可。

这道题就解决了！解决了？没有！
此题还有一坑点。输入这组数据试试：
2 0 0 4 4 0 4 4 8
答案是32？显然是24 啊！
需要注意，对边进行排序时，若横（纵）坐标相同，左边（讨论到时覆盖一层的边）应排在前面，否则就会尴尬。
具体见代码。

---

代码

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
struct idata{
    int x,y,q;
    bool p;
}h[2][10005];
bool cmp(idata a,idata b)
{return a.q==b.q?a.p>b.p:a.q<b.q;}
const int Q=800000;
int n,lazy[Q],ls[Q],rs[Q],co[Q],maxr=10000,maxl=-10000,t,tot;
void pd(int now)
{
    if(!ls[now])ls[now]=++tot;
    if(!rs[now])rs[now]=++tot;
    int v=lazy[now];
    co[ls[now]]+=v;
    lazy[ls[now]]+=v;
    co[rs[now]]+=v;
    lazy[rs[now]]+=v;
    lazy[now]=0;
}
void xiu(int now,int l,int r,int x,int y,int v)
{
    if(x<=l&&y>=r)
    {
        co[now]+=v;
        lazy[now]+=v;
        return;
    }
    if(lazy[now]!=0)pd(now);
    int mid=(l+r)>>1;
    if(x<mid){
        if(!ls[now])ls[now]=++tot;
        xiu(ls[now],l,mid,x,y,v);
    }
    if(y>mid){
        if(!rs[now])rs[now]=++tot;
        xiu(rs[now],mid,r,x,y,v);
    }
}
int gs(int now,int l,int r)
{
    if(r-l==1||((!ls[now])&&(!rs[now])))return (co[now]>0)*(r-l);
    if(lazy[now]!=0)pd(now);
    int mid=(l+r)>>1,temp=0;
    if(ls[now])temp+=gs(ls[now],l,mid);
    if(rs[now])temp+=gs(rs[now],mid,r);
    return temp;
}
int main()
{
    int i,ans=0,a,b,c,d,last,haha;
    scanf("%d",&n);
    for(i=1;i<=n;i++)
    {
        scanf("%d%d%d%d",&a,&b,&c,&d);
        h[0][2*i-1].p=h[1][2*i-1].p=false;
        h[0][2*i].p=h[1][2*i].p=true;
        h[0][2*i-1].x=h[0][2*i].x=a;
        h[0][2*i-1].y=h[0][2*i].y=c;
        h[0][2*i].q=b,h[0][2*i-1].q=d;
        h[1][2*i].x=h[1][2*i-1].x=b;
        h[1][2*i].y=h[1][2*i-1].y=d;
        h[1][2*i].q=a,h[1][2*i-1].q=c;
    }
    n*=2;
    sort(h[0]+1,h[0]+n+1,cmp);
    sort(h[1]+1,h[1]+n+1,cmp);
    for(t=0;t<2;t++)
    {
        memset(ls,0,sizeof(ls));
        memset(rs,0,sizeof(rs));
        memset(co,0,sizeof(co));
        memset(lazy,0,sizeof(lazy));
        tot=1;
        last=0;
        for(i=1;i<=n;i++)
        {
            if(h[t][i].p)xiu(1,maxl,maxr,h[t][i].x,h[t][i].y,1);
            else xiu(1,maxl,maxr,h[t][i].x,h[t][i].y,-1);
            haha=gs(1,maxl,maxr);
            ans+=abs(haha-last),last=haha;

        }
    }
    printf("%d",ans);
    return 0;
}

最后，感谢机智的何老板提供的反例数据。