概述:
RMQ(Range Minimum/Maximum Query),即区间最值查询,在一个序列中多次查询区间(l,r)中最大值,然而解决RMQ需要对ST算法有着深刻的理解
ST算法的实现
ST算法是一种解决静态询问区间最值问题的高效在线算法,算法上采用了倍增的思想,而实现上则是利用了动态规划。 时间复杂的为O(Nlog(N))
利用一个二维数组,而 f[i][j] 就代表这 [i,i+2^j) 范围上的最值。
然后在把 [i,i+2^j) 拆分成 [i,i+2 ^ (j-1))和[i+2 ^ (j-1),i+2 ^ j) 然后就可以推出动态转移方程:
F[ l ][ r ]= max/min { F[ l ][ r-1 ] , F[ l+( 1<<(r-1) ) ][ r-1 ] }
f 这个二维数组预处理已经结束,下面要解决怎样利用它。
设 d 为 [l,r] 的长度,所以 d =
log
2
(
r
−
l
+
1
)
\log{2}{(r-l+1)}
log2(r−l+1) ,ans = max/min{ f[l][d] ,f[r-2^d+1][d] } 这样做为了防止log没有正好是正数。
对应的区间[ l , r] 分成 [l ,l + 2 ^ d) 和 [r-2 ^ d +1, r+1) 自己带入一下就可以看出来了。
如果明白这个算法就去洛谷刷一下这道题的模板 题目链接
下面给出我的代码:
#include <bits/stdc++.h>
using namespace std;
inline int read() {
int ans=0,flag=1;
char ch=getchar();
while( (ch>'9' || ch<'0') && ch!='-' ) ch=getchar();
if(ch=='-') flag=-1,ch=getchar();
while(ch>='0' && ch<='9') ans=ans*10+ch-'0',ch=getchar();
return ans*flag;
}
const int maxn = 1e5+10;
int f[maxn][18]; //18是log2(maxn)大一点的数
int n,m;
int Query(int l, int r){
int d = log2(r-l+1);
return max(f[l][d],f[r-(1 << d) + 1][d]);
}
int main(){
n = read(),m = read();
for(int i = 1; i <= n; i++)f[i][0] = read();
for(int j = 1; j <= 17; j++){
for(int i = 1; i + (1 << j) - 1 <= n; i++)
f[i][j] = max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
}
for(int i = 1, l, r; i <= m; i++){
l = read(),r = read();
printf("%d\n",Query(l,r));
}
return 0;
}
ps: 自我感觉,能通过字面意思学会这个算法是很难的,毕竟牵扯到动态规划这种抽象的东西,所以看视频就显得尤为重要了。给出这个算法视频链接我感觉这个人讲的比较好。