关于二维矩阵的最大最小值的询问

最新推荐文章于 2024-05-17 01:35:36 发布
最新推荐文章于 2024-05-17 01:35:36 发布
阅读量924 收藏 1
点赞数 1
分类专栏： 个人总结 数据结构相关
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/BroDrinkWater/article/details/78334461
版权

今天考了这样子的一道题目，考试的时候yy出了二维RMQ的做法，调了半天，考完后，想了想二维线段树的做法,感觉很好想的，只是把一维线段树改一改。

二维RMQ

处理出 dp[i][j][l][r] 表示i往前 2l ，j往前 2r 的矩阵的最大值和最小值，首先是 dp[i][j][0][0] =原数值，然后每次枚举l，r然后枚举i，j就好了。像一维前缀和那样。询问的时候变成4个矩形就好了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>

using namespace std;

typedef long long LL;

#define REP(i,a,b) for(register int i = (a),i##_end_ = (b); i <= i##_end_; ++i)
#define DREP(i,a,b) for(register int i = (a),i##_end_ = (b); i >= i##_end_; --i)
#define mem(a,b) memset(a,b,sizeof(a))

template<typename T> inline bool chkmin(T &a,const T &b) { return a > b ? a = b , 1 : 0;}
template<typename T> inline bool chkmax(T &a,const T &b) { return a < b ? a = b , 1 : 0;}

int read()
{
    register int f = 1,s = 0;char c = getchar();
    while(!isdigit(c)) { if(c == '-')f = -1;c = getchar();}
    while(isdigit(c)) { s = s * 10 + c - '0';c = getchar();}
    return f * s;
}

const int maxn = 510;
int n,m;
int a[maxn][maxn],dp[maxn][maxn][10][10],s[maxn][maxn];
int f[maxn][maxn][10][10];
char p[20];

void init()
{
    for (int i = 0; (1 << i) <= n; i++)
    {
        for (int j = 0; (1 << j) <= m; j++)
        {
            if (i == 0 && j == 0) continue;
            for (int row = 1; row + (1 << i) - 1 <= n; row++)
                for (int col = 1; col + (1 << j) - 1 <= m; col++)
                {
                    if (i == 0)
         dp[row][col][i][j] = max(dp[row][col][i][j - 1], dp[row][col + (1 << (j - 1))][i][j - 1]);
                    else if (j == 0) 
         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]);
                }
        }
    }
    for (int i = 0; (1 << i) <= n; i++)
    {
        for (int j = 0; (1 << j) <= m; j++)
        {
            if (i == 0 && j == 0) continue;
            for (int row = 1; row + (1 << i) - 1 <= n; row++)
                for (int col = 1; col + (1 << j) - 1 <= m; col++)
                {
                    if (i == 0)
         f[row][col][i][j] = min(f[row][col][i][j - 1], f[row][col + (1 << (j - 1))][i][j - 1]);
                    else if (j == 0) 
         f[row][col][i][j] = min(f[row][col][i - 1][j], f[row + (1 << (i - 1))][col][i - 1][j]);
                    else
         f[row][col][i][j] = min(f[row][col][i][j - 1], f[row][col + (1 << (j - 1))][i][j - 1]);
                }
        }
    }

}

int Query_Max(int x1, int y1, int x2, int y2) {
    int kx = 0, ky = 0;
    while (x1 + (1 << (1 + kx)) - 1 <= x2) kx++;
    while (y1 + (1 << (1 + ky)) - 1 <= y2) ky++;
    int m1 = dp[x1][y1][kx][ky];
    int m2 = dp[x2 - (1 << kx) + 1][y1][kx][ky];
    int m3 = dp[x1][y2 - (1 << ky) + 1][kx][ky];
    int m4 = dp[x2 - (1 << kx) + 1][y2 - (1 << ky) + 1][kx][ky];

    return max(max(m1, m2), max(m3, m4));
}
int Query_Min(int x1, int y1, int x2, int y2) {
    int kx = 0, ky = 0;
    while (x1 + (1 << (1 + kx)) - 1 <= x2) kx++;
    while (y1 + (1 << (1 + ky)) - 1 <= y2) ky++;
    int m1 = f[x1][y1][kx][ky];
    int m2 = f[x2 - (1 << kx) + 1][y1][kx][ky];
    int m3 = f[x1][y2 - (1 << ky) + 1][kx][ky];
    int m4 = f[x2 - (1 << kx) + 1][y2 - (1 << ky) + 1][kx][ky];

    return min(min(m1, m2), min(m3, m4));
}

