预备知识

  1. st表(Sparse Table)

    主要用来解决区间最值问题(RMQ)以及维护区间的各种性质（比如维护一段区间的最大公约数）。

  2. 树状数组

    • 单点更新
    • 数组前缀和的查询

    拓展：原数组是差分数组时，可进行区间更新，单点查询，当然想区间查询也有办法（维护两个树状数组即可）。

区别

树状数组用来维护一个具有区间可减（加）性质的工具，所以可以用来维护区间前缀和。

区间最值不具有区间可减性，所以不能使用树状数组进行维护，而使用st表。

st表的思想

  • 初始化

    预处理数组, f[x][i] 表示区间[x, x + (1 << i) ) 的最值

    f[x][i] = min/max{f[x][i - 1], f[x + (1 << i - 1)][i - 1]}; // 1 << i - 1表示2的i-1次方

    其实这里就是一个动态规划的思想 :

    从x开始2的i次方的最值等于从 x 开始 (1 << i - 1) 范围内的最值加上从x + (1 << i - 1) 开始 (1 << i - 1)的最值

    类似于2 = 1 + 1
  • 询问区间内最值

    对于询问区间[l, r], 设 d = [log2(r - l + 1)]

    答案即为 : min/max{f[l][d], f[r - (1 << d) + 1][d]}

    至于为什么，很显然从左边开始(1 << d)区间内，要是够不到，就再从右开始

树状数组

数组模拟树形结构, 没啥好说的

题单

ST表

(Poj 3264 Balanced Lineup)

对应知识点

  • st表维护最值问题
想不出来查看思路 
  1.   利用st表开两个数组维护最大值和最小值

代码

