ST表(Sparse Table) -Updating

于 2023-07-14 19:25:05 发布
阅读量136 收藏
点赞数
文章标签： 算法 数据结构 c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_36549698/article/details/131730228
版权

ST 表是用于解决 可重复贡献问题 的数据结构，是一种离线算法。

可重复贡献问题：简单来说就是求max或min的顺序是不影响结果，每个值可以被多次包括、计算。

常见的可重复贡献问题有：区间最值、区间按位和、区间按位或、区间GCD等。

时间复杂度和空间复杂度

  • 时间复杂度：
    处理ST表：O（log2 n）
    查询ST表：O（1）
  • 空间复杂度：O（n log n）

优缺点

  • 优点： 时间复杂度低，代码量相对线段树少了不少。
  • 缺点： 维护信息有限，拓展较少，修改困难。

核心代码

处理ST表(储存区间最大值)
void initSparseTable(vector<int> a) {
  int la = a.size(), l2 = ceil(log2(la)) + 1;
  for (int i = 0; i < la; ++i) ST[i][0] = a[i];
  for (int j = 1; j < l2; ++j) {
    int pj = 1 << (j - 1);
    for (int i = 0; i + pj < la; ++i) {
      ST[i][j] = max(ST[i][j - 1], ST[i + pj][j - 1]);
    }
  }
}
查询ST表（储存区间最大值）
int Query(int l, int r) {
  int lt = r - l + 1;
  int t = lt == 1 ? 0 : ceil(log2(lt)) - 1;
  return max(ST[l][t], ST[r - (1 << t) + 1][t]);
}
反向ST表实现边建边查（储存最大值）
void UpdateSparseTable() {
  ST[n][0] = a[n];
  for (int j = 1; (1 << j) <= n; ++j) {
    ST[n][j] = max(ST[n][j - 1], ST[n - (1 << (j - 1))][j - 1]);
  }
}
int Query(int l, int r) {
  int lt = r - l + 1;
  int t = lt == 1 ? 0 : ceil(log2(lt)) - 1;
  return max(ST[r][t], ST[l + (1 << t) - 1][t]);
}

练习题目

1.质量检测

题意描述

求区间最小值。

解题思路

ST表

AC代码
/**
 * \link: https://www.luogu.com.cn/problem/P2251
 * \category: ST表 SparseTable 单调队列
 *
 * 
 * 
 **/
// #pragma GCC optimize("O2")

#include <algorithm>
#include <iostream>
#include <cmath>
#include <deque>

#define DBG(x) cout << #x << " = " << (x) << '\n'

using namespace std;

/** using SparseTable. */
const int kN = 1e6 + 233;
vector<vector<int> > ST(kN, vector<int>(21));
void initSparseTable(vector<int> &data) {
  int ld = data.size(), l2 = ceil(log2(ld)) + 1;
  for (int i = 0; i < ld; ++i) ST[i][0] = data[i];
  for (int j = 1; j < l2; ++j) {
    int pj = 1 << (j - 1);
    for (int i = 0; i + pj< ld; ++i) {
      ST[i][j] = min(ST[i][j - 1], ST[i + pj][j - 1]);
    }
  } 
}

int Query(int l, int r) {
  int lt = r - l + 1;
  int t = lt == 1 ? 0 : ceil(log2(lt)) - 1;
  return min(ST[l][t], ST[r - (1 << t) + 1][t]);
}

inline void solve() {
  int n, k; cin >> n >> k;
  vector<int> data(n);
  for (auto &it : data) cin >> it;
  initSparseTable(data);
  for (int i = k - 1; i < n; ++i) {
    cout << Query(i - k + 1, i) << '\n';
  }
  return;
}

