BZOJ原题链接

洛谷原题链接

显然可以用数据结构或\(ST\)表或单调队列来维护最值。

这里采用单调队列来维护。

先用单调队列维护每一行的最大值和最小值,区间长为正方形长度。

再用单调队列维护之前维护出的每行最值数组的每一列的最大值和最小值,区间同上。

最后维护出的数组其实就是以每个点为左上角的正方形中的最值,直接扫一遍求最小的差即可。

借用洛谷题解里大佬的图以更好说明:



(\(X,x\)分别是维护出原矩阵中行的最大、最小值,\(Y,y\)分别是维护\(X,x\)中列的最大、最小值)

#include<cstdio>
using namespace std;
const int N = 1010;
const int M = 2e6 + 10;
int a[N][N], ma_x[N][N], mi_x[N][N], ma_y[N][N], mi_y[N][N], qma[M], qmi[M];
inline int re()
{
int x = 0;
char c = getchar();
bool p = 0;
for (; c < '0' || c > '9'; c = getchar())
p |= c == '-';
for (; c >= '0' && c <= '9'; c = getchar())
x = x * 10 + c - '0';
return p ? -x : x;
}
inline int minn(int x, int y)
{
return x < y ? x : y;
}
int main()
{
int i, j, n, m, k, lmi, lma, rmi, rma, o, oo, mi = 1e9;
n = re();
m = re();
k = re();
for (i = 1; i <= n; i++)
for (j = 1; j <= m; j++)
a[i][j] = re();
for (i = 1; i <= n; i++)
{
lmi = lma = 1;
rmi = rma = 0;
for (j = 1; j < k; j++)
{
for (; lma <= rma && a[i][qma[rma]] <= a[i][j]; rma--);
for (; lmi <= rmi && a[i][qmi[rmi]] >= a[i][j]; rmi--);
qma[++rma] = j;
qmi[++rmi] = j;
}
for (; j <= m; j++)
{
for (; lma <= rma && j - qma[lma] + 1 > k; lma++);
for (; lmi <= rmi && j - qmi[lmi] + 1 > k; lmi++);
for (; lma <= rma && a[i][qma[rma]] <= a[i][j]; rma--);
for (; lmi <= rmi && a[i][qmi[rmi]] >= a[i][j]; rmi--);
qma[++rma] = j;
qmi[++rmi] = j;
ma_x[i][j - k + 1] = a[i][qma[lma]];
mi_x[i][j - k + 1] = a[i][qmi[lmi]];
}
}
for (i = 1, o = m - k + 1; i <= o; i++)
{
lmi = lma = 1;
rmi = rma = 0;
for (j = 1; j < k; j++)
{
for (; lma <= rma && ma_x[qma[rma]][i] <= ma_x[j][i]; rma--);
for (; lmi <= rmi && mi_x[qmi[rmi]][i] >= mi_x[j][i]; rmi--);
qma[++rma] = j;
qmi[++rmi] = j;
}
for (; j <= n; j++)
{
for (; lma <= rma && j - qma[lma] + 1 > k; lma++);
for (; lmi <= rmi && j - qmi[lmi] + 1 > k; lmi++);
for (; lma <= rma && ma_x[qma[rma]][i] <= ma_x[j][i]; rma--);
for (; lmi <= rmi && mi_x[qmi[rmi]][i] >= mi_x[j][i]; rmi--);
qma[++rma] = j;
qmi[++rmi] = j;
ma_y[j - k + 1][i] = ma_x[qma[lma]][i];
mi_y[j - k + 1][i] = mi_x[qmi[lmi]][i];
}
}
for (i = 1, o = n - k + 1, oo = m - k + 1; i <= o; i++)
for (j = 1; j <= oo; j++)
mi = minn(mi, ma_y[i][j] - mi_y[i][j]);
printf("%d", mi);
return 0;
}

