单调队列、RMQ--[HAOI2007]理想的正方形

传送门

啊一看到这个题就想到了那个矩形的dp
但是时间复杂度好像不太够？？
有一个叫RMQ的算法，是求区间最值问题的
要是会的话这题就是裸题嘛（然而我不会
RMQ好像是基于倍增的？？？
还是单调队列好···
用stl好像常数比较大，不如手写

---

求矩形区间最大值，可以先求出每行的区间最大值
将行缩成点，再用求出的行最大值求每列的最大值
定义rmx[i][j]是第i行j-n+1到j的区间最大值，rmn是最小值
lmx[i][j]是以i，j为右下角的区间最大值，lmn是最小值
每求出一列的都更新一下答案
好像复制粘贴就可以解决0.0

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#define maxn 1005
using namespace std;
int a,b,n,rmx[maxn][maxn],rmn[maxn][maxn],lmx[maxn],lmn[maxn];
int v[maxn][maxn],val[maxn],pos[maxn],ans;

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

inline void work_row(){
    int l,r;
    for(int i=1;i<=a;i++){
        l=r=1;
        for(int j=1;j<=b;j++){
            while(l<r && val[r-1]<=v[i][j]) r--;
            val[r]=v[i][j], pos[r]=j, r++;
            if(pos[l]==j-n) l++;
            if(j>=n) rmx[i][j]=val[l];
        }
        l=r=1;
        for(int j=1;j<=b;j++){
            while(l<r && val[r-1]>=v[i][j]) r--;
            val[r]=v[i][j], pos[r]=j, r++;
            if(pos[l]==j-n) l++;
            if(j>=n) rmn[i][j]=val[l];
        }
    }
}

void work_line(){
    ans=0x3f3f3f3f;
    int l,r;
    for(int i=n;i<=b;i++){
        l=r=1;
        for(int j=1;j<=a;j++){
            while(l<r && val[r-1]<=rmx[j][i]) r--;
            val[r]=rmx[j][i], pos[r]=j, r++;
            if(pos[l]==j-n) l++;
            if(j>=n) lmx[j]=val[l];
        }
        l=r=1;
        for(int j=1;j<=a;j++){
            while(l<r && val[r-1]>=rmn[j][i]) r--;
            val[r]=rmn[j][i], pos[r]=j, r++;
            if(pos[l]==j-n) l++;
            if(j>=n) lmn[j]=val[l];
        }
        for(int j=n;j<=a;j++) ans=min(ans,lmx[j]-lmn[j]);
    }
}

int main(){
    a=rd(); b=rd(); n=rd();
    for(int i=1;i<=a;i++)
        for(int j=1;j<=b;j++)
            v[i][j]=rd();
    work_row();
    work_line();
    printf("%d\n",ans);
    return 0;
}