[HAOI2007] 修筑绿化带

【问题描述】

为了增添公园的景致，现在需要在公园中修筑一个花坛，同时在画坛四周修建一片绿化带，让花坛被绿化带围起来。 
如果把公园看成一个M*N的矩形，那么花坛可以看成一个C*D的矩形，绿化带和花坛一起可以看成一个A*B的矩形。 
如果将花园中的每一块土地的“肥沃度”定义为该块土地上每一个小块肥沃度之和，那么， 
绿化带的肥沃度=A*B块的肥沃度-C*D块的肥沃度 
为了使得绿化带的生长得旺盛，我们希望绿化带的肥沃度最大。

【输入】:

第一行有6个正整数M,N,A,B,C,D 
接下来一个M*N的数字矩阵，其中矩阵的第i行j列元素为一个整数Xij,表示该花园的第i行第j列的土地“肥沃度”。

【输出】:

一个正整数，表示绿化带的最大肥沃程度。

【输入输出样例】
parterre.in
4 5 4 4 2 2 
20 19 18 17 16 
15 14 13 12 11 
10 9 8 7 6 
5 4 3 2 1 

parterre.out
132

【数据范围】

30%的数据，1<=M,N<=50 
100%的数据，1<=M,N<=1000,1<=A<=M,1<=B<=N,1<=C<=A-2,1<=D<=B-2，1<=“肥沃度”<=100

考虑答案为一个A*B矩阵权值和-其中最小C*D矩阵权值和

先拿单调队列滑出来第i,j-B+1 -> i,j格最小矩阵和

在拿单调队列维护上述最小矩阵和 那么可以处理出i-A+1,j-B+1 -> i,j格最小子矩阵 再更新答案

时间复杂度O(n*m)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <deque>
using namespace std;
const int maxn=1000+10;
int A[maxn][maxn];
int v[maxn][maxn];
int f[maxn][maxn];
deque<int>q;
int main(){
	freopen("parterre.in","r",stdin);
	freopen("parterre.out","w",stdout);
	int m,n,a,b,c,d;
	scanf("%d %d %d %d %d %d",&m,&n,&a,&b,&c,&d);
	for(int i=1;i<=m;i++)
		for(int j=1;j<=n;j++){
			scanf("%d",&A[i][j]);
			A[i][j]=A[i][j]+A[i-1][j]+A[i][j-1]-A[i-1][j-1];
		}
	memset(v,127/2,sizeof(v));
	for(int i=c;i<=m;i++){
		while(!q.empty())
			q.pop_back();
		for(int j=d;j<=n;j++){
			f[i][j]=A[i][j]+A[i-c][j-d]-A[i-c][j]-A[i][j-d];
			while(!q.empty()&&q.front()-d+1<=j-b+1)
				q.pop_front();
			if(!q.empty()){
				int u=q.front();
				v[i][j]=f[i][u];
			}
			while(!q.empty()&&f[i][q.back()]>=f[i][j])
				q.pop_back();
			q.push_back(j);
		}
	}
	int ans=0;
	for(int j=d;j<=n;j++){
		while(!q.empty())
			q.pop_back();
		for(int i=c;i<=m;i++){
			while(!q.empty()&&q.front()-c+1<=i-a+1)
				q.pop_front();
			if(i>=a&&j>=b&&!q.empty()){
				int summ=A[i][j]+A[i-a][j-b]-A[i-a][j]-A[i][j-b];
				ans=max(ans,summ-v[q.front()][j]);
			}
			while(!q.empty()&&v[q.back()][j]>=v[i][j])
				q.pop_back();
			q.push_back(i);
		}
	}
	printf("%d\n",ans);
	/*for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++)
			printf("%d ",A[i][j]);
		puts("");
	}*/
return 0;
}