题意:给定一个n*m大小的0-1矩阵,求所有元素都为1的最大子矩阵。
思路:把2维转化为一维(这一点有点像poj1050求最大子矩阵),枚举每一行。dp[i][j]的值表示第(i,j)个位置作为底,它上面的(包括自己)的连续1的个数。然后通过栈的单调性求每行的“最大值”(每行的情况如同poj2559)。可以用滚动数组优化。
#include <stdio.h>
#include <string.h>
#define min(a,b) ((a)<(b)?(a):(b))
#define max(a,b) ((a)>(b)?(a):(b))
#define N 2002
struct stack{
int id,num;
}stack[N];
int dp[N][N];
int n,m;
int test(int x[N]){
int top,i,temp,res=0;
x[m+1] = -1;
top = -1;
for(i = 1;i<=m+1;i++){
while(top>-1 && x[i]<stack[top].num){
temp = stack[top].num * (i-1 - (top? stack[top-1].id: 0));
res = max(res,temp);
top--;
}
stack[++top].id = i;
stack[top].num = x[i];
}
return res;
}
int main(){
freopen("a.txt","r",stdin);
while(scanf("%d %d",&n,&m)!=EOF){
int i,j,a,res=0;
memset(dp,0,sizeof(dp));
for(i = 1;i<=n;i++)
for(j = 1;j<=m;j++){
scanf("%d",&a);
if(a)
dp[i][j] = 1 + dp[i-1][j];
}
for(i = 1;i<=n;i++)
res = max(res,test(dp[i]));
printf("%d\n",res);
}
return 0;
}