线段树区间扫描线（hdu 1542）

题意描述：

给定一个矩阵，输入为左上角和右下角，现在要计算出所有矩阵重叠之后的面积总和。

解题思路：

根据x轴或者y轴建树，不过数据如果够大，必须先离散化，利用扫描线记录当前区间已经被覆盖的长度，那么现在覆盖的长度乘以下一条线段的距离，这个面积则是我们需要累加的面积。

其中的测试数据计算过程如下：

#include<iostream>
#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
#include<string.h>
using std::sort;
using std::unique;
const int N = 2222;
struct node
{
       double y,lx,rx;
       int s;
       bool operator <(const node &tmp)const{
           return y<tmp.y;
       }
}line[N<<4];
double X[N<<4],sum[N<<4];
int cnt[N<<4];
void update(int t,int l,int r,int L,int R,int val)
{
    if(L<=l&&r<=R)
    {

        cnt[t]+=val;
        if(cnt[t])sum[t]=X[r+1]-X[l];
        else
           if(l==r)sum[t]=0;
           else sum[t]=sum[t<<1]+sum[t<<1|1];
        return ;
    }
    int m=(l+r)>>1;
    if(L<=m)update(t<<1,l,m,L,R,val);
    if(R>m)update(t<<1|1,m+1,r,L,R,val);
    if(cnt[t])sum[t]=X[r+1]-X[l];
        else
           if(l==r)sum[t]=0;
           else sum[t]=sum[t<<1]+sum[t<<1|1];
}
int Fin(double k,int len)
{
    int l=0,r=len-1;
    while(l<r)
    {
         int m=(l+r)>>1;
         if(X[m]==k)return m;
         if(X[m]>k)r=m-1;
         else l=m+1;
    }
    return l;
}
int main()
{
    int n,T=0;
    double x1,x2,y1,y2;
    while(scanf("%d",&n)&&n)
    {
        int top=0,i,j;
        for(i=0,j=0;i<n;i++,j+=2)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            X[top++]=x1,X[top++]=x2;
            line[j].y=y1,line[j].lx=x1,line[j].rx=x2,line[j].s=1;
            line[j+1].y=y2,line[j+1].lx=x1,line[j+1].rx=x2,line[j+1].s=-1;
        }
        sort(X,X+top);
        top=unique(X,X+top)-X;
        sort(line,line+j);
        int l,r;
        double answer=0;
        memset(cnt,0,sizeof(cnt));
        memset(sum,0,sizeof(sum));
        for(int i=0;i<j-1;i++)
        {
            l=Fin(line[i].lx,top);
            r=Fin(line[i].rx,top)-1;
            if(r>=l)update(1,0,top-1,l,r,line[i].s);
            answer+=sum[1]*(line[i+1].y-line[i].y);
        }
        printf("Test case #%d\nTotal explored area: %.2lf\n\n",++T,answer);
    }
    return 0;
}