(2888)HDU-RMQ

最新推荐文章于 2019-04-22 11:29:27 发布

原创最新推荐文章于 2019-04-22 11:29:27 发布 · 386 阅读

0 ·

CC 4.0 BY-SA版权

ACM-ICPC 专栏收录该内容

166 篇文章

订阅专栏

本文深入探讨了二维RMQ（Range Minimum Query）模板题的解决策略，包括初始化和查询过程，通过实例演示如何高效地解决此类问题。

#include<iostream>
#include<cstdio>
#include<string.h>
#include<string>
#include<stack>
#include<set>
#include<algorithm>
#include<cmath>
#include<vector>
#include<map>

#define ll __int64
#define lll unsigned long long
#define MAX 1000009
#define eps 1e-8

using namespace std;
/*
二维RMQ模板题
同一维一样用dp[row][col][i][j]表示(row,col)到(row+2^i,col+2^j)矩形内的最大值
查询
*/

int mapp[309][309];
int dp[309][309][9][9];
int flag;

void RMQ_init2d(int m,int n)
{
for(int i=1; i<=m; i++)
{
for(int j = 1; j<=n; j++)
{
dp[i][j][0][0] = mapp[i][j];
}
}
int t = log((double)n) / log(2.0);

for(int i = 0; i<=t; i++)
{
for(int j = 0; j<=t; j++)
{
if(i==0&&j==0)
continue;
for(int row = 1; row+(1<<i)-1<= m; row++)
{
for(int col = 1; col+(1<<j)-1<= n; col++)
{
if(i)
dp[row][col][i][j] = max(dp[row][col][i-1][j],dp[row+(1<<(i-1))][col][i-1][j]);
else
dp[row][col][i][j] = max(dp[row][col][i][j-1],dp[row][col+(1<<(j-1))][i][j-1]);
}
}
}
}
}
int RMQ_2d(int x1,int y1,int x2,int y2)
{
int k1 = log(double(x2 - x1 + 1)) / log(2.0);
int k2 = log(double(y2 - y1 + 1)) / log(2.0);
int m1 = dp[x1][y1][k1][k2];
int m2 = dp[x2 - (1<<k1) + 1][y1][k1][k2];
int m3 = dp[x1][y2 - (1<<k2) + 1][k1][k2];
int m4 = dp[x2 - (1<<k1) + 1][y2 - (1<<k2) + 1 ][k1][k2];
int _max = max(max(m1,m2),max(m3,m4));
if(mapp[x1][y1]==_max||mapp[x1][y2]==_max||mapp[x2][y1]==_max||mapp[x2][y2]==_max)
flag = 1;
return _max;
}

int main()
{
int n,m,t;
int x1,x2,y1,y2;
while(~scanf("%d%d",&m,&n))
{
for(int i = 1; i<=m; i++)
{
for(int j = 1; j<=n; j++)
{
scanf("%d",&mapp[i][j]);
}
}
RMQ_init2d(m,n);
scanf("%d",&t);
while(t--)
{
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
flag = 0;
int _max = RMQ_2d(x1,y1,x2,y2);
if(flag == 1)
printf("%d yes\n",_max);
else
printf("%d no\n",_max);
}
}
return 0;
}