传送门:luoguP2219
题解
细节题
考虑对于一个右下角为 ( x , y ) (x,y) (x,y)的 C × D C\times D C×D的矩形,能围住它的 A × B A\times B A×B的矩形右下角坐标范围为 ( i , j ) ( x + 1 ≤ i ≤ x + A − C − 1 , y + 1 ≤ j ≤ y + B − D − 1 ) (i,j)(x+1\leq i\leq x+A-C-1,y+1\leq j\leq y+B-D-1) (i,j)(x+1≤i≤x+A−C−1,y+1≤j≤y+B−D−1)。
用单调栈求出这样的一个矩形范围内的最大值即可。
具体来说,设:
s
[
i
]
[
j
]
=
∑
x
′
=
1
i
∑
y
′
=
1
j
X
[
x
′
]
[
y
′
]
s[i][j]=\sum\limits_{x'=1}^i\sum\limits_{y'=1}^jX[x'][y']
s[i][j]=x′=1∑iy′=1∑jX[x′][y′]。
e [ i ] [ j ] e[i][j] e[i][j]表示右下角为 ( i , j ) (i,j) (i,j)的 A × B A\times B A×B的矩形中值的总和。
q [ i ] [ j ] q[i][j] q[i][j]表示右下角为 ( i , j ) (i,j) (i,j)的 C × D C\times D C×D的矩形中值的总和。
f [ i ] [ j ] = m a x ( e [ x ′ ] [ y ′ ] , i ≤ x ′ ≤ i + A − C − 2 , j ≤ y ′ ≤ j + B − D − 2 ) f[i][j]=max(e[x'][y'],i\leq x'\leq i+A-C-2,j\leq y'\leq j+B-D-2) f[i][j]=max(e[x′][y′],i≤x′≤i+A−C−2,j≤y′≤j+B−D−2)。
则 a n s = m a x ( f [ i + 1 ] [ j + 1 ] − q [ i ] [ j ] ) ans=max(f[i+1][j+1]-q[i][j]) ans=max(f[i+1][j+1]−q[i][j])
代码
#include<bits/stdc++.h>
using namespace std;
const int N=1010;
int s[N][N],e[N][N],f[N][N];
int n,m,a,b,c,d,ans;
struct P{
int pos,v;
P(int pos_=0,int v_=0):pos(pos_),v(v_){};
}stk[N];
int main(){
int i,j,x,y,z,l,r,res;
scanf("%d%d%d%d%d%d",&n,&m,&a,&b,&c,&d);
y=b-d-1;x=a-c-1;
for(i=1;i<=n;++i){
for(j=1;j<=m;++j){
scanf("%d",&z);
s[i][j]=s[i][j-1]-s[i-1][j-1]+s[i-1][j]+z;
if(i<a || j<b) e[i][j]=0;
else e[i][j]=s[i][j]-s[i-a][j]-s[i][j-b]+s[i-a][j-b];
}
}
for(i=1;i<=m;++i){
f[n][i]=e[n][i];
stk[(l=r=1)]=P(n,e[n][i]);
for(j=n-1;j>0;--j){
z=e[j][i];
for(;l<=r && stk[l].pos>=j+x;++l);
for(;l<=r && stk[r].v<=z;--r);
stk[++r]=P(j,z);f[j][i]=stk[l].v;
}
}
for(i=1;i<=n;++i){
stk[(l=r=1)]=P(m,f[i][m]);
for(j=m-1;j>0;--j){
z=f[i][j];
for(;l<=r && stk[l].pos>=j+y;++l);
for(;l<=r && stk[r].v<=z;--r);
stk[++r]=P(j,z);f[i][j]=stk[l].v;
}
}
for(i=c+1;i<n;++i)
for(j=d+1;j<m;++j){
res=s[i][j]-s[i-c][j]-s[i][j-d]+s[i-c][j-d];
ans=max(ans,f[i+1][j+1]-res);
}
printf("%d",ans);
return 0;
}