扫面线+线段树（hdu1542）

之前写过这个算法，时间长了就忘掉了，，现在不看书自己努力回想起来，对算法的理解，对线段树的理解感觉也更深了一点（可能心理作用，哈哈哈）

思路简单说一下吧

从做到右遍历每一条矩阵的边（左右边），看该边对右边的面积贡献是正还是负（矩阵左边为正，右边为负），在y轴上用线段树维护在y轴的贡献值，与x轴上该边与下一条边的差值相乘即可

// acm1.cpp: 定义控制台应用程序的入口点。
//
#include<iostream>
#include<cstdio>
#include<cmath>
#include<queue>
#include<vector>
#include<string.h>
#include<cstring>
#include<algorithm>
#include<set>
#include<map>
#include<fstream>
#include<cstdlib>
#include<ctime>
using namespace std;
typedef long long ll;
typedef long long lint;
const ll mod=1e9+7;
const int maxn=1e2+7;
const int maxm=1e5+7;
const double eps=1e-3;
int ar1[maxn],ar2[maxn],m,n;
bool bo[maxn];
struct line{
    double y1,y2,x;
    int ma;
}li[maxn*2];

struct node{
    int l,r;
    int mark;
    double length;
}no[maxn*8];
double xl[maxn],yl[maxn];
bool cmp(line l1,line l2){
    return l1.x<l2.x;
}
void build(int x,int y,int o){
    no[o].l=x;
    no[o].r=y;
    no[o].length=0;
    no[o].mark=0;
    if(y-x==1)return ;
    int mid=(x+y)>>1;
    if(mid>x)build(x,mid,o*2);
    if(mid<y)build(mid,y,o*2+1);
}
void down(int x){
    if(no[x].mark>0)no[x].length=yl[no[x].r]-yl[no[x].l];
    else if(no[x].r==no[x].l+1){
        no[x].length=0;
    }
    else{
        no[x].length=no[x*2].length+no[x*2+1].length;
    }
}
void update(int o,line & liMid){
    double func1=yl[no[o].l],func2=yl[no[o].r];
    if(liMid.y1<=func1&&liMid.y2>=func2){
        no[o].mark+=liMid.ma;
    }
    else{
        double funcMid=yl[(no[o].l+no[o].r)>>1];
        if(liMid.y1<funcMid)update(o*2,liMid);
        if(liMid.y2>funcMid)update(o*2+1,liMid);
    }
    down(o);
}
int main()
{
    int fir=0;
    while(~scanf("%d",&n)&&n){
        for(int i=0;i<n;i++){
            double x1,x2,y1,y2;
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            li[i*2].x=x1;
            li[i*2].y1=y1;
            li[i*2].y2=y2;
            li[i*2].ma=1;
            li[i*2+1].x=x2;
            li[i*2+1].y1=y1;
            li[i*2+1].y2=y2;
            li[i*2+1].ma=-1;
            yl[i*2]=y1;
            yl[i*2+1]=y2;
        }
        sort(li,li+2*n,cmp);
        sort(yl,yl+2*n);
        int len=unique(yl,yl+2*n)-yl;
        build(0,len-1,1);
        double ans=0;
        for(int i=0;i<2*n-1;i++){
            update(1,li[i]);
            ans+=no[1].length*(li[i+1].x-li[i].x);
        }
        printf("Test case #%d\nTotal explored area: %.2lf\n\n",++fir,ans);
    }
    return 0;
}