【线段树】[POJ Atlantis] 线段树扫描线

本文主要是介绍【线段树】[POJ Atlantis] 线段树扫描线，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

题目：

题目链接：[POJ Atlantis]
题解：
（这个题作为学习线段树的尾巴，，，结果还是学了好长时间）

先安利一下sdau_blue大佬的优秀博客，太有用了。

这个题就是有许多的矩形进行随意放置在二维平面内，重叠的面积只算一次，求总面积。

在这里就是离散化处理横坐标，对于纵坐标进行结构体处理一下，然后再以横坐标作为线段(区间)，对横坐标线段进行扫描，扫描的作用是每次更新下底边总长度和下底边个数,增加新面积即可。

代码：

#include<stdio.h>
#include<string.h>
#include<algorithm>
#define D double
using namespace std;
inline int read()
{int s=0,w=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}while(ch<='9'&&ch>='0')s=s*10+ch-'0',ch=getchar();return s*w;
}
const int sea=1e5+7;
int T,n,m,t,cnt[sea];
D ans,x[sea],sum[sea];
struct hit
{D l,r,h;int d;hit(){} hit(D l,D r,D h,int d):l(l),r(r),h(h),d(d) {}
}tr[sea*4];//结构体函数，，，，lrd的经典
bool cmp(hit x,hit y){return x.h<y.h;}
void up(int l,int r,int i)
{if(cnt[i]) sum[i]=x[r+1]-x[l];else sum[i]=sum[i<<1]+sum[i<<1|1];
}
void build(int x,int y,int l,int r,int i,int tt)
{if(x<=l&&y>=r){cnt[i]+=tt; up(l,r,i);return ;}int mid=(l+r)/2;if(x<=mid) build(x,y,l,mid,i<<1,tt);if(y>mid) build(x,y,mid+1,r,i<<1|1,tt);up(l,r,i);
}
int main()
{while(1){T=read(); if(!T) return 0;	n=m=0; ans=0;memset(sum,0,sizeof(sum));memset(cnt,0,sizeof(cnt));for(int i=1;i<=T;i++){D x1,y1,x2,y2;scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);x[++n]=x1; x[++n]=x2;tr[++m]=hit(x1,x2,y1,1);tr[++m]=hit(x1,x2,y2,-1);}sort(x+1,x+n+1); sort(tr+1,tr+m+1,cmp);int k=unique(x+1,x+n+1)-x-1;for(int i=1;i<m;i++){int l=lower_bound(x+1,x+k+1,tr[i].l)-x;int r=lower_bound(x+1,x+k+1,tr[i].r)-x; r--;//这里的r-1，，还不是很懂，应该是为了防止少算mid——mid-1之间少出现的小数所以上面要加1，下面减1应该是与其抵消的，，还是不很懂if(l<=r) build(l,r,1,k-1,1,tr[i].d);ans+=sum[1]*(tr[i+1].h-tr[i].h);}printf("Test case #%d\nTotal explored area: %.2f\n\n",++t,ans);}return 0;
}