ST表这一数据结构以前没有涉及过,今天在集训的时候听师兄讲了大有启发,深感其是个美妙的一种算法。
参考资料:ST 表 - OI Wiki
以及SCUT22级 吕浩天学长集训的 数据结构进阶 PPT
ST表用于解决可重复贡献问题,并且要满足结合律
可重复贡献问题 是指对于运算 ,满足
,则对应的区间询问就是一个可重复贡献问题。例如,最大值有
,gcd 有
,所以 RMQ 和区间 GCD 就是一个可重复贡献问题。像区间和就不具有这个性质,如果求区间和的时候采用的预处理区间重叠了,则会导致重叠部分被计算两次,这是我们所不愿意看到的。另外,
还必须满足结合律才能使用 ST 表求解。
ST表主要利用了 倍增 的思想,预处理时间 单次询问时间
,但不支持修改,下面以RMQ问题(区间最值问题)为例,解释ST表如何实现
构造递推关系(类似dp过程) i表示下标,j表示从
开始到
区间的最值
我们可以很容易得到关系:
这样我们便可以预处理完毕
如何利用这一关系求区间最值呢?
首先我们知道,ST表需要满足可重复贡献性质,而最大值时可重复贡献的,也就是说,两个相交区间的重合不影响答案
对于询问区间我们只需要找到两个交叉的区间,使着两个区间的并集覆盖
即可
通过数学知识,我们可以知道,也就是说,如果我们令指数值
,则从
与
可以覆盖完全
体现在代码上,则是返回
完整代码:
// Problem: P3865 【模板】ST 表
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P3865
// Memory Limit: 125 MB
// Time Limit: 800000 ms
//
// Powered by CP Editor (https://cpeditor.org)
#include <fstream>
#include <iostream>
#include <vector>
using namespace std;
const int N = 1e5 + 10;
inline int fr() {
int x = 0, f = 1;
char ch = getchar();
while (ch < '0' || ch > '9') {
if (ch == '-')
f = -1;
ch = getchar();
}
while (ch >= '0' && ch <= '9') {
x = x * 10 + ch - 48;
ch = getchar();
}
return x * f;
}
vector<int> A(N);
vector<vector<int>> f(N, vector<int>(31));
void Build_ST_table(int n) {
for (int k = 0; k <= 30; k++)
for (int i = 1; i + (1 << k) - 1 <= n; i++)
if (!k)
f[i][k] = A[i];
else
f[i][k] = max(f[i][k - 1], f[i + (1 << k - 1)][k - 1]);
return;
}
int Query(int l, int r) {
int k = __lg(r - l + 1);
return max(f[l][k], f[r - (1 << k) + 1][k]);
}
int main() {
int n, m;
n = fr();
m = fr();
for (int i = 1; i <= n; i++)
A[i] = fr();
Build_ST_table(n);
int l, r;
for (int i = 1; i <= m; i++) {
l = fr();
r = fr();
printf("%d\n", Query(l, r));
}
return 0;
}