主要这篇文章想说一些自己遇到的困惑
1. 线段树维护的内容是什么,以y轴建立扫描线为例 ,我们知道当遇到两条边的时候间距很好求出,然而高度却不好求出,并且随着出边与入边的改变,在y轴上被划分的小矩形块的一直在变化,因此,需要维护的就是这个长度。
2. 为什么这段代码当中 r要 -1 ,l 不用?
l = lower_bound(ally + 1,ally + 1 + cnt,e[i].a) - ally;
r = lower_bound(ally + 1,ally + 1 + cnt,e[i].b) - ally - 1;
试想一个问题,五个区间的间隔是多少? 显然是4,我们储存在线段树当中的 l r 并不是端点,而是一个个的区间 !当我们进行修改操作的时候,l代表的其实是从第几个区间开始,r代表第几个区间之前结束,因此 l = 1 ,r = 3 时,我们其实只经过了第 1 和第2个区间。
3.为什么下面代码是 r + 1 呢? 这个时候其实我们使用的是区间左右端点的y值,试想,第三个区间的终点,是不是对应的ally[4]呢(ally记录从小到大,离散后的y的值)?,
if(t[p].cnt) t[p].length = ally[t[p].r + 1] - ally[t[p].l];
4. 为什么cnt = 0 的时候有 t[p].length = t[p * 2].length + t[p * 2 + 1].length;呢?
我们要时刻维护线段树的正确性,当更新cnt值的时候,其子区间也得全部被更新为cnt = 0,因此,这样子写,先递归到最底层,然后逐层递归上来,每个子区间都被更新了。
下附代码:
#include<bits/stdc++.h>
using namespace std;
int n,m,cnt;
typedef long long ll;
const int N = 1001;
double ally[N];
struct node
{
double x,a,b;
int k;
node(){}
node(double a ,double b ,double c ,int d ):x(a),a(b),b(c),k(d) {}
}e[N];
bool cmp(node a ,node b)
{
return a.x < b.x;
}
struct SegmentTree
{
int l,r;
int cnt;
double length;
} t[N * 4];
void build(int p , int l ,int r)
{
t[p].l = l , t[p].r = r;
t[p].length = t[p].cnt = 0;
if(l==r){
// t[p].length = t[p].cnt = 0;
t[p].l = l , t[p].r = r;
return;
}
int mid = (l + r) / 2;
build(p * 2 ,l,mid);
build(p * 2 + 1,mid + 1,r);
}
void change(int p, int l,int r, int v )
{
if(l <= t[p].l && r >= t[p].r)
{
t[p].cnt += v;
if(t[p].cnt) t[p].length = ally[t[p].r + 1] - ally[t[p].l];
else t[p].length = t[p * 2].length + t[p * 2 + 1].length;
return ; //递归最底层
}
int mid =(t[p].l + t[p].r) / 2;
if(l <= mid) change(p*2, l,r, v);
if(r > mid) change(p * 2 + 1, l,r, v);
if(t[p].cnt) t[p].length = ally[t[p].r + 1] - ally[t[p].l];
else t[p].length = t[p * 2].length + t[p * 2 + 1].length;
}
int main(){
int k = 0;
while(~scanf("%d",&n) && n)
{
printf("Test case #%d\n",++k);
cnt = 0;
double x1,y1,x2,y2;
int num = 0;
int numy = 0;
for(int i = 1 ; i <= n ;i++)
{
scanf("%lf %lf %lf %lf",&x1,&y1,&x2,&y2);
e[++num] = node(x1,y1,y2,1);
e[++num] = node(x2,y1,y2,-1);
ally[++numy] = y1;
ally[++numy] = y2;
}
sort(e + 1,e + 1 + num,cmp);
sort(ally + 1 ,ally + 1 + numy);
cnt = unique(ally + 1 ,ally + 1 + numy) - ally - 1;
//for(int i = 1 ;i <= cnt ;i ++) cout <<ally[i] << endl;
//cout << cnt;
build(1,1,cnt - 1);
double ans = 0;
int l,r;
l = lower_bound(ally + 1,ally + 1 + cnt,e[1].a) - ally;
r = lower_bound(ally + 1,ally + 1 + cnt,e[1].b) - ally - 1;
change(1, l, r, e[1].k);
// cout << t[1].length <<endl;
// printf("%lf %lf -> %d %d\n",ally[l],ally[r],l,r);
for(int i = 2 ; i <= num; i++)
{
ans +=(e[i].x -e[i - 1].x)*(t[1].length);
// cout << t[1].length <<endl;
l = lower_bound(ally + 1,ally + 1 + cnt,e[i].a) - ally;
r = lower_bound(ally + 1,ally + 1 + cnt,e[i].b) - ally - 1;
// printf("%lf %lf -> %d %d\n",ally[l],ally[r],l,r);
change(1, l, r, e[i].k);
}
printf("Total explored area: %.2f\n\n",ans);
}
}
对于扫描线的介绍,这篇写的不错戳这里