ST表

Posted 2020-12-10 eqvpkbz

很好的用二进制来优化了(RMQ)中的操作

Body

定义一个数组(f[i][j])表示数列中([i,2 ^ j - 1])这一段区间的极值

可得

[f[i][j] = max(f[i][j-1],f[i+2^{j-1}][j-1]) ]

先说一下这个的正确性,因为我们知道

[f[i][j] = max{f[i,i+2^j-1]} ]

且

[[i, i + 2 ^ j - 1] = [i, i + 2 ^ {j - 1} - 1],[i + 2 ^ {j - 1}, i + 2 ^ j - 1] ]

其实也就是从中间分开，分成两个相等的区间
那么f[i][0]又是什么呢？
它就是每个数组对应的值
因为

[f[i][0] = max{f[i,i+2^j-1]} ]

[2^0 = 1 ]

所以

[f[i][0] = max{f[i,i]} ]

所以这个时候所有的f[i][0]都是当前节点的值
所以也就比较好理解了
然后还有一个比较重要的问题: 循环边界
首先我们知道(i+2^j-1 leqslant n)
所以因为是先枚举j，再枚举的i，所以这个时候就可以作为i的边界

那么(j)的边界呢？

[ecause i+2^j-1leqslant n ]

所以当(i)最小为(1)时

[2^j leqslant n ]

也就是
(jleqslant log_2{n})

好了，那么下面就是主程序

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
#define maxn 1000005
int f[maxn][21];
int ask(int l,int r)
{
    int k=log2(r-l+1); 
    return max(f[l][k],f[r-(1<<k)+1][k]);
}
int main()
{
    int n,m;
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&f[i][0]);
    int Max_j = log2(n);
    for(int j=1;j<=Max_j;j++)
        for(int i=1;i+(1<<j)-1<=n;i++)
            f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]); 
    for(int i=1;i<=m;i++)
    {
        int l,r;
        scanf("%d %d",&l,&r);
        printf("%d
",ask(l,r));
    }
    return 0;
}

Conclusion

ST表方便便捷，支持(O(1))查询，比较适合在不变的序列上查找最值

Update - 2018.9.15

比较重要的一点是循环的顺序，如果i和j的位置没有搞好的话那么将会出很大的锅

今天考试的时候才发现这里有循环顺序错误啥的