RMQ学习总结(附带简单例题)
RMQ (Range Minimum/Maximum Query)算法
状态转移方程是:
f [ i ] [ j ] = max( f [ i ] [ j - 1 ] , f [ i + 2^(j - 1) ] [ j - 1 ] )
模板:
找出区间内的最大值:
#include<bits/stdc++.h>
using namespace std;
int a[10002],f[10002][18],n;
void RMQ_prepare()
{
for(int i=1;i<=n;i++)
f[i][0]=a[i];//数据初始化,即从i开始向右走2的0次方的区间中的最大值
for(int j=1;(1<<j)<=n;j++)
for(int i=1;i+(1<<j)-1<=n;i++)//i+(1<<j)-1<=n是为了保证区间左端点不超出总数n
f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);//实质是动态规划
}
int query(int x,int y)
{
//int k=trunc(log2(y-x+1));//注意y-x要加一才为区间长度
int k = (int)(log(y - x + 1.0) / log(2.0));
return max(f[x][k],f[y-(1<<k)+1][k]);
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++){
scanf("%d",&a[i]);
}
int q;
scanf("%d",&q);
RMQ_prepare();//预处理
while(q--)
{
int x,y;
scanf("%d%d",&x,&y);
printf("%d\n",query(x,y));
}
return 0;
}
输出区间内最大值和最小值的差:
#include<bits/stdc++.h>
using namespace std;
int a[50005],b[50005],f_max[50005][20],f_min[50005][20],n;
void RMQ_prepare()
{
for(int i=1;i<=n;i++)
f_max[i][0]=a[i],f_min[i][0]=a[i];
for(int j=1;(1<<j)<=n;j++)
for(int i=1;i+(1<<j)-1<=n;i++){
f_max[i][j]=max(f_max[i][j-1],f_max[i+(1<<(j-1))][j-1]);
f_min[i][j]=min(f_min[i][j-1],f_min[i+(1<<(j-1))][j-1]);
}
}
int main()
{
int Q,x,y,aa,bb;
scanf("%d%d",&n,&Q);
for(int i=1;i<=n;i++){
scanf("%d",&a[i]);
}
RMQ_prepare();
for(int j=0;j<Q;j++){
scanf("%d%d",&x,&y);
int k = (int)(log(y - x + 1.0) / log(2.0));
aa=max(f_max[x][k],f_max[y-(1<<k)+1][k]);
bb=min(f_min[x][k],f_min[y-(1<<k)+1][k]);
printf("%d\n",aa-bb);
}
return 0;
}