N - Picture - poj 1177（扫描线求周长）

题意：求周长的，把矩形先进行融合后的周长，包括内周长

分析：刚看的时候感觉会跟棘手，让人无从下手，不过学过扫描线之后相信就很简单了吧（扫描线的模板- -），还是不说了，下面是一精确图，可以拿来调试数据

*****************************************************************************************************************

#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;

#define Lson r<<1
#define Rson r<<1|1

const int MAXN = 2e5+5;
const int oo = 1e9+7;

struct stgmentTree
{///len 保存区间包含的边的长度，cover 保存区间被覆盖的次数, sum保存有多少个不相连的区间段
    int L, R, len, cover, sum;
    bool lb, rb;///左边和右边是否被覆盖
    int mid(){return (R+L)>>1;}
}a[MAXN<<2];
///dir 等于 1 的时候代表左边，-1时候代表右边，右边会抵消左边
struct Point{int x, y1, y2, dir;}ege[MAXN];
int Hash[MAXN], nh;///保存离散化后的数据，nh表示元素个数

bool cmp(Point n1, Point n2)
{///把边按照x的值从左往右排序
    return n1.x < n2.x;
}
///查找一条边的长度, x1 < x2
int findSegLen(int x1, int x2)
{///用下标可以直接检索值
    return Hash[x2] - Hash[x1];
}
void buildTree(int r, int L, int R)
{
    a[r].L = L, a[r].R = R, a[r].cover=0;

if(L == R-1)return ;

buildTree(Lson, L, a[r].mid());
    buildTree(Rson, a[r].mid(), R);
}
void pushUp(int r)///合并操作
{///需要先注意本区间是否被覆盖
    if(a[r].cover != 0)
    {
        a[r].len= findSegLen( a[r].L, a[r].R );
        a[r].sum = a[r].lb = a[r].rb = 1;
    }
    else if(a[r].L == a[r].R - 1)
    {
        a[r].len = 0;
        a[r].sum = a[r].lb = a[r].rb = 0;
    }
    else
    {///如果本区间未被覆盖，并且不为叶子节点，那么就等于子区间的和
        a[r].len = a[Lson].len + a[Rson].len;

a[r].lb = a[Lson].lb, a[r].rb = a[Rson].rb;
        a[r].sum = a[Lson].sum + a[Rson].sum - (a[Lson].rb & a[Rson].lb);
    }

}
void upData(int r, int L, int R, int dir)
{
    if( a[r].L == L && a[r].R == R )
    {
        a[r].cover += dir;
        pushUp(r);

return ;
    }

if(R <= a[r].mid())
        upData(Lson, L, R, dir);
    else if(L >= a[r].mid())
        upData(Rson, L, R, dir);
    else
    {
        upData(Lson, L, a[r].mid(), dir);
        upData(Rson, a[r].mid(), R, dir);
    }

pushUp(r);
}

int main()
{
    int N;

while(scanf("%d", &N) != EOF)
    {
        int i, x1, x2, y1, y2, k=0, C=0, pre=0; nh = 0;

for(i=0; i<N; i++)
        {
            scanf("%d%d%d%d", &x1, &y1, &x2, &y2);
            ege[k].x=x1, ege[k].y1=y1, ege[k].y2=y2, ege[k++].dir=1;///左边y
            ege[k].x=x2, ege[k].y1=y1, ege[k].y2=y2, ege[k++].dir=-1;///右边y
            Hash[nh++] = y1, Hash[nh++] = y2;
        }

sort(Hash, Hash+nh);///排序去重复
        nh = unique(Hash, Hash+nh) - Hash;
        ///离散化后的数据是从0下标开始的
        buildTree(1, 0, nh-1);

///对边按照x排序
        sort(ege, ege+k, cmp);

for(i=0; i<k-1; i++)
        {
            int L = lower_bound(Hash, Hash+nh, ege[i].y1) - Hash;
            int R = lower_bound(Hash, Hash+nh, ege[i].y2) - Hash;

upData(1, L, R, ege[i].dir);

C += fabs(a[1].len - pre) + (ege[i+1].x - ege[i].x)*a[1].sum * 2;
            pre = a[1].len;
        }

printf("%d\n", C+pre);
    }

return 0;

}