其实一开始就觉得RMQ算法很实用,大一的时候好像只是记了个小模版,并没有深刻理解并运用到做题中去,现在也遇到了几种用RMQ算法在线询问区间最值问题,以前一直以为RMQ基于DP性质询问区间最值是nlogn,正解当然是O1了,下面详情看一下吧。
首先我们基于DP性质,令maxsum[i][j]为i的j次幂的最大值,从这里就可以观察出,我们所存储的区间长度都是2的幂次,而且对于一段询问(l,r),我们可以以l为开头,例如(3,17),长度为15,我们可以将其拆分为两端幂次和,所以所询问区间最值为【3】【10】与【10】【17】 ,两端长度为8的的区间的最值。复杂度为O1.DP性质的转移也没什么可说的,就是长度为16的可以由两个长度为8的转移过来呗。
#include<bits/stdc++.h>
using namespace std;
const int maxx=1e6+5;
int maxsum[maxx][21],minsum[maxx][21],a[maxx];
int n;
void rmq()
{
int i,j,k;
cin>>n;
for(i=1;i<=n;i++)
scanf("%d",&a[i]);
for(i=1;i<=n;i++)
{
maxsum[i][0]=a[i];
minsum[i][0]=a[i];
}
for(j=1;j<=23;j++)
for(i=1;i+(1<<j)-1<=n;i++)
{
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 query(int l,int r)
{
int k=log(r-l+1)/2-1;
return min(minsum[l][k],minsum[r-(1<<k)+1][k]);
// return max(maxsum[l][k],maxsum[y-(1<<k)+1][k]);
}
int main()
{
rmq();
return 0;
}