洛谷P2219 [HAOI2007]修筑绿化带（单调队列）

传送门

啧……明明以前做到过这种类型的题结果全忘了……

这种矩阵的，一般都是先枚举行，然后对列进行一遍单调队列，搞出右下角在每一行中合法位置时的最小权值

再枚举列，对行做一遍单调队列，用之前搞出来的最小权值再做一次单调队列，更新答案

感觉很难讲清楚啊……看代码好了……虽然代码可能更不清楚……

//minamoto
#include<bits/stdc++.h>
using namespace std;
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[1<<21],*p1=buf,*p2=buf;
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
inline int read(){
    #define num ch-'0'
    char ch;bool flag=0;int res;
    while(!isdigit(ch=getc()))
    (ch=='-')&&(flag=true);
    for(res=num;isdigit(ch=getc());res=res*10+num);
    (flag)&&(res=-res);
    #undef num
    return res;
}
const int N=1005;
int sum[N][N],sum1[N][N],sum2[N][N],q[N],mn[N][N];
int n,m,A,B,C,D,h,t,res;
void init(){
    for(int i=C+1;i<n;++i){
        h=1,t=0;
        for(int j=D+1;j<m;++j){
            while(h<=t&&sum2[i][q[t]]>=sum2[i][j]) --t;
            q[++t]=j;
            if(q[h]<=j-B+D+1) ++h;
            mn[i][j]=sum2[i][q[h]];
        }
    }
}
void solve(){
    for(int j=D+1;j<m;++j){
        h=1,t=0;
        for(int i=C+1;i<n;++i){
            while(h<=t&&mn[q[t]][j]>=mn[i][j]) --t;
            q[++t]=i;
            if(q[h]<=i-A+C+1) ++h;
            cmax(res,sum1[i+1][j+1]-mn[q[h]][j]);
        }
    }
}
int main(){
//    freopen("testdata.in","r",stdin);
    n=read(),m=read(),A=read(),B=read(),C=read(),D=read();
    if(n==2||m==2) return puts("0"),0;
    for(int i=1;i<=n;++i) for(int j=1;j<=m;++j)
    sum[i][j]=read()+sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1];
    for(int i=A;i<=n;++i) for(int j=B;j<=m;++j)
    sum1[i][j]=sum[i][j]-sum[i-A][j]-sum[i][j-B]+sum[i-A][j-B];
    for(int i=C;i<=n;++i) for(int j=D;j<=m;++j)
    sum2[i][j]=sum[i][j]-sum[i-C][j]-sum[i][j-D]+sum[i-C][j-D];
    init(),solve();
    printf("%d
",res);
    return 0;
}