ST表

最新推荐文章于 2021-10-13 18:54:52 发布

可可亚

最新推荐文章于 2021-10-13 18:54:52 发布

阅读量141

点赞数

分类专栏：数据结构文章标签：算法数据结构

本文链接：https://blog.csdn.net/weixin_44232130/article/details/113741448

版权

数据结构专栏收录该内容

2 篇文章 0 订阅

订阅专栏

本文介绍了如何使用ST表（Segment Tree）数据结构解决静态区间最大值查询问题，通过动态规划实现O(1)查询速度。适用于给定区间内的最大值查找，但适用于数据不更新的情况。核心包括预处理阶段的动态规划计算和查询时的区间覆盖策略。

摘要由CSDN通过智能技术生成

解决的问题

静态RMQ问题，区间查询最大值最小值
给定一个区间，多次让你查询 $l\sim r$ 之间的最大值
优点：

查询的时间复杂度为O(1)，查询速度极快
写起来比较简便，容易理解

缺点：

只能解决静态问题，无法修改区间内的值

思想

①预处理
主要思想为动态规划，用一个数组 $s t [i] [j]$ ，表示区间 $\sim i + 2^j - 1$ 内的最大值，意思就是从第 $i$ 个开始，之后 $2^j$ 个数中的最大值，因此，可以用一个数组就保存最后的值。显而易见，转移方程为 $st[i][j] = max(st[i][j-1],st[i+(2^{j-1})][j-1])$ 由一个范围前一半的最大值和后一半的最大值取一个最大值，即是总体最大值.
（用 $1 < < (j)$ 位运算表示 $2^j$ ）

void RMQ(){
    for (int i = 1; i <= n; i ++){
        st[i][0] = a[i];
    }
    //先枚举倍数，在枚举起点
    for (int j = 1; (1 << j) <= 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]);
        }
    }
}

②查询
那么现在能够查询从 $i$ 到之后的 $2^j$ 个数的最大值了，但是我们要求的是任意两个边界之间的最大值。
我们的方法是

预处理一个Log函数，求每一个值的log值，向下取整

void GetLog(){
    Log[1] = 0;
    for (int i = 2; i <= n; i ++){
        Log[i] = Log[i / 2] + 1;
    }
}

$\space k = log(r - l+1)$ ,那么 $2 * k$ 一定大于最后的总个数，我们分别取两个区间 $l\sim l + 2^k - 1$ 以及 $l-2^k+1 \sim r$ ，分别可以用 $s t [l] [k]$ 以及 $st[l-2^k+1][k]$ 表示，这两个区间一定覆盖整个 $l\sim r$ 的区间，求其最大值即可。

int query(int l, int r){
    int k = Log[r - l + 1]; //k为长度
    int res;
    res = max(st[l][k], st[r - (1 << k) + 1][k]);
    return res;
}

模版题

P3865 【模板】ST表

在这里插入图片描述

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define INF 0x7fffffff;
const int maxn = 1e5 + 5;
int a[maxn];
int Log[maxn];
int st[maxn][20];
int n, m;
void GetLog(){
    Log[1] = 0;
    for (int i = 2; i <= n; i ++){
        Log[i] = Log[i / 2] + 1;
    }
}
void RMQ(){
    for (int i = 1; i <= n; i ++){
        st[i][0] = a[i];
    }
    for (int j = 1; (1 << j) <= 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]);
        }
    }
}
int query(int l, int r){
    int k = Log[r - l + 1]; //k为长度
    int res;
    res = max(st[l][k], st[r - (1 << k) + 1][k]);
    return res;
}
int main()
{  
    cin >> n >> m;
    GetLog();
    for (int i = 1; i <= n; i ++) cin >> a[i];
    RMQ();
    for (int i = 1; i <= m; i ++){
        int l, r;
        scanf("%d%d", &l, &r);
        printf("%d\n", query(l, r));
    }
    return 0;
}