ヾ(•ω•`)o还是写不出来就点我吧 
  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <cmath>
  5. #include <algorithm>
  7. using namespace std;
  8. typedef long long LL;
  10. const int N = 5e4 + 10;
  11. int n, q; // n : 数的数量 q : 询问数量
  12. int dp_max[N][20]; // 存储区间最大值
  13. int dp_min[N][20]; // 存储区间最小值
  15. // st表初始化
  16. void st_init()
  17. {
  18.     for(int j = 1; j <= 21; j ++)
  19.         for(int i  = 1; i + (1 << j) - 1 <= n; i ++) {
  20.             dp_max[i][j] = max(dp_max[i][j - 1], dp_max[i + (1 << (j - 1))][j - 1]);
  21.             dp_min[i][j] = min(dp_min[i][j - 1], dp_min[i + (1 << (j - 1))][j - 1]);
  22.         }
  23. }
  25. // st表询问区间内最大值及其最小值,得到它们的差
  26. int query(int l, int r)
  27. {
  28.     // int k = log2(r - l + 1);
  29.     int k=(int)(log(double(r-l+1)) / log(2.0));
  30.     int x = max(dp_max[l][k], dp_max[r - (1 << k) + 1][k]);
  31.     int y = min(dp_min[l][k], dp_min[r - (1 << k) + 1][k]);
  32.     return x - y;
  33. }
  35. // 快读模板，过题用
  36. inline int read()
  37. {
  38. 	int x=0,f=1;char ch=getchar();
  39. 	while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
  40. 	while (ch>='0'&&ch<='9'){x=x*10+ch-48;ch=getchar();}
  41. 	return x*f;
  42. }
  44. int main()
  45. {
  46.     // ios::sync_with_stdio(false);
  47.     // cin >> n >> q;
  48.     n = read(), q = read();
  49.     for(int i = 1; i <= n; i ++)
  50.     {
  51.         int x;
  52.         x = read();
  54.         dp_max[i][0] = x;
  55.         dp_min[i][0] = x;
  56.     }
  57.     st_init();
  58.     for(int i = 1, l, r; i <= q; i ++)
  59.     {
  60.         l = read(), r = read();
  61.         printf("%d\n",query(l, r));
  62.     }
  64.     return 0;
  66. }

(CF 1547F Array Stabilization (GCD version))

对应知识点

  • st表维护区间内最大公约数 + 二分答案

题目大意

给定一个序列a，每次对每一对相邻的数求他们的最大公约数，并且保存为新的序列。Latex不会打，可能题意不清楚，可以移步去洛谷

解题核心

gcd(a,b,c)=gcd(gcd(a,b),c)

  • 断环成链
  • 二分查询区间最大长度
ヾ(•ω•`)o还是写不出来就点我吧 
  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <algorithm>
  5. #include <cmath>
  7. using namespace std;
  9. typedef long long LL;
  10. const int N = 400010;
  12. int T, a[N];
  13. int n;
  14. LL g;
  15. LL st[N][25];
  16. LL gcd(int a, int b)
  17. {
  18.     return b ? gcd(b, a % b) : a;
  19. }
  21. void st_init(){
  22.     for(int i = 1; i <= 2 * n; i ++) st[i][0] = a[i]; // 赋初值 + 断环成链
  24.     // st表维护区间最大公约数
  25.     for(int j = 1; (1 << j) <= 2 * n; j ++){
  26.         for(int i = 1; i + (1 << j)  <= 2 * n; i ++){
  27.             st[i][j] = gcd(st[i][j - 1], st[i + (1 << j - 1)][j - 1]);
  28.         }
  29.     }
  30. }
  32. LL query(int l, int r){
  33.     int k = (int)(log(((double)(r - l + 1)))/log(2.0));
  34.     return gcd(st[l][k], st[r - (1 << k) + 1][k]);
  35. }
  37. bool check(int x){
  38.     for(int i = 1; i <= n; i ++){
  39.         if(g != query(i, i + x)) return false;
  40.     }
  41.     return true;
  42. }
  43. int main()
  44. {
  45.     scanf("%d", &T);
  46.     while(T --)
  47.     {
  48.         scanf("%d", &n);
  49.         for(int i = 1; i <= n; i ++)
  50.         {
  51.             scanf("%d", &a[i]);
  52.             a[i + n] = a[i];
  53.             if(i == 1) g = a[i];
  54.             else g = gcd(g, a[i]);
  55.         }
  56.         st_init();
  57.         int l = 0, r = n - 1;
  58.         while(l < r){
  59.             int mid = (l + r) >> 1;
  60.             if(check(mid)) r = mid;
  61.             else l = mid + 1;
  62.         }
  63.         printf("%d\n", l);
  65.     }
  67.     return 0;
  68. }

(CF1549D Integers Have Friends)

知识点同上题

这题需要用差分消去余数，剩下还是用st表维护最大公约数

