RMQ(Range Minimum/Maximum Query),即区间最值查询,是指这样一个问题:对于长度为n的数列A,回答若干次询问RMQ(i,j),返回数列A中下标在区间[i,j]中的最小/大值。
本文介绍一种比较高效的ST算法解决这个问题。ST(Sparse Table)算法可以在O(nlogn)时间内进行预处理,然后在O(1)时间内回答每个查询。
其实RMQ的思路比较好理解,现在我们假设一个二维数组A[i][j]代表的是从顶点i开始走2^j步所到达的位置内的最大值或者是最小值。那么A[i][j]=max / min(A[i+2^(j-1)] [j-1],A[i][j-1])。RMQ的思想与LCA的树上倍增思想都是一样的,现在我们是不是可以预处理出A[i][j]数组,其实j只用开到20左右就可以了,因为2^20已经是10的9次方左右了(int)。
经过预处理之后,查询两点区间中的最大值或者是最小值都是O(1)的时间复杂度。看上图
A[1][2]=max(A[1][1],A[i+2^(j-1)][1]);那么这里的j是如何保证j>(end-star)/2,大一点也没有所谓的。因为可以可以是
求i,j之间的最大值=max(A[i][k],A[end-2^k][k]);
奈斯 贴膜版了;
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define MAXN 200010
int maxsum[MAXN][20],minsum[MAXN][20];//表示从第i个数起连续2^j个数中的最大值/最小值
void RMQ(int num)
{
for(int j=1;j<20;j++)
for(int i=1;i<=num;i++)
{
if(i+(1<<j)-1 <= num)
{
maxsum[i][j]=max(maxsum[i][j-1],maxsum[i+(1<<(j-1))][j-1]);
minsum[i][j]=min(minsum[i][j-1],minsum[i+(1<<(j-1))][j-1]);
}
}
}
int main()
{
int i,j,num,t,query;
int T,k=1;
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&num,&query);
for(i=1;i<=num;i++)
{
scanf("%d",&maxsum[i][0]);
minsum[i][0]=maxsum[i][0];
}
RMQ(num);
int st,en,maxl,minl;
printf("Scenario #%d:\n",k++);
while(query--)
{
scanf("%d%d",&st,&en);
int k=(int)((log(en-st+1))/log(2.0));
maxl=max(maxsum[st][k],maxsum[en-(1<<k)+1][k]);
minl=min(minsum[st][k],minsum[en-(1<<k)+1][k]);
printf("%d\n",minl);
}
}
return 0;
}
有一点注意的是RMQ函数中的两个循环是次数在外 i在里面,初始化—— scanf("%d",&maxsum[i][0]);
minsum[i][0]=maxsum[i][0];