【RMQ的总结】:
一般的区间最值问题,即一组数列,给出区间,问这个区间的最大值或最小值。
这个问题可以用线段树解决,但是如果只是查询最值,不涉及修改,RMQ的代码会更短小,容易调试。
这是一个离线算法,不支持值的更新。只能预处理后单纯的查询
RMQ通过预处理区间的最值,用O(n*logn)的时间把最值求完存数组,然后每次询问可以再O(1)的时间给出答案。
这样,这个问题的复杂度就降到了O(q + n*logn)
RMQ的实现。
设给出的数列存在数组a中
用一个二维数组dp,dp[i][j] 表示区间左端点为i,区间长度为2^j。
本算法的厉害之处就在第二维的j表示区间长度,并且是与2的倍数相关的数字,因此才能将复杂度降到log。
以最小值为例,求区间[ l , r ]的最小值:
当j=0时,毫无疑问 dp[i][0] = a[i]
当j>0时,区间长度为2^j ,即一个偶数。我们把这个区间平分,两段长为2^(j-1)的区间。
两段子区间为[ l , l+2^(j-1)-1] 和 [ l+2^(j-1) , r ]
那么当前区间的dp[i][j]就等于这两区间的最小值,把问题化为子问题,动态规划把dp预处理出来。
dp[i][j] = min { dp[l][j-1] , dp[l+2^(j-1)][j-1] }; //状态更新方程
由此得到dp数组,接下来询问区间最值时,只需查看dp数组。
关于查询。
虽然dp数组预处理了,但询问也有些技巧。例如询问区间[ l , r ]的最小值。
需要一个整数k,满足2^k < 区间长度。
如此一来,给定的区间又能分成两块(注意,此时分成两块,是可能有重叠部分的,因为给定的区间不一定偶数)
可以分成[ l , l+2^k]和 [ r-2^k+1 , r] (有重叠部分,但没有影响)
这两个区间在dp数组里是有的,所以直接求这两个的最小值即可
【例题】:NYOJ 119 士兵杀敌(三)
士兵杀敌(三)
-
描述
-
南将军统率着N个士兵,士兵分别编号为1~N,南将军经常爱拿某一段编号内杀敌数最高的人与杀敌数最低的人进行比较,计算出两个人的杀敌数差值,用这种方法一方面能鼓舞杀敌数高的人,另一方面也算是批评杀敌数低的人,起到了很好的效果。
所以,南将军经常问军师小工第i号士兵到第j号士兵中,杀敌数最高的人与杀敌数最低的人之间军功差值是多少。
现在,请你写一个程序,帮小工回答南将军每次的询问吧。
注意,南将军可能询问很多次。
-
输入
-
只有一组测试数据
第一行是两个整数N,Q,其中N表示士兵的总数。Q表示南将军询问的次数。(1<N<=100000,1<Q<=1000000)
随后的一行有N个整数Vi(0<=Vi<100000000),分别表示每个人的杀敌数。
再之后的Q行,每行有两个正正数m,n,表示南将军询问的是第m号士兵到第n号士兵。
输出
- 对于每次询问,输出第m号士兵到第n号士兵之间所有士兵杀敌数的最大值与最小值的差。 样例输入
-
5 2 1 2 6 9 3 1 2 2 4
样例输出
-
1 7
-
只有一组测试数据
【代码】:
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<string.h>
#include<iostream>
#include<algorithm>
using namespace std;
struct RMQ{
int Max[100010][20],Min[100010][20];
void init(int a[],int l,int r)
{
int top=log((r-l+1)*1.0)/log(2.0);
for(int i=l;i<=r;i++)
Min[i][0] = Max[i][0] = a[i];
for(int j=1;j<=top;j++)
for(int i=l;i+(1<<j)-1<=r;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 read(int l,int r)
{
int k=log(r*1.0-l+1)/log(2.0);
int u=max(Max[l][k],Max[r-(1<<k)+1][k]);
int v=min(Min[l][k],Min[r-(1<<k)+1][k]);
return u-v;
}
}R;
int a[100010];
int main()
{
int n,q,l,r;
scanf("%d%d",&n,&q);
for(int i=1;i<=n;i++)
scanf("%d",a+i);
R.init(a,1,n);
while(q--)
{
scanf("%d%d",&l,&r);
printf("%d\n",R.read(l,r));
}
}