Codeforce 375B 给定5000的布尔矩阵 求最大面积的全1子矩阵

题目链接：http://codeforces.com/problemset/problem/375/B

任意行变换后，求出最大的全1子矩阵面积

n^2的复杂度

用sum[i][j] 表示 i列到j列之间 有多少行是全1的

显然我们可以得到答案是 ans = max(sum[i][j] * (j-i +1))

 

对于sum数组， 我们首先得到 a[i][j] 表示 从(i,j) 开始向←连续1的个数

利用a数组 我们可以得到 sum[i][j] 表示 最大连续1的列 （是中空的）

然后sum数组做一次 左前缀和, 再做一次右前缀和

即可得到sum数组

 

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#define N 50
using namespace std;
int sum[N][N];//sum[i][j]表示 i列到j列都为1时 存在多少这样的行
int a[N][N];
int n, m;
int main(){
	int i, j, k;
	memset(a, 0, sizeof(a));

	while(~scanf("%d %d",&n,&m)){
		for(i = 1;i <= n;i++)
		{char s[N]; scanf("%s",s);
		for(j = 1;j <= m;j++)
			a[i][j] = s[j-1] == '1';
		memset(sum[i], 0, sizeof(sum[i]));
		}
		for(i = 1;i <= n; i++)//a[i][j]表示 i行j列的位置 →来连续1的个数
			for(j = 1;j <= m; j++)
				a[i][j] = a[i][j] ? (a[i][j-1] +1) : 0;
		for(i = 1;i <= n; i++)//这一步算出特定的sum[i][j]（i行j列之间有一行1（但i-1行j列这样包含的关系并未计算（就是中空的））
		{
			for(j = m;j >= 1; j--)if(a[i][j])//每行最大的连续1的列 计算一次
			{
				k = j - a[i][j] +1;
				sum[k][j] ++;
				j = k;//只算最大连续1列 且只算一次（中空的）
			}
		}
		for(i = 1;i <=m ;i++)//对于第i列
		{
			for(j = m-1; j>=i; j--)//从第i列到最后一列
				sum[i][j] += sum[i][j+1];
		}
		for(i = 1;i <= m; i++)//对于第i列
		{
			for(j = 1; j <= i; j++)//从第j列到第i列
				sum[j][i] += sum[j-1][i];
		}
		int ans = 0;
		for(i = 1;i <= m; i++)
			for(j = i;j <= m; j++)
				ans = max(ans, sum[i][j]*(j-i+1));
		printf("%d\n",ans);
	}
	return 0;
}
/*
1 1
1
2 2
10
11
4 3
100
011
000
101
11 16
0111110101100011
1000101100010000
0010110110010101
0110110010110010
0011101101110000
1001100011010111
0010011111111000
0100100100111110
1001000000100111
0110000011001000
1011111011010000

*/