hdu 1255 覆寫的面積（線段樹+離散化+掃描線）

覆寫的面積

Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)

Total Submission(s): 2993    Accepted Submission(s): 1473

Problem Description 給定平面上若幹矩形,求出被這些矩形覆寫過至少兩次的區域的面積.

Input 輸入資料的第一行是一個正整數T(1<=T<=100),代表測試資料的數量.每個測試資料的第一行是一個正整數N(1<=N<=1000),代表矩形的數量,然後是N行資料,每一行包含四個浮點數,代表平面上的一個矩形的左上角坐标和右下角坐标,矩形的上下邊和X軸平行,左右邊和Y軸平行.坐标的範圍從0到100000.

注意:本題的輸入資料較多,推薦使用scanf讀入資料.

Output 對于每組測試資料,請計算出被這些矩形覆寫過至少兩次的區域的面積.結果保留兩位小數.

Sample Input

2
5
1 1 4 2
1 3 3 7
2 1.5 5 4.5
3.5 1.25 7.5 4
6 3 10 7
3
0 0 1 1
1 0 2 1
2 0 3 1
        

Sample Output

7.63
0.00
        

Author Ignatius.L & weigang Lee  

Recommend Ignatius.L  

一天一棵樹...

題解：線段樹的離散化和掃描線，搞了好久...沒精力說了...看代碼吧...

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<map>
using namespace std;
map<int,double>mip;
map<double,int>mip2;
struct line{
    double x,y1,y2;
    int flag;
}seg[2008];
struct Tree{
    double one,zero,more;
    int cover;
}tree[8008];
int n,all,cou;
double c[5],hasy[2008];
int cmp(const void *a,const void *b)
{
    return *(double *)a>*(double *)b?1:-1;
}
int cmp2(const void *a,const void *b)
{
    struct line c=*(struct line *)a;
    struct line d=*(struct line *)b;
    return c.x>d.x?1:-1;
}
void build(int l,int r,int pos)
{
    int mid=(l+r)>>1;

    tree[pos].zero=mip[r]-mip[l];
    tree[pos].one=tree[pos].more=tree[pos].cover=0;
    if(l+1==r) return;
    build(l,mid,pos<<1);
    build(mid,r,pos<<1|1);
}
void jerry(int l,int r,int pos)
{
    if(l+1==r)
    {
        if(tree[pos].cover>1)
        {
            tree[pos].more=mip[r]-mip[l];
            tree[pos].one=tree[pos].zero=0;
        }
        else if(tree[pos].cover==1)
        {
            tree[pos].one=mip[r]-mip[l];
            tree[pos].more=tree[pos].zero=0;
        }
        else
        {
            tree[pos].more=tree[pos].one=0;
            tree[pos].zero=mip[r]-mip[l];
        }
        return;
    }
    if(tree[pos].cover==0)
    {
        tree[pos].one=tree[pos<<1].one+tree[pos<<1|1].one;
        tree[pos].more=tree[pos<<1].more+tree[pos<<1|1].more;
        tree[pos].zero=mip[r]-mip[l]-tree[pos].one-tree[pos].more;
    }
    else if(tree[pos].cover==1)
    {
        tree[pos].more=tree[pos<<1].one+tree[pos<<1|1].one;
        tree[pos].more+=tree[pos<<1].more+tree[pos<<1|1].more;
        tree[pos].one=tree[pos<<1].zero+tree[pos<<1|1].zero;
        tree[pos].zero=mip[r]-mip[l]-tree[pos].one-tree[pos].more;
    }
    else
    {
        tree[pos].more=mip[r]-mip[l];
        tree[pos].one=tree[pos].zero=0;
    }
}
void update(int l,int r,int pos,int templ,int tempr,int val)
{
    int mid=(l+r)>>1;

    if(templ<=l&&r<=tempr)
    {
        tree[pos].cover+=val;
        jerry(l,r,pos);
        return;
    }
    if(templ>=mid) update(mid,r,pos<<1|1,templ,tempr,val);
    else if(tempr<=mid) update(l,mid,pos<<1,templ,tempr,val);
    else
    {
        update(l,mid,pos<<1,templ,mid,val);
        update(mid,r,pos<<1|1,mid,tempr,val);
    }
    jerry(l,r,pos);
}
int main()
{
    int t,i,j;
    double res;

    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        mip.clear();
        for(cou=all=i=0;i<n;i++)
        {
            for(j=0;j<4;j++) scanf("%lf",c+j);

            hasy[cou]=seg[all].y1=seg[all+1].y1=c[1];
            cou++;
            hasy[cou]=seg[all].y2=seg[all+1].y2=c[3];
            cou++;
            seg[all].x=c[0];  seg[all].flag=1;
            seg[all+1].x=c[2];   seg[all+1].flag=-1;
            all+=2;
        }
        qsort(hasy,cou,sizeof(hasy[0]),cmp);
        qsort(seg,all,sizeof(seg[0]),cmp2);

        mip[1]=hasy[0];
        mip2[hasy[0]]=1;
        for(cou=i=1;i<all;i++)
        {
            if(hasy[i]==hasy[i-1]) continue;
            mip[++cou]=hasy[i];
            mip2[hasy[i]]=cou;
        }
        build(1,cou,1);
        seg[all].x=seg[all-1].x;
        for(res=i=0;i<all;i++)
        {
            update(1,cou,1,mip2[seg[i].y1],mip2[seg[i].y2],seg[i].flag);
            res+=tree[1].more*(seg[i+1].x-seg[i].x);
        }
        printf("%.2lf\n",res);
    }

    return 0;
}