【BZOJ 1047】 [HAOI2007]理想的正方形 单调队列（也可以dp水）

BZOJ上面是总时间，所以这一道题在BZOJ上是可以直接用滚动数组水过去的，但是还是不推荐，不过还是写出来吧

#include<cstdio>#include<cstring>#include<iostream>#define mmax(a,b,c,d) max(a,max(b,max(c,d)))#define mmin(a,b,c,d) min(a,min(b,min(c,d)))using namespace std;int f[1024][1024],m,n,K,g[1024][1024],mat[1024][1024];int read(){    int x=0;    char c=getchar();int flag=1;    for(;c>'9'||c<'0';c=getchar()){if(c=='-')flag=-1;}    for(;c>='0'&&c<='9';c=getchar())x=x*10+c-'0';    x=x*flag;    return x;} int main(){    scanf("%d%d%d",&n,&m,&K);    for(int i=1;i<=n;i++){        for(int j=1;j<=m;j++){            mat[i][j]=read();               f[i][j]=g[i][j]=mat[i][j];        }    }    for(int k=2;k<=K;k++){        for(int i=1;i+k-1<=n;i++){            for(int j=1;j+k-1<=m;j++){                f[i][j]=mmax(mat[i][j],f[i+1][j],f[i][j+1],f[i+1][j+1]);                g[i][j]=mmin(mat[i][j],g[i+1][j],g[i][j+1],g[i+1][j+1]);            }        }    }         int ans=1e9;    for(int i=1;i+K-1<=n;i++){        for(int j=1;j+K-1<=m;j++){            ans=min(ans,f[i][j]-g[i][j]);        }    }         printf("%d",ans);    return 0;}

好了，不要只是水题就满足了，正解应该是跑两次单调队列，一次计算一行连续k的最小最大，然后再利用跑出来的一列中连续k个元素的极值来跑行就可以计算出这一个矩阵了

#include<cstdio>#include<cstring>#include<iostream>#define maxn 1020using namespace std;int n,m,k,mat[maxn][maxn],f[3][maxn][maxn];pair<int,int> q[maxn*3];void front(int pos){//0是最小值1是最大值 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].first+1>k)l++;if(pos==0)while(l<=r&&q[r].second>=mat[i][j])r--;else while(l<=r&&q[r].second<=mat[i][j])r--;q[++r]=make_pair(j,mat[i][j]);f[pos][i][j]=q[l].second;}}for(int i=k;i<=m;i++){int l=1,r=0;for(int j=1;j<=n;j++){while(l<=r&&j-q[l].first+1>k)l++;if(pos==0)while(l<=r&&q[r].second>=f[pos][j][i])r--;else while(l<=r&&q[r].second<=f[pos][j][i])r--;q[++r]=make_pair(j,f[pos][j][i]);f[pos][j][i]=q[l].second;}}}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",&mat[i][j]);front(0),front(1);int ans=1e9;for(int i=k;i<=n;i++)for(int j=k;j<=m;j++)ans=min(ans,f[1][i][j]-f[0][i][j]);printf("%d",ans);return 0;}