二维RMQ 思想及模板

最新推荐文章于 2018-12-15 11:18:05 发布
最新推荐文章于 2018-12-15 11:18:05 发布
阅读量322 收藏
点赞数
分类专栏： RMQ 数据结构
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_40859951/article/details/83957716
版权

二维RMQ问题就是求一个矩阵N*M中的一个小块矩阵内的最值问题.其中dmin[i][j][ii][jj]=x表示以(i, j)为左上角,以(i+(1<<ii)-1, j+(1<<jj)-1 )为右下角的矩阵内的最小值.dmax的值类似.

         下面dmin[i][j][ii][jj]的值如何求呢?首先我们知道dmin[i][j][0][0]的值就是v[i][j],而假设dmin[i][j][ii][jj]中的ii不为0,那么dmin[i][j][ii][jj]= min(dmin[i][j][ii-1][jj], dmin[i+(1<<ii)][j][ii-1][jj] );如果ii为0,那么就按jj来求.

         其实上面的求法就是等于把二维问题转变为一维问题来求解.

         下面我们讨论如何查询结果.

         对于一个以(x1, y1)为左上角,以(x2, y2)为右下角的矩形,如何求它的最小值和最大值呢?下面假设我们求最小值:

         我们把(x1，y1)与(x2，y2)构成的矩形分成四小块，这四小块可能有重合部分，但是它们共同构成了目标矩形：

         dmin[x1][y1][ii][jj]

         dmin[x1][y2-(1<<jj)+1][ii][jj]

         dmin[x2-(1<<ii)+1][y1][ii][jj]

         dmin[x2-(1<<ii)+1][y2-(1<<jj)+1][ii][jj]

         (自己想象下上面4小块是怎么样的？)

         temp 1=min(dmin[x1][y1][ii][jj] , dmin[x1][y2-(1<<jj)+1][ii][jj])

         temp 2=min(dmin[x2-(1<<ii)+1][y1][ii][jj] ,dmin[x2-(1<<ii)+1][y2-(1<<jj)+1][ii][jj] )

         最终结果是min(temp1, temp2);

//POJ 2019
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAXN=300;
int val[MAXN][MAXN];
int dmin[MAXN][MAXN][10][10];
int dmax[MAXN][MAXN][10][10];
void initRMQ(int n,int m)//对n*m的矩阵初始化RMQ且矩阵下标从1开始
{
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            dmin[i][j][0][0]=dmax[i][j][0][0]=val[i][j];
    for(int ii=0;(1<<ii)<=n;ii++)
        for(int jj=0;(1<<jj)<=m;jj++)
        if(ii+jj)
            for(int i=1;i+(1<<ii)-1<=n;i++)
                for(int j=1;j+(1<<jj)-1<=m;j++)
                if(ii)
                {
                    dmin[i][j][ii][jj] = min(dmin[i][j][ii-1][jj] ,dmin[i+(1<<(ii-1))][j][ii-1][jj]);
                    dmax[i][j][ii][jj] = max(dmax[i][j][ii-1][jj] ,dmax[i+(1<<(ii-1))][j][ii-1][jj]);
                }
                else
                {
                    dmin[i][j][ii][jj] = min(dmin[i][j][ii][jj-1] , dmin[i][j+(1<<(jj-1))][ii][jj-1]);
                    dmax[i][j][ii][jj] = max(dmax[i][j][ii][jj-1] , dmax[i][j+(1<<(jj-1))][ii][jj-1]);
                }
}
int getMax(int x1,int y1,int x2,int y2)//RMQ查询
{
    int k1=0;
    while((1<<(k1+1))<=x2-x1+1)k1++;
    int k2=0;
    while((1<<(k2+1))<=y2-y1+1)k2++;
    x2 = x2 - (1<<k1)+1;
    y2 = y2 - (1<<k2)+1;
    return max(max(dmax[x1][y1][k1][k2],dmax[x1][y2][k1][k2]) ,max(dmax[x2][y1][k1][k2],dmax[x2][y2][k1][k2]) );
}
int getMin(int x1,int y1,int x2,int y2)//RMQ查询
{
    int k1=0;
    while((1<<(k1+1))<=x2-x1+1)k1++;
    int k2=0;
    while((1<<(k2+1))<=y2-y1+1)k2++;
    x2 = x2 - (1<<k1)+1;
    y2 = y2 - (1<<k2)+1;
    return min( min(dmin[x1][y1][k1][k2],dmin[x1][y2][k1][k2]) ,min(dmin[x2][y1][k1][k2],dmin[x2][y2][k1][k2]) );
}


例题：

题意：

        给你一个n*n的整数矩阵，然后有k个询问且给你一个常数b。对于每个询问(x,y)，你要输出以(x,y)为左上角，以(x+b-1，y+b-1)为左下角的矩形 内整数最大值-最小值的差。

