首先先介绍一下RMQ算法,RMQ算法是针对区间最值的一个算法,有些问题可以用线段树实现但是你不觉得 有些问题RMQ仅仅几行的代码机上超级快的查询速度是一件很美好的事情吗?!!!!
RMQ原理:定义一个二维数组dp[i][j],表示的是从第i个元素开始,长度(1<<j)内的最值,我们用动态规划的思想来实现
dp[i][j]=max{dp[i][j-1],dp[i+(1<<(j-1))][j-1];
dp[i][j]=min{dp[i][j-1],dp[i+(1<<(j-1))][j-1];
实现代码如下
void RMQ(int num) //预处理->O(nlogn)
{
for(int i=1;i<=num;i++)
dp[i][0]=a[i];
for(int j = 1; (1<<j) <= num; ++j)
for(int i = 1; i + (1 << j) - 1 <= num; ++i)
{
dp[i][j] = max(dp[i][j - 1], dp[i + (1 << (j - 1))][j - 1]);
}
}
之所以j在上i在下 我们得明白我们是先给第i个位置确定了一个初值,然后我们只能
从已经有了的值去推出值 所以只能慢慢来从短到长 由已知到未知 大家好好思考一下这个问题。
接着就是查询了,查询的思路很简单就是找到一个最大的k使(1<<k)<=r-l+1;
int query(int l,int r)
{
int k=0;
while(1<<(k+1)<=r-l+1)
k++;
return max(dp[l][k],dp[r-(1<<k)+1][k]);
}
查询的时间复杂度只有O(1) 。
最后我们来看一道例题 NYOJ119:http://acm.nyist.net/JudgeOnline/status.php?pid=119
代码如下 :
#include <iostream>
#include <cstdio>
using namespace std;
int n,m;
int maxsum[100000][20];
int minsum[100000][20];
void RMQ(int num) //预处理->O(nlogn)
{
for(int j = 1; (1<<j) <= num; ++j)
for(int i = 1; i + (1 << j) - 1 <= num; ++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=0;
while(1<<(k+1)<=r-l+1)
k++;
return max(maxsum[l][k],maxsum[r-(1<<k)+1][k])-min(minsum[l][k],minsum[r-(1<<k)+1][k]);
}
int main()
{
while(~scanf("%d%d",&n,&m))
{
int x;
for(int i=1;i<=n;i++)
{
scanf("%d",&x);
maxsum[i][0]=x;
minsum[i][0]=x;
}
RMQ(n);
int l,r;
while(m--)
{
scanf("%d%d",&l,&r);
printf("%d\n",query(l,r));
}
}
return 0;
}