POJ - 1177 矩形周长并

参考博客：http://www.cnblogs.com/scau20110726/archive/2013/04/13/3018687.html

方法：

把矩形分成横线和竖线去处理，可知是完全相同的操作，我们来讲下怎么算出横线部分，竖线部分就是照搬即可。
将横线保存在一个表中，按横线所处的竖直位置排序（升序），另外每条横线带一个标记值，原矩形的下线为1，上线为-1（对应过去就是插入线段和删除线段）
从低到高扫描横线,没扫到一条横线就能计算出一部分横线值。计算方法是算出现在总区间的被覆盖的长度，然后求出与上一次的总区间的覆盖长度的差（即相减求绝对值），因为每次添加了一条线段，如果没有没有使总区间覆盖长度发生变化，说明这条线段其实在多边形的内部，被覆盖掉了，不能计算，只要能引起总区间长度发生变化的，说明该线段不被覆盖不被包含，要计算。
而竖线部分的做法是一样的，把竖线保存在一个表中，按水平位置排序（升序），每条横线带一个标记值，原矩形的左线为1，右线为-1，然后同样地操作

我用的是第一种方法，思路和上面的博客一样，只有一个地方需要注意。

注意：扫描线的排序方法    return a.h==b.h?a.status>b.status:a.h<b.h;

为什么这么排序？样例 ：2 0 0 1 1 0 1 1 2    自己手动模拟一下这组样例就会明白了

代码：

///#include<bits/stdc++.h>
///#include<unordered_map>
///#include<unordered_set>
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<queue>
#include<bitset>
#include<set>
#include<stack>
#include<map>
#include<list>
#include<new>
#include<vector>
#define MT(a,b) memset(a,b,sizeof(a));
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double pai=acos(-1.0);
const double E=2.718281828459;
const int MAXN=1e5;
const ll mod=1e9+7;
const int INF=0x3f3f3f3f;

struct node
{
    ll l;
    ll r;
    ll h;           
    ll status;
    bool friend operator<(node a,node b)
    {
        ///重要的一步 样例 2 0 0 1 1 0 1 1 2
        return a.h==b.h?a.status>b.status:a.h<b.h;
    }
}line[10005];

ll x[10005];

struct xx
{
    ll s;
    ll sum;
}root[40005];

struct pp
{
    ll x1;
    ll y1;
    ll x2;
    ll y2;
}point[5005];

void pushup(int sign,int l,int r)
{
    if(root[sign].s)
        root[sign].sum=x[r+1]-x[l];
    else if(l==r)
        root[sign].sum=0;
    else
        root[sign].sum=root[sign<<1].sum+root[sign<<1|1].sum;
}

void updateroot(int sign,int l,int r,int a,int b,int f)
{
    if(l==a&&r==b)
    {
        root[sign].s+=f;
        pushup(sign,l,r);
        return ;
    }
    int mid=(l+r)>>1;
    if(b<=mid)
        updateroot(sign<<1,l,mid,a,b,f);
    else if(a>mid)
        updateroot(sign<<1|1,mid+1,r,a,b,f);
    else
    {
        updateroot(sign<<1,l,mid,a,mid,f);
        updateroot(sign<<1|1,mid+1,r,mid+1,b,f);
    }
    pushup(sign,l,r);
}

int main()
{
    ll n;
    ll x1,y1,x2,y2;
    while(scanf("%lld",&n)!=EOF)
    {
        memset(root,0,sizeof(root));
        int sign=0;
        for(int i=1;i<=n;i++)
        {
            scanf("%lld %lld %lld %lld",&x1,&y1,&x2,&y2);
            point[i]=pp{x1,y1,x2,y2};
            line[++sign]=node{x1,x2,y1,1};
            x[sign]=x1;
            line[++sign]=node{x1,x2,y2,-1};
            x[sign]=x2;
        }
        sort(line+1,line+1+sign);
        sort(x+1,x+1+sign);
        int d=unique(x+1,x+1+sign)-(x+1);
        ll len=abs(line[sign].r-line[sign].l);
        ll pre=0;
        for(int i=1;i<sign;i++)
        {
            int a=lower_bound(x+1,x+1+d,line[i].l)-x;
            int b=lower_bound(x+1,x+1+d,line[i].r)-x-1;
            updateroot(1,1,d-1,a,b,line[i].status);
            len+=abs(root[1].sum-pre);
            pre=root[1].sum;
        }
        memset(root,0,sizeof(root));
        sign=0;
        for(int i=1;i<=n;i++)
        {
            x1=point[i].x1;
            y1=point[i].y1;
            x2=point[i].x2;
            y2=point[i].y2;
            line[++sign]=node{y1,y2,x1,1};
            x[sign]=y1;
            line[++sign]=node{y1,y2,x2,-1};
            x[sign]=y2;
        }
        sort(x+1,x+1+sign);
        sort(line+1,line+1+sign);
        d=unique(x+1,x+1+sign)-(x+1);
        len+=abs(line[sign].r-line[sign].l);
        pre=0;
        for(int i=1;i<sign;i++)
        {
            int a=lower_bound(x+1,x+1+d,line[i].l)-x;
            int b=lower_bound(x+1,x+1+d,line[i].r)-x-1;
            updateroot(1,1,d-1,a,b,line[i].status);
            len+=abs(root[1].sum-pre);
            pre=root[1].sum;
        }
        printf("%lld\n",len);
    }
    return 0;
}