【HDU 1542】Atlantis（线段树+离散化，矩形面积并)

求矩形面积并，离散化加线段树。

扫描线法：

用平行x轴的直线扫,每次ans+=（下一个高度-当前高度）*当前覆盖的宽度。

#include<algorithm>
#include<cstdio>
#include<cstring>
#define dd double
#define ll long long
#define N 201
using namespace std;
struct P{dd s,e,h;int f;}p[N];
struct Tree{dd sum;int c;}t[N<<];
dd a,b,c,d,x[N],ans;
int n,m,num;
int cmp(const P &a,const P &b){
    return a.h<b.h;
}
void pushUp(ll rt,ll l,ll r){
    if(t[rt].c)t[rt].sum=x[r+]-x[l];//r+1是因为节点[l,r]表示区间[x[l],x[r+1]]。
    else if(l==r)t[rt].sum=;
    else t[rt].sum=t[rt<<].sum+t[rt<<|].sum;
}
void update(ll s,ll e,ll rt,ll l,ll r,ll v){
    if(s<=l&&r<=e) t[rt].c+=v;
    else {
        if(l>e||r<s)return;
        ll m=l+r>>;
        update(s,e,rt<<,l,m,v);
        update(s,e,rt<<|,m+,r,v);
    }
    pushUp(rt,l,r);
}
int main()
{
    while(scanf("%d",&n),n){
        int k=;
        for(int i=;i<n;i++){
            scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
            x[++k]=a,p[k]=(P){a,c,b,};
            x[++k]=c,p[k]=(P){a,c,d,-};
        }
        sort(x+,x++k);
        sort(p+,p++k,cmp);
        m=;
        for(int i=;i<=k;i++)
            if(x[i]!=x[i-])x[m++]=x[i];
        ans=;
        for(int i=;i<=k;i++){//共k条线段,每次计算p[i].h到p[i+1].h之间的面积，第k次相当于清空所有，酱就不用初始化线段树了。
            int l=lower_bound(x,x+m,p[i].s)-x;
            int r=lower_bound(x,x+m,p[i].e)-x-;//r-1同上原因
            update(l,r,,,m-,p[i].f);
            ans+=t[].sum*(p[i+].h-p[i].h);
        }
        printf("Test case #%d\nTotal explored area: %.2lf\n\n",++num,ans);
    }
    return ;
}

另一种方法还是线段树，这里扫描线用的是平行y轴的直线，每次增加的面积是当前扫描的竖线所在的高度区间的最后一次的x与当前x的差值乘上区间的高度。所以每次增加的不一定是一个矩形，而是多个矩形并。

　

#include<cstdio>
#include<algorithm>
#define dd double
using namespace std;
#define N 201
struct P{dd x,y1,y2;int f;}p[N];
struct TREE{dd y1,y2,x;int c,f;}tree[N<<];
dd x1,y1,x2,y2,y[N];
int n,k,num;
int cmp(const P &a,const P &b){
    return a.x<b.x;
}
void build(int i,int l,int r){
    tree[i].c=tree[i].f=;
    tree[i].y1=y[l];//直接将线段树节点代表的区间存在线段树里
    tree[i].y2=y[r];
    if(l+==r){
        tree[i].f=;
        return;
    }
    int mk=(l+r)>>;
    build(i<<,l,mk);
    build(i<<|,mk,r);
}
dd insert(int i,dd x,dd l,dd r,int flag){
    if(r<=tree[i].y1||l>=tree[i].y2)
        return ;
    if(tree[i].f){//离散后的一个最小区间，叶子节点
        dd ans;
        if(tree[i].c>)//全覆盖
            ans=(x-tree[i].x)*(tree[i].y2-tree[i].y1);//(当前x-该区间最后的x)*区间高度
        else
            ans=;
        tree[i].x=x;//该区间最新的x
        tree[i].c+=flag;//更新覆盖
        return ans;
    }
    return insert(i<<,x,l,r,flag)+insert(i<<|,x,l,r,flag);
}
int main(){
    while(scanf("%d",&n),n){
        k=;
        for(int i=;i<=n;i++){
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            y[++k]=y1,p[k]=(P){x1,y1,y2,};
            y[++k]=y2,p[k]=(P){x2,y1,y2,-};
        }
        sort(y+,y+k+);
        sort(p+,p+k+,cmp);
        //没有去重，其实数量少，没必要。
        build(,,k);
        dd ans=;
        for(int i=;i<=k;i++)
            ans+=insert(,p[i].x,p[i].y1,p[i].y2,p[i].f);
        printf("Test case #%d\nTotal explored area: %.2f\n\n",++num,ans);
    }
    return ;
}