RMQ (模板)

RMQ是查找区间最小或者最大值的问题

一般有三种

• 1、暴力
• 2、线段树进行区间维护
• 3、就是主要说一下dp的解法

这种操作可以说是一种预处理，我们要求一个区间
[L , R ] 的最大值，就把这个区间从中间劈开，求两个区间的最大值然后比较，就这样一直分解，直到区间长度为1
那么也就是这个区间最大的是它本身。
太懒了，所以我只记录一下模板，具体的请看传送门

poj 3264

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <queue>
#include <vector>
#include <string>
#include <cmath>
#include <cstdlib>

using namespace std;

const int MAXN = 50000 + 10;
const int Size = 20;

int dp_max[MAXN][Size] , dp_min[MAXN][Size] , arr[MAXN];

//预处理，以2 ^k为单位
inline void RMQ(int n)
{
	for(int j = 1; 1<<j <= n; j++){
		for(int i = 1; i + (1<<j) - 1 <= n; i++){
			dp_max[i][j] = max(dp_max[i][j - 1] , dp_max[i + (1<<(j - 1))][j - 1]);
			dp_min[i][j] = min(dp_min[i][j - 1] , dp_min[i + (1<<(j - 1))][j - 1]);
		}
	}
}

inline int Query(int L , int R , int flag  = 0)
{
	//把区间【l , r】分成两个长度为2^k的区间（可能有部分重叠）公式
	int k = (int)(log( R - L + 1 ) * 1.0 / log( 2.0 ));
	if(flag){
		return max(dp_max[L][k] , dp_max[R - (1<<k) + 1][k]);
	}
	return  min(dp_min[L][k] , dp_min[R - (1<<k) + 1][k]);
}

int main()
{
	int m,n;
	scanf("%d%d",&n,&m);
	for(int i = 1; i <= n; i++){
		scanf("%d",arr + i);
		dp_max[i][0] = arr[i];
		dp_min[i][0] = arr[i];
	}
	RMQ(n);
	int l,r;
	while( m-- ){
		scanf("%d%d",&l,&r);
		int ans_max = Query( l , r , 1);
		int ans_min = Query( l , r );
		cout<<ans_max - ans_min<<endl;
	}
	return 0;
}

poj 3368

• 题意：

给一个单调不递减的数列，求某一区间内的连续重复最多数字的次数。
与模板不同的是我们要开一个num数组，去记录以这个位置的点为尾时，这个数字连续出现了几次，我们要找给定区间内num数组的最大值，但是要注意，左边界要进行特判
因为会出现这种情况
-1 ，1 ，-1 ，1， 1 【1， 1， 10， 10， 10】
我要求括号范围内重复次数最多的数字，如果我直接计算这个区间内的num的最大值，会把左侧不属于这个区间的1算进去，毕竟以区间内第二个结尾的重复1确实是最大的，但是并不是这个区间内重复次数最多的，所以对左边界进行特判，从这个区间内第一个不等于左边界的数到右区间内查找num最大值再和区间内前面和左边界连续且重复的数的数量进行比较
注意：这里可能会出现整个区间都与左边界相同，这样在求RMQ时会越出边界，需要特判一下，如果L > R就返回0

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <queue>
#include <vector>
#include <string>
#include <cmath>
#include <cstdlib>

using namespace std;

const int INF  = 0x3f3f3f3f;
const int MAXN = 100100;
const int Size = 20;

int dp[MAXN][Size] ,num[MAXN] , arr[MAXN];

inline void RMQ(int n)
{
	for(int i = 1; i <= n; i++){
		dp[i][0] = num[i];
	} 
	for(int j = 1; 1<<j <= n; j++){
		for(int i = 1; i + (1<<j) - 1 <= n; i++){
			dp[i][j] = max(dp[i][j - 1] , dp[i + (1<<(j - 1))][j - 1]);
		}
	}
}

inline int Query(int L , int R)
{
	//这里要注意！
	if(L > R){
		return 0;
	}
	int k = (int)(log( R - L + 1 ) * 1.0 / log( 2.0 ));
	return  max(dp[L][k] , dp[R - (1<<k) + 1][k]);
}

int main()
{
	int n;
	arr[0] = INF;
	while(scanf("%d",&n),n){
		int m;
		scanf("%d",&m);
		//memset(dp, 0, sizeof dp);
		memset(num,  0, sizeof num);
		for(int i  = 1; i <= n; i++){
			scanf("%d",arr + i);
			if(arr[i] == arr[i - 1]){
				num[i] = num[i - 1] + 1;
				continue;
			}
			num[i] = 1;
		}
		RMQ(n);
		while( m-- ){
			int l , r;
			scanf("%d%d",&l,&r);
			int k = l;
			//左起第一个数的次数不能RMQ查询，因为查询区间之外可能有这个数
			while(k <= r && arr[l] == arr[k]){
				k++;
			}
			int ans = max(k - l , Query( k , r ));
			printf("%d\n",ans);
		}
	}
	return 0;
}