bzoj1047/luogu2216 理想的正方形 (单调队列)

开b组单调队列，分别维护此时某一列中的最大/最小值

然后我每次把它们的头取出来，塞到维护行的单调队列里，就是n*n的最大/最小值

#include<bits/stdc++.h>
#define pa pair<int,int>
#define CLR(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
const int maxn=1e3+10;

inline ll rd(){
    ll x=0;char c=getchar();int neg=1;
    while(c<'0'||c>'9'){if(c=='-') neg=-1;c=getchar();}
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*neg;
}

struct Q{
    int q[maxn][2],h,t;
    inline void psh(int x,int y,bool b){
        while(h&&h>=t&&(b?q[h][0]<=x:q[h][0]>=x)) h--;
        q[++h][0]=x,q[h][1]=y;
        if(!t) t=h;
    }
    inline int pop(int lim){
        while(h>t&&q[t][1]<=lim) t++;
        return q[t][0];
    }
    inline void clr(){h=t=0;}
}ama[maxn],ami[maxn],bma,bmi;
int N,A,B,arr[maxn][maxn];

int main(){
    //freopen("","r",stdin);
    int i,j,k;
    A=rd(),B=rd();N=rd();
    for(i=1;i<=A;i++){
        for(j=1;j<=B;j++)
            arr[i][j]=rd();
    }
    int ans=2e9;
    for(i=1;i<=A;i++){
        for(j=1;j<=B;j++){
            ama[j].psh(arr[i][j],i,1);
            ami[j].psh(arr[i][j],i,0);
        }
        if(i>=N){
            bma.clr(),bmi.clr();
            for(j=1;j<=B;j++){
                int a=ama[j].pop(i-N),b=ami[j].pop(i-N);
                bma.psh(a,j,1);
                bmi.psh(b,j,0);
                if(j>=N) ans=min(ans,bma.pop(j-N)-bmi.pop(j-N));
            }
        }
    }
    printf("%d
",ans);
    return 0;
}