一般的暴力做法的时间复杂度为O(a*b*n*n),一定会超时。其中a*b是省不掉的,所以我们可以在对n*n的区间寻找最大值和最小值内进行优化。
https://www.luogu.org/blog/denial/solution-p2216
以上是本题的解题思路,但具体如何实现这个思路,我们就需要用到单调队列了。对于X(行内最大值)和Y(列最大值)我们维护一个从大到小的单调队列,然后对于读进来的点,我们判断一下如果它比队尾元素大,那就让队尾元素出队,然后将该元素放到队尾。然后因为长度必须为k,所以我们判断一下当前的位置与队首元素所在的坐标之差如果等于k(之所以不是大于等于是因为之前已经更新过了),那就让队首元素出队,而且能保证下一个队首元素一定在范围之内,然后我们标记X数组为该元素坐标所对应的位置。对于最小值也是类似的判断,只不过队列是反着来的。更新完Y数组之后,我们再枚举,得到区间最小值即可。
#include<cstdio>
#include<iostream>
#include<queue>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=1010;
int X[maxn][maxn],x[maxn][maxn];
int Y[maxn][maxn],y[maxn][maxn];
int map[maxn][maxn];
int n,m,k;
int main()
{
scanf("%d%d%d",&n,&m,&k);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
scanf("%d",&map[i][j]);
}
}
for(int i=1;i<=n;i++)
{
deque<int> Q,q;
for(int j=1;j<=m;j++)
{
while(!Q.empty()&&map[i][Q.back()]<map[i][j]) Q.pop_back();
while(!q.empty()&&map[i][q.back()]>map[i][j]) q.pop_back();
Q.push_back(j);
q.push_back(j);
if(j-Q.front()==k) Q.pop_front();
if(j-q.front()==k) q.pop_front();
if(j>=k) X[i][j-k+1]=map[i][Q.front()],x[i][j-k+1]=map[i][q.front()];
}
}
for(int j=1;j<=m-k+1;j++)
{
deque<int> Q,q;
for(int i=1;i<=n;i++)
{
while(!Q.empty()&&X[Q.back()][j]<X[i][j]) Q.pop_back();
while(!q.empty()&&x[q.back()][j]>x[i][j]) q.pop_back();
Q.push_back(i);
q.push_back(i);
if(i-Q.front()==k) Q.pop_front();
if(i-q.front()==k) q.pop_front();
if(i>=k) Y[i][j]=X[Q.front()][j],y[i][j]=x[q.front()][j];
}
}
int ans=1e9+10;
for(int i=k;i<=n;i++)
{
for(int j=1;j<=m-k+1;j++)
{
ans=min(ans,Y[i][j]-y[i][j]);
}
}
cout<<ans<<endl;
return 0;
}