我的扫面线第一题,一开始看网上讲的都好抽象,最后还是研究别人代码整明白的,所以我要写一个直观的,哈哈哈!!希望大家都能看懂
如图虚线将整个图型分成三个矩形,我们现将每个点的x进行排序,也就是(10,15,20,25.5)
这样就知道矩形的长了,还差宽,然后我们就要开始扫描线啦
先每个点的y值离散话,因为我们要求的是映射在x轴的线段,然后建树,共有t哥点,样立t-1等于4
下面的线段树我每个离散花后的点对应的y值我都标记上了
将一号边加入线段,先加入一号边,二号边三号边四号边具体过程建议自己模拟一下
我每个线段树的点用s标记这条边被覆盖还是没被都盖,每条边的flag表示这天边是入度还是出度
注意啊!!我这个图x,y轴画反啦
#include<iostream> #include<cstdio> #include<algorithm> using namespace std; #define N 300 struct Node{ double x;double y1;double y2; int flag; }node[N]; bool cmp(Node a,Node b){ return a.x-b.x<0.0000001; } double y[N]; struct node{ int l;int r;double ml;double mr;int s;double len; }a[N*3]; void build(int i,int left,int right){ a[i].l=left; a[i].r=right; a[i].ml=y[left]; a[i].mr=y[right]; a[i].s=0; a[i].len=0; if(a[i].l+1==a[i].r){ return ; } int mid=(left+right)>>1; build(i*2,left,mid); build(i*2+1,mid,right);//建树时注意这里不是mid+1,因为做相减的时候如果mid+1这么建回到值左孩子的右边与有孩子的左边无法进行运算 } void callen(int i){ if(a[i].s>0){//注意这里不是所有边都是左孩子的长度加上右孩子的长度,他存在一个覆盖问题 a[i].len=a[i].mr-a[i].ml; }else if(a[i].r-a[i].l==1){ a[i].len=0; }else{ a[i].len=a[i*2].len+a[i*2+1].len; } return ; } void updata(int i,Node b){ if(a[i].ml==b.y1&&a[i].mr==b.y2){ a[i].s+=b.flag; callen(i); return ; } if(b.y2<=a[i*2].mr) updata(i*2,b); else if(b.y1>=a[i*2+1].ml) updata(i*2+1,b); else{ Node temp=b; temp.y2=a[i*2].mr; updata(i*2,temp); temp=b; temp.y1=a[i*2+1].ml; updata(i*2+1,temp); } callen(i); return ; } int main(){ // freopen("in.txt","r",stdin); int n,t,p=1,te; double x1,x2,y1,y2; while(scanf("%d",&n),n){ t=1; for(int i=0;i<n;i++){ scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2); node[t].x=x1; node[t].y1=y1; node[t].y2=y2; node[t].flag=1;//入边 y[t++]=y1; node[t].x=x2; node[t].y1=y1; node[t].y2=y2; node[t].flag=-1;//出边 y[t++]=y2; } sort(node+1,node+t,cmp); sort(y+1,y+t); build(1,1,t-1); updata(1,node[1]); double sum=0; for(int i=2;i<t;i++){ sum+=a[1].len*(node[i].x-node[i-1].x); updata(1,node[i]); } printf("Test case #%d\n",p++); printf("Total explored area: %.2lf\n\n",sum); } return 0; }