程序博客网 > 苹果下软件赚钱的软件

线段树入门&lazy思想

来源：互联网 发布：苹果下软件赚钱的软件 编辑：程序博客网 时间：2024/06/13 22:40

转自：http://blog.csdn.net/u012469987/article/details/41357377


线段树将区间分成若干个子区间，子区间又继续分，直到区间为一个点（区间左值等于右值）

对于父区间[a,b]，其子区间为[a,(a+b)/2]和[(a+b)/2+1,b] 

用于求区间的值，如区间最值、区间的和等。

代码实现中，约定结点下标从1开始，所以某结点下标为x，那么左儿子下标为2x，右儿子下标为2x+1，父结点下标为x/2。



常用符号符号等价意义rt<<1rt*2左子树的编号
rt<<1|1
rt*2+1右子树的编号
 (l+r)>>1
(l+r)/2区间长度的一半


常用宏定义

[cpp] view plain copy
print?
  1. #define Mid ((l+r)>>1)        //注意括号  
  2. #define lson rt<<1,l,Mid  //左结点  
  3. #define rson rt<<1|1,Mid+1,r  //右结点  


建树

建树丛根结点开始，递归建立左右子树，直到叶子结点，然后反向赋值，父结点的值 = F(左结点的值，右结点的值），这个F是依据题意变的，如果是区间最大则为max()

[cpp] view plain copy
print?
  1. void build(int rt,int l,int r)      //建编号为rt的区间为[l,r]的树，主函数传进来的固定是(1,1,n)  
  2. {  
  3.     if(l==r){            //叶子结点赋初值，注意下标，Max的是编号，val原数组的是l，看图可以理解  
  4.         Max[rt] = val[l];  
  5.     }else{                      //建左右子树  
  6.         build(lson);  
  7.         build(rson);  
  8.         Max[rt] = max( Max[rt<<1], Max[rt<<1|1]);   //父结点Max值为Max(左子结点，右子结点）  
  9.     }  
  10. }  


查询

查询为区间查询（只是查询某个点的话不需要线段树），即在区间里查询某个特性值，每次查询都是从跟结点开始往下，根据查询区间和当前区间的区间位置判断是要去左右子区间查询，还是直接返回。如果被查询区间是查询区间的子区间则直接返回子区间的值，如在[1,6]里查询[1,12]就返回[1,6]的值，不再往下查询。
[cpp] view plain copy
print?
  1. void query(int rt,int l,int r,int L,int R)  //在[l,r]里查询[L,R]的值，[L,R]一直不变，[l,r]变  
  2. {  
  3.     if(L <= l && r <= R){     
  4.         ans1 = max(ans1,Max[rt]);  
  5.         ans2 = min(ans2,Min[rt]);  
  6.     }else{  
  7.         if( L <= Mid) //查询区间在当前区间的左半区间有内容，如在[1,6]里查询[2,3]  
  8.             query(lson,L,R);      
  9.         if( R > Mid)  //同理去右子区间，注意不能有else，因为有横跨左右的情况，如[1,6]里查询[2,5]  
  10.             query(rson,L,R);  
  11.     }  
  12. }  


更新

更新分为单点更新和区间更新，区间更新等会在下面讲述，而单点更新跟普通区间查询差不多

[cpp] view plain copy
print?
  1. void update(int rt,int l,int r,int pos,int num)  
  2. {  
  3.     if(l == r && r == pos){     //到对应的叶结点  
  4.         Max[rt] = num;  
  5.     }else{  
  6.         if( pos <= Mid)  
  7.             update(lson,pos,num);  
  8.         if( pos > Mid)          //或者直接else，点不可能同时在两个区间里  
  9.             update(rson,pos,num);  
  10.         Max[rt] = max( Max[rt<<1], Max[rt<<1|1]);  
  11.     }  
  12. }  

数组大小要为原数据范围的4倍，证明点这里

题目

POJ 3264 Balanced Lineup 求区间的最大值-最小值

[cpp] view plain copy
print?
  1. #include <cstdio>  
  2. #include <cmath>  
  3. #include <iostream>  
  4. using namespace std;  
  5. #define FOR( i , a , b ) for ( int i = a ; i <= b ; ++ i )  
  6. #define RE freopen("1.in","r",stdin);  
  7. #define bug(x) cout<<#x<<":"<<(x)<<endl;  
  8.   
  9. #define Mid ((l+r)>>1)      //括号！  
  10. #define lson rt<<1,l,Mid  
  11. #define rson rt<<1|1,Mid+1,r  
  12.   
  13. const int maxn= 50005;  
  14. const int inf=0x3f3f3f3f;  
  15. int Max[maxn<<2],Min[maxn<<2];  
  16. int ans1,ans2;  
  17.   
  18. void build(int rt,int l,int r)  
  19. {  
  20.     if(l==r){  
  21.         scanf("%d",&Max[rt]);  
  22.         Min[rt] = Max[rt];  
  23.     }else{  
  24.         build(lson);  
  25.         build(rson);  
  26.         Max[rt] = max( Max[rt<<1], Max[rt<<1|1]);  
  27.         Min[rt] = min( Min[rt<<1], Min[rt<<1|1]);  
  28.     }  
  29. }  
  30.   
  31. void query(int rt,int l,int r,int L,int R)  
  32. {  
  33.     if(L <= l && r <= R){  
  34.         ans1 = max(ans1,Max[rt]);  
  35.         ans2 = min(ans2,Min[rt]);  
  36.     }else{  
  37.         if( L <= Mid)  
  38.             query(lson,L,R);  
  39.         if( R > Mid)  
  40.             query(rson,L,R);  
  41.     }  
  42. }  
  43.   
  44. int main()  
  45. {  
  46.     int n,m,L,R;  
  47.     while(scanf("%d%d",&n,&m)!=EOF){  
  48.         build(1,1,n);  
  49.         while(m--){  
  50.             ans1 = -inf,ans2 = inf;  
  51.             scanf("%d%d",&L,&R);  
  52.             query(1,1,n,L,R);  
  53.             printf("%d\n", ans1-ans2);  
  54.         }  
  55.     }  
  56.     return 0;  
  57. }  


HDU 1754  I Hate It单点更新，区间查询最大

[cpp] view plain copy
print?
  1. #include <cstdio>  
  2. #include <cmath>  
  3. #include <iostream>  
  4. using namespace std;  
  5. #define FOR( i , a , b ) for ( int i = a ; i <= b ; ++ i )  
  6. #define RE freopen("1.in","r",stdin);  
  7. #define bug(x) cout<<#x<<":"<<(x)<<endl;  
  8.   
  9. #define Mid ((l+r)>>1)       
  10. #define lson rt<<1,l,Mid  
  11. #define rson rt<<1|1,Mid+1,r  
  12.   
  13. const int maxn= 200000+5;  
  14. const int inf=0x3f3f3f3f;  
  15. int Max[maxn<<2];  
  16. int ans1,ans2;  
  17.   
  18. void build(int rt,int l,int r)  
  19. {  
  20.     if(l==r){  
  21.         scanf("%d",&Max[rt]);  
  22.     }else{  
  23.         build(lson);  
  24.         build(rson);  
  25.         Max[rt] = max( Max[rt<<1], Max[rt<<1|1]);  
  26.     }  
  27. }  
  28. void update(int rt,int l,int r,int pos,int num)  
  29. {  
  30.     if(l == r && r == pos){  
  31.         Max[rt] = num;  
  32.     }else{  
  33.         if( pos <= Mid)  
  34.             update(lson,pos,num);  
  35.         if( pos > Mid)  
  36.             update(rson,pos,num);  
  37.         Max[rt] = max( Max[rt<<1], Max[rt<<1|1]);  
  38.     }  
  39. }  
  40. int query(int rt,int l,int r,int L,int R)  
  41. {  
  42.     if(L <= l && r <= R){  
  43.         return Max[rt];  
  44.     }else{  
  45.         int tmp = -1;  
  46.         if( L <= Mid)  
  47.             tmp = max(tmp,query(lson,L,R));  
  48.         if( R > Mid)  
  49.             tmp = max(tmp,query(rson,L,R));  
  50.         return tmp;  
  51.     }  
  52. }  
  53.   
  54. int main()  
  55. {  
  56. //    RE  
  57.     int n,m,L,R;  
  58.     char op;  
  59.     while(scanf("%d%d",&n,&m)!=EOF){  
  60.         build(1,1,n);  
  61.         getchar();  
  62.         while(m--){  
  63.             scanf("%c%d%d%*c",&op,&L,&R);  
  64.             if(op=='Q')  
  65.                 printf("%d\n",query(1,1,n,L,R));  
  66.             else  
  67.                 update(1,1,n,L,R);  
  68.         }  
  69.     }  
  70.     return 0;  
  71. }  


HDU 1166 敌兵布阵

单点更新，查询区间和【此题树状数组解法转见→】

[cpp] view plain copy
print?
  1. #include <cstdio>  
  2. #include <cmath>  
  3. #include <iostream>  
  4. using namespace std;  
  5. #define FOR( i , a , b ) for ( int i = a ; i <= b ; ++ i )  
  6. #define RE freopen("1.in","r",stdin);  
  7. #define bug(x) cout<<#x<<":"<<(x)<<endl;  
  8.   
  9. #define Mid ((l+r)>>1)  
  10. #define lson rt<<1,l,Mid  
  11. #define rson rt<<1|1,Mid+1,r  
  12.   
  13. const int maxn= 50000+5;  
  14. const int inf=0x3f3f3f3f;  
  15. int sum[maxn<<2];  
  16. int ans1,ans2;  
  17.   
  18. void build(int rt,int l,int r)  
  19. {  
  20.     if(l==r){  
  21.         scanf("%d",&sum[rt]);  
  22.     }else{  
  23.         build(lson);  
  24.         build(rson);  
  25.         sum[rt] = sum[rt<<1] + sum[rt<<1|1];  
  26.     }  
  27. }  
  28. void update(int rt,int l,int r,int pos,int num)  
  29. {  
  30.     if(l == r && r == pos){       
  31.         sum[rt] += num;  
  32.     }else{  
  33.         if( pos <= Mid)  
  34.             update(lson,pos,num);  
  35.         else          
  36.             update(rson,pos,num);  
  37.         sum[rt] = sum[rt<<1] + sum[rt<<1|1];  
  38.     }  
  39. }  
  40.   
  41. int query(int rt,int l,int r,int L,int R)  
  42. {  
  43.     if(L <= l && r <= R){  
  44.         return sum[rt];  
  45.     }else{  
  46.         int tmp = 0;  
  47.         if( L <= Mid)  
  48.             tmp += query(lson,L,R);  
  49.         if( R > Mid)  
  50.             tmp += query(rson,L,R);  
  51.         return tmp;  
  52.     }  
  53. }  
  54.   
  55. int main()  
  56. {  
  57. //    RE  
  58.     int n,m,L,R,t;  
  59.     char op[10];  
  60.     scanf("%d",&t);  
  61.     for(int cas=1;cas<=t;cas++){  
  62.         scanf("%d",&n);  
  63.         build(1,1,n);  
  64.         printf("Case %d:\n",cas);  
  65.         while(scanf("%s",op),op[0]!='E'){  
  66.             scanf("%d%d",&L,&R);  
  67.             if(op[0]=='Q')  
  68.                 printf("%d\n", query(1,1,n,L,R));  
  69.             else if(op[0]=='A')  
  70.                 update(1,1,n,L,R);  
  71.             else  
  72.                 update(1,1,n,L,-R);  
  73.         }  
  74.     }  
  75.     return 0;  
  76. }  




Lazy


区间成段更新

Lazy：正常来说，区间改值，当更改某个区间的值的时候，子区间也该跟着更改，这样容易TLE。

Lazy思想就是更新到某个区间的时候，就先给这个区间打上标记，标记内容是需要更新的值，并把子区间的值改为子区间对应的值，清除该区间的lazy标记；然后return，不去更新子区间。当下一次更新或查询等需要访问该区间的子区间的时候再把该区间的lazy和其他信息送回子区间。


举个简单粗暴的例子：

对应下面的那个图，假如目的是求和，现在要给[1,6] 的值都加2，那么我们从[1,12]->[1,6]，然后[1,6]的sum值加上区间长度[ (6-1+1)*2 ]，再把[1,6]的add[i]设置为2，就不再往下更新了【这里极大提高效率】。下一次更新/查询[1,6]的子区间时，我们将[1,6]原存的add值下传给[1,6]的两个直接子区间，再往下更新。假设在这种情况下，我们再更新[1,6]加3，则[1,6]的add值为2+3=5，然后我们查询[1,3]，则从上往下经过[1,6]时把[1,6]的add值给了子区间[1,3]和[4,6]，同时把sum[子区间]跟着子区间长度和add[父结点]改动，清除add[父节点]。【如果是查询间接子区间，则连续传递add值，也就是连续pushDown】


详细例子：假设update()是区间改值，query()是求和，所有叶子区间的和都为1，则[7,8]和[7,9]在build()的时候就附上了值（图中绿色字体）。假设此时我们更新[7,9]的值，改为2，则线段树从[1,12]->[7,12]->[7,9]，然后把[7,9]打上值为2的标记，求和（求和直接用区间长度*此时更新的值），然后不去更新[7,8]和[9,9]了，他们值仍然是2和1，lazy值为0。


然后我们查询[7,8]，当遍历经过[7,9]时

[cpp] view plain copy
print?
  1. if(add[i])  
  2.        pushDown(i);  

成立，把[7,9]的lazy标记2传给子区间[7,8]和[9,9],分别求这2个子区间的和，把[7,9]的lazy标记去掉，然后继续遍历，到[7,8]的时候直接返回答案。



HDU 1698 Just a Hook 区间改值，求和

[cpp] view plain copy
print?
  1. #include <cstdio>  
  2. #include <cmath>  
  3. #include <iostream>  
  4. #include <algorithm>  
  5. using namespace std;  
  6. #define ll long long  
  7. #define FOR( i , a , b ) for ( int i = a ; i <= b ; ++ i )  
  8. #define RE freopen("1.in","r",stdin);  
  9. #define WE freopen("1.out","w",stdout);  
  10. #define MOD 10009  
  11. #define bug(x) cout<<#x<<":"<<(x)<<endl;  
  12.   
  13. #define Mid ((l+r)>>1)  
  14. #define lson rt<<1,l,Mid  
  15. #define rson rt<<1|1,Mid+1,r  
  16. const int maxn = 100010;  
  17. int sum[maxn<<2],add[maxn<<2];  
  18.   
  19. void build(int rt,int l,int r)  
  20. {  
  21.     add[rt] = 0;  
  22.     if(l == r){  
  23.         sum[rt] = 1;  
  24.     }else{  
  25.         build(lson);  
  26.         build(rson);  
  27.         sum[rt]=sum[rt<<1]+sum[rt<<1|1];  
  28.     }  
  29. }  
  30.   
  31. void pushDown(int rt,int len)  
  32. {  
  33.     add[rt<<1] = add[rt<<1|1] = add[rt];  
  34.     sum[rt<<1] = (len-(len>>1))*add[rt];  
  35.     sum[rt<<1|1] = (len>>1)*add[rt];  
  36.     add[rt] = 0;  
  37. }  
  38.   
  39. void update(int rt,int l,int r,int L,int R,int z)  
  40. {  
  41.     if(L <= l && r <= R){  
  42.         add[rt] = z;  
  43.         sum[rt] = (r-l+1)*z;  
  44.     }else{  
  45.         if(add[rt])  
  46.             pushDown(rt,r-l+1);  
  47.         if(L <= Mid)  
  48.             update(lson,L,R,z);  
  49.         if(R > Mid)  
  50.             update(rson,L,R,z);  
  51.         sum[rt] = sum[rt<<1] + sum[rt<<1|1];  
  52.     }  
  53. }  
  54.   
  55. int main()  
  56. {  
  57. //    RE;  
  58.     int t,n,q,x,y,z;  
  59.     int cnt=1;  
  60.     scanf("%d",&t);  
  61.     while(t--){  
  62.         scanf("%d%d",&n,&q);  
  63.         build(1,1,n);  
  64.         while(q--){  
  65.             scanf("%d%d%d",&x,&y,&z);  
  66.             update(1,1,n,x,y,z);  
  67.         }  
  68.         printf("Case %d: The total value of the hook is %d.\n", cnt++,sum[1]);  
  69.     }  
  70.     return 0;  
  71. }  

相关链接：线段树题目合集                  二维线段树                       线段树模板


阅读全文
1 0
苹果下软件赚钱的软件 苹果下软件赚钱的软件
原创粉丝点击
热门IT博客
按键精灵时间限制源码按键精灵时间限制源码
淘宝店铺二维码生成器淘宝店铺二维码生成器
php ip
词典 知乎
网络电话短信软件
尺八制作数据
环球网 知乎
阿里云申请域名
java程序员考试
知豆山东生产基地
mac内存不足怎么清理
cms内容发布系统
在数据库中的文件类型
淘宝联盟推广赚钱吗
大石家具城怎么样 知乎
淘宝店铺定时发布
redis php使用
宣传片制作软件
淘宝1cd出坐骑
微信服务号 域名 端口
热门问题 老师的惩罚 人脸识别 我在镇武司摸鱼那些年 重生之率土为王 我在大康的咸鱼生活 盘龙之生命进化 天生仙种 凡人之先天五行 春回大明朝 姑娘不必设防，我是瞎子 红萝卜炒鸡蛋有毒吗 为什么说土豆炒鸡蛋有毒 食用鹅膏菌 鹅膏菌 可食用鹅膏菌图片 毒黄瓜 灭鼠药 驱鼠药 刀塔致命毒液 金钻有毒吗 金钻有毒为什么还这么多人养 雨花毓秀 毓秀园 毓秀园二手房 毓秀园小区 骆老师 三比一提升 第一章爷爷比爸爸厉害 一张课桌比一把椅子贵10元 为什么机票改签比买一张还贵 等比中项 等比中项公式 生活中的比 我比想象中更爱你 等比中项的公式 生活中的比有哪些 生活中哪里还有黄金分割比 生活中的黄金分割比 生活中有哪些比 等比中项的公式是什么 空气中氧气和氮气的体积比 生活中的比教学设计 比丘国剧情 比丘尼 比康股份股吧 比亚乔 比亚乔摩托车 比亚乔部落 比亚乔150 比亚乔medley 比亚乔踏板

程序博客网，程序员的互联网技术博客家园。csdn论坛精品 msdn技术资料都在这里