int Sum(int x1,int y1,int x2,int y2)
{
    return s[x2][y2]-s[x2][y1-1]-s[x1-1][y2]+s[x1-1][y1-1];
}

int main()
{
    freopen("phalanx.in","r",stdin);
    freopen("phalanx.out","w",stdout);
    n = read(),m = read();
    REP(i,1,n)REP(j,1,m)
    {
        a[i][j] = read();
        dp[i][j][0][0] = f[i][j][0][0] = a[i][j];
        s[i][j] = s[i-1][j] + s[i][j-1] - s[i-1][j-1] + a[i][j];
    }
    init();
    int T = read();
    while(T--)
    {
        cin>>p;
        int x1 = read()+1,y1 = read()+1,x2 = read()+1,y2 = read()+1;
        if(p[0] == 'S')printf("%d\n",Sum(x1,y1,x2,y2));
        else if(p[1] == 'A')printf("%d\n",Query_Max(x1,y1,x2,y2));
        else printf("%d\n",Query_Min(x1,y1,x2,y2));
    }
    return 0;
}

二维线段树

其实跟一维线段树差不多，我们首先在每一行维护一个线段树，然后再在列上维护一个线段树，然后每次合并就好了。询问的时候先查对应的行，然后就是，递归找。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>

using namespace std;

typedef long long LL;

#define REP(i,a,b) for(register int i = (a),i##_end_ = (b); i <= i##_end_; ++i)
#define DREP(i,a,b) for(register int i = (a),i##_end_ = (b); i >= i##_end_; --i)
#define mem(a,b) memset(a,b,sizeof(a))

template<typename T> inline bool chkmin(T &a,const T &b) { return a > b ? a = b , 1 : 0;}
template<typename T> inline bool chkmax(T &a,const T &b) { return a < b ? a = b , 1 : 0;}

int read()
{
    register int f = 1,s = 0;char c = getchar();
    while(!isdigit(c)) { if(c == '-')f = -1;c = getchar();}
    while(isdigit(c)) { s = s * 10 + c - '0';c = getchar();}
    return f * s;
}

const int maxn = 810;
int n,m;
int a[maxn][maxn];
int s[maxn][maxn];

int Max[maxn<<2][maxn<<2],Min[maxn<<2][maxn<<2];
//[]第一维表示第一颗线段树，[]第二维表示第二棵线段树的根节点编号。

void Build(int H,int row,int h,int l,int r)
{
//递归建第二维。
    if(l == r)
    {
        Max[H][h] = a[row][l];
        Min[H][h] = a[row][l];
        return ;
    }
    int mid = (l + r) >> 1;
    Build(H,row,h<<1,l,mid);Build(H,row,h<<1|1,mid+1,r);
    Max[H][h] = max(Max[H][h<<1],Max[H][h<<1|1]);
    Min[H][h] = min(Min[H][h<<1],Min[H][h<<1|1]);
}

void merge(int H,int h1,int h2,int h,int l,int r)
{
//合并两个儿子，把第二维合并。
    if(l == r)
    {
        Max[H][h] = max(Max[h1][h],Max[h2][h]);
        Min[H][h] = min(Min[h1][h],Min[h2][h]);
        return ;
    }
    int mid = (l + r) >> 1;
    merge(H,h1,h2,h<<1,l,mid);
    merge(H,h1,h2,h<<1|1,mid+1,r);
    Max[H][h] = max(Max[h1][h],Max[h2][h]);
    Min[H][h] = min(Min[h1][h],Min[h2][h]);
}

void build(int h,int l,int r)
{
//普通线段树建树，到叶子节点的时候再递归建树。回溯的时候，merge两个儿子。
    if(l == r)
    {
        Build(h,l,1,1,m);
        return ;
    }
    int mid = (l + r) >> 1;
    build(h<<1,l,mid);build(h<<1|1,mid+1,r);
    merge(h,h<<1,h<<1|1,1,1,m);
}