/** using monotonic queue. */
inline void solve2() {
  int n, k; cin >> n >> k;
  vector<int> data(n);
  for (auto &it : data) cin >> it;
  vector<int> ans;
  deque<int> que;
  for (int i = 0; i < n; ++i) {
    while (que.size() && data.at(i) <= data.at(que.back())) que.pop_back(); 
    while (que.size() && que.front() <= i - k) que.pop_front();
    if (i >= k - 1) {
      if (que.size()) ans.push_back(data.at(que.front()));
      else ans.push_back(data.at(i));
    }
    que.push_back(i);
  }
  for (auto item : ans) cout << item << '\n';
}

bool rt = false;
signed main() {
  ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
#ifndef ONLINE_JUDGE
  freopen("/Users/Shared/test.in", "r", stdin);
#endif
  if (rt) { int T; cin >> T; while (T--) solve(); }
  else solve();
  return (0 ^ 0);
}

2. Balanced Lineup G

题意描述

求区间最大值减最小值

解题思路

双倍ST表

AC代码
/**
 * \problem: P2880 [USACO07JAN] Balanced Lineup G
 * \link: https://www.luogu.com.cn/problem/P2880
 * \category: SparseTable 树状树组
 * \date: Fri Jul 14 09:43:28 CST 2023 - Fri Jul 14 10:29:30 CST 2023
 * 
 **/
#include <iostream>
#include <vector>
#include <cmath>

using namespace std;

const int kN = 2e5 + 233;
vector<vector<int> > ST1(kN, vector<int>(31));
vector<vector<int> > ST2(kN, vector<int>(31));

void initSparseTable(vector<int> data) {
  int ld = data.size(), l2 = (int)ceil(log2(ld)) + 1;
  for (int i = 0; i < ld; ++i) {
    ST1[i][0] = data.at(i);
    ST2[i][0] = data.at(i);
  }
  for (int j = 1; j < l2; ++j) {
    int pj = 1 << (j - 1);
    for (int i = 0; i + pj < ld; ++i) {
      ST1[i][j] = max(ST1[i][j - 1], ST1[i + pj][j - 1]);
      ST2[i][j] = min(ST2[i][j - 1], ST2[i + pj][j - 1]);
    }
  }
}

int Query(int l, int r, bool isMax) {
  int lt = r - l + 1;
  int t = lt == 1 ? 0 : (int) ceil(log2(lt)) - 1;
  if (isMax)  return max(ST1[l][t], ST1[r - (1 << t) + 1][t]);
  else        return min(ST2[l][t], ST2[r - (1 << t) + 1][t]);
}

inline void solve() {
  int n, q; cin >> n >> q;
  vector<int> data(n);
  for (auto &it : data) cin >> it;
  initSparseTable(data);
  while (q--) {
    int l, r; cin >> l >> r;
    int maxx = Query(l - 1, r - 1, true);
    int minn = Query(l - 1, r - 1, false);
//    cout << "test:" << maxx <<  ' ' << minn << '\n';
    cout << maxx - minn << '\n';
  }
}

bool rt = false;

signed main() {
  ios::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr);
#ifndef ONLINE_JUDGE
  freopen("test.in", "r", stdin);
#endif
  if (rt) { int T; cin >> T; while (T--) solve(); }
  else solve();
  return (0 ^ 0);
}

3.最大数

题目描述

向末尾添加元素的同时求区间最大值

解题思路

反向ST表实现

γυ

AC代码
/**
 * \link: https://www.luogu.com.cn/problem/P1198
 * \category: SparseTable
 *
 * \date:
 **/
// #pragma GCC optimize("O2")

#include <algorithm>
#include <iostream>
#include <cmath>
#define DBG(x) cout << #x << " = " << (x) << '\n'
#define int long long

using namespace std;

const int kN = 2e5 + 233;
int n = 0;
vector<int> a(kN);
vector<vector<int> > ST(kN, vector<int>(31));

void UpdateSparseTable() {
  ST[n][0] = a[n];
  for (int j = 1; (1 << j) <= n; ++j) {
    ST[n][j] = max(ST[n][j - 1], ST[n - (1 << (j - 1))][j - 1]);
  }
}

int Query(int l, int r) {
  int lt = r - l + 1;
  int t = lt == 1 ? 0 : ceil(log2(lt)) - 1;
  return max(ST[r][t], ST[l + (1 << t) - 1][t]);
}

