hdu1542（线段树——矩形面积并）

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

分析：离散化+扫描线+线段树

#pragma comment(linker,"/STACK:102400000,102400000")
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>
#include <cstdlib>
#include <stack>
#include <vector>
#include <set>
#include <map>
#define LL long long
#define mod 10007
#define inf 0x3f3f3f3f
#define N 2015
#define FILL(a,b) (memset(a,b,sizeof(a)))
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
using namespace std;
struct line
{
    double l,r,h;
    int d;
    line(){}
    line(double l,double r,double h,int d):l(l),r(r),h(h),d(d){}
    bool operator<(const line &a)const
    {
        return h<a.h;
    }
}s[N];
double sum[N<<],has[N];
int col[N<<];
void Pushup(int l,int r,int rt)
{
    if(col[rt])sum[rt]=has[r+]-has[l];//每个点表示的是has[i]~has[i+1]的长度，故has[r+1]-has[l]
    else if(l==r)sum[rt]=;//节点没有子节点加上来，故清零
    else sum[rt]=sum[rt<<]+sum[rt<<|];//加上子孙节点被覆盖着的线段长度
}
void update(int L,int R,int d,int l,int r,int rt)
{
    if(L<=l&&r<=R)
    {
        col[rt]+=d;
        Pushup(l,r,rt);
        return;
    }
    int m=(l+r)>>;
    if(L<=m)update(L,R,d,lson);
    if(m<R)update(L,R,d,rson);
    Pushup(l,r,rt);
}
int bin(double key,double a[],int n)
{
    int l=,r=n-;
    while(l<=r)
    {
        int m=(l+r)>>;
        if(a[m]==key)return m;
        if(a[m]>key)r=m-;
        else l=m+;
    }
    return -;
}
int main()
{
    int n,cas=;
    double x1,y1,x2,y2;
    while(scanf("%d",&n)&&n)
    {
        int k=;
        for(int i=;i<n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            has[k]=x1;
            s[k++]=line(x1,x2,y1,);
            has[k]=x2;
            s[k++]=line(x1,x2,y2,-);
        }
        sort(s,s+k);
        sort(has,has+k);
        int m=;
        for(int i=;i<k;i++)//去重
            if(has[i]!=has[i-])has[m++]=has[i];
        double ans=;
        for(int i=;i<k;i++)
        {
            int L=bin(s[i].l,has,m);
            int R=bin(s[i].r,has,m)-;
            update(L,R,s[i].d,,m-,);//区间更新
            ans+=sum[]*(s[i+].h-s[i].h);//sum[1]表示被覆盖着的总长度
        }
        printf("Test case #%d\nTotal explored area: %.2lf\n\n",cas++,ans);
    }
}