poj 1177 又是一道十分恶心到线段树。 主要是用线段树到区间操作。可以用来节省线性查找的时间。 离散化也十分好的一道题。

这道题写来将近半天。

题意

就是求所覆盖到矩形外面的周长（不包括里面的空心矩形）。

 

算法： 线段树+离散化

 

思路： 把x和y轴作成离散化（就是用超元离散化   98陈宏论文有讲）。然后把 y轴作成，线段树，然后把x轴 扫描一遍。

特别注意的就是，线段树与以往有所不同。 这个线段树到中点的处理.  线段树左右代表的是点.而我们要用是区间.所以 要保证r-l>1

 

下面两篇讲的还不错.

http://angels1.0.blog.163.com/blog/static/84580504200893171819117/

http://hi.baidu.com/soonyu/blog/item/a66a6dedfa9a2a2c63d09f28.html

 

 

下面我的代码

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <string.h>
#include <cmath>
#include <vector>
using namespace std;

#define e 10010
struct node{
    int len,pl,pr;
    int l,r;
    int cover;
    int lennum;
    int sum;
    }tree[e];


int cnt=-1;
struct Nodee{
        int sta,end;
        int x;

        bool flag;
        bool operator < (const Nodee & a)const{
            return x<a.x;
            }
}yline[e];

int Index[e];

void build(int nod,int l,int r)
{
    if(l<r)
    {
    tree[nod].l=l;
    tree[nod].r=r;
    tree[nod].len=Index[r]-Index[l];
    tree[nod].pl=tree[nod].pr=0;
    tree[nod].sum=0;
    tree[nod].lennum=0;
    tree[nod].cover=0;

    if(r-l<=1)    return ;
    int mid= ( l+r )/2;

    build(nod<<1,l,mid);
    build(nod<<1|1,mid,r);
    }
}

void updatey(int nod)
{
        if(tree[nod].cover>0)
            tree[nod].sum=tree[nod].len;
        else if(tree[nod].r -tree[nod].l>1)
            tree[nod].sum=tree[nod<<1].sum+tree[nod<<1|1].sum;
        else
            tree[nod].sum=0;
}
void updatex(int nod)
{
        if(tree[nod].cover>0)
        {
            tree[nod].lennum=1;
            tree[nod].pl=tree[nod].pr=1;
        }
        else if(tree[nod].r -tree[nod].l >1)
            {
                   tree[nod].pr=tree[nod<<1|1].pr;
                    tree[nod].pl=tree[nod<<1].pl;
                    tree[nod].lennum= tree[nod<<1].lennum+tree[nod<<1|1].lennum- tree[nod<<1].pr * tree[nod<<1|1].pl;
            }
        else
        {
                tree[nod].pr=tree[nod].pl=0;
                tree[nod].lennum=0;
        }
}

void insert(int nod,int l,int r,int num)
{

   // printf("%d %d %d %d/n",l,r,tree[nod].l,tree[nod].r);

    if(l<=tree[nod].l  && tree[nod].r <=r)
        tree[nod].cover+=num;
   else
   {
        int mid=(tree[nod].l+tree[nod].r)/2;
        if(l<mid)
            insert(nod<<1,l,r,num);
        if(r>mid)
            insert(nod<<1|1,l,r,num);
   }


    updatey(nod);
    updatex(nod);
}

int findd(int w){
        return lower_bound( Index , Index+cnt,w ) - Index;
    }

int main()
{
    #ifndef ONLINE_JUDGE
     freopen("//home//mabodx//桌面//in.txt","r",stdin);

#endif
    int i,j,n,m,k,l,t;
    scanf("%d",&t);
    int x1,x2,y1,y2;
    for(i=0; i<t; i++)
    {
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            Index[i<<1]=y1;
            Index[i<<1|1]=y2;

            yline[i<<1].x=x1;
            yline[i<<1].sta=y1;
            yline[i<<1].end=y2;
            yline[i<<1].flag=1;

            yline[i<<1|1].x=x2;
            yline[i<<1|1].sta=y1;
            yline[i<<1|1].end=y2;
            yline[i<<1|1].flag=0;
    }


    sort(Index,Index+t*2);


    sort(yline,yline+t*2);
     cnt=-1;

    for(i=1; i<=2*t-1; i++)
    {
            if(Index[i]!=Index[i-1])
                 Index[++cnt]=Index[i-1];
    }

    Index[++cnt]=Index[2*t-1];

    build(1,0,cnt);


    int ans=0;
    int total=0;

    for(i=0; i<2*t-1; i++)
    {

        if(yline[i].flag)
         insert(1, findd(yline[i].sta),findd(yline[i].end),1);
        else
            insert(1,findd(yline[i].sta), findd(yline[i].end),-1);

        ans+=(tree[1].lennum * 2*(yline[i+1].x-yline[i].x ));

        ans+=abs(tree[1].sum-total);

        total=tree[1].sum;
    }
    insert(1,findd(yline[i].sta),findd(yline[i].end),-1);
    ans+=abs(tree[1].sum- total);
    printf("%d/n", ans);
}