inline void solve() {
  int q, MOD; cin >> q >> MOD;
  int t = 0;
  while (q--) {
    char ch; cin >> ch;
    switch(ch) {
      case 'A': {
        int x; cin >> x;
        a[++n] = (x + t) % MOD;
        UpdateSparseTable();
        break;
      }
      case 'Q' : {
        int x; cin >> x;
        cout << (t = Query(n - x + 1, n)) << '\n';
        break;
      }
    }
  }
}

bool rt = false;
signed main() {
  ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
#ifndef ONLINE_JUDGE
  freopen("test.in", "r", stdin);
#endif
  if (rt) { int T; cin >> T; while (T--) solve(); }
  else solve();
  return (0 ^ 0);
}
YaeSaraki
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
ST表
唐大帅
04-15 1241
P3865 【模板】ST表 题目背景 这是一道ST表经典题——静态区间最大值 请注意最大数据时限只有0.8s，数据强度不低，请务必保证你的每次查询复杂度为O(1) 题目描述 给定一个长度为N的数列，和M次询问，求出每一次询问的区间内数字的最大值。 输入输出格式 输入格式： 第一行包含两个整数N, M 分别表示数列的长度和询问的个数。 第二行包含N个整数（记为 ...
RMQ问题之ST表
最新发布
PengFly123的博客
04-29 521
上面方法2中，f[i][j]表示区间[i,j]的最大值，现在重新定义f[i][j]：以i为起点，长度为2^j的区间，如此一来，f数组的第二维就很小了，比如方法二中的f数组第二维应该是n(n<=10^5)，由于是二维数组，所以整个数组的大小为10^10，这个数组占多少空间呢？4、第二个区间为：[y - 2^j + 1, y]，这个区间的最大值就是f[y - 2^j + 1][j]，所以区间[x, y]的最大值为：max(f[x][j], f[y - 2^j + 1][j])。并向下取整转为int；
Multimodal Sparse Representation-Based Classification for Lung Needle Biopsy Images
02-21
In this study, a novel method, multimodal sparse representation-based classification (mSRC), is proposed for classifying lung needle biopsy images. In the data acquisition procedure of our method, ...
sparse-to-dense:ICRA 2018“稀疏密密麻麻
04-30
稀疏到密集 此存储库由麻省理工学院的 和进行了深度回归神经网络的训练和测试，以用于 。 一个视频演示。 此仓库提供了中该论文的原始实现。 版本可以在找到。 此仓库可用于培训和测试 基于RGB（或灰度图像）的深度...
leetcode2sumc-sparse-table:稀疏表
07-06
稀疏表 前缀和的限制 运算必须是可逆的，即可以通过在更大范围和更小范围上执行运算来计算 它仅适用于范围总和查询 不能用于获取max , min , gcd , bitwise and , bitwise or query 例如，如果给定范围0....4和0.......
Sparse_Table_algorithm.rar_RMQ_Sparse table_Table
09-19
**稀疏表（Sparse Table）算法详解** 稀疏表（Sparse Table）是一种数据结构，用于高效地解决区间最值查询问题（Range Minimum Query, RMQ）。在计算机科学，特别是算法竞赛和离线算法中，这种数据结构具有重要的...
tf.SparseTensor - 知乎1
08-03
tf.SparseTensor tf.SparseTensor 是 TensorFlow 中的一种数据结构，用于表示稀疏矩阵。稀疏矩阵是一种特殊的矩阵，其中大部分元素都是零，只有少部分元素是非零的。 在 TensorFlow 中，稀疏矩阵是通过 tf.Sparse...
sparse r-cnn网络详细
04-28
Sparse R-CNN是一个用于目标检测的深度学习网络，它在Dense R-CNN的基础上进行了改进，通过稀疏卷积实现了更高效的计算。以下是Sparse R-CNN网络的详细介绍： 1. Backbone网络：Sparse R-CNN使用ResNet作为其主干...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值