hdu1255 覆盖的面积 扫描线+线段树+离散化

覆盖的面积

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

Total Submission(s): 3464    Accepted Submission(s): 1689

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   |   We have carefully selected several similar problems for you:   1542  1540  1828  1823  3016 

 

题意：如题

思路：同poj1151，只是要维护两个变量len1(该段全覆盖1次的长度)，len2(该段全覆盖>=2次的长度)就行了，详见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define L(x) (x<<1)
#define R(x) (x<<1|1)
const int MAXN=2000+100;
int n,edge_cnt;
double y[MAXN];
struct node
{
    int l,r,lazy;
    double lf,rf,len1,len2;
}segtree[MAXN<<2];
struct Line
{
    int flag;
    double x,y1,y2;
}line[MAXN];
bool cmp(Line a,Line b)
{
    return a.x<b.x;
}
void build(int rt,int l,int r)
{
    segtree[rt].l=l; segtree[rt].r=r;
    segtree[rt].lazy=segtree[rt].len1=segtree[rt].len2=0;
    segtree[rt].lf=y[l]; segtree[rt].rf=y[r];
    if(l+1 == r)
        return ;
    int mid=(l+r)>>1;
    build(L(rt),l,mid); build(R(rt),mid,r);
}
void cal(int rt)
{
    if(segtree[rt].lazy <= -2){
        segtree[rt].len2=segtree[rt].rf-segtree[rt].lf;
        return ;
    }
    if(segtree[rt].lazy == -1){
        segtree[rt].len1=segtree[rt].rf-segtree[rt].lf;
        if(segtree[rt].l+1 == segtree[rt].r )
            segtree[rt].len2=0;
        else
            segtree[rt].len2=segtree[L(rt)].len1+segtree[R(rt)].len1; // 注意本身lazy=-1说明这段已经全覆盖，子段lazy=-1对于本段就是有两次覆盖了
    }
    else {
        if(segtree[rt].l+1 == segtree[rt].r )
            segtree[rt].len1=segtree[rt].len2=0;
        else {
            segtree[rt].len2=segtree[L(rt)].len2+segtree[R(rt)].len2;
            segtree[rt].len1=segtree[L(rt)].len1+segtree[R(rt)].len1;
        }
    }
}
void update(int rt,Line e)
{
    if( segtree[rt].lf == e.y1 && segtree[rt].rf == e.y2){
        segtree[rt].lazy+=e.flag; cal(rt);
        return ;
    }
    if(e.y1>=segtree[R(rt)].lf)
        update(R(rt),e);
    else if(e.y2<=segtree[L(rt)].rf)
        update(L(rt),e);
    else{
        Line temp=e;
        temp.y2=segtree[L(rt)].rf;
        update(L(rt),temp);
        temp=e;
        temp.y1=segtree[R(rt)].lf;
        update(R(rt),temp);
    }
    cal(rt);
}
int main()
{
    //freopen("text.txt","r",stdin);
    int T;
    double x1,y1,x2,y2;
    scanf("%d",&T);
    while(T--)
    {
        edge_cnt=0;
        scanf("%d",&n);
        for(int i=1;i<=n;i++){
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            line[++edge_cnt].x=x1; line[edge_cnt].y1=y1;
            line[edge_cnt].y2=y2; line[edge_cnt].flag=-1; y[edge_cnt]=y1;
            line[++edge_cnt].x=x2; line[edge_cnt].y1=y1;
            line[edge_cnt].y2=y2; line[edge_cnt].flag=1; y[edge_cnt]=y2;
        }
        sort(line+1,line+edge_cnt+1,cmp);
        sort(y+1,y+edge_cnt+1);
        build(1,1,edge_cnt);
        update(1,line[1]);
        double ans=0;
        for(int i=2;i<=edge_cnt;i++){
            ans+=segtree[1].len2*(line[i].x-line[i-1].x);
            update(1,line[i]);
        }
        printf("%.2lf\n",ans);
    }
    return 0;
}