RMQ/ST表

应用倍增的思想,主要用来解决区间最值问题,可以做到O(NlogN)预处理,O(1)查询,相比于线段树代码更短,但是不支持修改,是静态数据结构,本质就是一个动态规划。

f(i,j)表示起点为i,区间大小为2j的最大值,即区间[i,i+2j1]里的最大值,那么边界就是f(i,0),即w[i]

递推的时候让子区间成倍增长,得递推式f(i,j)=max(f(i,j1),f(i+2j1,j1),如下图。
image

代码

void init() {
    for (int j = 0; j < M; j++) {
        for (int i = 1; i + (1 << j) - 1 <= n; i++) {
            if (!j) f[i][j] = w[i];
            else f[i][j] = max(f[i][j - 1], f[i + (1 << j - 1)][j - 1]);
        }
    }
}

查询任意区间[l,r]时,先计算一个k,使得2k是小于区间长度得一个最大得k,那么从l开始得2k个数和以r为结尾得2k个数就覆盖了整个区间[l,r],如下图。
image

即使有重叠也没有关系,那么区间最大值就是max(f(l,k),f(r2k+1,k))
代码

int query(int l, int r) {
    int len = r - l + 1;
    int k = log(len) / log(2);
    return max(f[l][k], f[r - (1 << k) + 1][k]);
}

例题AcWing 1273. 天才的记忆

#include <bits/stdc++.h>

using namespace std;

const int N = 2e5 + 10, M = 18;
int f[N][M], w[N]; //f[i][j]表示从i开始长度为2^j的区间最大值是多少
int n, m;

void init() {
    for (int j = 0; j < M; j++) {
        for (int i = 1; i + (1 << j) - 1 <= n; i++) {
            if (!j) f[i][j] = w[i];
            else f[i][j] = max(f[i][j - 1], f[i + (1 << j - 1)][j - 1]);
        }
    }
}

int query(int l, int r) {
    int len = r - l + 1;
    int k = log(len) / log(2);
    return max(f[l][k], f[r - (1 << k) + 1][k]);
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    
    cin >> n;
    for (int i = 1; i <= n; i++) cin >> w[i];
    
    init();
    
    cin >> m;
    while (m--) {
        int l, r;
        cin >> l >> r;
        cout << query(l, r) << endl;
    }
    
    return 0;
}
posted @   Xxaj5  阅读(54)  评论(0编辑  收藏  举报
编辑推荐:
· AI与.NET技术实操系列:基于图像分类模型对图像进行分类
· go语言实现终端里的倒计时
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
阅读排行:
· 25岁的心里话
· 闲置电脑爆改个人服务器(超详细) #公网映射 #Vmware虚拟网络编辑器
· 零经验选手,Compose 一天开发一款小游戏!
· 因为Apifox不支持离线,我果断选择了Apipost!
· 通过 API 将Deepseek响应流式内容输出到前端
点击右上角即可分享
微信分享提示