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;
}