ST表的简单应用

1. ST表简介

ST表是基于倍增思想，用于解决重复贡献问题的数据结构。
ST表最广泛的应用就是解决RMQ问题（区间最值问题），ST表可以通过O(n logn)的时间进行预处理，O(1)的时间进行问题查询。但是ST表不支持单点修改和区间修改。

2.ST表的原理

ST表利用了重复贡献问题，例如：$Max(n,n)=n,Min(n,n)=n，GCD(n,n)=n$ ，求解时区间有重叠部分并不会影响我们对结果的计算，所以 RMQ 和区间 GCD 就是一个可重复贡献问题。像区间和就不具有这个性质，如果求区间和的时候采用的预处理区间重叠了，则会导致重叠部分被计算两次，这是我们所不愿意看到的。

我们用$f[i][j]$来表示从第 $i$ 个数起连续 $2^j$ 个数中的最小值，也就是区间$[i,i+2^j-1]$的最值。

显然 $f[i][0]=x_i，$我们可以得到状态转移方程：

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

所以在预处理的时候，我们不妨先枚举区间长度 $j$，再枚举区间左端点 $i$，只要能保证$i+(1<<j)-1\le n$，那么我们便可以将各种区间给预处理出来。

看不懂的可以自己跟着代码试着走一下

void st(int n)
{
    for (int j = 1; j <= 20; j++) //从j开始枚举区间长度
        for (int i = 1; i <= n; i++)
            if (i + (1 << j)-1 <= n)
                f[i][j] = min(f[i][j - 1], f[i + (1 << (j - 1))][j - 1]);
}

以上就是预处理部分。而对于查询，可以简单实现 如下：

我们可以找距离区间长度$（R-L+1）$最近的一个 $2^n$ 次方数 $k$ ，那么我们就可以用 $L$ 为左端点 长度为$k$的区间 $and$ 以 $R$ 为右端点 区间长度为 $k$ 的两个区间 来更新所求的区间。

不懂的可以画画图！

也就是说对于每个询问区间$[l，r]$，我们把它分成两个部分：$[l,l+2^s-1]$与$[r-2^s+1，r]$，其中，$s=\lfloor lo{g}_2(r-l+1)\rfloor$（下取整），两部分结果的最大值就是回答。

3.例题

洛谷 P1816 忠诚

题目描述

老管家是一个聪明能干的人。他为财主工作了整整 $10$ 年。财主为了让自已账目更加清楚，要求管家每天记 $k$ 次账。由于管家聪明能干，因而管家总是让财主十分满意。但是由于一些人的挑拨，财主还是对管家产生了怀疑。于是他决定用一种特别的方法来判断管家的忠诚，他把每次的账目按 $1,2,3\dots$ 编号，然后不定时的问管家问题，问题是这样的：在 $a$ 到 $b$ 号账中最少的一笔是多少？为了让管家没时间作假，他总是一次问多个问题。

输入格式

输入中第一行有两个数 $m,n$，表示有 $m$ 笔账 $n$ 表示有 $n$ 个问题。

第二行为 $m$ 个数，分别是账目的钱数。

后面 $n$ 行分别是 $n$ 个问题，每行有 $2$ 个数字说明开始结束的账目编号。

输出格式

在一行中输出每个问题的答案，以一个空格分割。

样例 #1

样例输入 #1

10 3
1 2 3 4 5 6 7 8 9 10
2 7
3 9
1 10

样例输出 #1

2 3 1

蒟蒻代码

#include <bits/stdc++.h>
using namespace std;

#define int long long
#define xx first
#define yy second
typedef pair<int, int> PII;
const int inf = 0x3f3f3f3f;
const int N = 1e5 + 10;

int v[N], n, m, k, x, t, T;
int f[N][20];

void ClearFloat()
{
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
}

void st(int n)
{
    for (int j = 1; j <= 20; j++)
        for (int i = 1; i <= n; i++)
            if (i + (1 << j)-1 <= n)
                f[i][j] = min(f[i][j - 1], f[i + (1 << (j - 1))][j - 1]);
}

signed main()
{
    ClearFloat();
    cin >> n >> m;
    for (int i = 1; i <= n; i++)
        cin >> f[i][0];
    st(n);
    while (m--)
    {
        int l, r;
        cin >> l >> r;

        int k = log2(r - l + 1);
        cout << min(f[l][k], f[r - (1 << k) + 1][k]) << " ";
    }
    
    return 0;
}

记录一下第一次用MarkDown语法写的文章 😁