hihocoder 1580 Matrix [dp+枚举+rmq]

题意：给你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);
	}
	
}