RMQ算法之一——ST算法

一.RMQ问题

对于长度为n的数列A，回答若干如后询问：RMQ（A，i，j）（i，j<=n）——返回数列A中下标在i->j的最小（大）值，即：求区间最值

因此，求解RMQ问题的算法可视为RMQ算法，这里主要介绍ST（Sparse_Table）算法。

二.ST算法

ST算法只适用于静态区间求最值，动态区间的话用线段树比较好。
ST算法本质相当于动态规划，以求最大值为例
我们用f[i][j]来表示以i为起点，连续2^j个数的区间的最大值（例如：f[2][2]表示第2到第5个数的最值）
本算法共分为两部分，预处理与查询
（需要对两个数组（一个二维f[i][j]，一个一维log[i]）进行预处理）
（1）预处理

f[i][j]预处理

设a[n]为原序列，因为f[i][0]为从第i个数开始连续2^0=1个数的最值，即第i个数本身，所以f数组初始化为f[i][0]=a[i]；
然后，将每个要求的区间分为两部分，分别求两部分的最大值，而整个区间的最大值为两个最大值中较大的那个，由此可列出转移方程：f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
(注：(i<<j)为位运算，表示2^j)

void RMQ()
{
    int i,j;
    for(i=1;i<=n;i++)
        f[i][0]=a[i];
    for(j=1;(1<<j)<=n;j++)
        for(i=1;i+(1<<(j-1))<=n;i++)
            f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
    (此处为求最大值，若求最小值，只需把max换成min即可)
}

f[i][j]预处理注意i与j的双重循环，j为第一重，i为第二重，不能倒过来

如果i在外，j在内的话，我们更新的顺序就变成了从1开始的前1、2、4、8…个元素。

当j等于3，即求一个长为2³=8的区间的时候:
dp[1][3]=max(max(ans[0],ans[1],ans[2],ans[3]),max(ans[4],ans[5],ans[6],ans[7])))

但此时我们并没有计算max(ans[4],ans[5],ans[6],ans[7])，所以这样的方法是错误的。

log数组预处理

注意：
1.log数组取值向下取整
2.log[i]=log[i/2]+1;
(此处与下方查询部分对照之后比较好理解)

void get_log()
{
    int i;
    log[1]=0;
    for(i=2;i<=n+1;i++)
        log[i]=log[i/2]+1;
}

(2)查询

设查询区间为[l,r]，长度为L，则：L=r-l+1；
设2^k=L，但L不一定刚好等于2^k，又因为log向下取整，所以2^k可能小于L；
解决办法：从l向右取2^k个（即：f[i][k]），从r向右取2^k个数（即：f[r-(1<<k)+1][k]），那么最终结果:answer=max(f[l][k],f[r-(1<<k)+1][k]);

int search(int l,int r)
{
    if(l>r)
        return 0;
    int k=log[r-l+1];
    return max(f[l][k],f[r-(1<<k)+1][k]);
}
int main()
{
    int i,l,r,m,ans;
    ......
    scanf("%d",&m);
    for(i=1;i<=m;i++)
    {
        scanf("%d %d",&l,&r);
        ans=search(l,r);
        printf("%d\n",ans);
    }
    ......
    return 0;
}

(先这样，后续应该大概可能也许会补充…=.=||)