[离散化+线段树+扫描线]POJ_1151_Atlantis

题意：求矩形面积并

前天直接用“矩形切割”（我也不知道我那方法是不是矩形切割，姑且打上引号吧），过掉了这个题，那个猥琐的算法达到了O(n^3)，如果题目数据很大的话，还是不行的，所以这两天还是继续学习扫描线，发现网上的很多文章讲的都不算很清楚，于是翻除了1999年陈宏的论文，开始研究，最终总算搞定了扫描线这个玩意儿。

注意点：要理解论文中提到的测度（即线段并的长度）和连续段的含义（求矩形周长并需要用到），这两点很关键，然后建树的时候用段树比较好一点，用点树我还没想到这么处理，建段树的话，数组稍微开大点（最好是5倍左右），否则会悲剧的。。

各种蛋疼：最后忘记去掉freopen..WA数次，蛋碎。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAXN = 5555;
struct Line{
    double l,r,h;
    int s;
    Line(){}
    Line(double ll,double rr,double hh,int ss):l(ll),r(rr),h(hh),s(ss){}
    bool operator < (const Line &l) const {
        return h<l.h;
    }
}L[MAXN<<2];
double X[MAXN<<2],sum[MAXN<<2];
int cnt[MAXN<<2],num_x,n,line_num;
#define lson l,m,rt<<1
#define rson m,r,rt<<1|1
void pushUP(int rt,int l,int r){
    if(cnt[rt]){
        sum[rt] = X[r]-X[l];
    }else{
        sum[rt] = sum[rt<<1]+sum[rt<<1|1];
    }
}
void update(int L,int R,int s,int l,int r,int rt){
    if(L<=l&&R>=r){
        cnt[rt]+=s;
        pushUP(rt,l,r);
        return;
    }
    //cout<<L<<" "<<R<<" "<<l<<" "<<r<<endl;system("pause");
    int m = (l+r)>>1;
    if(m>L)update(L,R,s,lson);
    if(m<R)update(L,R,s,rson);
    pushUP(rt,l,r);
}
int search(double a[],double k){
    int l=0,r=num_x-1,m;
    while(l<=r){
        m = (l+r)>>1;
        if(a[m]==k)return m;
        else if(k<a[m])r=m-1;
        else l=m+1;
    }
    return -1;
}
int main(){
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int cas=0;
    while(scanf("%d",&n),n){
        num_x = line_num = 0;
        for(int i=0;i<n;i++){
            double a,b,c,d;
            scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
            X[num_x++] = a;
            X[num_x++] = c;
            L[line_num++] = Line(a,c,b,1);
            L[line_num++] = Line(a,c,d,-1);
        }
        sort(X,X+num_x);
        sort(L,L+line_num);
        int size = 1;
        for(int i=1;i<num_x;i++)if(X[i]!=X[i-1])X[size++]=X[i];
        num_x = size;
        memset(cnt,0,sizeof(cnt));
        memset(sum,0,sizeof(sum));
        double ans = 0;
        int l,r;
        l = search(X,L[0].l);
        r = search(X,L[0].r);
        update(l,r,L[0].s,0,num_x-1,1);
        for(int i=1;i<line_num;i++){
            ans+=sum[1]*(L[i].h-L[i-1].h);
            l = search(X,L[i].l);
            r = search(X,L[i].r);
            update(l,r,L[i].s,0,num_x-1,1);
        }
        printf("Test case #%d\nTotal explored area: %.2lf\n\n",++cas,ans);
    }
    return 0;
}