hdu1542线段树+离散化+扫描线详解

用

struct tree{

double l,r,h;

int d;

tree(){}

tree(double x1,double x2,double y,int c):l(x1),r(x2),h(y),d(c){}

bool operator <(const tree &a)const{

return h<a.h;

}

}s[MAX];

s表示一条线段,分别有线段的左右端点l,r,线段的高h,线段是矩形的上底边还是还是下底边d:-1/1,

假设给出的矩形如上图,红色字表示相应边的高,

共9条线段

经从小到大排序后的边的高度h是:2,3,3,4,5,6,7,7,8,9

                    相应的d是:1,1,1,-1,1,1,-1,-1,-1,-1

1,对第一条线段进行扫描:

得到的总下底边长sum[1]为黑色线部分,同时mark标记某个区间的下底边情况,1表示该区间(线段)全部可作为下底边

的到的面积ans+=sum[1]*(s[i+1].h-s[i].h)即为蓝色部分

2,对第2条线段继续扫描

sum为1.,2两部分合起来的总长度,此时s[3].h-s[2].h

=3-3=0,所以增加的面积是0

对第3条线段继续扫描

得到的总底边长度不变,新增加的面积为新增加的蓝色区域,忘记说了经过扫描后,得到2线段和3线段之间的mark=1，因为该区域下底边比上底边多就只有一条,1线段和2,3线段交集部分mark,因为上下底边差为2条

对第4条线段扫面,即高为4的那条线段,由于这条线段为上底边,所以d=-1，所以扫描后2线段和3线段之间的mark为0(1+-(1)=0),其余部分为1,则重新得到的总下底边sum[1]为线段2和线段3的长度

新增加的面积为蓝色面积

后面同理继续扫面线段增加新面积............

总结下,mark用来记录区间内上下底边差的个数,作用是为了统计总的下底边的长度sum时不出错,使得不会重复计算某个面积

其实这样扫描相当于把给定的矩形通过边右分成若干个举行进行计算,而统计总下底边进行计算相当于同时计算多个矩形的面积，达到效率增快

原图分成的若干矩形,即任何高度不想同的两条线段之间都可以分成矩形来计算

 

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<queue>
#include<algorithm>
#include<map>
#include<iomanip>
#define INF 99999999
using namespace std;

const int MAX=200+10;
int mark[MAX<<2];//记录某个区间的下底边个数
double sum[MAX<<2];//记录某个区间的下底边总长度
double hash[MAX];//对x进行离散化,否则x为浮点数且很大无法进行线段树 

//以横坐标作为线段(区间),对横坐标线段进行扫描
//扫描的作用是每次更新下底边总长度和下底边个数,增加新面积 
struct seg{//线段 
	double l,r,h;
	int d;
	seg(){}
	seg(double x1,double x2,double H,int c):l(x1),r(x2),h(H),d(c){}
	bool operator<(const seg &a)const{
		return h<a.h;
	}
}s[MAX];

void Upfather(int n,int left,int right){
	if(mark[n])sum[n]=hash[right+1]-hash[left];//表示该区间整个线段长度可以作为底边 
	else if(left == right)sum[n]=0;//叶子结点则底边长度为0(区间内线段长度为0) 
	else sum[n]=sum[n<<1]+sum[n<<1|1];
}

void Update(int L,int R,int d,int n,int left,int right){
	if(L<=left && right<=R){//该区间是当前扫描线段的一部分,则该区间下底边总长以及上下底边个数差更新 
		mark[n]+=d;//更新底边相差差个数 
		Upfather(n,left,right);//更新底边长 
		return;
	}
	int mid=left+right>>1;
	if(L<=mid)Update(L,R,d,n<<1,left,mid);
	if(R>mid)Update(L,R,d,n<<1|1,mid+1,right);
	Upfather(n,left,right);
}

int search(double key,double* x,int n){
	int left=0,right=n-1;
	while(left<=right){
		int mid=left+right>>1;
		if(x[mid] == key)return mid;
		if(x[mid]>key)right=mid-1;
		else left=mid+1;
	}
	return -1;
}

int main(){
	int n,num=0;
	double x1,x2,y1,y2;
	while(cin>>n,n){
		int k=0;
		for(int i=0;i<n;++i){
			cin>>x1>>y1>>x2>>y2;
			hash[k]=x1;
			s[k++]=seg(x1,x2,y1,1);
			hash[k]=x2;
			s[k++]=seg(x1,x2,y2,-1);
		}
		sort(hash,hash+k);
		sort(s,s+k);
		int m=1;
		for(int i=1;i<k;++i)//去重复端点 
		    if(hash[i] != hash[i-1])hash[m++]=hash[i];
        double ans=0;
        //memset(mark,0,sizeof mark);
        //memset(sum,0,sizeof sum);如果下面是i<k-1则要初始化,因为如果对第k-1条线段扫描时会使得mark,sum为0才不用初始化的 
		for(int i=0;i<k;++i){//扫描线段 
			int L=search(s[i].l,hash,m);
			int R=search(s[i].r,hash,m)-1;
			Update(L,R,s[i].d,1,0,m-1);//扫描线段时更新底边长度和底边相差个数
			ans+=sum[1]*(s[i+1].h-s[i].h);//新增加面积 
		}
		printf("Test case #%d\nTotal explored area: %.2lf\n\n",++num,ans);
	}
	return 0;
}
/*
这里注意下扫描线段时r-1:int R=search(s[i].l,hash,m)-1;
计算底边长时r+1:if(mark[n])sum[n]=hash[right+1]-hash[left];
解释:假设现在有一个线段左端点是l=0,右端点是r=m-1
则我们去更新的时候,会算到sum[1]=hash[mid]-hash[left]+hash[right]-hash[mid+1]
这样的到的底边长sum是错误的,why?因为少算了mid~mid+1的距离,由于我们这利用了
离散化且区间表示线段,所以mid~mid+1之间是有长度的,比如hash[3]=1.2,hash[4]=5.6,mid=3
所以这里用r-1,r+1就很好理解了 
*/