RMQ学习记录

RMQ学习记录
RMQ即区间最值算法，它可以在O(nlogn)时间内进行预处理，然后在O(1)时间内回答每个查询。
具体怎么实现呢？
解决方法：


（一）首先是预处理，用动态规划（DP）解决。


设a[i]是要求区间最值的数列，maxsum[i][j]表示从第i个数起连续2^j个数中的最大值。
即maxsum[i][j]=max(i，i+1......i+2^j-1)
例如：a数列为3 2 4 5 6 8 1 2 9 7
maxsum[1][0]表示第1个数起，长度为2^0=1的最大值，a[1]。
同理 maxsum[1][1]=max(a[1],a[2])=max(3,2)=3.
maxsum[1][2]=max(a[1],a[2],a[3],a[4])=max(3,2,4,5)=5
以此类推...
我们易知maxsum[i][0]=a[i]，即dp的初始化。


（二）具体代码


我们把F[i，j]平均分成两段（因为f[i，j]一定是偶数个数字），从i到i+2^(j-1)-1为一段，i+2^(j-1)到i+2^j-1为一段。当i=1,j=3时就是3,2,4,5 和 6,8,1,2这两段。maxsum[i][j]就是这两段各自最大值中的最大值。
于是我们得到了状态转移方程maxsum[i][j]=max（maxsum[i][j-1], maxsum[i+2^(j-1)][j-1]）。
核心代码如下:
for(int j=1;(1<<j)<=n;j++)  
      for(int i=1;i+(1<<j)-1<=n;i++)  

        maxsum[i][j]=max(maxsum[i][j-1],maxsum[i+(1<<(j-1))][j-1]);  

注意：i和j要倒一倒

（三）查询


假如我们需要查询的区间为(i,j)，那么我们需要找到覆盖这个闭区间(左边界取i，右边界取j)的最小幂（可以重复，比如查询5，6，7，8，9，我们可以查询5678和6789）。
因为这个区间的长度为j - i + 1,所以我们可以取k=log2( j-i+1)，则有：i到j的最大值为max{maxsum[i]maxsum[ j - 2 ^ k + 1, k]}。

举例说明，要求区间[2，8]的最大值，k = log2（8 - 2 + 1）= 2，即求max(F[2, 2]，F[8 - 2 ^ 2 + 1, 2]) = max(F[2, 2]，F[5, 2])；

洛谷p3865

#include<bits/stdc++.h>
using namespace std;
int maxsum[1000001][101],a[1000001];
int main()
{
int n,m;
scanf("%d%d",&n,&m); 
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
}
for(int i=1;i<=n;i++)
{
maxsum[i][0]=a[i];
}
    for(int j=1;(1<<j)<=n;j++)  
      for(int i=1;i+(1<<j)-1<=n;i++)  
        maxsum[i][j]=max(maxsum[i][j-1],maxsum[i+(1<<(j-1))][j-1]);  
for(int i=1;i<=m;i++)
{
int x,y;
scanf("%d%d",&x,&y);
int k=trunc(log2(y-x+1));
printf("%d\n",max(maxsum[x][k],maxsum[y-(1<<k)+1][k]));
}
return 0;
}