hdu 1542

链接：http://acm.hdu.edu.cn/showproblem.php?pid=1542

题意：

求所给矩形的覆盖面积

题解：

利用扫描线的思想，先将坐标离散化，之后以y轴分成多个矩形求解，可以让下边界+1上边界-1

问题就转化为了：求区间中有多少个非0数，要求支持区间+1 -1操作

我们可以通过维护区间最小值以及最小值的个数来完成这件事情

代码：

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include <iostream>
#define maxn 5000
#define INF 59999999
#define eps 1e-6
#define mid (h+t)/2
using namespace std;
struct re
{
    double a;int b;
}a[maxn],b[maxn];
struct ree
{
    int h,t,x,lazy;
    double sum,tot;
}p[maxn*];
double c[maxn];
int a1[maxn];
bool cmp(re x,re y)
{
  if (x.a<y.a) return(true); else return(false);
}
void build(int x,int h,int t)
{
    p[x].h=h; p[x].t=t; p[x].x=;
    if (h==t)
    {
      p[x].sum=p[x].tot=c[h+]-c[h];
      return;
    }
    build(x*,h,mid);
    build(x*+,mid+,t);
    p[x].sum=p[x].tot=p[x*].sum+p[x*+].sum;
}
void updata(int x)
{
        p[x].x=min(p[x*].x,p[x*+].x);
        if (p[x*].x==p[x*+].x)
        {
            p[x].sum=p[x*].sum+p[x*+].sum;
        }
        else
        {
            if (p[x*].x<p[x*+].x) p[x].sum=p[x*].sum;
            else p[x].sum=p[x*+].sum;
        }
        return;
}
void down(int x)
{
    if (p[x].lazy==) return;
    p[x].x+=p[x].lazy;
    p[x*].lazy+=p[x].lazy;
    p[x*+].lazy+=p[x].lazy;
    p[x].lazy=;
}
void insert(int x,int h,int t,int sum)
{
    down(x);
    if (p[x].h>t|| p[x].t<h) return;
    if (h<=p[x].h &&p[x].t<=t)
    {
      p[x].lazy+=sum; down(x); return;
    }
    insert(x*,h,t,sum); insert(x*+,h,t,sum);
    updata(x);
}
double query(int x,int h,int t)
{
    down(x);
    if (p[x].h>t||p[x].t<h) return();
    if (h<=p[x].t && p[x].t<=t)
    {
        if (p[x].x==) return(p[x].tot-p[x].sum); else return(p[x].tot);
    }
    return(query(x*,h,t)+query(x*+,h,t));
}
int main()
{int n,o=;
while (cin>>n&&n!=)
{
    o++;
  memset(p,,sizeof(p));
  for (int i=;i<=*n;i++) a[i].b=i,b[i].b=i;
  for (int i=;i<=n;i++)
  {
      cin>>a[*i-].a>>b[*i-].a>>a[*i].a>>b[*i].a;
  }
  sort(a+,a++*n,cmp);
  sort(b+,b++*n,cmp);
  int ll=; a[].a=INF;
  for (int i=;i<=*n;i++)
  {
      if (abs(a[i].a-a[i-].a)>eps) ll++;
      a1[a[i].b]=ll;
      c[ll]=a[i].a;
  }
  double ans=;
    build(,,ll-);
  for (int i=;i<*n;i++)
  {
      int pp,tmp=b[i].b;
      if (tmp%==) pp= ;else pp=-;
    insert(,a1[(tmp+)/*-],a1[(tmp+)/*]-,pp);
    ans+=(b[i+].a-b[i].a)*query(,,ll-);
  }
  cout<<"Test case #"<<o<<endl<<"Total explored area: " ;
  printf("%.2f\n\n",ans);
  }
}