什么是ST表

引入

• 最朴素的问题是这样的，给定$n$个数，现在想求出$k_1$到$k_2$的区间最大值，其中$$1\le k_1\le k_2\le n$$
• 这样的问题其实是普遍的，最常见的例子就是高考，所有考生的成绩已经都知道了，如何求某个区间段内的最高分，便于报考

分析

• 此类问题都属于RMQ（Range Minimum Query）问题，也就是区间范围最值问题，朴素的方法是$for$一圈，寻找区间最值，这种方法的问题在于对于多个查询，比如$10^6$次查询，如果区间长度也是$10^6$，那么时间复杂度的数量级是$10^{12}$

解法

模板题

DP

• 朴素算法的最大问题就是每一次都需要一次$for$来寻找最优解，那么可以考虑能否采用备忘录思想，将最优解存储下来，以便于$O(1)$的求出每一次查询的区间最值
• 设$dp[i][j]$表示从第$i$个元素到第$j$个元素的区间最大值，设$Data[i]$为第$i$个元素，那么显然问题具有最优子结构性质，也就是$i$到$j$的区间最大值可以由$i$到$j-1$的区间最大值推出，所以很容易列出状态转移方程为$$dp[i][j]=max(dp[i][j-1],Data[i])$$这时候，查询时间复杂度是$O(1)$，但这带来了一个新的问题，如果要查询第$10^6$到$10^7$元素之间的最大值，那么空间显然无法达到，也就是算法空间复杂度太高，同时，从时间复杂度上看，我们需要两圈$for$循环，时间复杂度是$O(n^2)$，也不能够达到要求，需要考虑优化

倍增优化

• 考虑到上面算法的第二维设置成了区间终点，那么换一种思路，考虑设置为和区间长度相关的量，如果直接设置为区间长度，空间优化也不大，所以采用倍增思想，令$dp[i][j]$表示为起点为$i$，区间长度为$2^j$的区间的最大值，那么应该如何进行状态转移呢？
• 设需要查询的区间长度为$len$，$size=log2(len)$，区间左端点为$l$，右端点为$r$，$dp[i][j]$表示为起点为$i$，区间长度为$2^j$的区间的最大值那么有$len=r-l+1$，考虑将$dp[i][j]$划分为这样两个区间，$dp[i][j-1]$和$dp[i+2^{j-1}][j-1]$，对应于查询区间为$[i,i+2^j]$划分为$[i,i+2^{j-1}]$和$[i+2^{j-1}][i+2^j]$，取这两个区间的最大值就是$dp[i][j]$也就是状态转移方程为$$dp[i][j]=max(dp[i][j-1],dp[i+2^{j-1}][j-1])$$当然，初值应该为$dp[i][0]=Data[i]，$这样$dp$数组就填充完了
• 接下来考虑如何使用$dp$数组，如果问题是求区间$[l,r]$的最大值，需要在$[l,r]$内进行一次划分，划分为$[l,l+2^{size}]$和$[r-2^{size}+1,r]$这样两个集合之间必然有交点，考虑集合两端点的大小可证明得到，这里略去，那么取这两部分的最大值也就得到了$dp[i][j]$，这样查询的时间复杂度就优化到了$O(1)$
• 这个$dp$数组就是ST表，将$dp$数组转换为ST数组如下

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;
const int MAXN = 1e6 + 100;
const double eps = 1e-6;
int ST[MAXN][64];
int main(){
    //freopen("input.txt", "r", stdin);
    //freopen("output.txt", "w", stdout);
    int n, m, l, r;
    scanf("%d%d", &n, &m);
    for(int i=1;i<=n;i++) scanf("%d", &ST[i][0]);
    for(int j=1;j<=log2(n);j++){
        for(int i=1;i+(1 << (j - 1)) <= n;i++){
            ST[i][j] = max(ST[i][j - 1], ST[i + (1 << (j - 1))][j - 1]);
        }
    }
    while(m--){
        scanf("%d%d", &l, &r);
        int len = r - l + 1;
        int sz = log2(len);
        printf("%d\n", max(ST[l][sz], ST[r - (1 << sz) + 1][sz]));
    }
    return 0;
}

练习

洛谷

• 求区间长度为m的区间最大值，起点已知

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;
const int MAXN = 1e6 + 100;
const double eps = 1e-6;
int ST[MAXN][64];
int main(){
    //freopen("input.txt", "r", stdin);
    //freopen("output.txt", "w", stdout);
    int n, m;
    scanf("%d%d", &n, &m);
    for(int i=1;i<=n;i++) scanf("%d", &ST[i][0]);
    for(int j=1;j<=log2(n);j++){
        for(int i=1;i+(1 << (j - 1)) <= n;i++){
            ST[i][j] = min(ST[i][j - 1], ST[i + (1 << (j - 1))][j - 1]);
        }
    }
    for(int i=1;i<=n-m+1;i++){
        int len = m;
        int j = i + len - 1;
        int sz = log2(len);
        printf("%d\n", min(ST[i][sz], ST[j - (1 << sz) + 1][sz]));
    }
    return 0;
}