扫面线+线段树（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+;
 const int maxn=1e2+;
 const int maxm=1e5+;
 const double eps=1e-;
 int ar1[maxn],ar2[maxn],m,n;
 bool bo[maxn];
 struct line{
     double y1,y2,x;
     int ma;
 }li[maxn*];

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