ヾ(•ω•`)o还是写不出来就点我吧 
  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <algorithm>
  5. #include <cmath>
  7. using namespace std;
  8. typedef long long LL;
  9. const LL N = 200010;
  11. LL a[N], b[N];
  12. LL st[N][25];
  13. LL n, T;
  14. LL gcd(LL a, LL b){
  15.     return b ? gcd(b, a % b) : a;
  16. }
  18. LL query(LL l, LL r){
  19.     LL k = (LL)(log((double)(r - l  + 1)) / log(2.0));
  20.     return gcd(st[l][k], st[r - (1 << k) + 1][k]);
  21. }
  22. int main()
  23. {
  24.     scanf("%d", &T);
  25.     while(T --)
  26.     {
  27.         scanf("%d", &n);
  28.         for(LL i = 1; i <= n; i ++) scanf("%lld", &a[i]);
  29.         if(n == 1)
  30.         {
  31.             printf("1\n");
  32.             continue;
  33.         }
  34.         for(LL i = 1; i < n; i ++){
  35.             b[i] = a[i] - a[i + 1];
  36.             if(b[i] < 0) b[i] = - b[i];
  37.         }
  38.         n--;
  39.         for(LL i = 1; i <= n; i ++) st[i][0] = b[i];
  40.         for(LL j = 1; (1 << j) <= n; j ++)
  41.             for(LL i = 1; i + (1 << j) - 1 <= n; i ++)
  42.             st[i][j] = gcd(st[i][j - 1], st[i + (1 << j - 1)][j - 1]);
  43.         LL ans = 0;
  44.         for(LL i = 1, l, r; i <= n; i ++){
  45.             l = i, r = n;
  46.             while(l < r){
  47.                 LL mid = l + r + 1 >> 1;
  48.                 if(query(i, mid) == 1) r = mid - 1;
  49.                 else l = mid;
  50.             }
  51.             if(b[i] != 1) ans = max(ans, l - i + 1);
  52.         }
  54.         printf("%lld\n", ans + 1);
  55.     }
  57.     return 0;
  58. }

树状数组

(HDU 1166 敌兵布阵)

对应知识点

  • 树状数组裸题（模板的使用）
ヾ(•ω•`)o还是写不出来就点我吧 
  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <cmath>
  5. #include <algorithm>
  7. using namespace std;
  9. const int N = 50010;
  10. int T, n;
  11. char s[8];
  12. int l, r;
  13. int tr[N], a[N];
  14. inline int read()
  15. {
  16. 	int x=0,f=1;char ch=getchar();
  17. 	while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
  18. 	while (ch>='0'&&ch<='9'){x=x*10+ch-48;ch=getchar();}
  19. 	return x*f;
  20. }
  22. int lowbit(int x){
  23.     return x & -x;
  24. }
  25. void add(int k, int c){
  26.     while(k <= n) tr[k] += c, k += lowbit(k);
  27. }
  29. int ask(int x){
  30.     int res = 0;
  31.     while(x) res += tr[x], x -= lowbit(x);
  32.     return res;
  33. }
  34. int main()
  35. {
  36.     T = read();
  37.     int cnt = 1;
  38.     while(T -- )
  39.     {
  40.         n = read();
  41.         printf("Case %d:\n", cnt ++);
  42.         for(int i = 1; i <= n; i ++){
  43.             a[i] = read();
  44.             add(i, a[i]);
  45.         }
  47.         while(1){
  48.             scanf("%s", s);
  49.             if(s[0] == 'Q'){
  50.                 l = read(), r = read();
  51.                 printf("%d\n", ask(r) - ask(l - 1));
  52.             }else if(s[0] == 'A'){
  53.                 l = read(), r = read();
  54.                 add(l, r);
  55.             }else if(s[0] == 'S'){
  56.                 l = read(), r = read();
  57.                 add(l, -r);
  58.             }else if(s[0] == 'E'){
  59.                 break;
  60.             }
  61.         }
  62.         memset(tr, 0, sizeof(tr));
  63.         memset(a, 0, sizeof(a));
  64.     }
  66.     return 0;
  67. }

(HDU 1754 I Hate It)

对应知识点

  • 树状数组维护最大值(本来准备用st表的，但是有修改操作)

    具体做法就是把树状数组维护的区间和变成区间最大值，一开始实在没有想到这么写

    如果看不懂可以把这题的操作和树状数组最基本的两个操作代码对比起来看
