HDU1542 Atlantis 線段樹掃描面積

第一次做面積掃描的題，首先進行數據的離散化排序，建立線段樹，把整個圖形從下到上一次求值，

如圖，我們每次更新最下面的邊的長度，首先是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;
}