Notonlysuccess 线段树

很早前写的那篇线段树专辑至今一直是本博客阅读点击量最大的一片文章,当时觉得挺自豪的,还去pku打广告,但是现在我自己都不太好意思去看那篇文章了,觉得当时的代码风格实在是太丑了,很多线段树的初学者可能就是看着这篇文章来练习的,如果不小心被我培养出了这么糟糕的风格,实在是过意不去,正好过几天又要给集训队讲解线段树,所以决定把这些题目重新写一遍,顺便把近年我接触到的一些新题更新上去~;并且学习了splay等更高级的数据结构后对线段树的体会有更深了一层,线段树的写法也就比以前飘逸,简洁且方便多了.
在代码前先介绍一些我的线段树风格:

maxn是题目给的最大区间,而节点数要开4倍,确切的来说节点数要开大于maxn的最小2x的两倍
lson和rson分辨表示结点的左儿子和右儿子,由于每次传参数的时候都固定是这几个变量,所以可以用预定于比较方便的表示
以前的写法是另外开两个个数组记录每个结点所表示的区间,其实这个区间不必保存,一边算一边传下去就行,只需要写函数的时候多两个参数,结合lson和rson的预定义可以很方便
PushUP(int rt)是把当前结点的信息更新到父结点
PushDown(int rt)是把当前结点的信息更新给儿子结点
rt表示当前子树的根(root),也就是当前所在的结点
整理这些题目后我觉得线段树的题目整体上可以分成以下四个部分:

单点更新:最最基础的线段树,只更新叶子节点,然后把信息用PushUP(int r)这个函数更新上来
hdu1166 敌兵布阵 http://acm.hdu.edu.cn/showproblem.php?pid=1166
题意:O(-1)
思路:O(-1)
线段树功能:update:单点增减 query:区间求和

[cpp]  view plain copy print ?

1. #include <cstdio>  
2. #define lson l , m , rt << 1  
3. #define rson m + 1 , r , rt << 1 | 1  
4. const int maxn = 55555;  
5. int sum[maxn<<2];  
6. void PushUP(int rt) {  
7.         sum[rt] = sum[rt<<1] + sum[rt<<1|1];  
8. }  
9. void build(int l,int r,int rt) {  
10.         if (l == r) {  
11.                 scanf("%d",&sum[rt]);  
12.                 return ;  
13.         }  
14.         int m = (l + r) >> 1;  
15.         build(lson);  
16.         build(rson);  
17.         PushUP(rt);  
18. }  
19. void update(int p,int add,int l,int r,int rt) {  
20.         if (l == r) {  
21.                 sum[rt] += add;  
22.                 return ;  
23.         }  
24.         int m = (l + r) >> 1;  
25.         if (p <= m) update(p , add , lson);  
26.         else update(p , add , rson);  
27.         PushUP(rt);  
28. }  
29. int query(int L,int R,int l,int r,int rt) {  
30.         if (L <= l && r <= R) {  
31.                 return sum[rt];  
32.         }  
33.         int m = (l + r) >> 1;  
34.         int ret = 0;  
35.         if (L <= m) ret += query(L , R , lson);  
36.         if (R > m) ret += query(L , R , rson);  
37.         return ret;  
38. }  
39. int main() {  
40.         int T , n;  
41.         scanf("%d",&T);  
42.         for (int cas = 1 ; cas <= T ; cas ++) {  
43.                 printf("Case %d:\n",cas);  
44.                 scanf("%d",&n);  
45.                 build(1 , n , 1);  
46.                 char op[10];  
47.                 while (scanf("%s",op)) {  
48.                         if (op[0] == 'E') break;  
49.                         int a , b;  
50.                         scanf("%d%d",&a,&b);  
51.                         if (op[0] == 'Q') printf("%d\n",query(a , b , 1 , n , 1));  
52.                         else if (op[0] == 'S') update(a , -b , 1 , n , 1);  
53.                         else update(a , b , 1 , n , 1);  
54.                 }  
55.         }  
56.         return 0;  
57. }  

hdu1754 I Hate It
http://acm.hdu.edu.cn/showproblem.php?pid=1754
题意:O(-1)
思路:O(-1)
线段树功能:update:单点替换 query:区间最值

[cpp]  view plain copy print ?

1.   

[cpp]  view plain copy print ?

1. #include <cstdio>  
2. #include <algorithm>  
3. using namespace std;  
4.    
5. #define lson l , m , rt << 1  
6. #define rson m + 1 , r , rt << 1 | 1  
7. const int maxn = 222222;  
8. int MAX[maxn<<2];  
9. void PushUP(int rt) {  
10.         MAX[rt] = max(MAX[rt<<1] , MAX[rt<<1|1]);  
11. }  
12. void build(int l,int r,int rt) {  
13.         if (l == r) {  
14.                 scanf("%d",&MAX[rt]);  
15.                 return ;  
16.         }  
17.         int m = (l + r) >> 1;  
18.         build(lson);  
19.         build(rson);  
20.         PushUP(rt);  
21. }  
22. void update(int p,int sc,int l,int r,int rt) {  
23.         if (l == r) {  
24.                 MAX[rt] = sc;  
25.                 return ;  
26.         }  
27.         int m = (l + r) >> 1;  
28.         if (p <= m) update(p , sc , lson);  
29.         else update(p , sc , rson);  
30.         PushUP(rt);  
31. }  
32. int query(int L,int R,int l,int r,int rt) {  
33.         if (L <= l && r <= R) {  
34.                 return MAX[rt];  
35.         }  
36.         int m = (l + r) >> 1;  
37.         int ret = 0;  
38.         if (L <= m) ret = max(ret , query(L , R , lson));  
39.         if (R > m) ret = max(ret , query(L , R , rson));  
40.         return ret;  
41. }  
42. int main() {  
43.         int n , m;  
44.         while (~scanf("%d%d",&n,&m)) {  
45.                 build(1 , n , 1);  
46.                 while (m --) {  
47.                         char op[2];  
48.                         int a , b;  
49.                         scanf("%s%d%d",op,&a,&b);  
50.                         if (op[0] == 'Q') printf("%d\n",query(a , b , 1 , n , 1));  
51.                         else update(a , b , 1 , n , 1);  
52.                 }  
53.         }  
54.         return 0;  
55. }  

