bzoj 1047（单调队列）

传送门
题解：
设读入的矩阵为v[i][j]。
先求出每行len列矩阵的最大/小值g[i][j]=max/min(v[i][j]~v[i+len-1][j])，再求出每len行len列矩阵的最大/小值ret[0/1][i][j]=max/min(g[i][j]~g[i+len-1][j])
最后O(n^2)统计答案。
扯了这么多其实我自己也没弄太懂，改天精神好点再看。。。

#include<bits/stdc++.h>using namespace std;const int MAXN=1e3+4;int n,m,len,ans=0x3f3f3f3f;int v[MAXN][MAXN],g[MAXN][MAXN],ret[2][MAXN][MAXN];struct NODE {    int id,val;}q[MAXN];inline int read() {    int x=0;char c=getchar();    while (c<'0'||c>'9') c=getchar();    while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();    return x;}inline void solve(int type) {    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].id>=len) ++l;            if (!type) while (l<=r&&v[i][j]>=q[r].val) --r;            else while (l<=r&&v[i][j]<=q[r].val) --r;            q[++r]=(NODE){j,v[i][j]};            g[i][j]=q[l].val;        }    }    for (int j=len;j<=m;++j) {        int l=1,r=0;        for (int i=1;i<=n;++i) {            while (l<=r&&i-q[l].id>=len) ++l;            if (!type) while (l<=r&&g[i][j]>=q[r].val) --r;            else while (l<=r&&g[i][j]<=q[r].val) --r;            q[++r]=(NODE){i,g[i][j]};            ret[type][i][j]=q[l].val;        }    }}int main() {//  freopen("bzoj 1047.in","r",stdin);    n=read(),m=read(),len=read();    for (int i=1;i<=n;++i)        for (int j=1;j<=m;++j)            v[i][j]=read();    solve(0);    solve(1);    for (int i=len;i<=n;++i)        for (int j=len;j<=m;++j)            ans=min(ans,ret[0][i][j]-ret[1][i][j]);    printf("%d\n",ans);    return 0;}