ヾ(•ω•`)o还是写不出来就点我吧 
  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cmath>
  4. #include <algorithm>
  6. using namespace std;
  8. const int N = 2e6 + 10;
  10. int lowbit(int x){
  11.     return x & -x;
  12. }
  13. inline int read()
  14. {
  15. 	int x=0,f=1;char ch=getchar();
  16. 	while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
  17. 	while (ch>='0'&&ch<='9'){x=x*10+ch-48;ch=getchar();}
  18. 	return x*f;
  19. }
  20. int h[N], a[N];
  21. int n, m;
  22. // 树状数组维护最大值 时间复杂度O(logn^2)
  23. void update(int x)
  24. {
  25.     while(x <= n){
  26.         h[x] = a[x];
  27.         for(int i = 1; i < lowbit(x);i <<= 1)
  28.         h[x] = max(h[x], h[x - i]);
  29.         x += lowbit(x);
  30.     }
  31. }
  33. int query(int x, int y)
  34. {
  35. 	int ans = 0;
  37.     // 如果y-lowbit(y) >= x即直接取已经维护的h[y]最大值
  38.     // 否则逐个逼近最大值
  39. 	while (y >= x)
  40. 	{
  41. 		ans = max(a[y], ans);
  42. 		y --; // 已经取了第一个值
  43. 		for (; y-lowbit(y) >= x; y -= lowbit(y))
  44. 			ans = max(h[y], ans);
  45. 	}
  46. 	return ans;
  47. }
  49. char ch[3];
  50. int main()
  51. {
  52.     while(~scanf("%d%d", &n, &m)){
  53.         for(int i = 1; i <= n; i ++) a[i] = read(), update(i);
  54.         for(int i = 1, l, r; i <= m; i ++){
  55.             scanf("%s", ch);
  56.             l = read(), r = read();
  57.             if(ch[0] == 'Q')
  58.             {
  59.                 printf("%d\n", query(l, r));
  60.             }else{
  61.                 a[l] = r;
  62.                 update(l);
  63.             }
  64.         }
  65.     }
  66.     return 0;
  67. }

(Poj 3468 A Simple Problem with Integers)

对应知识点

  • 树状数组的应用

    这题较模板而言的不同之处是要进行区间修改以及区间和

    对于区间修改我们简单的用树状数组维护一个差分数组

    对于区间和，我们可以维护两个树状数组，b[i] 和 i * b[i] 利用这两个数组去取得区间和

    (x + 1) * sum(Tr1, x) - sum(tr2, x)
ヾ(•ω•`)o还是写不出来就点我吧 
  1. #include <cstdio>
  2. #include <cstring>
  3. #include <iostream>
  4. #include <algorithm>
  5. #include <cmath>
  7. using namespace std;
  8. const int N = 100010;
  9. typedef long long LL;
  11. int a[N];
  12. LL Tr1[N], tr2[N];
  13. int n, m;
  14. int lowbit(int x){
  15.     return x & -x;
  16. }
  17. void add(LL tr[], int x, LL c){
  18.     while(x <= n) tr[x] += c, x += lowbit(x);
  19. }
  21. LL sum(LL tr[], int x){
  22.     LL res = 0;
  23.     while(x) res += tr[x], x-=lowbit(x);
  24.     return res;
  25. }
  27. LL prefix_sum(int x){
  28.     return (x + 1) * sum(Tr1, x) - sum(tr2, x);
  29. }
  30. int main()
  31. {
  32.     scanf("%d%d", &n, &m);
  33.     for(int i = 1; i <= n; i ++) scanf("%d", &a[i]);
  34.     for(int i = 1; i <= n; i ++){
  35.         int b = a[i] - a[i - 1];
  36.         add(Tr1, i, b);
  37.         add(tr2, i, (LL)i * b);
  38.     }
  40.     while(m --){
  41.         char op[2];
  42.         int l, r, d;
  43.         scanf("%s%d%d", op, &l, &r);
  44.         if(*op == 'Q'){
  45.             printf("%lld\n", prefix_sum(r) - prefix_sum(l - 1));
  46.         }else{
  47.             scanf("%d", &d);
  48.             add(Tr1, l, d),add(tr2, l, (LL)l * d);
  49.             add(Tr1, r + 1, -d),add(tr2, r + 1, (LL)(r + 1)*(-d));
  50.         }
  51.     }
  53.     return 0;
  54. }

(Poj 2182 Lost cows)

对应知识点

  • 树状数组 + 二分

思路

tr数组存每个序号是否还有牛，故初始化为1即可。从后往前，牛的序号为前面比他小的再加上1

二分得到答案，二分的是比他小的牛的数量和

