RMQ问题的解析和高效率算法--uva11235

RMQ--范围最小值问题 ， 给出一个n个元素的数组A1.......An ， 设计一个数据结构，支持查询Query()： min{Al ......Ar}

对于这个问题，如果我们每次用一个循环来查找，那就太慢了 ， 即使我们用前缀和的思想我是不能提高效率的。

我们另d(i , j)表示从 i 开始的， 长度为 2^j ， 的一段元素的最小值 ， 则可以通过递推的方法计算 d(i , j)：d(i , j) = min{d(i , j-1) , d(i+2^(j-1) , j-1)} ， 注意 2^j <= n；

则因此d数组中的元素不会超过 nlog(n) ， 而每一项都可以在常数时间计算完 ， 固时间效率为O(nlong(n)).

初始化操作代码：

void rmq_init()
{
	int i , j;
	for(i = 1; i <= n; i++) d[i][0] = A[i];
	for(i = 1; (1<<i) <= n; i++)
	{
		for(j = 1; j+(1<<i)-1 <= n; j++)
		{
			d[j][i] = min(d[j][i-1] , d[j+(1<<(j-1))][j-1]);
		}
	}
}

查询：

查询操作很简单 ， 假设查询区间为：( r  , l ) ， 我们令 k 满足：2^k <= r - l + 1 ， 那么这时我们只需要查询 ， 以 l 开头 和以 r 结尾的长度为 2^k 的区间 ，  由于我们是找最小值 ， 因此有元素重复是没问题的。

代码：

void rmq_find(int l , int r)
{
	int k = 0;
	while(1<<(k+1) <= r-l+1) k++;
	return min(d[l][k] , d[r-(1<<k)+1][k]);
}

练习题目 uva 11235：

这个题目关键在于：这个数组是非降序的 ， 所以我们在一开始就能把所有相等的元素聚集在一起，这样就可以把整个数组进行游程编码，这样我们只需要用3个数组 ， 记录：每个段的（相等的元素为一个段）开始处、结尾处、元素的数量。

这样我们就能把每个段当成一个元素 ， 然后用RMQ算法来初始化 ， 只不过这时记录的是最大值 ， 而不是最小值。

我们每次查询（L , R）分3 部分结果 的最大值，从L 到L 所在段结尾的处的元素 ， 从 R开始段到R的元素，处于L , R之间完整段的最大值。

代码：

#include <iostream>
#include <stdio.h>
#include <string.h>
using namespace std;

#define max(x , y)  (x)>(y)?(x):(y)

int a[100010] , b[100010][20] , c[100010];
int e[100010] , f[100010];
int n;

void rmq_init(int k)
{
	int i , j;
	for(i = 1; (1<<i)<=k; i++)
	{
		for(j = 1; j+(1<<i)-1<=k; j++)
		{
			b[j][i] = max(b[j+(1<<(i-1))][i-1] , b[j][i-1]);
		}
	}

}

int rmq(int l , int r)
{
	int k = 0;
	while((1<<(k+1)) <= r-l+1)  k++;
	return max(b[l][k] , b[r-(1<<k)+1][k]);
}

int main()
{
	int q;
	while(scanf("%d" , &n) && n)
	{
		scanf("%d" , &q);
		memset(b , 0 , sizeof(b));
		int i , j;
		for(i = 1; i <= n; i++)
		{
			scanf("%d" , &a[i]);
		}
		int k = 1;
		e[k] = 1;
		c[1] = 1;
		for(i = 2; i <= n; i++)
		{
			c[i] = k;
			if(a[i] != a[i-1])
			{
				b[k][0] = i-e[k];
				f[k] = i-1;
				k += 1;
				e[k] = i;
				c[i] = k;
			}
		}
		b[k][0] = i-e[k];
		f[k] = i-1;
		rmq_init(k);
		int x , y ;
		int g , h;
		for(i = 1; i <= q; i++)
		{
			scanf("%d %d" , &x , &y);
			g = c[x] , h = c[y];
			if(g == h)  printf("%d\n" , y-x+1);
			else if((h-g) == 1)
				printf("%d\n" , max(f[g]-x+1 , y-e[h]+1));
			else
			{
				k = max(f[g]-x+1 , rmq(g+1 , h-1));
	
				printf("%d\n" , max(k , y-e[h]+1));
			}
		}
	}
	return 0;
}