RMQ应用场景
区间最值查询,用来解决当查询区间范围很大且查询次数较多时暴力求解造成的超时问题
- 例如:求区间最大值、区间最小值
- 求区间最大差值
- 区间最大公约数等
时间复杂度
O(nlogn)预处理,O(1)查询最值
核心思想:倍增思想
参考链接RMQ算法讲解-CSDN博客 ST算法解决RMQ问题详解(图文并茂,保证看懂)_rmq代码c++思路及模板-CSDN博客
RMQ模版
我们可以利用二分的思想,将查询的区间分为长度相等的两部分,取两个区间的最值作为答案。
dp[i][j] 储存第i个位置开始,往后2^j个元素中的最值。即下标[i,i+2^j-1]的区间最值。
我们将区间一分为二,将2^j个元素拆解为2个 2^(j-1)个元素,也就是将[i,i+2^j-1]的区间分为[i,i+2^(j-1)-1]和[i+2^(j-1),i+2^j -1]
那显然问题就很好解决了,我们先初始化每个长度唯一的区间最值(dp[i][0]就表示第i个数字本身),然后可以得到
转移方程dp[i][j] = min(dp [i][j - 1], dp [i + (1 << j - 1)][j - 1])
注:i<<(j-1)==i<<j-1
RMQ预处理模版
void rmq_init()
{
for(int i=1;i<=N;i++)
dp[i][0]=arr[i];//初始化
for(int j=1;(1<<j)<=N;j++)
for(int i=1;i+(1<<j)-1<=N;i++)
dp[i][j]=min(dp[i][j-1],dp[i+(1<<j-1)][j-1]);
}
RMQ查询模版
要注意一下我们利用位运算求出来的每dp【i】【j】的j只是2的次幂 也就是2,4,8,16……并不是每个位置都求出来了所以不能直接查询
例如如果我需要查询【2,13】 那我们并不可以直接求出,而是需要取两个相交区间的最值
int rmq(int l,int r)
{
int k=log2(r-l+1);
return min(dp[l][k],dp[r-(1<<k)+1][k]);
}