ヾ(•ω•`)o还是写不出来就点我吧 
  1. #include <iostream>
  2. #define lowbit(x) x & (-x)
  4. using namespace std;
  6. const int N = 1e5;
  7. int tr[N], pre[N], ans[N];
  8. int n;
  10. void add(int x, int k)
  11. {
  12.     for(int i = x; i <= n; i += lowbit(i)) tr[i] += k;
  13. }
  15. int ask(int x)
  16. {
  17.     int sum = 0;
  18.     for(int i = x; i; i -= lowbit(i)) sum += tr[i];
  19.     return sum;
  20. }
  22. int findx(int x)
  23. {
  24.     int l = 1, r = n;
  25.     while(l < r)
  26.     {
  27.         int mid = l + r >> 1;
  28.         if(ask(mid) < x)
  29.             l = mid + 1;
  30.         else
  31.             r = mid;
  32.     }
  34.     return l;
  35. }
  36. int main()
  37. {
  38.     cin >> n;
  39.     pre[1] = 0;
  40.     for(int i = 2; i <= n; i ++) cin >> pre[i];
  42.     for(int i = 1; i <= n; i ++){
  43.         tr[i] = lowbit(i);
  44.         // cout << tr[i] << ' ';
  45.     }
  47.     for(int i = n; i > 0; i --){
  48.         int x = findx(pre[i] + 1);
  49.         add(x, - 1);
  50.         ans[i] = x;
  51.     }
  53.     for(int i = 1; i <= n; i ++) cout << ans[i] << endl;
  54.     return 0;
  55. }

(Poj 2299 Ultra-QuickSort)

过题

直接用归并排序

树状数组的后面学一手再补上来

ヾ(•ω•`)o还是写不出来就点我吧 
  1. #include <cstdio>
  3. typedef long long LL;
  5. const int N = 500010;
  7. int n;
  8. LL q[N], w[N];
  10. LL merge_sort(int l, int r)
  11. {
  12.     if (l == r) return 0;
  14.     int mid = l + r >> 1;
  15.     LL res = merge_sort(l, mid) + merge_sort(mid + 1, r);
  16.     int i = l, j = mid + 1, k = 0;
  17.     while (i <= mid && j <= r)
  18.         if (q[i] <= q[j]) w[k ++ ] = q[i ++ ];
  19.         else
  20.         {
  21.             res += mid - i + 1;
  22.             w[k ++ ] = q[j ++ ];
  23.         }
  24.     while (i <= mid) w[k ++ ] = q[i ++ ];
  25.     while (j <= r) w[k ++ ] = q[j ++ ];
  27.     for (i = l, j = 0; i <= r; i ++, j ++ ) q[i] = w[j];
  29.     return res;
  30. }
  32. int main()
  33. {
  34.     while (scanf("%d", &n), n)
  35.     {
  36.         for (int i = 0; i < n; i ++ ) scanf("%d", &q[i]);
  38.         printf("%lld\n", merge_sort(0, n - 1));
  39.     }
  41.     return 0;
  42. }

写到最后

第一次写博客，如果有不足请提出，我很乐意听取意见。

现在处于一个大量知识的学习时间段，希望博客能够让我所学积淀下来。

感谢观看(￣︶￣*)) 完结撒花

