RMQ (Range Minimum/Maximum Query)

    RMQ (Range Minimum/Maximum Query)问题是指：对于长度为n的数列A，回答若干询问RMQ(A,i,j)(i,j<=n)，返回数列A中下标在i,j里的最小(大）值，也就是说，RMQ问题是指求区间最值的问题。

①朴素做法：我们遍历区间[l,r]，找出最大值/最小值，一次遍历复杂度是O(n)
如果n次查询，最坏的复杂度能到达O(n*n)，这是我们无法接受的！

②对于区间最大值/最小值，我们能使用很多优秀的算法，比如线段树/树状数组。一次查询的时间复杂度在logN，n次查询，复杂度是nlogn，这个我们能接受！

③今天的主角不是线段树，也不是树状数组。是一个打表的做法，很多人称它是ST表，它的预处理复杂度是n*logn，惊人的是查询复杂度是O（1）

     下面将详细的为你讲解，如何求出ST表，以及如何查询ST表。

①首先求ST表是一个动态规划(dp)的思想，f[i][j]表示从i位置开始，后面连续2^j个数字的最大值/最小值。下面所有例子为了方便，我们以最小值为例子，最大值也同理的。

②动态规划的问题是将一个问题划分成若干的子问题，然后依次解决！

③由于上面f[2][2]区间长度为4，我们是不是可以把区间长度分成2和2，这样规模是不是减小了！

④推动态方程！！！！把特殊的一般化！！！

⑤动态方程有了，可以实现代码了！！！ 首先我们先确定下j的范围，由于j表示的是2^j个连续数字。所以，如果有一个长度为n的数组，j的范围[0,log(n)/log(2)] 。f[i][0]也是可以预处理的，0表示2^0=1个。

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <string>
#include <vector>
#include <map>
#include <queue>
#include <deque>
#include <set>
#include <stack>
#include <cctype>
#include <cmath>
#include <cassert>

using namespace std;

typedef long long ll;

const ll mod=1000000007;
ll powmod(ll a,ll b) {ll res=1;a%=mod; assert(b>=0); for(;b;b>>=1){if(b&1)res=res*a%mod;a=a*a%mod;}return res;}
ll gcd(ll a,ll b) { return b?gcd(b,a%b):a;}

int f[100][100];//具体问题具体分析数组大小,这里的100,100毫无意义，为了方便处理 
int a[]={0,5,4,1,2,7,8,6};//下面代码下标是1开始，这里第一个0只是站位,无实际意义 

void st()//计算st表 
{
	int n=7;
	
	for(int i=1;i<=n;i++)f[i][0]=a[i];//j=0就是当前位置的a[i]，这个我们先处理了
	
	int m=(int)(log((double)n)/log(2.0)); //计算出j最大范围
	
	for(int j=1;j<=m;j++)//外层循环j 
	{//至于为什么要j外层循环，i内层循环。是因为我们递推的过程,比如f[2][2]=min(f[2][1],f[4][1]),要知道j=2的f数组,必须知道j=1的数组。 
		for(int i=1;i+(1<<j)-1<=n;i++)//内层循环i
		{//i到i+(1<<j)-1是整个区间,那么i+(1<<j)-1是右边界,当然不能超过n
 			f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]);//转移方程 
		}
	} 
	
	//打印st表
	printf("a[i]:");
	for(int i=1;i<=n;i++)
	{
		printf("%d%c",a[i]," \n"[i==n]);
	}
	printf("-------------------------------\n");
	for(int j=1;j<=m;j++)
	{
		for(int i=1;i+(1<<j)-1<=n;i++)
		{
			printf("从第%d位开始连续%d个数字的最小值是%d\n",i,(1<<j),f[i][j]); 
		}
	} 
} 

int main()
{
	st();
	return 0;
}

结果如下：

⑥st表已经出来了，接下来就是最后一步，如何O（1）的查询了！
比如我们要查询区间[1,5]的最小值
首先我们要算出j能覆盖的范围

k=（int）(log((double)(r-l+1))/log(2.0))

区间【1,5】-> k=2

f[1][2]只覆盖了4个，还差一个5没覆盖，我们还需要一个f值

最后的公式就是 ： min(f[1][2],f[2][2])这样就把整个区间覆盖了！！
min(f[l][k], f[r-(1<<k)+1][k]);

查询部分代码：

int rmp_find(int l, int r)  //求区间 l 到 r 之间的最值
{
    int k = (int)(log(1.0 * (r-l+1)) / log(2.0));
    return min(f[l][k], f[r-(1<<k)+1][k]);
}

最后我们来看看查询所有区间的结果：

愿，自己不再辜负自己，奇迹不再辜负奇迹！