Problem Description
输入一个n,给你n个矩阵的左下角坐标,和右上角坐标。让你求所有面积并的和。
思路:
扫描线参考博客:http://www.cnblogs.com/scau20110726/archive/2013/04/12/3016765.html
博客是从下到上扫描,个人习惯从左到右扫描。我们左边扫描过来的有效长度设为len,那么线段树就是维护区间最长len,还有cover值(用来标记该区间的整段线段是否有效,大于0的时候有效,等于0的时候无效,根据扫描线博客)。
#include<bits/stdc++.h>
using namespace std;
#define lson root<<1
#define rson root<<1|1
#define N 105
#define MID int mid = (l + r) / 2
struct node
{
int cover;
double len;
};
node tree[N<<3];//线段树
double yy[N<<2];//用来存所有点的纵坐标
struct Node
{
int cover;//左边是1,右边是-1.
double x;//该点的横坐标
double yd, yu;//下纵,上纵坐标
bool operator < (const Node &b) const{
return x < b.x;
}
};
Node a[N<<2];//存各点的数据
void build(int root, int l, int r)
{
tree[root].cover = tree[root].len = 0;//该区间线段无效所以0,有效长度0
if(l + 1 == r)
return;
MID;
build(lson, l, mid);
build(rson, mid, r);
}
void get_len(int root, int l, int r)
{
if(tree[root].cover)//该区间整段有效
tree[root].len = yy[r] - yy[l];
else if(l + 1 == r)//叶子节点,同时区间无效
tree[root].len = 0;
else tree[root].len = tree[lson].len + tree[rson].len;//左右儿子len的和
}
void updata(int root, int l, int r, int ul, int ur, int cover)
{
if(ul <= l && r <= ur)//更新cover,len
{
tree[root].cover += cover;
get_len(root, l, r);
return;
}
MID;
if(ul < mid) updata(lson, l, mid, ul, ur, cover);
if(ur > mid) updata(rson, mid, r, ul, ur, cover);
get_len(root, l, r);//归并,向上更新
}
int main()
{
int n, cnt, i, cas = 1;
double x1, x2, y1, y2;
while(~scanf("%d", &n) && n)
{
cnt = 0;
for(i = 0; i < n; i++)
{
scanf("%lf %lf %lf %lf", &x1, &y1, &x2, &y2);
yy[cnt] = y1;
a[cnt++] = ((Node){1, x1, y1, y2});//左下角
yy[cnt] = y2;
a[cnt++] = ((Node){-1, x2, y1, y2});//右上角
}
sort(yy, yy + cnt);//离散化一下
n = unique(yy, yy + cnt) - yy;
build(1, 0, n - 1);//初始化建树,0-(n-1)是离散化后对应的最小yy到最大yy的范围
sort(a, a + cnt);//按x从小到大排序
int ul, ur;
double ans = 0;
for(i = 0; i < cnt - 1; i++)
{
ul = lower_bound(yy, yy + n, a[i].yd) - yy;
ur = lower_bound(yy, yy + n, a[i].yu) - yy;
updata(1, 0, n - 1, ul, ur, a[i].cover);//每次更新有效长度
ans += tree[1].len * (a[i + 1].x - a[i].x);//有效长度是底边,有效长度×高是面积。
}
printf("Test case #%d\n", cas++);
printf("Total explored area: %.2lf\n\n", ans);
}
}