//递归询问.
int query_max(int H,int h,int l,int r,int q,int w)
{
    if(l == q && r == w)return Max[H][h];
    int mid = (l + r) >> 1;
    if(w <= mid)return query_max(H,h<<1,l,mid,q,w);
    else if(q > mid)return query_max(H,h<<1|1,mid+1,r,q,w);
    else return max(query_max(H,h<<1,l,mid,q,mid),query_max(H,h<<1|1,mid+1,r,mid+1,w));
}

int Query_max(int h,int l,int r,int q,int w,int ll,int rr)
{
    if(l == q && r == w)return query_max(h,1,1,m,ll,rr);
    int mid = (l + r) >> 1;
    if(w <= mid)return Query_max(h<<1,l,mid,q,w,ll,rr);
    else if(q > mid)return Query_max(h<<1|1,mid+1,r,q,w,ll,rr);
    else return max(Query_max(h<<1,l,mid,q,mid,ll,rr),Query_max(h<<1|1,mid+1,r,mid+1,w,ll,rr));
}

int query_min(int H,int h,int l,int r,int q,int w)
{
    if(l == q && r == w)return Min[H][h];
    int mid = (l + r) >> 1;
    if(w <= mid)return query_min(H,h<<1,l,mid,q,w);
    else if(q > mid)return query_min(H,h<<1|1,mid+1,r,q,w);
    else return min(query_min(H,h<<1,l,mid,q,mid),query_min(H,h<<1|1,mid+1,r,mid+1,w));
}

int Query_min(int h,int l,int r,int q,int w,int ll,int rr)
{
    if(l == q && r == w)return query_min(h,1,1,m,ll,rr);
    int mid = (l + r) >> 1;
    if(w <= mid)return Query_min(h<<1,l,mid,q,w,ll,rr);
    else if(q > mid)return Query_min(h<<1|1,mid+1,r,q,w,ll,rr);
    else return min(Query_min(h<<1,l,mid,q,mid,ll,rr),Query_min(h<<1|1,mid+1,r,mid+1,w,ll,rr));
}

int Q;
char ss[20];

int Sum(int x1,int y1,int x2,int y2)
{
    return s[x2][y2] - s[x2][y1-1] - s[x1-1][y2] + s[x1-1][y1-1];
}

