编程之美-2.15 子数组之和的最大值（二维）

给一个二维数组，求子数组的和的最大值，也就是求子矩阵的最大和。

这道题的解题思路与2.14求数组的子数组之和的最大值有类似之处，具体步骤如下：

设原数组为A，行数为row,列数为column.

1. 先求矩阵P(row+1,column+1), P[i,j]表示在A中以A[i-1,j-1]到A[0,0]为对角线的子矩阵的元素之和。

初始化：P[0,i]=0,i=0,1,...colum; P[i,0]=0,i=0,1,...row;

迭代求P[i][j]=P[i - 1][j] + P[i][j - 1] - P[i - 1][j - 1] + A[i-1][j-1]; //注意，因为res比v多了一行一列，因此res中的第[i,j]元素对应于v中第[i-1,j-1]的元素。

2. 将问题从二维转化为一维，假设已将确定了子矩阵的上下边界，分别为第r_min, r_max行，现在想要确定左右边界。在r_min, r_max的区域里，把每一列的元素看成一个整体，BC[i]即表示在此区域的第i列的和，BC[i]=A[r_min,i]+A[r_min+1,i]+...+A[r_max,i];因此问题为先求BC数组的每一个元素（i=0,1,...,column）,计算公式为BC[i]=P[r_max][i] - P[r_min - 1][i] - P[r_max][i- 1] + P[r_min - 1][i- 1]（注意这里也有1中出现的P的行列数比A大1，编程时需要转换一下的小问题具体可参考程序） ；然后转换为题2.14求数组BC的子数组之和的最大值。

3. 令r_min从0-(row-1), r_max从r_min-(row-1), 迭代所有可能的上下区域，并利用步骤2求该区域的最大子数组之和，如果发现更大的和则更新最大值。

#include <vector>
using namespace std;
vector<vector<int>> calculateSum(const vector<vector<int>> &v, int row, int column)
{
	vector<vector<int>>res(row+1);
	for (int i = 0; i <= row; i++)
		res[i].resize(column+1);
	for (int i = 0; i <= row; i++)
		res[i][0] = 0;
	for (int i = 0; i <= column; i++)
		res[0][i] = 0;
	for (int i = 1; i <= row; i++)
		for (int j = 1; j <= column; j++)
			res[i][j] = res[i - 1][j] + res[i][j - 1] - res[i - 1][j - 1] + v[i-1][j-1];
	//注意，因为res比v多了一行一列，因此res中的第[i,j]元素对应于v中第[i-1,j-1]的元素。
	return res;
}

int ithColMaxsum(const vector<vector<int>>&P,int r_min, int r_max, int iColumn)
{
	//r_min,r_max,iColumn加1是因为P是在原数据矩阵v的最上边和最左边各加了一行和一列，均设为0
	//因此MaxSum中调用ithColMaxsum求v[i,j]对应元素和在P矩阵里相当于求p[i+1,j+1]对应的元素和。
	r_min ++;
	r_max++;
	iColumn++;
	return P[r_max][iColumn] - P[r_min - 1][iColumn] - P[r_max][iColumn - 1] + P[r_min - 1][iColumn - 1];
}
int MaxSum(vector<vector<int>>v)
{
	int row = v.size();
	int column = v[1].size();
	int maxSum = v[0][0];
	vector<vector<int>>p = calculateSum(v, row, column);
	for (int i = 0; i < row; i++)
	{
		for (int j = i; j < row; j++)
		{
			int startNum = ithColMaxsum(p,i,j,0);
			int allNum = ithColMaxsum(p, i, j, 0);
			for (int k = 1; k < column; k++)
			{
				if (startNum<0)
					startNum = ithColMaxsum(p, i, j, k);
				else
					startNum += ithColMaxsum(p, i, j, k);
				if (startNum > allNum)
					allNum = startNum;
			}
			if (allNum > maxSum)
				maxSum = allNum;
		}
	}
	return maxSum;
}

int main(void)
{
	vector<vector<int>> v = { { 3, 2, -1 }, { -5, 5, 5 }, { -6, 4, -3 }, { 5, 6, -4 } };
	int res=MaxSum(v);
	return 0;
}

注意：用vector定义二维数组v（n*m），这里n为行数，m为列数，则可通过如下方法：

vector<vector<int> > v(n); //行数n一定不能少，它表示v有n个元素
for(int i=0;i<n;i++) 
v[i].resize(m);