令f[i][j][k][x]表示第i行,从j~k列状态为x的时的最大值。x=0表示第一个矩形;x=1表示第二个矩形;x=2表示第三个矩形,然后构造两个辅助数组转移即可。
AC代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
#define N 205
using namespace std;
int m,n,s[N][N],dp[2][N][N][3],g[N][N][2];
int main(){
scanf("%d%d",&m,&n); int i,j,k;
for (i=1; i<=m; i++)
for (j=1; j<=n; j++){
scanf("%d",&k); s[i][j]=s[i][j-1]+k;
}
memset(dp[0],192,sizeof(dp[0])); memset(g,192,sizeof(g));
int last,now=0,ans=0,x;
for (i=1; i<=m; i++){
last=now; now^=1;
memset(dp[now],200,sizeof(dp[now]));
for (j=1; j<=n; j++)
for (k=j; k<=n; k++) if (s[i][j-1]==s[i][k]){
x=k-j+1;
dp[now][j][k][0]=max(dp[last][j][k][0],0)+x;
dp[now][j][k][1]=max(dp[last][j][k][1],g[j][k][0])+x;
ans=max(ans,dp[now][j][k][2]=max(dp[last][j][k][2],g[j][k][1])+x);
}
if (i==m) break;
memset(g,192,sizeof(g));
for (j=2; j<n; j++)
for (k=n-1; k>=j; k--)
g[j][k][0]=max(max(g[j-1][k][0],g[j][k+1][0]),dp[now][j-1][k+1][0]);
for (j=n-2; j; j--)
for (k=j+2; k<=n; k++)
g[j][k][1]=max(max(g[j+1][k][1],g[j][k-1][1]),dp[now][j+1][k-1][1]);
}
printf("%d\n",ans);
return 0;
}
by lych
2016.4.9