题意:给你n*m的矩阵和一个p,必须用p替换掉一个位置的值,问最大的子矩阵和是多少。
题解:我们枚举矩阵的上边界和下边界,定义dp[305][2]表示以第i列为结尾的子矩阵被替换一次或者没被替换的最大和。这样我们就可以通过预处理某一列上边界到下边界的和,因为对于当前列我们必定替换最小的值,所以我们还需要rmq预处理某一列的区间最值。复杂度O(n^3)。
注意:dp的时候第一个位置不能用1状态转移= =!,还有假如我们取整个矩阵的时候只能取1的情况,不能取0的情况。
AC代码:
#include<stdio.h>
#include<algorithm>
#include<string.h>
#define inf 1000000000
using namespace std;
typedef long long ll;
const int MAXN = 305;
int rmq[MAXN][MAXN][10];
int sum[MAXN][MAXN];
int a[MAXN][MAXN];
int mm[MAXN];//最大的小于等于i的2^mm[i]
void initRMQ(int pos,int n)
{
mm[0]=-1;
for(int i=1;i<=n;i++)
{
mm[i]=((i&(i-1))==0)?mm[i-1]+1:mm[i-1];
rmq[pos][i][0]=a[i][pos];
}
for(int j=1;j<=mm[n];j++)
for(int i=1;i+(1<<j)-1<=n;i++)
rmq[pos][i][j]=min(rmq[pos][i][j-1],rmq[pos][i+(1<<(j-1))][j-1]);
}
int RMQ(int pos,int x,int y)
{
int k=mm[y-x+1];
return min(rmq[pos][x][k],rmq[pos][y-(1<<k)+1][k]);
}
int dp[MAXN][2];
int main()
{
int n,m,p;
while(~scanf("%d%d%d",&n,&m,&p))
{
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&a[i][j]);
for(int i=1;i<=m;i++)
{
initRMQ(i,n);
for(int j=1;j<=n;j++)
sum[i][j]=sum[i][j-1]+a[j][i];
}
int ans=p;
for(int i=1;i<=n;i++)
for(int j=i;j<=n;j++)
{
for(int i=0;i<MAXN;i++)
dp[i][0]=dp[i][1]=-inf;
int haha=0;
for(int k=1;k<=m;k++)
{
haha+=sum[k][j]-sum[k][i-1];
dp[k][0]=sum[k][j]-sum[k][i-1];
dp[k][1]=sum[k][j]-sum[k][i-1]-RMQ(k,i,j)+p;
dp[k][0]=max(dp[k][0],dp[k-1][0]+sum[k][j]-sum[k][i-1]);
dp[k][1]=max(dp[k][1],dp[k-1][1]+sum[k][j]-sum[k][i-1]);
dp[k][1]=max(dp[k][1],dp[k-1][0]+sum[k][j]-sum[k][i-1]-RMQ(k,i,j)+p);
ans=max(ans,dp[k][1]);
if(!(k==m&&haha==dp[k][0]&&(i==1&&j==n)))ans=max(ans,dp[k][0]);
}
}
printf("%d\n",ans);
}
}