【POJ1177】Picture-线段树+离散化+扫描线

测试地址：Picture
题目大意：平面上有 n <script type="math/tex" id="MathJax-Element-180">n</script>个边和坐标轴平行的矩形，求这些矩形组成的图形的轮廓长度。
做法：本题需要用到线段树+离散化+扫描线。
求多矩形轮廓周长也是线段树的经典应用之一，由于我不会插图，在这里讲不太清楚，具体请看这篇文章。
处理的方法和POJ1151（求矩形面积并）比较相似，这里就不赘述了。
以下是本人代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
int n,y[5010][2],x1[5010],y1[5010],x2[5010],y2[5010];
struct forsort
{
    int val;
    int id,idd;
}f[10010];
struct oper
{
    int x,y1,y2,add;
}p[10010];
struct segnode
{
    int cov,sum,num,lb,rb;
}seg[40010];

bool cmp(forsort a,forsort b)
{
    return a.val<b.val;
}

bool cmp0(oper a,oper b)
{
    return a.x<b.x;
}

void pushup(int no,int l,int r)
{
    if (seg[no].cov) seg[no].sum=f[r+1].val-f[l].val;
    else if (l==r) seg[no].sum=0.0;
         else seg[no].sum=seg[no<<1].sum+seg[no<<1|1].sum;
    if (seg[no].cov) seg[no].lb=seg[no].rb=seg[no].num=1;
    else if (l==r) seg[no].lb=seg[no].rb=seg[no].num=0;
         else
         {
            if (seg[no<<1].rb&&seg[no<<1|1].lb) seg[no].num=seg[no<<1].num+seg[no<<1|1].num-1;
            else seg[no].num=seg[no<<1].num+seg[no<<1|1].num;
            seg[no].lb=seg[no<<1].lb,seg[no].rb=seg[no<<1|1].rb;
         }
}

void buildtree(int no,int l,int r)
{
    seg[no].cov=seg[no].sum=seg[no].num=seg[no].lb=seg[no].rb=0;
    if (l==r) return;
    int mid=(l+r)>>1;
    buildtree(no<<1,l,mid);
    buildtree(no<<1|1,mid+1,r); 
}

void modify(int no,int l,int r,int s,int t,int d)
{
    if (l>=s&&r<=t)
    {
        seg[no].cov+=d;
        pushup(no,l,r);
        return;
    }
    int mid=(l+r)>>1;
    if (s<=mid) modify(no<<1,l,mid,s,t,d);
    if (t>mid) modify(no<<1|1,mid+1,r,s,t,d);
    pushup(no,l,r);
}

int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d%d%d",&x1[i],&y1[i],&x2[i],&y2[i]);
        f[2*i-1].val=y1[i],f[2*i-1].id=i,f[2*i-1].idd=0;
        f[2*i].val=y2[i],f[2*i].id=i,f[2*i].idd=1;
    }
    sort(f+1,f+2*n+1,cmp);
    for(int i=1;i<=2*n;i++)
        y[f[i].id][f[i].idd]=i;
    for(int i=1;i<=n;i++)
    {
        p[2*i-1].x=x1[i],p[2*i].x=x2[i];
        p[2*i-1].y1=p[2*i].y1=y[i][0];
        p[2*i-1].y2=p[2*i].y2=y[i][1];
        p[2*i-1].add=1,p[2*i].add=-1;
    }
    sort(p+1,p+2*n+1,cmp0);

    buildtree(1,1,2*n);
    int ans=0,last=0;
    for(int i=1;i<=2*n;i++)
    {
        int num=seg[1].num;
        modify(1,1,2*n,p[i].y1,p[i].y2-1,p[i].add);
        ans+=abs(seg[1].sum-last);
        if (i>1) ans+=num*2*(p[i].x-p[i-1].x);
        last=seg[1].sum;
    }
    printf("%d",ans);

    return 0; 
}