作用:
st表,也称RMQ算法,其用处就是求出给定区间的最值问题。
复杂度:
预处理:O(logn)
查询:O(1)
算法简介: ST表主要运用了dp的思想来转移值,不同的是ST表运用了倍增的思想来保证其处理的复杂度,ST表每次增加2^j个区间长度。
假设我们使用二维数组rmq来存储ST表,那么该数组的第一维度表示的是该段区间起始的位置,第二维表示的是该区间单调长度。例如,数组rmq[i, j],其中i表示的就是该区间的起点,j表示该区间长度为2^j,即此时该区间范围是[i, i + 2^j - 1]。
在维护ST表的时候,我们通常会将一段区间划分为两份,然后分别求出这两个区间的最值,再取其中的最值。例如区间[i, j],我们会分别求出[i, i + 2^(k - 1) - 1]和区间[i + 2^(k - 1), j]的最值(其中k表示区间长度的2的k次幂),然后取两者的最值。ST表一般用二维数组存储,如rmq [i] [j],i表示该区间的起点是i,长度为2^j。当j = 0时,区间长度就为2的0次方,即为1,这时候我们要对数组进行初始化,即为rmq[i] [0]赋初始值,表示一个单个点的值。
维护操作完成后,我们需要进行查询操作,对于区间范围[l, r] (其中 l 表示区间左端点,r表示区间右端点),我们先要求出其区间长度len,一般len = r - l + 1。然后我们需要求出len对应的2的幂k,这里可以通过math自带的log()函数求出,如:k = log(len) / log(2);也可以事先线性处理一下前n个数对应的值,这样会比使用log()函数快上一些。查询时我们要将一段区间分为两个相等 小区间来查询,这时就有区间没办法被完整的分为两个相等大小的小区间,这里允许对一段区间进行重复查询,这不会影响我们查询的结果。
如:待查询的区间为[3,11],实际上查询的两个子区间是[3, 10]和[4 ,11],应为其区间长度为k = log(11 - 3 + 1)) = 3。
代码模板
//这里以求区间最大值为例
# include<iostream>
# include<algorithm>
# include<cmath>
using namespace std;
const int maxn = 2e5 + 5;
const int N = 18;
int n, m;
int arr[maxn];
int rmq[maxn][N];
int query(int l, int r){
int len = r - l + 1;
int k = log(len) / log(2);
return max(rmq[l][k], rmq[r - (1 << k) + 1][k]);//查询并返回两个子区间中的最值
}
void init(){
for(int j = 0; j < N; j++){
for(int i = 1; i + (1 << j) - 1 <= maxn; i++){
if(j == 0)
rmq[i][j] = arr[i];//初始化赋值
else
rmq[i][j] = max(rmq[i][j - 1], rmq[i + (1 << j - 1)][j - 1]);//每次使用两个子区间的最值来更新该区间的最值,j表示该区间的长度为2^j
}
}
}
int main(){
ios::sync_with_stdio(false);
cin.tie(0);cout.tie(0);
cin >> n;
for(int i = 1; i <= n; i++)
cin >> arr[i];
init();
cin >> m;
while(m--){
int l, r;
cin >> l >> r;
cout << query(l, r) << endl;
}
return 0;
}