分析:

        基本的二维RMQ查询.

        类似于二维树状数组问题,二维RMQ问题就是求一个矩阵N*M中的一个小块矩阵内的最值问题.其中dmin[i][j][ii][jj]=x表示以(i , j)为左上角,以( i+(1<<ii)-1, j+(1<<jj)-1 )为右下角的矩阵内的最小值.dmax的值类似.

       下面dmin[i][j][ii][jj]的值如何求呢?首先我们知道dmin[i][j][0][0]的值就是v[i][j],而假设dmin[i][j][ii][jj]中的ii不为0,那么dmin[i][j][ii][jj] =min(dmin[i][j][ii-1][jj], dmin[i+(1<<ii)][j][ii-1][jj] );如果ii为0,那么就按jj来求.

       其实上面的求法就是等于把二维问题转变为一维问题来求解.

       下面我们讨论如何查询结果.

        对于一个以(x1 , y1)为左上角,以(x2 , y2)为右下角的矩形,如何求它的最小值和最大值呢?下面假设我们求最小值:

        我们把(x1，y1)与(x2，y2)构成的矩形分成四小块，这四小块可能有重合部分，但是它们共同构成了目标矩形：

        dmin[x1][y1][ii][jj]

        dmin[x1][y2-(1<<jj)+1][ii][jj]

        dmin[x2-(1<<ii)+1][y1][ii][jj]

        dmin[x2-(1<<ii)+1][y2-(1<<jj)+1][ii][jj]

        (自己想象下上面4小块是怎么样的？)

        temp 1=min( dmin[x1][y1][ii][jj] , dmin[x1][y2-(1<<jj)+1][ii][jj])

        temp 2=min( dmin[x2-(1<<ii)+1][y1][ii][jj] ,dmin[x2-(1<<ii)+1][y2-(1<<jj)+1][ii][jj] )

        最终结果是min(temp1, temp2)；
代码：

<span style="font-size:18px;">//POJ 2019
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAXN=255;
int val[MAXN][MAXN];
int dmin[MAXN][MAXN][8][8];
int dmax[MAXN][MAXN][8][8];
void initRMQ(int n,int m)
{
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            dmin[i][j][0][0]=dmax[i][j][0][0]=val[i][j];
    for(int ii=0;(1<<ii)<=n;ii++)
        for(int jj=0;(1<<jj)<=m;jj++)
        if(ii+jj)
            for(int i=1;i+(1<<ii)-1<=n;i++)
                for(int j=1;j+(1<<jj)-1<=m;j++)
                if(ii)
                {
                    dmin[i][j][ii][jj] = min(dmin[i][j][ii-1][jj] ,dmin[i+(1<<(ii-1))][j][ii-1][jj]);
                    dmax[i][j][ii][jj] = max(dmax[i][j][ii-1][jj] ,dmax[i+(1<<(ii-1))][j][ii-1][jj]);
                }
                else
                {
                    dmin[i][j][ii][jj] = min(dmin[i][j][ii][jj-1] , dmin[i][j+(1<<(jj-1))][ii][jj-1]);
                    dmax[i][j][ii][jj] = max(dmax[i][j][ii][jj-1] , dmax[i][j+(1<<(jj-1))][ii][jj-1]);
                }
}
int getMax(int x1,int y1,int x2,int y2)
{
    int k1=0;
    while((1<<(k1+1))<=x2-x1+1)k1++;
    int k2=0;
    while((1<<(k2+1))<=y2-y1+1)k2++;
    x2 = x2 - (1<<k1)+1;
    y2 = y2 - (1<<k2)+1;
    return max(max(dmax[x1][y1][k1][k2],dmax[x1][y2][k1][k2]) ,max(dmax[x2][y1][k1][k2],dmax[x2][y2][k1][k2]) );
}
int getMin(int x1,int y1,int x2,int y2)
{
    int k1=0;
    while((1<<(k1+1))<=x2-x1+1)k1++;
    int k2=0;
    while((1<<(k2+1))<=y2-y1+1)k2++;
    x2 = x2 - (1<<k1)+1;
    y2 = y2 - (1<<k2)+1;
    return min( min(dmin[x1][y1][k1][k2],dmin[x1][y2][k1][k2]) ,min(dmin[x2][y1][k1][k2],dmin[x2][y2][k1][k2]) );
}
 
int main()
{
    int n,b,k;
    while(scanf("%d%d%d",&n,&b,&k)==3&&n&&b&&k)
    {
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                scanf("%d",&val[i][j]);
        initRMQ(n,n);
        while(k--)
        {
            int x1,y1;
            scanf("%d%d",&x1,&y1);
            int ans=getMax(x1,y1,x1+b-1,y1+b-1)-getMin(x1,y1,x1+b-1,y1+b-1);
            printf("%d\n",ans);
        }
    }
    return 0;
}
</span>

 