hdu1394 Minimum Inversion Number
http://acm.hdu.edu.cn/showproblem.php?pid=1394
题意:求Inversion后的最小逆序数
思路:用O(nlogn)复杂度求出最初逆序数后,就可以用O(1)的复杂度分别递推出其他解
线段树功能:update:单点增减 query:区间求和

[cpp]  view plain copy print ?

1. #include <cstdio>  
2. #include <algorithm>  
3. using namespace std;  
4.    
5. #define lson l , m , rt << 1  
6. #define rson m + 1 , r , rt << 1 | 1  
7. const int maxn = 5555;  
8. int sum[maxn<<2];  
9. void PushUP(int rt) {  
10.         sum[rt] = sum[rt<<1] + sum[rt<<1|1];  
11. }  
12. void build(int l,int r,int rt) {  
13.         sum[rt] = 0;  
14.         if (l == r) return ;  
15.         int m = (l + r) >> 1;  
16.         build(lson);  
17.         build(rson);  
18. }  
19. void update(int p,int l,int r,int rt) {  
20.         if (l == r) {  
21.                 sum[rt] ++;  
22.                 return ;  
23.         }  
24.         int m = (l + r) >> 1;  
25.         if (p <= m) update(p , lson);  
26.         else update(p , rson);  
27.         PushUP(rt);  
28. }  
29. int query(int L,int R,int l,int r,int rt) {  
30.         if (L <= l && r <= R) {  
31.                 return sum[rt];  
32.         }  
33.         int m = (l + r) >> 1;  
34.         int ret = 0;  
35.         if (L <= m) ret += query(L , R , lson);  
36.         if (R > m) ret += query(L , R , rson);  
37.         return ret;  
38. }  
39. int x[maxn];  
40. int main() {  
41.         int n;  
42.         while (~scanf("%d",&n)) {  
43.                 build(0 , n - 1 , 1);  
44.                 int sum = 0;  
45.                 for (int i = 0 ; i < n ; i ++) {  
46.                         scanf("%d",&x[i]);  
47.                         sum += query(x[i] , n - 1 , 0 , n - 1 , 1);  
48.                         update(x[i] , 0 , n - 1 , 1);  
49.                 }  
50.                 int ret = sum;  
51.                 for (int i = 0 ; i < n ; i ++) {  
52.                         sum += n - x[i] - x[i] - 1;  
53.                         ret = min(ret , sum);  
54.                 }  
55.                 printf("%d\n",ret);  
56.         }  
57.         return 0;  
58. }  

      成段更新
(通常这对初学者来说是一道坎),需要用到延迟标记(或者说懒惰标记),简单来说就是每次更新的时候不要更新到底,用延迟标记使得更新延迟到下次需要更新or询问到的时候
hdu1698 Just a Hook
http://acm.hdu.edu.cn/showproblem.php?pid=1698
题意:O(-1)
思路:O(-1)
线段树功能:update:成段替换 (由于只query一次总区间,所以可以直接输出1结点的信息

[cpp]  view plain copy print ?

1. #include <cstdio>  
2. #include <algorithm>  
3. using namespace std;  
4.    
5. #define lson l , m , rt << 1  
6. #define rson m + 1 , r , rt << 1 | 1  
7. const int maxn = 111111;  
8. int h , w , n;  
9. int col[maxn<<2];  
10. int sum[maxn<<2];  
11. void PushUp(int rt) {  
12.         sum[rt] = sum[rt<<1] + sum[rt<<1|1];  
13. }  
14. void PushDown(int rt,int m) {  
15.         if (col[rt]) {  
16.                 col[rt<<1] = col[rt<<1|1] = col[rt];  
17.                 sum[rt<<1] = (m - (m >> 1)) * col[rt];  
18.                 sum[rt<<1|1] = (m >> 1) * col[rt];  
19.                 col[rt] = 0;  
20.         }  
21. }  
22. void build(int l,int r,int rt) {  
23.         col[rt] = 0;  
24.         sum[rt] = 1;  
25.         if (l == r) return ;  
26.         int m = (l + r) >> 1;  
27.         build(lson);  
28.         build(rson);  
29.         PushUp(rt);  
30. }  
31. void update(int L,int R,int c,int l,int r,int rt) {  
32.         if (L <= l && r <= R) {  
33.                 col[rt] = c;  
34.                 sum[rt] = c * (r - l + 1);  
35.                 return ;  
36.         }  
37.         PushDown(rt , r - l + 1);  
38.         int m = (l + r) >> 1;  
39.         if (L <= m) update(L , R , c , lson);  
40.         if (R > m) update(L , R , c , rson);  
41.         PushUp(rt);  
42. }  
43. int main() {  
44.         int T , n , m;  
45.         scanf("%d",&T);  
46.         for (int cas = 1 ; cas <= T ; cas ++) {  
47.                 scanf("%d%d",&n,&m);  
48.                 build(1 , n , 1);  
49.                 while (m --) {  
50.                         int a , b , c;  
51.                         scanf("%d%d%d",&a,&b,&c);  
52.                         update(a , b , c , 1 , n , 1);  
53.                 }  
54.                 printf("Case %d: The total value of the hook is %d.\n",cas , sum[1]);  
55.         }  
56.         return 0;  
57. }