BZOJ1047或洛谷2216 [HAOI2007]理想的正方形的更多相关文章

  1. 洛谷 2216 [HAOI2007]理想的正方形

    题目戳这里 一句话题意 给你一个a×b的矩形,求一个n×n的子矩阵,矩阵里面的最大值和最小值之差最小. Solution 这个题目许多大佬都是单调队列,但是我不是很会,只好用了比较傻逼的方法: 首先我 ...

  2. 洛谷 P2216 [HAOI2007]理想的正方形

    P2216 [HAOI2007]理想的正方形 题目描述 有一个a*b的整数组成的矩阵,现请你从中找出一个n*n的正方形区域,使得该区域所有数中的最大值和最小值的差最小. 输入输出格式 输入格式: 第一 ...

  3. 洛谷P2216: [HAOI2007]理想的正方形 单调队列优化DP

    洛谷P2216 )逼着自己写DP 题意: 给定一个带有数字的矩阵,找出一个大小为n*n的矩阵,这个矩阵中最大值减最小值最小. 思路: 先处理出每一行每个格子到前面n个格子中的最大值和最小值.然后对每一 ...

  4. 洛谷P2216 HAOI2007 理想的正方形 (单调队列)

    题目就是要求在n*m的矩形中找出一个k*k的正方形(理想正方形),使得这个正方形内最值之差最小(就是要维护最大值和最小值),显然我们可以用单调队列维护. 但是二维平面上单调队列怎么用? 我们先对行处理 ...

  5. 【DP】【单调队列】洛谷 P2216 [HAOI2007]理想的正方形 题解

        算是单调队列的复习吧,不是很难 题目描述 有一个$a\times b$的整数组成的矩阵,现请你从中找出一个$n\times n$的正方形区域,使得该区域所有数中的最大值和最小值的差最小. 输入 ...

  6. [洛谷P2216][HAOI2007]理想的正方形

    题目大意:有一个$a\times b$的矩阵,求一个$n\times n$的矩阵,使该区域中的极差最小. 题解:二维$ST$表,每一个点试一下是不是左上角就行了 卡点:1.用了一份考试时候写的二维$S ...

  7. 洛谷 P2216 [HAOI2007]理想的正方形 || 二维RMQ的单调队列

    题目 这个题的算法核心就是求出以i,j为左上角,边长为n的矩阵中最小值和最大值.最小和最大值的求法类似. 单调队列做法: 以最小值为例: q1[i][j]表示第i行上,从j列开始的n列的最小值.$q1 ...

  8. [Luogu 2216] [HAOI2007]理想的正方形

    [Luogu 2216] [HAOI2007]理想的正方形 题目描述 有一个a*b的整数组成的矩阵,现请你从中找出一个n*n的正方形区域,使得该区域所有数中的最大值和最小值的差最小. 输入输出格式 输 ...

  9. 洛谷 P2216 [HAOI2007]理想正方形

    洛谷 巨说这是一道单调队列好题,但是我并不是用单调队列做的诶. 如果往最暴力的方向去想,肯定是\(n^3\)的\(dp\)了. \(f[i][j][k]\)代表当前正方形的左上角定点是\((i,j)\ ...

随机推荐

  1. git 分支合并develop 重新拉取

    1.本地包 zb-local 上传到对应的线上分支 zb-online cd 进入本地所在文件夹,然后执行下面三项 git add . git commit -m '修改' git push   或者 ...

  2. vue watch详细用法

    在vue中,使用watch来响应数据的变化.watch的用法大致有三种.下面代码是watch的一种简单的用法 html: <input type="text" v-model ...

  3. 使用AJAX实现文件上传时Illegal invocation错误

    在参数里面加 processData:false, 就行了

  4. spring 代理注解 <aop:aspectj-autoproxy />

    spring默认使用jdk的代理方式,使用jdk的代理方式我们知道,代理的类需要实现一个接口,若果没有就会报,java.lang.NoSuchMethodException: com.sun.prox ...

  5. Collection集合总结,List和set集合的用法,HashSet和LinkedHashSetde用法

    首先:Collection是List集合和Set集合的父类,同时,它们三个都为接口,不能直接实例化,所以需要,List和set的子类来实例化. List的子类: ArrayList类和LinkedLi ...

  6. Displacement Mapping

    [Displacement Mapping] Displacement Mapping(移位贴图). Normal maps的另一个应用是displacement mapping,在这个应用中,细节的 ...

  7. Mybatis多个in查询

    Map<String, Object> params = null; List<Map<String, Object>> list=new ArrayList(); ...

  8. 解题(MiGong--迷宫问题(深度搜索))

    题目描述 定义一个二维数组N*M(其中2<=N<=10;2<=M<=10),如5 × 5数组下所示: int maze[5][5] = { 0, 1, 0, 0, 0, 0,  ...

  9. Excel学习之图表创建

    前言:由于工作的需要,现在需要恶补Excel知识,温故而知新,不记录下来自己的学习的点点滴滴怎么能行 Excel版本:2016 一.插入如下图所示的数据源: 1.勾选图表区域范围 2.选中区域范围后直 ...

  10. jquery获取焦点位于的元素

    参考 https://zhidao.baidu.com/question/497311181016643684.html if ($(":focus").length > 0 ...