RMQ 问题是指:对于长度为 n 的数列 A,回答若干询问 RMQ (A , i , j ) ( i , j ≤ n),返回数列A中下标在 i , j 里的最小(大)值。
ST算法只适用于静态区间求最值问题。
3 2 4 5 6 8 1 2 9 7
预处理,用DP解决。
设a[i]是要求区间最值的数列,dp[i, j]表示从第i个数起连续2^j个数中的最大值。dp[1,0]表示第1个数起,长度为2^0=1的最大值,其实就是3这个数。 dp[1,2]=5,dp[1,3]=8,dp[2,0]=2,dp[2,1]=4……可以看出dp[i,0]其实就等于a[i],这样就得到dp的初值。
状态转移方程。我们把dp[i,j]平均分成两段(因为f[i,j]一定是偶数个数字),从i到i+2^(j-1)-1为一段,i+2^(j-1)到i+2^j-1为一段(长度都为2^(j-1))。举例说明,当i=1,j=3时就是3,2,4,5 和 6,8,1,2这两段。dp[i,j]就是这两段的最值中的最值。于是得到动态规划方程dp[i, j]=max(dp[i,j-1], dp[i + 2^(j-1),j-1])。
查询
取k=[log2(j-i+1)],则有:rmq(i , j)=max{dp[i,k],dp[j-2^k+1,k]}。 举例说明,要求区间[2,8]的最大值,就要把它分成[2,5]和[5,8]两个区间,因为这两个区间的最大值我们可以直接由dp[2,2]和dp[5,2]得到。
洛谷模板题
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e5+10;
int dp[maxn][30],a[maxn];
void rmq(int n)
{
for(int i=1;i<=n;i++)
dp[i][0]=a[i];
for(int j=1;(1<<j)<=n;j++) //注意j和i的范围
for(int i=1;i+(1<<(j-1))<=n;i++)
dp[i][j]=max(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
}
int query(int l,int r)
{
int k=(int)(log(r-l+1)/log(2)); //k向下取整
return max(dp[l][k],dp[r-(1<<k)+1][k]);
}
int main()
{
int n,m,l,r;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)scanf("%d",&a[i]);
rmq(n);
while(m--){
scanf("%d%d",&l,&r);
printf("%d\n",query(l,r));
}
}