【算法】ST表
想学习一下LCA倍增,先 水一个黄题 学一下ST表
ST表
介绍:
这是一个运用倍增思想,通过动态规划来计算区间最值的算法
算法步骤:
求出区间最值
回答询问
求出区间最值:
用f[i][j]来存储从第 j 个点开始,向后 2 ^ i - 1 个点中的最值(包括本身)
利用二分法的思想,将区间 [ j,j +(2 ^ i)- 1 ] 平均(大概)分成两半
可以算出,区间 [ j,j +(2 ^ i)- 1 ] 的长度为 2 ^ i
所以一半的长度为 2 ^ i - 1
那么分成的两个区间就为 [ j,j +(2 ^(i - 1)- 1 ] 和 [ j +(2 ^ i - 1 ),j +(2 ^ i)- 1 ] 。
毫无疑问,
f[i][j] = max(f[i - 1][j],f[i - 1][j +(1 << i - 1)])
这样递推式就出来了
现在来想一下:
f[0][j]就是从 j 开始向后数第 2 ^ 0 - 1 个点的最值,区间为 [ j,j ]
不用考虑,f[0][j]就是第 j 个数本身
好了,现在边界也得出来了,可以开始 dp 了
上代码:
void prew() {
F1(i, 1, n) f[0][i] = a[i]; // 给边界赋值,a[i] 存的是数列的第 i 个数
int kf = log2(n); // 得出数列最多可以向后跳几个 2 的幂,n 为数列长度
F1(i, 1, kf) { // 枚举区间的长度 2 ^ i
for (int j = 1; j + (1 << i) - 1 <= n; j++) // 枚举起点
f[i][j] = max(f[i - 1][j], f[i - 1][j + (1 << i - 1)]); // 递推式
}
}
回答询问:
由于 log2 运算可能会出现实数,而我们又用整数类型来存储,所以可能会出现两个区间不能完全覆盖整个区间的情况,得出的 f[i][j]不够精准
最简单的方法就是用两个区间覆盖
反正又没要求两个区间不能重叠~~
可以选用f[k][l]和f[k][r-(1<<k)+1]来覆盖f[l][r]
所以f[l][r] = max(f[k][l],f[k][r -(1 << k)+ 1])(k 为区间 [l,r] 的长度的 log2)
再上代码:
int ask(int l, int r) {
int k = log2(r - l + 1); // 求出区间最大的 log2 值
return max(f[k][l], f[k][r - (1 << k) + 1]); // 返回区间 [l,r] 的最大值
}
完整代码:
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iomanip>
#include <cmath>
#include <algorithm> // 妈妈再也不怕我的头文件不够使啦~~
#define MAXN 100100
#define INF 0x3f3f3f3f
#define LL long long
#define F1(i, a, b) for (LL i = a; i <= b; ++i) // 懒人必备神器
#define F2(i, a, b) for (LL i = a; i >= b; --i)
using namespace std;
int f[31][MAXN], a[MAXN];
//f[i][j]表示从 j 往后 2 ^ i - 1 个数中的最大值
int n, m;
inline int read() { // 快读
int sto = 0, fg = 1;
char ch = getchar();
while (ch < '0' || ch > '9') {
if (ch == '-') fg = -1;
ch = getchar();
}
while (ch >= '0' && ch <= '9') {
sto = (sto << 1) + (sto << 3) + (ch ^ 48);
ch = getchar();
}
return sto * fg;
}
void prew() { // 预处理 dp
F1(i, 1, n) f[0][i] = a[i];
int kf = log2(n);
F1(i, 1, kf) {
for (int j = 1; j + (1 << i) - 1 <= n; j++)
f[i][j] = max(f[i - 1][j], f[i - 1][j + (1 << i - 1)]);
}
}
int ask(int l, int r) { // 回答询问
int k = log2(r - l + 1);
return max(f[k][l], f[k][r - (1 << k) + 1]);
}
int main()
{
int l, r, ans;
n = read(); m = read();
F1(i, 1, n) a[i] = read();
prew();
F1(i, 1, m) {
l = read(); r = read();
ans = ask(l, r);
printf("%d\n", ans);
}
return 0;
}
模板题:
洛谷P3865
【算法】ST表的更多相关文章
- 【基础算法-ST表】入门 -C++
前言 学了树状数组看到ST表模板跃跃欲试的时候发现完全没思路,因为给出的查询的时间实在太短了!几乎是需要完成O(1)查询.所以ST表到底是什么神仙算法能够做到这么快的查询? ST表 ST表是一个用来解 ...
- 线段树(two value)与树状数组(RMQ算法st表)
士兵杀敌(三) 时间限制:2000 ms | 内存限制:65535 KB 难度:5 描述 南将军统率着N个士兵,士兵分别编号为1~N,南将军经常爱拿某一段编号内杀敌数最高的人与杀敌数最低的人进行比 ...
- 浅谈ST表
发现自己学的一直都是假的ST表QWQ. ST表 ST表的功能很简单 它是解决RMQ问题(区间最值问题)的一种强有力的工具 它可以做到$O(nlogn)$预处理,$O(1)$查询最值 算法 ST表是利用 ...
- ST表学习笔记
ST表是一种利用DP思想求解最值的倍增算法 ST表常用于解决RMQ问题,即求解区间最值问题 接下来以求最大值为例分步讲解一下ST表的建立过程: 1.定义 f[i][j]表示[i,i+2j-1]这个长度 ...
- 浅谈 倍增/ST表
命题描述 给定一个长度为 \(n\) 的序列,\(m\) 次询问区间最大值 分析 上面的问题肯定可以暴力对吧. 但暴力肯定不是最优对吧,所以我们直接就不考虑了... 于是引入:倍增 首先,倍增是个什么 ...
- LCA 算法(一)ST表
介绍一种解决最近公共祖先的在线算法,st表,它是建立在线性中的rmq问题之上. 代码: //LCA: DFS+ST(RMQ) #include<cstdio> #include&l ...
- 算法学习 - ST表 - 稀疏表 - 解决RMQ问题
2017-08-26 21:44:45 writer:pprp RMQ问题就是区间最大最小值查询问题: 这个SparseTable算法构造一个表,F[i][j] 表示 区间[i, i + 2 ^ j ...
- [算法模板]ST表
[算法模板]ST表 ST表和线段树一样,都能解决RMQ问题(范围最值查询-Range Minimum Query). 我们开一个数组数组\(f[maxn][maxn\log_2]\)来储存数据. 定义 ...
- 【算法学习笔记】RMQ问题与ST表
\(0.\) RMQ问题 P1816 人话翻译 给定一个长度为\(n\)的数列\(a\),然后有\(m\)组询问,每次询问一个区间\([l,r]\)的最小值. 其中\(m,n\leq10^5\) \( ...
- 【JZOJ5064】【GDOI2017第二轮模拟day2】友好城市 Kosarajo算法+bitset+ST表+分块
题面 在Byteland 一共有n 座城市,编号依次为1 到n,这些城市之间通过m 条单向公路连接. 对于两座不同的城市a 和b,如果a 能通过这些单向道路直接或间接到达b,且b 也能如此到达a,那么 ...
随机推荐
- 看完这篇,保证让你真正明白:分布式系统的CAP理论、CAP如何三选二
引言 CAP 理论,相信很多人都听过,它是指: 一个分布式系统最多只能同时满足一致性(Consistency).可用性(Availability)和分区容错性(Partition tolerance) ...
- 【译】对Rust中的std::io::Error的研究
原文标题:Study of std::io::Error 原文链接:https://matklad.github.io/2020/10/15/study-of-std-io-error.html 公众 ...
- [leetcode]82. Remove Duplicates from Sorted List
第一题:遍历链表,遇到重复节点就连接到下一个. public ListNode deleteDuplicates(ListNode head) { if (head==null||head.next= ...
- 96. Unique Binary Search Trees1和2
/* 这道题的关键是:动态表尽量的选取,知道二叉搜索树中左子树的点都比根节点小,右子树的点都比根节点大 所以当i为根节点,左子树有i-1个点,右子树有n-i个点,左右子树就可以开始递归构建,过程和一开 ...
- JAVA数据结构(十一)—— 堆及堆排序
堆 堆基本介绍 堆排序是利用堆这种数据结构而设计的一种排序算法,堆排序是一种选择排序,最坏,最好,平均时间复杂度都是O(nlogn),不稳定的排序 堆是具有以下性质的完全二叉树:每个节点的值都大于或等 ...
- centos7安装手册
一:centos7下载 下载centos可以访问阿里云完成centos 7.7下载(建议使用迅雷下载):http://mirrors.aliyun.com/centos/7.9.2009/isos/x ...
- 冒泡排序算法JAVA实现版
/***关于冒泡排序,从性能最低版本实现到性能最优版本实现*/public class BubbleSortDemo { public static void sort(int array[]) { ...
- layer做阻塞式弹出层的方法
今天遇到一个问题: layer弹出一个confirm提示窗,然后confirm还没有点击对应的按钮的时候,就已经执行了后续代码,我这里做出的判断是,是否需要进行后续操作,但是因为layer.confi ...
- 解析SwiftUI布局细节(二)循环轮播+复杂布局
前言 上一篇我们总结的主要是VStack里面的东西,由他延伸到 @ViewBuilder, 接着我们上一篇总结的我们这篇内容主要说的是下面的几点,在这些东西说完后我准备解析一下苹果在SiwftUI文档 ...
- 坐标转换成SVG的path路径
大家好,我是一个刚入职的前端小白,入职后一直做关于svg 的东西,我将自以为很方便的方法提供给大家. function svgPathCurv(a,b,curv) { /* * 弯曲函数. * a:a ...