码迷,mamicode.com
首页 > 其他好文 > 详细

ST表

时间:2020-06-12 00:53:41      阅读:60      评论:0      收藏:0      [点我收藏+]

标签:ant   程序   play   names   isp   数组   比较   是什么   今天   

很好的用二进制来优化了\(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\)的边界呢?

\[\because i+2^j-1\leqslant n \]

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

\[2^j \leqslant n \]

也就是
\(j\leqslant \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\n",ask(l,r));
    }
    return 0;
}

Conclusion

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

Update - 2018.9.15

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

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

ST表

标签:ant   程序   play   names   isp   数组   比较   是什么   今天   

原文地址:https://www.cnblogs.com/eqvpkbz/p/13096805.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!