int main()
{
    freopen("phalanx.in","r",stdin);
    freopen("phalanx.out","w",stdout);
    n = read(),m = read();
    REP(i,1,n) REP(j,1,m)
    {
        a[i][j] = read();
        s[i][j] = s[i-1][j] + s[i][j-1] - s[i-1][j-1] + a[i][j];
    }
    build(1,1,n);   
    Q = read();
    while(Q--)
    {
        cin>>ss;
        int x1 = read()+1,y1 = read()+1,x2 = read()+1,y2 = read()+1;
        if(ss[0] == 'S')printf("%d\n",Sum(x1,y1,x2,y2));
        else if(ss[1] == 'A')printf("%d\n",Query_max(1,1,n,x1,x2,y1,y2));
        else printf("%d\n",Query_min(1,1,n,x1,x2,y1,y2));
    }
    return 0;
}
BroDrinkWater
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
Python使用min、max函数查找二维数据矩阵中最小、最大值的方法
09-20
主要介绍了Python使用min、max函数查找二维数据矩阵中最小、最大值的方法,涉及Python矩阵生成、遍历、查找等相关操作技巧,需要的朋友可以参考下
二维RMQ（查找区域最大值）
FIVE
12-15 483
int dp[505][505][10][10]; int val[505][505]; int num[505]; void rmq(int n, int m){ int nn = log2(n), mm = log2(m); for(int i=1; i&lt;=n; i++) for(int j=1; j&lt;=m; j++) dp[i][j][0][0] =...
二维线段树——区间最值系列
mumei的博客
01-27 2251
最近这几天学习了一下二维线段树，二维线段树主要有两种写法，四分树和树套树，暂时还没写过四分树，因为这个东西确实不常用，而且不好写也不好调。 树套树的写法思路其实不难，首先我们知道，我们通常的线段树的操作都是在线段上进行的，即一维的，当推广到二维上，即矩阵时，我们就把它叫做二维线段树。 二维线段树其实就是我们再操作时先找到对应x轴的区间（即第一维），之后再找到对应y轴的区间（即第二维），进行相...
C语言经典算法实例1：求二维数组最大最小值_c二维数组最小值
最新发布
2401_85012088的博客
05-17 283
定义宏和声明数组 代码如下所示定义了MAXN ，代表了ＭＡＸＮ为常数20。声明了数组a。
二维倍增RMQ求区间最大值
CindyMarshall的博客
06-18 929
n,m&lt;=300,询问数&lt;=1e6给定n,m,一个n,m的矩阵，然后q组询问，每次询问给定一个子矩阵的左上坐标，右下坐标，求这个子矩阵的最大值。蒟蒻先RE60，后发现边界处理有误，改完才AC#include&lt;cstdio&gt; #include&lt;iostream&gt; using namespace std; int n,m,a[400][400]; int f[400...
hdu4819 正宗二维线段树 区间求最值，单点跟新
旺旺的博客
08-20 543
题意：给你一个矩阵，随后m行，x，y，l，从x行y列，扩大l/2个格子，不能出界，表示要查询的矩阵大小，查询最小值与最大值，然后将他们的和除 2输出来，然后将他们的和除2更新给 x，y。 思路：就是构造二维线段树，然后就是单点跟新，区间最值了。 /* 题意：给你一个矩阵，随后m行，x，y，l，从x行y列，扩大l/2个格子，不能出界，表示要查询的矩阵大小， 查询最小值与最大值，然后
二维矩阵求极小值
09-02 1090
问题描述 给定一个不包含重复元素的N行M列二维矩阵，求矩阵的极小值。极小值的定义如下：若一个值是极小值，四连通方向上的的不越界的邻居都比它大。 思路 二维比较难想，我们从一维说起。 对于一维数组来说，如何寻找极小值？ 二分法，找到中间位置x，如果x-1处的元素比x+1处的元素小，那么0~x-1之间必然存在一个极小值，从而缩小了范围。 升到二维，这个问题等价于。在平原上，寻找一个坑。这其实...
java 求解二维数组列最小值
08-30
在Java编程中，处理二维数组是一项常见的任务，尤其是在数据处理和矩阵运算的场景。本问题主要探讨如何在Java中找到二维数组每一列的最小值，并将这些最小值组成一个新的一维数组。这个问题的关键在于理解二维数组的...
Python获取二维矩阵每列最大值的方法
12-23
因为做项目中间有一个很小的环节需要这个功能，... 得到矩阵中每一列最大的值 ''' res_list=[] for j in range(len(martix[0])): one_list=[] for i in range(len(martix)): one_list.append(int(martix[i][j]))
基于数组指针实现二维数组中最小值所在行的查找与显示程序
11-16
用于比较指针数组和数组指针实现算法的差异，利用函数返回指针实现。 可作为研究C语言指针使用的范例程序。 有源代码和实现结果图。
第03章 方法与数组 06 二维数组与最大最小值算法
09-13
本章节“第03章 方法与数组 06 二维数组与最大最小值算法”将深入探讨这两个核心概念。 首先，我们来看二维数组的定义和创建。在Java中，二维数组可以用以下方式声明： ```java int[][] array = new int[行数][列...
矩阵（二维数组）局部极大/小值-python实现
银河系搬运工的博客
07-23 1560
题目来源： 某为面试/算法第四版：Algs4-1.4.19矩阵的局部最小元素 参考思路：传送 CODE： import numpy as np ''' def findMin(): arr=np.random.rand(10,10) index_arr=np.zeros((10,10)) for i in range(arr.shape[0]): for j in range(arr.shape[1]): if i>0 and i&
渐近思想——寻找局部最小值
weixin_43012409的博客
05-30 356
渐近思想——寻找局部最小值 数组的局部最小元素 给定一个长度为N的数组a，找到一个局部最小元素：满足a[i-1]>a[i]<a[i+1]的索引i；要求时间复杂度为O(lgN)。 前置条件：数组中不含有重复的元素。 思路 ​ 首先，要注意到如果没有前置条件，那么就不可能在O(lgN)的限制内解决该问题。因为最起码需要将数组遍历一遍来确定是不是真的有可能存在局部最小值（比如全是1的情况下就不存在）。 ​ 接着，我们将数组拓展为函数，那么对于某个点A来说，存在另一个点B比它的值小，那么沿着AB方向
二维固定区间最值
2301_79199219的博客
03-11 181
然后我们得到了max_m和min_m,接下来只要知道i - m1 + 1到 i行的max_m和min_m的最值,相当于把各个max_m和min_m拼接到一起了。朴素的想法是直接暴力i到i+n1-1,一般的题也够用,如果要优化的话,就再来一次单调队列.解法有,单调队列,优先队列,滑动窗口.(优先队速度稍微慢点,多两logn)两次单调队列写法,只要把朴素后半部分改成如下的就行了。假如我们求n,m区间,n1,m1大小的区间的极差。问题处理的核心思想是,把二维的东西化成一维。相当于处理二维化一维,一维求区间最大。
9.3.3 二维数组最大值 (20分)
MortyC99430的博客
06-11 4835
编程求出一个m行n列矩阵中最大值，以及其所在的行号和列号（从0开始），输入的元素各不相同。 输入格式: 输入在第1行中给出m,n（1≤m≤20,（1≤n≤20），在第2~m+1行中每行中给出n个整数，数字间以空格分隔。 输出格式: 第一行最大值 第二行最大值的行号 第三行最大值的列号 输入样例: 3 5 2 7 9 8 1 5 4 3 19 12 17 22 11 4 6 输出样例: 22 2 1 #include<stdio.h> int main() { int m,n,f,
最大二维子区间和 （贪心+前缀和）
07-22 256
【题目描述】 　　　　已知矩阵的大小定义为矩阵中所有元素的和。给定一个矩阵，你的任务是找到最大的非空(大小至少是1 * 1)子矩阵。 【题目链接】 　　　　http://noi.openjudge.cn/ch0406/1768/ 【算法】 　　　　众所周知一维子区间最大值要如何求，那么为求二维子区间最大和可以枚举任意两行，这两行间的同一列的数之和可以看成是一维子区间的一个值，问...
线性代数笔记29——正定矩阵和最小值
我是8位的
11-28 3502
　　原文链接 |https://mp.weixin.qq.com/s/wX6wmVSqJUTgbmk8Z1r2_w 判断正定矩阵 　　给出一个矩阵： 　　有4个途径可以判定该矩阵是否是正定矩阵（注意这个矩阵的4个元素中有2个b，这是因为正定矩阵是对称矩阵，如果A的次对角线的元素不相等，A就不是对称的，也就没有必要进一步判定是否是正定的）： 所有特征值大于0，λ1>0，λ2&gt...
6-3 求二维数组中的最大值
热门推荐
weixin_52555162的博客
10-28 1万+
本题要求实现一个函数，可计算二维数组中的最大值。 函数接口定义： int fun(int a[][N]); 在这里解释接口参数。例如：其中N和D都是用户传入的参数。N的值不超过int的范围；D是[0, 9]区间内的个位数。函数须返回N中D出现的次数。 裁判测试程序样例： #include <stdio.h> #define M 3 #define N 4 int fun(int a[][N]); int main() { int x[M][N...
hdu1823二维线段树，单点跟新，区间求最大值。
旺旺的博客
08-16 463
/* 思路：以身高为x，活泼度为y建立二维线段树，然后其余就是二维线段树的区间最值和单点跟新. 然后有 升高和活泼度相同，但是缘分不同的数据。 还有输入查询数据的时候 有可能 区间大的输入在区间小的的前面 */ #include #include #include #include using namespace std; int m; int sum[102*4+10][1002*4+10];
matalb 统计二维数组中最小值
05-29
你可以使用Matlab中的min函数来计算二维数组中的最小值。如果你有一个矩阵A，你可以使用以下代码来获取最小值： ``` min_value = min(min(A)); ``` 这个代码会先计算A中每一列的最小值，然后再从这些最小值中选出...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值