hdu1542线段树+离散化+扫描线详代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<queue>
#include<algorithm>
#include<map>
#include<iomanip>
using namespace std;
const int MAXN=2000;
int mark[MAXN*2];//记录某个区间的下底边个数
double sum[MAXN*2];//记录某个区间的下底边总长度
double hash1[MAXN];//对x进行离散化,否则x为浮点数且很大无法进行线段树

//以横坐标作为线段(区间),对横坐标线段进行扫描
//扫描的作用是每次更新下底边总长度和下底边个数,增加新面积

struct node{//线段树 
   double l;
   double r;
   double h;
   int d;
   node(){}
   node(double l1,double r1,double h1,int d1):l(l1),r(r1),h(h1),d(d1){}
   bool operator<(const node &a)const
   {
      return h<a.h;
   }
}q[MAXN];
void upfather(int n,int l,int r)
{
   if(mark[n])
      sum[n]=hash1[r+1]-hash1[l];//表示该区间整个线段长度可以作为底边
   else if(l==r)
      sum[n]=0;//叶子结点则底边长度为0(区间内线段长度为0)
   else
      sum[n]=sum[n*2]+sum[n*2+1];
}
void updata(int L,int R,int d,int n,int l,int r)
{
   if(L<=l&&R>=r)//该区间是当前扫描线段的一部分,则该区间下底边总长以及上下底边个数差更新
   {
      mark[n]+=d;;//更新底边相差差个数
      upfather(n,l,r);//更新底边长
      return;
   }
   int mid=(l+r)/2;
   if(L<=mid)
      updata(L,R,d,n*2,l,mid);
   if(R>mid)
      updata(L,R,d,n*2+1,mid+1,r);
   upfather(n,l,r);
}
int find(double ss,double *k,int kk)
{
   int l=0,r=kk-1;
   while(l<=r)
   {
      int mid=(l+r)/2;
      if(k[mid]==ss)
      return mid;
      if(k[mid]<=ss)
         l=mid+1;
      else
         r=mid-1;
   }
   return -1;
}
int main(void)
{
   int n,num=1,i,j;
   while(cin>>n,n)
   {
      double x1,x2,y1,y2,s;
      int k=0,m=1;
      for(i=0;i<n;i++)
      {
         cin>>x1>>y1>>x2>>y2;
         hash1[k]=x1;
         q[k++]=node(x1,x2,y1,1);
         hash1[k]=x2;
         q[k++]=node(x1,x2,y2,-1);
      }
      sort(hash1,hash1+k);
      sort(q,q+k);
      for(i=1;i<k;i++)//去重复端点
         if(hash1[i]!=hash1[i-1])
            hash1[m++]=hash1[i];
      s=0;
      for(i=0;i<k;i++)//扫描线段
      {
         int L=find(q[i].l,hash1,m);
         int R=find(q[i].r,hash1,m)-1;
         updata(L,R,q[i].d,1,0,m-1);//扫描线段时更新底边长度和底边相差个数
         s=s+sum[1]*(q[i+1].h-q[i].h);//新增加面积
      }
      printf("Test case #%d\nTotal explored area: %.2lf\n\n",num++,s);
   }
   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就很好理解了
*/