本文主要是介绍POJ 1177 Picture (线段树扫描线),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
题意:
给定n个矩形(0 <= n < 5000)的左下角坐标(x1,y1)和右上角坐标(x2,y2) (-10000 <= x1,x2,y1,y2 <= 10000)
求所有矩形重合后的图形的周长,如下图(图片来自POJ 1177):
做法:线段树扫描线。
由于值域不大,所以不需要离散化,直接将Y值向正方向平移10001个单位,然后用线段树直接做。
扫描线就是用垂直于x轴的线来从左到右扫描整个图形,然后用线段树来记录这条线上有多少个单位被覆盖。
struct Node{
int cover;//完全覆盖层数
int lines;//分成多少个线段
bool L,R;//左右端点是否被覆盖
int CoverLength;//覆盖长度
int Length;//总长度
};
为了计算与x轴平行的边,需要记录扫描线上有多少条线段,用lines 表示,然后 L,R表示该区间的左右端点是否被覆盖,方便合并区间的时候维护 lines。
Length是本区间的总长度,CoverLength是本区间被覆盖的长度,cover是指本区间被完全覆盖的次数,非完全覆盖不改变cover值。
然后,每个矩形,遇到左侧的边就将这条边加入扫描线,遇到右侧的边就把这条边从扫描线删除。
每次操作后更新周长。
具体见代码:
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#define inf 0x1fffffff
#define out(i) <<#i<<"="<<(i)<<" "
#define OUT1(a1) cout out(a1) <<endl
#define OUT2(a1,a2) cout out(a1) out(a2) <<endl
#define OUT3(a1,a2,a3) cout out(a1) out(a2) out(a3)<<endl
#define ls l,m,rt<<1
#define rs m+1,r,rt<<1|1
#define maxn 20007 //直接按y值建树
using namespace std;struct Node{int cover;//完全覆盖层数int lines;//分成多少个线段bool L,R;//左右端点是否被覆盖int CoverLength;//覆盖长度 int Length;//总长度 Node(){}Node(int cover,int lines,bool L,bool R,int CoverLength):cover(cover),lines(lines),L(L),R(R),CoverLength(CoverLength){}Node operator +(const Node &B){//连续区间的合并 Node C;C.cover=0;C.lines=lines+B.lines-(R&&B.L);C.CoverLength=CoverLength+B.CoverLength;C.L=L;C.R=B.R;C.Length=Length+B.Length;return C;}void Show(){printf("cover=%d lines=%d L=%d R=%d CoverLength=%d Length=%d\n",cover,lines,L,R,CoverLength,Length);}
}K[maxn<<2];
void PushUp(int rt){//更新非叶节点 if(K[rt].cover){K[rt].CoverLength=K[rt].Length;K[rt].L=K[rt].R=K[rt].lines=1;}else{K[rt]=K[rt<<1]+K[rt<<1|1];}
}
void Build(int l,int r,int rt){//建树 if(l==r){K[rt].L=K[rt].R=K[rt].cover=K[rt].lines=K[rt].CoverLength=0;K[rt].Length=1;return;}int m=(l+r)>>1;Build(ls);Build(rs);K[rt].cover=0;PushUp(rt);
}
void Cover(int L,int R,int l,int r,int rt){//覆盖 if(L <= l && r <= R){if(!K[rt].cover){//如果覆盖前为空,则更新节点内容 K[rt].L=K[rt].R=true;K[rt].lines=1;K[rt].CoverLength=K[rt].Length; }K[rt].cover++;return;}int m=(l+r)>>1;if(L <= m) Cover(L,R,ls);if(R > m) Cover(L,R,rs);PushUp(rt);
}
void Uncover(int L,int R,int l,int r,int rt){//移除覆盖 if(L <= l && r <= R){//如果移除后为空,则更新节点 --K[rt].cover;if(!K[rt].cover) {if(l==r){//特判叶节点 K[rt].L=K[rt].R=K[rt].cover=K[rt].lines=K[rt].CoverLength=0;}else K[rt]=K[rt<<1]+K[rt<<1|1];//非叶节点用重载的+号 }return;}int m=(l+r)>>1;if(L <= m) Uncover(L,R,ls);if(R > m) Uncover(L,R,rs);PushUp(rt);
}int n;
struct Lines{//线段 int x,y1,y2,c;//x:线的横坐标 y1,y2:线的纵坐标端点 c:加入1 移除0 Lines(){}Lines(int x,int y1,int y2,int c):x(x),y1(y1),y2(y2),c(c){}bool operator < (const Lines &B)const{return x < B.x;} //重载小于号方便排序
}L[10000];
int main(void)
{ while(~scanf("%d",&n)){int x1,y1,x2,y2;for(int i=0;i<n;++i){scanf("%d%d%d%d",&x1,&y1,&x2,&y2);L[i<<1]=Lines(x1,y1+10001,y2+10001,1);L[i<<1|1]=Lines(x2,y1+10001,y2+10001,0);}int nL=n<<1;sort(L,L+nL);//按x值排序 int PreX=L[0].x;//前X坐标 int ANS=0;//目前累计答案 int PreLength=0;//前线段总长int PreLines=0;//前线段数量 Build(1,20001,1);for(int i=0;i<nL;++i){//操作 if(L[i].c) Cover(L[i].y1,L[i].y2-1,1,20001,1);else Uncover(L[i].y1,L[i].y2-1,1,20001,1);//更新横向的边界 ANS+=2*PreLines*(L[i].x-PreX);PreLines=K[1].lines;PreX=L[i].x;//更新纵向边界 ANS+=abs(K[1].CoverLength-PreLength);PreLength=K[1].CoverLength;}//输出答案 printf("%d\n",ANS);}
return 0;
}
这篇关于POJ 1177 Picture (线段树扫描线)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!