hdu 1255 覆盖的面积(矩形面积并,多次覆盖)

覆盖的面积

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


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
 

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

题意:就是求覆盖两次以上的图形的面积

分析:只要在求矩形面积并的基础上加上一个统计覆盖超过两次的长度,如果要求覆盖三次,那就再开一个,貌似有点土,不过好一点的方法没想到。。。

代码:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define ls rt<<1
#define rs rt<<1|1
#define lson l,m,ls
#define rson m,r,rs
using namespace std;
const int mm=2222;
const int mn=mm<<2;
struct seg
{
    double x,y1,y2;
    int val;
}g[mm];
double len[mn],ret[mn],y[mm],L,R;
int t[mn],val;
void build(int n)
{
    while(n--)ret[n]=len[n]=t[n]=0;
}
void updata(int l,int r,int rt)
{
    if(L<=y[l]&&R>=y[r])t[rt]+=val;
    else
    {
        int m=(l+r)>>1;
        if(L<y[m])updata(lson);
        if(R>y[m])updata(rson);
    }
    if(t[rt]>1)ret[rt]=len[rt]=y[r]-y[l];
    else if(t[rt]>0)len[rt]=y[r]-y[l],ret[rt]=len[ls]+len[rs];
    else if(l>=r)ret[rt]=len[rt]=0;
    else len[rt]=len[ls]+len[rs],ret[rt]=ret[ls]+ret[rs];
}
bool cmp(seg a,seg b)
{
    return a.x<b.x;
}
int main()
{
    int i,m,n,cs;
    double ans;
    scanf("%d",&cs);
    while(cs--)
    {
        scanf("%d",&n);
        for(i=0;i<n;++i)
        {
            scanf("%lf%lf%lf%lf",&g[i].x,&y[i],&g[i+n].x,&y[i+n]);
            g[i].y1=y[i],g[i].y2=y[i+n],g[i].val=1;
            g[i+n].y1=y[i],g[i+n].y2=y[i+n],g[i+n].val=-1;
        }
        sort(y,y+n+n);
        sort(g,g+n+n,cmp);
        for(m=i=0;i<n+n;++i)
            if(y[m]<y[i])y[++m]=y[i];
        build(m<<2);
        for(ans=i=0;i<n+n;++i)
        {
            L=g[i].y1,R=g[i].y2,val=g[i].val;
            updata(0,m,1);
            if(g[i].x<g[i+1].x)ans+=(g[i+1].x-g[i].x)*ret[1];
        }
        printf("%.2lf\n",ans);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值