1、思路:
先用求区间最大连续子序列之和的方法求出未修改时题中最大子矩阵之是多少,并求出最大子矩阵的坐标
(先确定x坐标为L,R,再寻找y方向上的最大连续子序列之和,从而求出最大子矩阵及其坐标)。
如果修改的坐标不在子矩阵里,就对最大子矩阵没有影响;
如果修改的坐标在最大子矩阵里,就遍历子矩阵中的每个坐标,判断修改一个坐标对子矩阵产生的影响,
比较修改后的最大子矩阵的值和要修改的点的上下左右四个方向上的最大子矩阵的大小,得到新的最大子矩阵,
然后求出这些最大子矩阵中的最小值。
2、代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn = 205;
const int INF = 10000500;
int a[maxn][maxn],r1[maxn],r2[maxn],c1[maxn],c2[maxn];
//r1[i]表示1~i区间内的最大子矩阵(i上边的最大子矩阵)
//r2[i]表示i~n区间内的最大子矩阵(i下边的最大子矩阵)
//c1[j]表示1~j区间内的最大子矩阵(j左边的最大子矩阵)
//c2[j]表示j~n区间内的最大子矩阵(j右边的最大子矩阵)
int n,m,p,f[maxn],ma,xx1,yy1,xx2,yy2;
int main(void){
while(~scanf("%d%d%d",&n,&m,&p)){
r1[0] = r2[n+1] = -INF;
c1[0] = c2[m+1] = -INF;
for(int i=1;i<=n;i++){
r1[i] = r2[i] = -INF;
for(int j=1;j<=m;j++){
c1[j] = c2[j] = -INF;
scanf("%d",&a[i][j]);
}
}
ma = -INF;
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++) f[j] = 0;
for(int j=i;j<=n;j++){
for(int k=1;k<=m;k++) f[k] += a[j][k];
r1[j] = max(r1[j],r1[j-1]);
int sum = 0,st = 1;
for(int k=1;k<=m;k++){
if(sum<=0) sum = f[k],st = k;
else sum += f[k];
r1[j] = max(r1[j],sum);
if(ma<sum){
ma = sum;
xx1 = i;yy1 = st;
xx2 = j;yy2 = k;
}
}
}
}
//
for(int i=n;i>=1;i--){
for(int j=1;j<=m;j++) f[j] = 0;
for(int j=i;j>=1;j--){
for(int k=1;k<=m;k++) f[k] += a[j][k];
r2[j] = max(r2[j],r2[j+1]);
int sum = 0;
for(int k=1;k<=m;k++){
if(sum<=0) sum = f[k];
else sum += f[k];
r2[j] = max(r2[j],sum);
}
}
}
//
for(int i=1;i<=m;i++){
for(int j=1;j<=n;j++) f[j] = 0;
for(int j=i;j<=m;j++){
for(int k=1;k<=n;k++) f[k] += a[k][j];
c1[j] = max(c1[j],c1[j-1]);
int sum = 0;
for(int k=1;k<=n;k++){
if(sum<=0) sum = f[k];
else sum += f[k];
c1[j] = max(c1[j],sum);
}
}
}
//
for(int i=m;i>=1;i--){
for(int j=1;j<=n;j++) f[j] = 0;
for(int j=i;j>=1;j--){
for(int k=1;k<=n;k++) f[k] += a[k][j];
c2[j] = max(c2[j],c2[j+1]);
int sum = 0;
for(int k=1;k<=n;k++){
if(sum<=0) sum = f[k];
else sum += f[k];
c2[j] = max(c2[j],sum);
}
}
}
int tp = ma;
for(int i=xx1;i<=xx2;i++){
for(int j=yy1;j<=yy2;j++)
if(a[i][j]>p){
int mx = tp-a[i][j]+p;
if(i-1>=1) mx = max(mx,r1[i-1]);
if(i+1<=n) mx = max(mx,r2[i+1]);
if(j-1>=1) mx = max(mx,c1[j-1]);
if(j+1<=m) mx = max(mx,c2[j+1]);
ma = min(ma,mx);
}
}
printf("%d\n",ma);
}
return 0;
}