题目链接:http://poj.org/problem?id=1151
大致题意:给出n个矩形,求出n个矩形一共覆盖的面积
做法:现将坐标离散化,然后用线段树维护
即让线段按照x坐标的大小从小到大排列,y坐标按照从小到大排列去掉重复的值
然后通过一个扫描线来求扫描线覆盖的y的长度。
线段的扫描按照x的大小从小到大扫描,求出当前扫描线覆盖的矩阵竖线的长度,然后乘以下条线段的跨度,则为这个区域矩阵覆盖的面积
其实,还是我柴的板硬。。。。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define MAXN 220
struct Node
{
int l, r;
int c;
double cnt, lf, rf;
} segTree[MAXN*3];
struct Line
{
double x, y1, y2;
int f;
} line[MAXN];
bool cmp(Line a, Line b)
{
return a.x < b.x;
}
double y[MAXN];
void Build(int t, int l, int r)
{
segTree[t].l = l;
segTree[t].r = r;
segTree[t].cnt = segTree[t].c = 0;
segTree[t].lf = y[l];
segTree[t].rf = y[r];
if(l + 1 == r) return ;
int mid = (l+r) >> 1;
Build(t<<1, l, mid);
Build(t<<1|1, mid, r);
}
void calen(int t)
{
if(segTree[t].c > 0)
{
segTree[t].cnt = segTree[t].rf - segTree[t].lf;
return ;
}
if(segTree[t].l+1 == segTree[t].r)
segTree[t].cnt = 0;
else
segTree[t].cnt = segTree[t<<1].cnt + segTree[t<<1|1].cnt;
}
void update(int t, Line e)
{
if(e.y1 == segTree[t].lf && e.y2 == segTree[t].rf)
{
segTree[t].c += e.f;
calen(t);
return ;
}
if(e.y2 <= segTree[t<<1].rf)
update(t<<1, e);
else if(e.y1 >= segTree[t<<1|1].lf)
update(t<<1|1, e);
else
{
Line tmp = e;
tmp.y2 = segTree[t<<1].rf;
update(t<<1, tmp);
tmp = e;
tmp.y1 = segTree[t<<1|1].lf;
update(t<<1|1, tmp);
}
calen(t);
}
int main()
{
int i, n, t, cas = 0;
double x1, x2, y1, y2;
while(cin>>n && n)
{
t = 1;
for(int i=1; i<=n; i++)
{
cin>>x1>>y1>>x2>>y2;
line[t].x = x1, line[t].y1 = y1;
line[t].y2 = y2, line[t].f = 1;
y[t++] = y1;
line[t].x = x2, line[t].y1 = y1;
line[t].y2 = y2, line[t].f = -1;
y[t++] = y2;
}
sort(line+1, line+t, cmp);
sort(y+1, y+t);
Build(1,1,t-1);
update(1, line[1]);
double res = 0;
for(int i=2; i<t; i++)
{
res += segTree[1].cnt * (line[i].x - line[i-1].x);
update(1, line[i]);
}
printf("Test case #%d\nTotal explored area: %.2f\n\n",++cas, res);
}
return 0;
}