倍增RMQ算法浅谈

题目：

求某一区间的最大值。

 

f[a][b]中a代表的是当前的位置，b代表的是以a为起点往后移动的区间长度2^b。

void bz(int n)
{
    for(int i=1;i<=n;i++)
        f[i][0]=a[i];//先定义长度1为自己本身。
    for(int j=1;(1<<j)<=n;j++)//以区间逐步移长。//tip1
    {
        for(int i=1;i+(1<<j)-1<=n;i++)//点的位置逐步往后移。//tip2
        　　f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
    }
}

其实tip1与tip2是利用以小的区间的最大值来逐步往上推更大的区间的最值。具体怎么推呢就需要利用到倍增思想（其实跟二分很相似）。

我们来模拟一下吧！

f[1][1]=max(f[1][0],f[2][0]);f[1][2]=max(f[1][1],f[3][1]);

f[1][1]={1,2},f[1][0]={1}，f[2][0]={2},f[1][2]={1,2,3,4},f[3][1]={3,4};//大括号里面代表的是位置。

因此可以看出要求f[1][1]就直接看1位置与2位置的最值，f[1][2]是由f[1][1]与f[3][1]比较得出的。在这里也体现了由小区来推广到大区间，因为所求的大区间都可以通过二进制关系分为两个小区间问题（若重叠也是没关系的），而最小的问题也就是最初定义f[a][0]为a位置本身的问题；

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=1e3+7;
int num[maxn],f[1000][1000];
int main()
{
    int n,a,b,c,d,q,l,r;
    while(~scanf("%d",&n))
    {
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&num[i]);
        }
        for(int i=1;i<=n;i++)
        {
            f[i][0]=num[i];
        }
        for(int i=1;i<20;i++)
        {
            for(int j=1;j+(1<<(i-1))-1<=n;j++)
            {
                f[j][i]=max(f[j][i-1],f[j+(1<<(i-1))][i-1]);
            }
        }
        cin>>q;
        for(int i=1;i<=q;i++)
        {
            cin>>l>>r;
            int k=(int)(log(r-l+1.0)/log(2.0));
            printf("%d\n",max(f[l][k],f[r-(1<<k)+1][k]));
        }
    }
}