单调队列优化dp(八)——#10182. 「一本通 5.5 练习 3」理想的正方形

题目链接：https://loj.ac/problem/10182
解题思路

如果是一维的话，就是一个单调队列就能解决，现在变成二维了，其实是一样的，只不过，要减少的是一列，增加的是一列，每次找下一列的最小/最大值加入进去，其它按照一般单调队列维护即可。
AC代码

#include <iostream>
#include <stdio.h>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const int N=1e3+5;
const ll inf=0x3f3f3f3f;
int a,b,n,l_mi,r_mi,l_ma,r_ma;
int get_int()
{
    int sign=1;
    char c;
    while((c=getchar())>'9'&&c<'0')
    if(c=='-')
    sign=-1;
    int res=c-'0';
    while((c=getchar())<='9'&&c>='0')
    res=res*10+c-'0';
    return res*sign;
}
ll get()
{
    ll sign=1;
    char c;
    while((c=getchar())>'9'&&c<'0')
    if(c=='-')
    sign=-1;
    ll res=c-'0';
    while((c=getchar())<='9'&&c>='0')
    res=res*10+c-'0';
    return res*sign;
}
ll ans;
ll mp[N][N];
struct node
{
    int pos;
    ll val;
}q_ma[N*N],q_mi[N*N];
int main()
{
    scanf("%d%d%d",&a,&b,&n);
    ans=inf;
    for(int i=1;i<=a;++i)
    for(int j=1;j<=b;++j)
    scanf("%lld",&mp[i][j]);
    for(int i=n;i<=a;++i)
    {
        l_mi=r_mi=l_ma=r_ma=0;
        for(int k=1;k<=n;++k)
        {
            ll mi=inf;
            ll ma=0;
            for(int j=i-n+1;j<=i;++j)
            {
                mi=min(mi,mp[j][k]);
                ma=max(ma,mp[j][k]);
            }
            while(l_mi<r_mi&&q_mi[r_mi-1].val>=mi)
            r_mi--;
            while(l_ma<r_ma&&q_ma[r_ma-1].val<=ma)
            r_ma--;
            q_mi[r_mi].val=mi;
            q_mi[r_mi++].pos=k;
            q_ma[r_ma].val=ma;
            q_ma[r_ma++].pos=k;
        }
        //cout<<q_ma[l_ma].val<<" "<<q_mi[l_mi].val<<endl;
        ans=min(ans,q_ma[l_ma].val-q_mi[l_mi].val);
        for(int j=n+1;j<=b;++j)
        {
            while(l_mi<r_mi&&(j-q_mi[l_mi].pos)>=n)
            l_mi++;
            while(l_ma<r_ma&&(j-q_ma[l_ma].pos)>=n)
            l_ma++;
            ll mi=inf;
            ll ma=0;
            for(int k=i-n+1;k<=i;++k)
            {
                mi=min(mi,mp[k][j]);
                ma=max(ma,mp[k][j]);
            }
            while(l_mi<r_mi&&q_mi[r_mi-1].val>=mi)
            r_mi--;
            while(l_ma<r_ma&&q_ma[r_ma-1].val<=ma)
            r_ma--;
            q_mi[r_mi].pos=j;
            q_mi[r_mi++].val=mi;
            q_ma[r_ma].pos=j;
            q_ma[r_ma++].val=ma;
            ans=min(ans,q_ma[l_ma].val-q_mi[l_mi].val);
            //cout<<q_ma[l_ma].val<<" "<<q_mi[l_mi].val<<endl;
        }
    }
    printf("%lld\n",ans);
    return 0;
}