[note]ST表学习笔记

本文介绍了ST表在解决可重复贡献问题(如RMQ和区间GCD)中的应用,利用倍增思想进行预处理,满足结合律,以区间最值问题为例详细解释了ST表的构造和查询方法。
摘要由CSDN通过智能技术生成

ST表这一数据结构以前没有涉及过,今天在集训的时候听师兄讲了大有启发,深感其是个美妙的一种算法。

参考资料:ST 表 - OI Wiki

以及SCUT22级 吕浩天学长集训的 数据结构进阶 PPT

ST表用于解决可重复贡献问题,并且要满足结合律

可重复贡献问题 是指对于运算 opt,满足x opt x = x ,则对应的区间询问就是一个可重复贡献问题。例如,最大值有 max(x,x) = x,gcd 有gcd(x,x) = x,所以 RMQ 和区间 GCD 就是一个可重复贡献问题。像区间和就不具有这个性质,如果求区间和的时候采用的预处理区间重叠了,则会导致重叠部分被计算两次,这是我们所不愿意看到的。另外,opt还必须满足结合律才能使用 ST 表求解。

ST表主要利用了 倍增 的思想,预处理时间O(NlogN)   单次询问时间O(1),但不支持修改,下面以RMQ问题(区间最值问题)为例,解释ST表如何实现

例题:【模板】ST 表 - 洛谷

构造递推关系(类似dp过程)F[i][j]   i表示下标,j表示从 i 开始到 i+2^j-1 区间的最值

我们可以很容易得到关系:

F[i][0] = A[i]

F[i][j] = max(F[i][j - 1], F[i + (1 << (j - 1))][j - 1])

 这样我们便可以预处理完毕

如何利用这一关系求区间最值呢?

首先我们知道,ST表需要满足可重复贡献性质,而最大值时可重复贡献的,也就是说,两个相交区间的重合不影响答案

对于询问区间[l,r]我们只需要找到两个交叉的区间,使着两个区间的并集覆盖[l,r]即可

通过数学知识,我们可以知道\frac{x}{2}< 2^{\left \lfloor \log_{2} {x} \right \rfloor}<=x,也就是说,如果我们令指数值k = log_{2}{r-l+1},则从[l,l+2^k-1][r-2^k+1,r]可以覆盖完全

体现在代码上,则是返回  max(F[l][k], F[r - (1 << k) + 1][k])

完整代码:

// Problem: P3865 【模板】ST 表
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P3865
// Memory Limit: 125 MB
// Time Limit: 800000 ms
//
// Powered by CP Editor (https://cpeditor.org)

#include <fstream>
#include <iostream>
#include <vector>
using namespace std;

const int N = 1e5 + 10;

inline int fr() {
    int x = 0, f = 1;
    char ch = getchar();
    while (ch < '0' || ch > '9') {
        if (ch == '-')
            f = -1;
        ch = getchar();
    }
    while (ch >= '0' && ch <= '9') {
        x = x * 10 + ch - 48;
        ch = getchar();
    }
    return x * f;
}

vector<int> A(N);
vector<vector<int>> f(N, vector<int>(31));

void Build_ST_table(int n) {
    for (int k = 0; k <= 30; k++)
        for (int i = 1; i + (1 << k) - 1 <= n; i++)
            if (!k)
                f[i][k] = A[i];
            else
                f[i][k] = max(f[i][k - 1], f[i + (1 << k - 1)][k - 1]);
    return;
}

int Query(int l, int r) {
    int k = __lg(r - l + 1);
    return max(f[l][k], f[r - (1 << k) + 1][k]);
}

int main() {
    int n, m;
    n = fr();
    m = fr();
    for (int i = 1; i <= n; i++)
        A[i] = fr();

    Build_ST_table(n);

    int l, r;
    for (int i = 1; i <= m; i++) {
        l = fr();
        r = fr();
        printf("%d\n", Query(l, r));
    }

    return 0;
}

  • 9
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值