RMQ(Range Minimum/Maximum Query):区间最值查询
(最近才知道原来RMQ是一大类问题,就是区间最值查询的问题,可以用线段树解决)
另外就是用下面写的ST(Sparse Table[稀疏表])算法【ST算法用了倍增的思想】:
ST是一种在线算法,预处理部分O(nlogn),查询部分O(1)
预处理部分:
max_[i][j]=max(max_[i][j-1],max_[i + (1 << j-1)][j-1]);
min_[i][j]=min(min_[i][j-1],min_[i + (1 << j-1)][j-1]);
预处理部分的复杂度:O(nlogn)
释意:
- max_[ i ][ j ],从区间的第 i 的位置开始,长度为2 ^ j 的区间上的最大值。
- 我们要求区间的最大值,可以将区间从中间分成两个子区间,然后取两个子区间的最大值中的较大者。区间长度一分为二,后边第二个下标变为 j-1 是很好理解的。
- max_[ i ] [ j-1 ]表示的区间是[ i , i+2 ^ (j-1)-1 ]
- 而max_[ i + (1 << j-1) ][ j-1 ])表示的区间是[ i+2 ^ (j-1) , i+2 ^ j-1]
- (可能有些人跟我一样不理解这个区间,想个例子其实超级好理解:比如从区间的1位置开始,长度为8,分成两个区间其实是[ 1 , 4 ]和[ 5 , 8 ],和上面的式子联系一下就明白了)
查询部分:
int query_max(int l,int r)
{
int k=log2(r-l+1);//向下取整
return max(max_[l][k],max_[r-(1<<k)+1][k]);
}
查询部分复杂度:O(1)
释意:
- 我们取 k 为区间长度以2为底的指数(向下取整)
- 如果区间长度刚好是2的次幂,max_[ l ][ k ]和max_[ r-(1<<k)+1 ][ k ]两者表示的是一样的区间,并且已经预处理过了,我们是不是可以直接返回已经预处理好的值?
- 如果不是呢,我们就向下取整log2( r-l+1 ),然后将区间分成有重合部分的长度都为 2 ^ k 的两个区间:[ l , l+2 ^ k-1 ] [ r-2^k+1 , r ]这时候我们手动分的这两个区间是已经预处理好的,是不是可以直接返回其中的较大值即可?
区间最小值与上所述同理
例题+代码:
- 题意:区间最小值
- 思路:模板
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <cmath>
#include <cstring>
#include <string>
#include <vector>
#include <set>
#include <stack>
#include <list>
#include <map>
#define INF 0x3f3f3f3f
#define P(x) x>0 ? x : 0
#define MID l + r >> 1
#define lsn rt << 1
#define rsn rt << 1 | 1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
using namespace std;
typedef long long ll;
typedef vector<int>:: iterator VITer;
const int maxN=5e4+5;
int max_[maxN][20],min_[maxN][20];
int N,Q;
void RMQ_pre()
{
for(int j=1;(1<<j)<=N;j++)
{
for(int i=1;i + (1 << j-1)<=N;i++)
{
max_[i][j]=max(max_[i][j-1],max_[i + (1 << j-1)][j-1]);
min_[i][j]=min(min_[i][j-1],min_[i + (1 << j-1)][j-1]);
}
}
}
int query_min(int l,int r)
{
int k=log2(r-l+1);//向下取整
return min(min_[l][k],min_[r-(1<<k)+1][k]);
}
int query_max(int l,int r)
{
int k=log2(r-l+1);
return max(max_[l][k],max_[r-(1<<k)+1][k]);
}
int main()
{
while(~scanf("%d%d", &N, &Q))
{
for(int i=1;i<=N;i++)
{
scanf("%d",&max_[i][0]);
min_[i][0]=max_[i][0];
}
RMQ_pre();
while(Q--)
{
int a,b;
scanf("%d%d", &a, &b);
printf("%d\n",query_max(a,b)-query_min(a,b));
}
}
return 0;
}