HDU1255-覆盖的面积

覆盖的面积

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

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

这题和HDU1542几乎完全一样，链接地址：http://blog.csdn.net/a664607530/article/details/52554335有详细讲，这个题目唯一的不同就是这是个求重合的面积，而HDU1542是求并面积，只要将cover>0改成cover>1就可以求出覆盖面积了，我两题的扫描方式不同

#include <iostream>
#include <algorithm>
#include <cstdio>

using namespace std;

const int maxn=1100;

struct LINE
{
    double  x,y_down,y_up;
    int  flag;
    friend bool operator<(LINE a,LINE b)
    {
        return  a.x<b.x;
    }
} line[2*maxn];

struct TREE
{
    double y_down,y_up;
    double x;
    int cover;
    bool flag;
} tree[1000*maxn];

int n;
double x1,y1,x2,y2;
int sum;
double y[2*maxn];

void build(int k,int l,int r)
{
    tree[k].x=-1;
    tree[k].cover=0;
    tree[k].y_down=y[l];
    tree[k].y_up=y[r];
    tree[k].flag=false;
    if(l+1==r)
    {
        tree[k].flag=true;
        return;
    }
    int mid=(l+r)>>1;
    build(2*k,l,mid);
    build(2*k+1,mid,r);
}

double query(int k,double x,double l,double r,int flag)
{
    if(r<=tree[k].y_down||l>=tree[k].y_up)
        return 0;
    if(tree[k].flag)
    {
        if(tree[k].cover>1)
        {
            double temp=tree[k].x;
            double ans=(x-temp)*(tree[k].y_up-tree[k].y_down);
            tree[k].x=x;
            tree[k].cover+=flag;
            return ans;
        }
        else
        {
            tree[k].cover+=flag;
            tree[k].x=x;
            return 0;
        }
    }
    double ans1,ans2;
    ans1=query(2*k,x,l,r,flag);
    ans2=query(2*k+1,x,l,r,flag);
    return ans1+ans2;
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        sum=1;
        for(int i=1; i<=n; i++)
        {
            scanf("%lf %lf %lf %lf",&x1,&y1,&x2,&y2);
            y[sum]=y1;
            line[sum].x=x1;
            line[sum].y_down=y1;
            line[sum].y_up=y2;
            line[sum].flag=1;
            sum++;
            y[sum]=y2;
            line[sum].x=x2;
            line[sum].y_down=y1;
            line[sum].y_up=y2;
            line[sum].flag=-1;
            sum++;
        }
        sort(y+1,y+sum);
        sort(line+1,line+sum);
        build(1,1,sum-1);
        double ans=0;
        for(int i=1; i<sum; i++)
            ans+=query(1,line[i].x,line[i].y_down,line[i].y_up,line[i].flag);
        printf("%.2lf\n",ans);
    }
    return 0;
}