第一次做面积扫描的题,首先进行数据的离散化排序,建立线段树,把整个图形从下到上一次求值,
如图,我们每次更新最下面的边的长度,首先是AB长度,之后可以AB*AP,这是第一部分面积,接着继续更新向上面的一层,即PB'长度,之后PB'*CP加在之前的面积上,依次类推,我们需要把每个矩形的上下边区分一下,用来更新上一层的长度,如图再向上的长度不是CP',而是Q'P',这个可以结合代码来理解,在离散化数据的时候,我们不能直接把线段的右端点的下标直接带入线段树的更新右端点值,因为我们的离散化做完后是没有下标相同的时候,对于线段树来说不是很规范(其实我也不是太懂。。。),线段树的每个区间代表的是该区域内横向线段长度,如果是上边和下边数相等就可以取零,左端点和右端点相等时也取零,然后就瞎写一通,乱搞了出来(上面的说明可以结合代码来理解)。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<ctime>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<set>
#include<bitset>
#include<map>
#include<stack>
#include<queue>
#include<vector>
#include<utility>
#define INF 0x3f3f3f3f
#define inf 2*0x3f3f3f3f
#define llinf 1000000000000000000
#define pi acos(-1)
#define mod 1000000007
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
struct node
{
double l,r,h;
int s;
node(double l=0,double r=0,double h=0,double s=0):l(l),r(r),h(h),s(s){}
bool operator<(const node &a)const
{
return h<a.h;
}
}point[2005];
double xl[2005<<2],X[2005];
int cnt[2005<<2];
void pushup(int l,int r,int rt)
{
if(cnt[rt])xl[rt]=X[r+1]-X[l];
else if(l==r)xl[rt]=0;
else xl[rt]=xl[rt<<1]+xl[rt<<1|1];
}
void update(int L,int R,int c,int l,int r,int rt)
{
if(L<=l&&r<=R)
{
cnt[rt]+=c;
pushup(l,r,rt);
return ;
}
int m=(l+r)>>1;
if(L<=m)update(L,R,c,lson);
if(m<R)update(L,R,c,rson);
pushup(l,r,rt);
}
int Bin(double key,int n,double X[])
{
int l=1,r=n;
while(l<=r)
{
int m=(l+r)>>1;
if(X[m]==key)return m;
if(X[m]<key)l=m+1;
else r=m-1;
}
return -1;
}
int main()
{
int n,cas,num,xnum;
double x1,y1,x2,y2;
for(int cas=1;;cas++)
{
scanf("%d",&n);
if(n==0)break;
num=1;xnum=2;
memset(cnt,0,sizeof(cnt));
memset(xl,0,sizeof(xl));
for(int i=0;i<n;i++)
{
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
X[num]=x1;point[num++]=node(x1,x2,y1,1);
X[num]=x2;point[num++]=node(x1,x2,y2,-1);
}
sort(X+1,X+num);sort(point+1,point+num);
for(int i=2;i<num;i++)
{
if(X[i]!=X[i-1])X[xnum++]=X[i];
}
double ret=0;xnum--;
for(int i=1;i<num-1;i++)
{
int L=Bin(point[i].l,xnum,X),R=Bin(point[i].r,xnum,X)-1;
if(L<=R)update(L,R,point[i].s,1,xnum,1);
ret+=xl[1]*(point[i+1].h-point[i].h);
}
printf("Test case #%d\nTotal explored area: %.2lf\n\n",cas,ret);
}
return 0;
}