最大子矩阵和

##子矩阵和的解释：：

####一个MN的矩阵，找到此矩阵的一个子矩阵，并且这个子矩阵的元素的和是最大的，输出这个最大的值。
例如：33的矩阵：
-1 3 -1
2 -1 3
-3 1 2
和最大的子矩阵是：
3 -1
-1 3
1 2

例题如下：

输入
第1行：M和N，中间用空格隔开（2 <= M,N <= 500)。
第2 - N + 1行：矩阵中的元素，每行M个数，中间用空格隔开。(-10^9 <= M[i] <= 10^9)
输出
输出和的最大值。如果所有数都是负数，就输出0。
输入示例
3 3
-1 3 -1
2 -1 3
-3 1 2
输出示例
7

我们把每一列第i行到第j行之间的和求出来，形成一个数组dp,于是一个第i行到第j行之间的最大子矩阵和对应于这个和数组dp的最大子段和。

代码如下：

#include<bits/stdc++.h>
using namespace std;
int main(){
    int m,n,r[510][510],dp[510];
    scanf("%d%d",&m,&n);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            scanf("%d",&r[i][j]);//r用来存储矩阵
        }
    }
    int maxn=-1;//maxn用来存最大的子矩阵和

    for(int i=1;i<=n;i++){//第一维度代表从第i行为起始行
        for(int j=i;j<=n;j++){//第二维度代表从第i行到第j行的矩阵每一列的矩阵和
            int sum=0;//用来存子矩阵和
            for(int k=1;k<=m;k++){
                dp[k]=(j==i)?r[i][k]:(dp[k]+r[j][k]);//如果j==i代表起始点等于dp[k]=r[i][k],否则直接求该列的子矩阵和
                if(sum<=0){
                    sum=0;
                }
                sum+=dp[k];//sum用来求矩阵和
                if(sum>maxn){
                    maxn=sum;
                }
                //printf("%d\n",dp[k]);
            }
        }
    }
    printf("%d\n",maxn);
    return 0;
}