HDU 1255 覆盖的面积（线段树+扫描线+离散化）

覆盖的面积

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

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

 

题意：

给你几个矩形，求它们的共同面积。

POINT：

cnt[x]即代表x区间的线段存在了几次，>=2即是重合了一次或多次。 用法和 HDU 1542相同。

然后每次计算就询问一次。乘以高度就行。

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <map>
#include <string.h>
#include <algorithm>
#include <vector>
#include <fstream>
using namespace std;
#define lt 2*x
#define rt 2*x+1
#define LL long long
const double eps = 1e-9;
const int N = 1001*10;
struct node
{
    double l,r,h;
    int f;
}len[N];
double sum[N];
double Hash[N];
int addcnt[N];
int cnt[N],xnum,lnum;
void init()
{
    lnum=0,xnum=0;
    memset(sum,0,sizeof sum);
    memset(len,0,sizeof len);
    memset(addcnt,0,sizeof addcnt);
    memset(cnt,0,sizeof cnt);
    memset(Hash,0,sizeof Hash);
}
bool cmd(node a,node b)
{
    return a.h<b.h;
}

void pushdown(int x,int l,int r)
{

    addcnt[lt]+=addcnt[x];
    addcnt[rt]+=addcnt[x];
    cnt[lt]+=addcnt[x];
    cnt[rt]+=addcnt[x];
    addcnt[x]=0;

    
}
void add(int x,int l,int r,int ll,int rr,int f)
{
    if(addcnt[x]) pushdown(x,l,r);
    if(ll<=l&&rr>=r)
    {
        cnt[x]+=f;
        addcnt[x]+=f;
    }
    else
    {
        int mid=(l+r)>>1;
        if(ll<=mid) add(lt,l,mid,ll,rr,f);
        if(mid<rr) add(rt,mid+1,r,ll,rr,f);
    }
}
double query(int x,int l,int r)
{
    if(addcnt[x]) pushdown(x,l,r);
    double ans=0;
    if(cnt[x]>=2)   ans+=Hash[r+1]-Hash[l];
    else if(l==r)
    {
        ans+=0;
    }
    else
    {
        int mid=(l+r)>>1;
        ans+=query(lt,l,mid);
        ans+=query(rt,mid+1,r);
    }
    return ans;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        init();
        int n;
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            double x1,x2,y1,y2;
            scanf("%lf %lf %lf %lf",&x1,&y1,&x2,&y2);
            Hash[++xnum]=x1,Hash[++xnum]=x2;
            len[++lnum].f=1,len[lnum].h=y1,len[lnum].l=x1,len[lnum].r=x2;
            len[++lnum].f=-1,len[lnum].h=y2,len[lnum].l=x1,len[lnum].r=x2;
        }
        sort(len+1,len+1+lnum,cmd);
        sort(Hash+1,Hash+1+xnum);
        int m=1;
        for(int i=2;i<=xnum;i++)
        {
            if(Hash[i-1]!=Hash[i]) Hash[++m]=Hash[i];
        }
        double ans=0;
        for(int i=1;i<=lnum;i++)
        {
            double now;
            int l=lower_bound(Hash+1,Hash+1+m,len[i].l)-Hash;
            int r=lower_bound(Hash+1,Hash+1+m,len[i].r)-Hash-1;
            add(1,1,m,l,r,len[i].f);
            now=query(1,1,m);
            ans+=now*(len[i+1].h-len[i].h);
        }
        printf("%.2lf\n",ans+eps);
    }

    return 0;
}