求N*M矩阵中的最大子矩阵和(DP)

题目来自2018南京大学夏令营机试第一题。
计算过程如下图所示：就是固定矩阵的上边界，每一次计算的时候将相同列值相加，构成一个新的数组，然后对这个数组求最大连续和。时间复杂度为nnm
值得指出的是，原博客中有两个问题：
①最大值的初始值应该设置为无穷小，而不是0，因为可能子矩阵之和最大值为负数。
②maxsub(arr, m)函数传入的参数应该是矩阵列数m，而不是行数n

参考博客：
https://blog.csdn.net/qq_40421671/article/details/81252273

#include <iostream>
#include <cstring>
#include<algorithm>
using namespace std;

int maxsub(int a[], int n)
{
	int i, maxsum = 0x80000000, b = 0;
	for (i = 0; i < n; i++)
	{
		if (b > 0)
			b += a[i];
		else   //前面的连续子串之和小于等于0就重置b
			b = a[i];
		maxsum = max(maxsum, b);
	}
	return maxsum;
}

int main()
{
	int n, m, i, j, k, maxsubrec, maxsubarr;
	int dp[101][101], arr[101];
	while (cin >> n >> m)
	{    //输入矩阵 
		for (i = 0; i < n; i++)
			for (j = 0; j < m; j++)
				cin >> dp[i][j];
		//初始化最大和为int的最小值 
		maxsubrec = 0x80000000;
		//从上到下，依次计算 
		for (i = 0; i < n; i++)
		{  //上界每下移一行，就要重新开始计算每一列的和,就要重新初始化统计每列之和的arr 
			memset(arr, 0, sizeof(arr));
			//从上界开始，逐行增加 
			for (j = i; j < n; j++)
			{
				for (k = 0; k < m; k++)
					arr[k] += dp[j][k];//arr[k]即计算第k列的数字之和 
				maxsubarr = maxsub(arr, m);//计算每次最大子段和 
				if (maxsubarr > maxsubrec) maxsubrec = maxsubarr;//取子段和最大值 
			}
		}
		cout << maxsubrec << endl;
	}
}