RMQ 问题 2 道 PKU 3264 && PKU 2019

今天学习了下RMQ问题的ST算法；

下面是引用以为大牛的解析，讲解的非常透彻，明了。。。。ORZ。。。。

RMQ问题ST算法

ST算法O(nlogn)预处理,O(1)的查询指定区间的最值(以最小值为例)

基本上是把待求区间[l,r]分为两段长为len的区间

左边一段为[l,l+len-1],右边一段为[r-len+1,r]

len必须使得两段区间覆盖待求区间

设所求数组为w

那么,所求最小值就是两个区间的最小值间的最小值

即min(min{w[i],l<=i<=l+len-1},min{w[j],r-len+1<=j<=r})

若都在预先处理中先求得两个区间的最小值

则每次查询的复杂度都是O(1)

---

 

对len做一个限制:只能为2的幂

在预处理中求出所有mi[b][t] : 以b为起点,长为2^t的区间的最小值.

则求解min(min{w[i],l<=i<=l+len-1},min{w[j],r-len+1<=j<=r})

就变成min(mi[l][t],mi[r-2^t+1][r]),其中t可以由此得出,以保证两段区间可以覆盖待求区间:

t=ln(r-l+1)/ln(2)

---

 

可以看到mi[b][t]=min(mi[b][t-1],mi[b+2^(t-1)-1][t-1])

特别地对于所有mi[i][0],其值都是w[i];

由此自底向上推出所有的mi[b][t]

mi大小为n*logn,预处理时间复杂度为O(nlogn),查询时间复杂度为O(1)

 

个人理解RMQ问题的ST算法就是构建多颗二叉树，只不过多颗二叉树没有共同的跟结点，只有共同的叶子结点。。。

 

题意：给定一个序列，求：

给定一个下标区间，求在该区间里面最大值-最小值；

下面是该题的实现代码：

PKU 3264

一道一维的rmq；

#include<iostream> #include<cmath> using namespace std ; #define N 50010 int maxl[N][16], minl[N][16]; int n, m, a[N]; void S_table() { int l = int(log((double)n)/log(2.0)) ; for(int j=1; j <= l; j ++) { for(int i=1; i + (1 << (j - 1)) -1 <= n;i++) { maxl[i][j] = max(maxl[i][j-1], maxl[i + (1<<(j - 1))][j-1]) ; minl[i][j] = min(minl[i][j-1], minl[i + (1<<(j - 1))][j-1]) ; } } } int rmq(int l, int r) { int k = int(log((double)(r - l + 1))/log(2.0)) ; int a1 = max(maxl[l][k], maxl[r - (1 << k) + 1][k]) ; int a2 = min(minl[l][k], minl[r - (1 << k) + 1][k]) ; return a1 - a2 ; } int main() { while(~scanf("%d%d",&n,&m)) { for(int i=1;i <= n;i++) { scanf("%d",&a[i]) ; maxl[i][0] = minl[i][0] = a[i] ; } S_table() ; int a, b ; while(m --) { scanf("%d%d",&a,&b) ; printf("%d/n",rmq(a,b)); } } return 0 ; }

 

2019 二维的rmq；但是这里b是确定的 并且矩阵的范围很小（250 * 250）， 直接暴力也可以解决。。。所以这里没能体现rmq的优势，但是PKU有些题还有有那么些些小缺点，数据范围小，弱了点，有一两次遇到这样的问题了吧，题目和pku上的一样，但是数据扩大的很多倍。这题要是数据范围扩大一些，暴力是行不通的，我用的是ST算法搞的，要是这里b不是确定的，暴力应该也不能过。。。

代码：

#include<iostream> #include<cmath> using namespace std ; #define N 255 int maxl[N][N][8], minl[N][N][8]; int n, b, k, a, c; void S_table() { int l = int(log(double(n)) / log(2.0)) ; for(int k = 1; k <= n ; k ++) { for(int j= 1;j <= l;j ++) { for(int i=1;i + (1 << (j-1))-1 <= n;i ++) { maxl[k][i][j] = max(maxl[k][i][j-1], maxl[k][i + (1 << (j -1 ))][j - 1]) ; minl[k][i][j] = min(minl[k][i][j-1], minl[k][i + (1 << (j -1 ))][j - 1]) ; } } } } int rmq(int x, int y) { int l = int(log(double(b))/log(2.0)) ; int maxt = -1, mint = 300 ; for(int k = x; k < x+b ; k++) { maxt = max(maxt, max(maxl[k][y][l], maxl[k][y+b - (1 << l)][l])); mint = min(mint, min(minl[k][y][l], minl[k][y+b - (1 << l)][l])); } return maxt - mint ; } int main() { while(EOF != scanf("%d%d%d", &n, &b, &k)) { for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { scanf("%d",&a); maxl[i][j][0] = minl[i][j][0] = a ; } } S_table() ; while(k --) { scanf("%d%d", &a, &c) ; printf("%d/n", rmq(a,c)) ; } } return 0; }