凌晨四点的洛杉矶fly
关注
专栏目录
二维RMQ
qq_42569807的博客
04-19 1351
二维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> #incl...
ACM比赛经验、刷题记录及模板库总结（更新中）
lordofadventure的博客
03-08 1911
文章目录前言第一部分 经验分享及感受第二部分 刷题记录一、基础算法&程序语言二、综合题（现场题）三、动态规划四、图论五、字符串算法六、数据结构七、数学八、网络流九、博弈十、随机算法十一、其他第三部分 专项题库第四部分 模板库附录A 题目索引附录B 常用平台参考文献 前言 本文所提及的部分题目代码，可以在https://github.com/volactina/ACM上找到 第一部分 经验分...
POJ 2019 Cornfields(简单二维RMQ)
code
03-24 1648
POJ2019 Cornfields(简单二维RMQ) 分析:基本的二维RMQ查询. 类似于二维树状数组问题,二维RMQ问题就是求一个矩阵N*M中的一个小块矩阵内的最值问题.其中dmin[i][j][ii][jj]=x表示以(i , j)为左上角,以( i+(1为右下角的矩阵内的最小值.dmax的值类似.        下面dmin[i][j][ii][jj]的值如何求呢?首先我们知道
HDU 2888 Check Corners(简单二维RMQ)
code
03-24 1079
HDU 2888 Check Corners(简单二维RMQ) 分析:基本的二维RMQ应用. 注意题目中存的矩阵是无符号整数的矩阵. AC代码:3625ms #include #include #include using namespace std; const int MAXN =305; unsigned int val[MAXN][MAXN]; unsigned int
二维RMQ问题
weixin_30794851的博客
11-03 252
前置知识 一维RMQ及其拓展 问题引入 题目地址IN 对于一个n×mn\times mn×m的矩阵，每个格子有一个值，有QQQ个询问，每次询问你一个子矩阵中的最大值。 1≤n,m≤500,Q≤1061\leq n,m\leq 500,Q\leq 10^61≤n,m≤500,Q≤106 暴力 每次花子矩阵大小的复杂度去查询。 复杂度最坏O(...
【模板】二维rmq
sadnohappy
10-14 501
二维rmq模板
数据结构：一维与二维RMQ详解
"WarBlood 模板.md - ACM竞赛相关数据结构知识：一维和二维RMQ（Range Maximum Query）" 本文档介绍了在ACM竞赛中常见的数据结构问题，特别是关于RMQ（Range Maximum Query，范围最大值查询）的实现。RMQ允许我们在...
ACM 算法 常用模板
11-23
#### RMQ-ST 一维及二维算法 **知识点解析：** 区间查询问题在算法竞赛中也非常常见，特别是求解区间最小值或最大值的问题。Sparse Table（稀疏表）算法提供了一种高效的解决方案，能够在线性时间内预处理数据，在...
ACM编程题模板和各种经典算法数据结构实现代码
03-06
- 二维树状数组可以应用于二维平面上的区间查询问题。 4. **TRIE树(K叉)** - TRIE树是一种用于存储字符串集合的树形数据结构，可以用于前缀搜索。 5. **TRIE树(左儿子右兄弟)** - 左儿子右兄弟表示法是另一种...
吉林大学ACM算法模板全面指南
- **高级数据结构**：左偏树合并（时间复杂度O(LOGN)）、树状数组、二维树状数组、Trie树（K叉和特殊结构）、后缀数组的两种时间复杂度实现、离线范围查询(RMQ)算法。 这些知识点旨在帮助学生掌握基本的算法思想、...
二维RMQ模板
maze_illusion的博客
08-30 250
int val[maxn][maxn],dp[maxn][maxn][9][9],mm[maxn]; void Rmq_init( int n , int m ) { mm[0] = -1; for ( int i=1 ; i&lt;maxn ; i++ ) mm[i] = ((i&amp;(i-1))==0)?mm[i-1]+1:mm[i-1]; for...
二维rmq
huyuncong的专栏
05-24 890
rmq推广
Monitor（二维RMQ）
FIVE
12-15 246
Recently Luba bought a monitor. Monitor is a rectangular matrix of size n × m. But then she started to notice that some pixels cease to work properly. Luba thinks that the monitor will become broken t...
二维RMQ(Range Minimum Query)
未完成
08-24 1057
二维RMQ(Range Minimum Query)先来看一维情形：给定一个数组x以及大量的查询，查询内容是任意区间内的最小值，例如数组长度为n，则查询为从x[a]到x[b]中的最小值，0问题是如何进行预处理，使得每次查询时间为O(1)。解法：计算出所有长度为2^i的区间的最小值，则任[a,b]区间可以拆为两个长度为2^i的区间的并集，i=log(b-a)。例如查询1到100间的最小值，[
一维RMQ和二维RMQ模板以及用法
inf_force
07-23 3086
RMQ可以用来查找区间最值（最大或最小）
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值