hdu1542 Atlantis

题目

矩形面积并基础题。

考虑把每一个矩形按照边界分成小矩形，用扫描线就好了，然后就是统计另一个方向上的长度，这个用线段树维护一下就好了。

比如这张图，大概就是这样分割的，对于x轴的一段，对应了y轴上的一部分。

这样，就可以x轴扫面线，然后y轴线段树区间修改就好了。这个线段树也比较好写，不需要pushdown操作，因为每次我们只需要整个的长度。

还有一点，这道题的坐标都是double类型的，就先离散化一下，至此，这道题就被完美解决了。

//矩形面积并 
#include<bits/stdc++.h>
#define N 10000
using namespace std;
double x1,x2,y,y2,tmp[N+5],l,r;
struct seg{
    double l,r,h;int ind;
    seg(){}
    seg(double l,double r,double h,int ind):l(l),r(r),h(h),ind(ind){}
    bool operator < (const seg &a)const
    {
        return h<a.h;
    }
};seg A[N+5];
struct Tree{
    int l,r,cnt;
    double sum;
};Tree T[4*N+5];
int n,id,cnt;
double f;
inline void build(int rt,int l,int r)
{
    T[rt].sum=T[rt].cnt=0;
    T[rt].l=l,T[rt].r=r;
    if(l==r){return;}
    int mid=(T[rt].l+T[rt].r)/2;
    build(rt*2,l,mid),build(rt*2+1,mid+1,r);
}
inline void pushup(int rt)
{
    if(T[rt].cnt)T[rt].sum=tmp[T[rt].r+1]-tmp[T[rt].l];
    else if(T[rt].l==T[rt].r)T[rt].sum=0;
    else T[rt].sum=T[rt*2].sum+T[rt*2+1].sum;
}
inline void modify(int rt,int l,int r,int x)
{
    if(T[rt].l==l&&T[rt].r==r)
    {
        T[rt].cnt+=x;
        pushup(rt);
        return;
    }
    int mid=(T[rt].l+T[rt].r)/2;
    if(r<=mid)modify(rt*2,l,r,x);
    else if(l>mid)modify(rt*2+1,l,r,x);
    else modify(rt*2,l,mid,x),modify(rt*2+1,mid+1,r,x);
    pushup(rt);
}
int main()
{
    freopen("in.txt","r",stdin);
    while(scanf("%d",&n)!=EOF)
    {
        if(n==0)return 0;
        id++;f=0.0;
        printf("Test case #%d\n",id);
        printf("Total explored area: ");
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y,&x2,&y2);
            A[i]=seg(x1,x2,y,1);A[i+n]=seg(x1,x2,y2,-1);
            tmp[i]=x1,tmp[i+n]=x2;
        }
        n<<=1;sort(A+1,A+n+1);
        sort(tmp+1,tmp+n+1);
        cnt=unique(tmp+1,tmp+n+1)-tmp-1;
        build(1,1,cnt);
        for(int i=1;i<n;i++)
        {
            l=lower_bound(tmp+1,tmp+cnt+1,A[i].l)-tmp;
            r=lower_bound(tmp+1,tmp+cnt+1,A[i].r)-tmp;
            if(l<r)modify(1,l,r-1,A[i].ind);
            //cout<<T[1].sum<<endl; 
            f+=T[1].sum*(A[i+1].h-A[i].h);
        }
        printf("%.2f\n",f);
        printf("\n");
    }
    return 0;
}