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);

}