二维RMQ
两种做法:
1) O(n*m*log(m))-O(n)
2) O(n*m*log(n)*log(m)-O(1)
第一种方法:
把每一行都当成一维RMQ处理
代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <set>
#include <vector>
#include <map>
#include <queue>
#include <set>
#include <algorithm>
using namespace std;
#define MAXN 250
int dp[MAXN][MAXN][20];
int dp1[MAXN][MAXN][20];
int a[MAXN][MAXN];
int n,m;
//RMQ预处理;
void st()
{
for(int i=1;i<=n;i++)
{
for(int k=0;(1<<k)<=m;k++)
{
for(int j=1;j+(1<<k)-1<=m;j++)
{
if(k==0)//第i行的一维RMQ中,从第j个数开始,长度为1的范围中,最大值最小值都是他自己;相当于一维数组中的maxx[i][0] = minn[i][0] = x;
{
maxx[i][j][k] = minn[i][j][k] = mapp[i][j];
}
else //第i行的一维RMQ;
{
maxx[i][j][k] = max(maxx[i][j][k-1], maxx[i][j+(1<<(k-1))][k-1]);
minn[i][j][k] = min(minn[i][j][k-1], minn[i][j+(1<<(k-1))][k-1]);
}
}
}
}
}
int rmq2dmax(int x,int y,int x1,int y1)
{
int k = log2(y1-y+1);
int mm = max(maxx[x][y][k], maxx[x][y1-(1<<k)+1][k]);
for(int i=x+1;i<=x1;i++)
mm = max(mm,max(maxx[i][y][k],maxx[i][y1-(1<<k)+1][k]));
return mm;
}
int rmq2dmin(int x,int y,int x1,int y1)
{
int k = log2(y1-y+1);
int mm = min(minn[x][y][k], minn[x][y1-(1<<k)+1][k]);
for(int i=x+1;i<=x1;i++)
mm = min(mm,min(minn[i][y][k], minn[i][y1-(1<<k)+1][k]));
return mm;
}
第二种方法:
首先,dp[ i ][ j ][ k ][ l ]表示从第一个点也就是左下角(i,j)起,横坐标长度为2^k, 纵坐标长度为2^l,即右上角(i+2^k-1,j+2^k-1)中的最大值。
一维RMQ是把一个区间划分成两个区间求最值,同理,在二维RMQ中,我们可以把一个大的矩阵划分成四个小矩阵求最值
但是要考虑到一种特殊情况:这个矩阵行为1或列为1时,需要特判断。
一共四种情况:
一个点:k=0,l=0(从点(i, j)横纵坐标长度都为1,所包含的点的最大值,即只有它本身) maxx[ i ][ j ][ k ][ l ] = mapp[ i ][ j ];
一列: k=0,l!=0(从点(i, j)横坐标长度为1,送坐标长度>1所包含的点的最大值,即