最大子矩阵全1的个数

/*
方法是：

1、先将0/1矩阵读入x，对每一个非零元素x[i][j]，将其更新为：在本行，它前面的连续的1的个数+1（+1表示算入自身）

　　比如，若某一行为0 1 1 0 1 1 0，则更新为0 1 2 0 1 2 0

2、对每一个非零元素x[i][j]，在第j列向上和向下扫描，直到遇到比自身小的数，若扫描了y行，则得到一个大小为x[i][j]*(y+1)的全1子矩阵（+1表示算入自身所在行）

　　比如，若某一列为[0 3 4 3 5 2 1]'(方便起见，这里将列表示成一个列向量)，我们处理这一列的第3个元素，也就是4，它向上可以扫描。。个元素，向下可以扫描。。个元素，于是得到一个4×。。的全1子矩阵。

3、在这些数值中取一个最大的
 */
#include<iostream>
using namespace std;
void CalMax(int a[][100],int m,int n)
{
	int maxsum=0;
	for (int i = 0; i < m; ++i)
	{
		for (int j = 0; j < n; ++j)
		{
			int tmp=0;
			int value=a[i][j];
			int count=0;
			for (int up = i-1; up >=0; --up)
			{
				if (a[up][j]>=value&&value!=0)
				{
					count++;
				}
				else
					break;
			}
			for (int down = i+1; down<m ; ++down)
			{
				if (a[down][j]>=value&&value!=0)
				{
					count++;
				}
				else
					break;
			}
			if (count!=0)
			{
				tmp=a[i][j]*(count+1);
			}
			if (tmp>maxsum)
			{
				maxsum=tmp;
			}
		}
	}
	cout<<maxsum<<endl;
}
void CalSum(int a[][100],int m,int n)
{
	for (int i = 0; i < m; ++i)
	{
		for (int j = 1; j <n ; ++j)
		{
			if (a[i][j]!=0)
			{
				a[i][j]+=a[i][j-1];
			}
		}
	}
	cout<<endl;
	CalMax(a,m,n);
}
int main()
{
	int a[100][100];
	int m,n;
	cin>>m>>n;
	int data;
	for (int i = 0; i < m; ++i)
	{
		for (int j = 0; j < n; ++j)
		{
			cin>>data;
			a[i][j]=data;
		}
	}
	CalSum(a,m,n);
	for (int i1 = 0; i1 < m; ++i1)
	{
		for (int j1 = 0; j1 < n; ++j1)
		{
			cout<<a[i1][j1]<<" ";
		}
		cout<<endl;
	}
	cout<<endl;
     return 0;
}

1 1 1 1 1 0 1
0 0 1 1 0 1 1 
0 0 1 1 1 0 1
1 1 0 0 1 0 0
得到
1 2 3 4 5 0 1
0 0 1 2 0 1 1
0 0 1 2 3 0 1
1 2 0 0 1 0 0