hdu 1255 覆盖的面积

覆盖的面积

http://acm.hdu.edu.cn/showproblem.php?pid=1255

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

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 1698 1166 1540 1394 

 

题意：求矩形交的面积

线段树维护扫描线，记录区间完全覆盖1次的长度sum1、完全覆盖2次的长度sum2

如果区间覆盖次数>=2，sum1=sum2=区间长度

如果区间只覆盖一次，sum1=区间长度，sum2=左右子区间的sum1之和（自己本身完全覆盖了一次，左右子区间只需覆盖一次即可）

如果区间没有覆盖，sum1，sum2=左右子区间的和

 

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1001
using namespace std;
int n,opl,opr,w;
double ans,x1,x2,y3,y4,has[N<<1],sum1[N<<4],sum2[N<<4];
int f[N<<4];
struct node
{
    double l,r,h;
    int f;
    bool operator < (node p)const
    {
        return h<p.h;
    } 
}e[N<<1];
void up(int k,int l,int r)
{
    if(f[k]>=2) sum2[k]=has[r+1]-has[l];
    else if(f[k]==1) 
    {
        sum1[k]=has[r+1]-has[l];
        sum2[k]=sum1[k<<1]+sum1[k<<1|1];
    }
    else
    {
        sum1[k]=sum1[k<<1]+sum1[k<<1|1];
        sum2[k]=sum2[k<<1]+sum2[k<<1|1];
    }
}
void change(int k,int l,int r)
{
    if(opl<=l && r<=opr)
    {
        f[k]+=w;
        up(k,l,r);
        return;
    }
    int mid=l+r>>1;
    if(opl<=mid) change(k<<1,l,mid);
    if(opr>mid) change(k<<1|1,mid+1,r);
    up(k,l,r);
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y3,&x2,&y4);
            has[i*2-1]=x1; has[i*2]=x2;
            e[i*2-1].l=x1; e[i*2-1].r=x2; e[i*2-1].f=-1; e[i*2-1].h=y4;
            e[i*2].l=x1;   e[i*2].r=x2;   e[i*2].f=1;    e[i*2].h=y3;
        }
        sort(e+1,e+n*2+1);
        sort(has+1,has+n*2+1);
        int cnt=unique(has+1,has+n*2+1)-(has+1);
        for(int i=1;i<=n*2;i++)
        {
            opl=lower_bound(has,has+cnt+1,e[i].l)-has;
            opr=lower_bound(has,has+cnt+1,e[i].r)-has-1;
            w=e[i].f;
            change(1,1,cnt);
            ans+=sum2[1]*(e[i+1].h-e[i].h);
        }
        printf("%.2lf
",ans);
        memset(f,0,sizeof(f));
        memset(sum1,0,sizeof(sum1));
        memset(sum2,0,sizeof(sum2));
        ans=0;
    }
}