st表树状数组入门题单的更多相关文章

  1. POJ 2299 Ultra-QuickSort 求逆序数 （归并或者数状数组）此题为树状数组入门题！！！

    Ultra-QuickSort Time Limit: 7000MS   Memory Limit: 65536K Total Submissions: 70674   Accepted: 26538 ...

  2. HDU 1166 敌兵布阵（线段树/树状数组模板题）

    敌兵布阵 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Submi ...

  3. bzoj1103树状数组水题

    (卧槽,居然规定了修改的两点直接相连,亏我想半天) 非常水的题,用dfs序(而且不用重复,应该是直接规模为n的dfs序)+树状数组可以轻松水 收获:树状数组一遍A(没啥好骄傲的,那么简单的东西) #i ...

  4. UESTC 1584 Washi与Sonochi的约定【树状数组裸题+排序】

    题目链接:UESTC 1584 Washi与Sonochi的约定 题意:在二维平面上,某个点的ranked被定义为x坐标不大于其x坐标,且y坐标不大于其y坐标的怪物的数量.(不含其自身),要求输出n行 ...

  5. 敌兵布阵 HDU - 1166 （树状数组模板题，线段树模板题）

    思路:就是树状数组的模板题,利用的就是单点更新和区间求和是树状数组的强项时间复杂度为m*log(n) 没想到自己以前把这道题当线段树的单点更新刷了. 树状数组: #include<iostrea ...

  6. poj2299树状数组入门，求逆序对

    今天入门了树状数组 习题链接 https://blog.csdn.net/liuqiyao_01/article/details/26963913 离散化数据:用一个数组来记录每个值在数列中的排名,不 ...

  7. 树状数组训练题1：弱弱的战壕（vijos1066）

    题目链接:弱弱的战壕 这道题似乎是vijos上能找到的最简单的树状数组题了. 原来,我有一个错误的思想,我的设计是维护两个树状数组,一个是横坐标,一个是纵坐标,然后读入每个点的坐标,扔进对应的树状数组 ...

  8. POJ 2352 stars （树状数组入门经典！！！）

    Stars Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 54352   Accepted: 23386 Descripti ...

  9. 树状数组 简单题 cf 961E

    题目链接 : https://codeforces.com/problemset/problem/961/E One day Polycarp decided to rewatch his absol ...

随机推荐

  1. IOS 集成 Bilibili IJKPlayer播放器，播放rtmp视频流

    因为公司项目需要,我一个连iPhone都没用过的人竟然跑去开发iOS APP.近一段时间一直忙于赶项目,到今天差不多了,所以记录一下当时遇到的各种坑,先从ios 集成 ijkplayer播放器说起! ...

  2. 从需求去理解 Linux dbus与基于dbus协议的无agent软件管理

    What is IPC IPC [Inter-Process Communication] 进程间通信,指至少两个进程或线程间传送数据或信号的一些技术或方法.在Linux/Unix中,提供了许多IPC ...

  3. 事务保存点savepoint

    一.

  4. Input 只能输入数字，数字和字母等的正则表达式

    JS只能输入数字,数字和字母等的正则表达式 1.文本框只能输入数字代码(小数点也不能输入) <input onkeyup="this.value=this.value.replace( ...

  5. IDE集成管理Tomcat的基本原理

    知道IDE是怎样控制Tomcat的,对更清晰地理解Java Web的执行过程有帮助.在此以IntelliJ IDEA为例,简要描述一下IDE集成管理Tomcat的基本原理. 首先是两个重要的环境变量: ...

  6. 密码学系列之:bcrypt加密算法详解

    目录 简介 bcrypt的工作原理 bcrypt算法实现 bcrypt hash的结构 hash的历史 简介 今天要给大家介绍的一种加密算法叫做bcrypt, bcrypt是由Niels Provos ...

  7. 关联数组VS索引数组

    关联数组和常规说的数组类似,它包含标量抄数据,可用索引值来单独选择这些数据,和常规数组不同的是, 关联数组的索引值不是非负的整数而是任意的标量袭.这些标量称为百Keys,可以在以后用于检索数组中的数值 ...

  8. php curl发送数据和文件

    function mycurl($file, $url, $aid) { // 如果文件名是中文名,将中文字符编码转换一下 $file=iconv("UTF-8","gb ...

  9. lumen-phpunit 单元测试

    lumen-框架5.8为例 1,把vendor下的bin目录放到环境变量里面: 2,设置路由 $router->get('syn', ['uses' => 'syn\syn@diction ...

  10. git 要求密码的解决方法：【生成gitLab公钥】：以及如何配置GitLab中的SSH key

    参考链接: https://www.cnblogs.com/yjlch1016/p/9692840.html https://blog.csdn.net/u011925641/article/deta ...