HDU 1255 覆盖的面积 ( 扫描线 + 离散 求矩阵大于k次面积并 )

覆盖的面积

Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 4040    Accepted Submission(s): 1995

Problem Description

给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积.

 

Input

输
入数据的第一行是一个正整数T(1<=T<=100),代表测试数据的数量.每个测试数据的第一行是一个正整数N(1<=N&
lt;=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

 

 

交G++迷之WA。。无语

C++可过。

 

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>
#include <map>
using namespace std;
#define root 1,(n<<1)+10,1
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define lr rt<<1
#define rr rt<<1|1

const int N =  ;

struct Point {
    double x , y1 , y2 ;
    int v ;
    bool operator < ( const Point &a ) const {
        return x < a.x ;
    }
}p[N];
int n , tot , tt ;
vector<double>e;
double c[N<<] ;

void addpoint( double x , double y1 , double y2 , int v ) {
    p[tot].x = x , p[tot].y1 = y1 , p[tot].y2 = y2 , p[tot].v = v , tot++ ;
}
map<double,int>mp;

void lisan() {
    mp.clear();
    sort( e.begin() , e.end() ) ;
    tt =  ;
    for( int i =  ; i < e.size() ; ++i ) {
        if( e[i] != e[i-] ) e[tt++] = e[i] ;
    }
    for( int i =  ; i < tt ; ++i ) {
        c[i+] = e[i] ;
        mp[e[i]] = i+ ;
    }
}

int cnt[N<<] , lazy[N<<] ;

double date[N<<] ;

void build( int l , int r , int rt ) {
    cnt[rt] = lazy[rt] =  ;
    if( l == r ) {
        date[rt] = c[l+] - c[l];
        return ;
    }
    int mid = (l+r)>>;
    build(lson) , build(rson) ;
    date[rt] = date[lr] + date[rr];
}

void Down( int rt ) {
    if( lazy[rt] ) {
        cnt[lr] += lazy[rt] ;
        cnt[rr] += lazy[rt] ;
        lazy[lr] += lazy[rt] ;
        lazy[rr] += lazy[rt] ;
        lazy[rt] =  ;
    }
}

void Up( int rt ) {
    cnt[rt] = min( cnt[lr] , cnt[rr] );
}

void update( int l , int r , int rt , int L , int R , int v ) {
    if( l == L && r == R ) {
        cnt[rt] += v ; lazy[rt] += v ;
        return ;
    }
    int mid = (l+r)>> ;
    Down(rt) ;
    if( R <= mid ) update(lson,L,R,v) ;
    else if( L > mid ) update( rson,L,R,v);
    else update(lson,L,mid,v),update(rson,mid+,R,v);
    Up(rt);
}

double query( int l , int r , int rt , int L , int R ) {
    if( cnt[rt] >  ) return date[rt] ;
    if( l == r ) return  ;
    Down(rt);
    int mid = (l+r)>>;
    if( R <= mid ) return  query(lson,L,R);
    else if( L > mid ) return query(rson,L,R);
    else return query(lson,L,mid) + query(rson,mid+,R);
}

int main() {
    #ifdef LOCAL
        freopen("in.txt","r",stdin);
    #endif // LOCAL
    int _ ; scanf("%d",&_);
    while( _-- ) {
        scanf("%d",&n);
        e.clear();
        tot =  ;
        for( int i =  ; i < n ; ++i ) {
            double x1 , y1 , x2 , y2 ;
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            if( y1 > y2 ) swap( y1 , y2 ) ;
            addpoint( x1 , y1 , y2 ,  );
            addpoint( x2 , y1 , y2 , - );
            e.push_back(y1);
            e.push_back(y2);
        }
        lisan();
        sort( p , p + tot ) ;
        build(,tt-,);
        double ans =  ;
        update(,tt-,,mp[p[].y1],mp[p[].y2]-,p[].v);
        for( int i =  ; i < tot ; ++i ) {
            ans += ( p[i].x - p[i-].x ) * query(,tt-,,,tt-) ;
            update(,tt-,,mp[p[i].y1],mp[p[i].y2]-,p[i].v);
        }
        printf("%.2lf\n",ans);
    }
    return  ;
}