传送门
题解:
设读入的矩阵为v[i][j]。
先求出每行len列矩阵的最大/小值g[i][j]=max/min(v[i][j]~v[i+len-1][j]),再求出每len行len列矩阵的最大/小值ret[0/1][i][j]=max/min(g[i][j]~g[i+len-1][j])
最后O(n^2)统计答案。
扯了这么多其实我自己也没弄太懂,改天精神好点再看。。。
#include<bits/stdc++.h>
using namespace std;
const int MAXN=1e3+4;
int n,m,len,ans=0x3f3f3f3f;
int v[MAXN][MAXN],g[MAXN][MAXN],ret[2][MAXN][MAXN];
struct NODE {
int id,val;
}q[MAXN];
inline int read() {
int x=0;char c=getchar();
while (c<'0'||c>'9') c=getchar();
while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
return x;
}
inline void solve(int type) {
for (int i=1;i<=n;++i) {
int l=1,r=0;
for (int j=1;j<=m;++j) {
while (l<=r&&j-q[l].id>=len) ++l;
if (!type) while (l<=r&&v[i][j]>=q[r].val) --r;
else while (l<=r&&v[i][j]<=q[r].val) --r;
q[++r]=(NODE){j,v[i][j]};
g[i][j]=q[l].val;
}
}
for (int j=len;j<=m;++j) {
int l=1,r=0;
for (int i=1;i<=n;++i) {
while (l<=r&&i-q[l].id>=len) ++l;
if (!type) while (l<=r&&g[i][j]>=q[r].val) --r;
else while (l<=r&&g[i][j]<=q[r].val) --r;
q[++r]=(NODE){i,g[i][j]};
ret[type][i][j]=q[l].val;
}
}
}
int main() {
// freopen("bzoj 1047.in","r",stdin);
n=read(),m=read(),len=read();
for (int i=1;i<=n;++i)
for (int j=1;j<=m;++j)
v[i][j]=read();
solve(0);
solve(1);
for (int i=len;i<=n;++i)
for (int j=len;j<=m;++j)
ans=min(ans,ret[0][i][j]-ret[1][i][j]);
printf("%d\n",ans);
return 0;
}