poj1177Picture（线段树-周长并）

链接

神奇的扫描线啊

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

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

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

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

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

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

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

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

 

  1 #include <iostream>
  2 #include<cstdio>
  3 #include<cstring>
  4 #include<algorithm>
  5 #include<stdlib.h>
  6 using namespace std;
  7 #define N 10010
  8 #define LL long long
  9 struct node
 10 {
 11     int lx,rx,y,d;
 12     node(){}
 13     node(int a,int b,int c,int d):lx(a),rx(b),y(c),d(d){}
 14     bool operator < (const node &S) const
 15     {
 16         if(y == S.y) return d > S.d;
 17         return y<S.y;
 18     }
 19 }mat[N];
 20 int que[N],sum[N<<2],lg[N<<2],rg[N<<2],seg[N<<2];
 21 int cnt[N<<2];
 22 LL ans=0;
 23 int bin(int k,int low,int high)
 24 {
 25     int m;
 26     while(low<=high)
 27     {
 28         m  = (low+high)/2;
 29         if(que[m]==k)
 30         return m;
 31         if(que[m]>k)
 32         high = m-1;
 33         else
 34         low = m+1;
 35     }
 36     return m;
 37 }
 38 void up(int l,int r,int w)
 39 {
 40     if(cnt[w])//入边就更新这段的值 及区间数 （1个区间 2条Y边）
 41     {
 42         sum[w] = que[r+1]-que[l];
 43         seg[w] = 2;
 44         lg[w] = 1;
 45         rg[w] = 1;
 46     }
 47     else if(l==r)
 48     {
 49         sum[w] = seg[w] = lg[w] = rg[w] = 0;//出边的话 消除这条边
 50     }
 51     else//非出入边 向上更新 与平常线段树一样 等于子孩子相加
 52     {
 53         sum[w] = sum[w<<1]+sum[w<<1|1];
 54         seg[w] = seg[w<<1]+seg[w<<1|1];
 55         lg[w] = lg[w<<1];
 56         rg[w] = rg[w<<1|1];
 57         if(lg[w<<1|1]&&rg[w<<1])//左右相连 减去重复的
 58         seg[w]-=2;
 59     }
 60 }
 61 void update(int a,int b,int d,int l,int r,int w)
 62 {
 63     if(a<=l&&b>=r)
 64     {
 65         cnt[w]+=d;
 66         up(l,r,w);
 67         return ;
 68     }
 69     int m = (l+r)>>1;
 70     if(a<=m)
 71     update(a,b,d,l,m,w<<1);
 72     if(b>m)
 73     update(a,b,d,m+1,r,w<<1|1);
 74     up(l,r,w);
 75 }
 76 int main()
 77 {
 78     int i,n,a,b,c,d;
 79     while(scanf("%d",&n)!=EOF)
 80     {
 81         int num=0;
 82         for(i = 1; i <= n ; i++)
 83         {
 84             scanf("%d%d%d%d",&a,&b,&c,&d);
 85             que[num] = a;
 86             mat[num++] = node(a,c,b,1);
 87             que[num] = c;
 88             mat[num++] = node(a,c,d,-1);
 89         }
 90         sort(que,que+num);
 91         sort(mat,mat+num);
 92         int k = 1;
 93         for(i = 1 ; i < num ; i++)
 94         {
 95             if(que[i]!=que[i-1])
 96             que[k++] = que[i];
 97         }
 98         memset(cnt,0,sizeof(cnt));
 99         memset(sum,0,sizeof(sum));
100         memset(lg,0,sizeof(lg));
101         memset(rg,0,sizeof(rg));
102         memset(seg,0,sizeof(seg));
103         int o=0;
104         ans=0;
105         mat[num].y = 0;
106         for(i = 0 ; i < num ; i++)
107         {
108             int l = bin(mat[i].lx,0,k-1);
109             int r = bin(mat[i].rx,0,k-1)-1;
110             update(l,r,mat[i].d,0,k-1,1);
111             ans+=seg[1]*(mat[i+1].y-mat[i].y);
112             ans+=abs(sum[1]-o);
113             o = sum[1];
114         }
115         printf("%lld\n",ans);
116     }
117     return 0;
118 }
View Code

 

posted @ 2013-09-13 19:38  _雨  阅读(290)  评论(0)  编辑  收藏  举报