poj1177Picture（线段树-周长并）

链接

神奇的扫描线啊

估计之前刷面积并的时候太急了 没来得及理解 。。

有一大段代码是与面积并一模一样的 都是离散化 更新

面积并是扫描的x  这次也一样

因为周长只算外围的 所以扫描到一条x边时 要减去上一次扫描到的 对于Y方向上 就是与面积并不大一样的地方了

看篇带图的讲解吧 网上的题解都大同小异 估计出自一人之手

http://blog.sina.com.cn/s/blog_691ce2b7010177dz.html

它是扫描的Y  其实方法是一样的 横过来就是

对于我的代码 自己感觉不好理解在于求区间数

 #include <iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<stdlib.h>
 using namespace std;
 #define N 10010
 #define LL long long
 struct node
 {
     int lx,rx,y,d;
     node(){}
     node(int a,int b,int c,int d):lx(a),rx(b),y(c),d(d){}
     bool operator < (const node &S) const
     {
         if(y == S.y) return d > S.d;
         return y<S.y;
     }
 }mat[N];
 int que[N],sum[N<<],lg[N<<],rg[N<<],seg[N<<];
 int cnt[N<<];
 LL ans=;
 int bin(int k,int low,int high)
 {
     int m;
     while(low<=high)
     {
         m  = (low+high)/;
         if(que[m]==k)
         return m;
         if(que[m]>k)
         high = m-;
         else
         low = m+;
     }
     return m;
 }
 void up(int l,int r,int w)
 {
     if(cnt[w])//入边就更新这段的值 及区间数 （1个区间 2条Y边）
     {
         sum[w] = que[r+]-que[l];
         seg[w] = ;
         lg[w] = ;
         rg[w] = ;
     }
     else if(l==r)
     {
         sum[w] = seg[w] = lg[w] = rg[w] = ;//出边的话 消除这条边
     }
     else//非出入边 向上更新 与平常线段树一样 等于子孩子相加
     {
         sum[w] = sum[w<<]+sum[w<<|];
         seg[w] = seg[w<<]+seg[w<<|];
         lg[w] = lg[w<<];
         rg[w] = rg[w<<|];
         if(lg[w<<|]&&rg[w<<])//左右相连 减去重复的
         seg[w]-=;
     }
 }
 void update(int a,int b,int d,int l,int r,int w)
 {
     if(a<=l&&b>=r)
     {
         cnt[w]+=d;
         up(l,r,w);
         return ;
     }
     int m = (l+r)>>;
     if(a<=m)
     update(a,b,d,l,m,w<<);
     if(b>m)
     update(a,b,d,m+,r,w<<|);
     up(l,r,w);
 }
 int main()
 {
     int i,n,a,b,c,d;
     while(scanf("%d",&n)!=EOF)
     {
         int num=;
         for(i = ; i <= n ; i++)
         {
             scanf("%d%d%d%d",&a,&b,&c,&d);
             que[num] = a;
             mat[num++] = node(a,c,b,);
             que[num] = c;
             mat[num++] = node(a,c,d,-);
         }
         sort(que,que+num);
         sort(mat,mat+num);
         int k = ;
         for(i =  ; i < num ; i++)
         {
             if(que[i]!=que[i-])
             que[k++] = que[i];
         }
         memset(cnt,,sizeof(cnt));
         memset(sum,,sizeof(sum));
         memset(lg,,sizeof(lg));
         memset(rg,,sizeof(rg));
         memset(seg,,sizeof(seg));
         int o=;
         ans=;
         mat[num].y = ;
         for(i =  ; i < num ; i++)
         {
             int l = bin(mat[i].lx,,k-);
             int r = bin(mat[i].rx,,k-)-;
             update(l,r,mat[i].d,,k-,);
             ans+=seg[]*(mat[i+].y-mat[i].y);
             ans+=abs(sum[]-o);
             o = sum[];
         }
         printf("%lld\n",ans);
     }
     return ;
 }