RMQ(Range Minimum/Maximum Query),即区间最值查询,是指这样一个问题:对于长度为n的数列A,回答若干询问RMQ(A,i,j)(i,j<=n),返回数列A中下标在i,j之间的最小/大值。
st算法就是动态规划,
首先弄清楚dp[i][j]的含义,即表示从i开始的2的j次方个数的最值,将区间分成两部分,每部分都有2的j-1次方个数,一部分是从i到i+2的j-1次方,另一部分是从i+2(j-1)到i+2j,可以得到状态转移方程为dp[i, j]=max(dp[i,j-1], dp[i + 2^(j-1),j-1])。
对于答案的询问
j的长度从1到log2(n),i的遍历从1到n-(2j)
先引入一个k值,每个区间均为2^(k-1)个数字,所以两部分区间为[i, 2^(k-1) + i -1]和[j - 2^(k-1) + 1,j]。显然必须要满足2个子区间能够完全覆盖[i , j],即要求(2^(k-1) + i -1) + 1 >= (j - 2^(k-1) +1),易得2^k >= (j - i + 1),取k的最小值,使2个子区间尽可能的短。所以k = log2(j - i + 1);(向下取整)
所以对于每次询问易得ans = max(dp[i][k],dp[j - 2^k + 1][k])。
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
using namespace std;
const int maxn=1e5+100;
int n,m,x,y;
int dp[maxn][30],a[maxn];
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
dp[i][0]=a[i];
}
for(int j=1;j<=log2(n);j++)
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]);
}
while(m--)
{
scanf("%d%d",&x,&y);
int k=log2(y-x+1);
printf("%d\n",max(dp[x][k],dp[y-(1<<k)+1][k]));
}
return 0;
}