学习了最大子矩阵,网上找了各种帖子,还是那组经典的数据。
Example:
0 -2 -7 0
9 2 -6 2
-4 1 -4 1
-1 8 0 -2
其中左上角的子矩阵:
9 2
-4 1
-1 8
此子矩阵的值为9+2+(-4)+1+(-1)+8=15。
求最大子矩阵,可以转化为求最大子段和问题(三种方法)。对于一个给定的矩阵N*N,假设它的最大子矩阵为:第x行到第y行,第i列到第j列。那么就把每一行的相同列加起来啦,并存到一个数组中。求这个数组的最大子段和,就是这个矩阵的最大子矩阵。
/**
**author :Or_me **
╭︿︿︿╮
{/ A C /}
( (OO) )
︶︶︶
** **
** **
** 2014 年 7月 日**
**/
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <cstdio>
#include <string>
#include <cctype>
#include <climits>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
int maxsub(int n,int *a)//求最大字段和
{
int ans=-INT_MAX,sum=0;
for (int i=0;i<n;i++)
{
if (sum>0)
sum+=a[i];
else
sum=a[i];
if (sum>ans)
ans=sum;
}
return ans;
}
int main ()
{
int n;
while (~scanf("%d",&n))
{
int matrix[110][110],dp[110],sum=-INT_MAX,max;
for (int i=0;i<n;i++)
for (int j=0;j<n;j++)
scanf("%d",&matrix[i][j]);
for (int i=0;i<n;i++)
{
for (int x=0;x<110;x++) //每次使用前数组清零
dp[x]=0;
for (int j=i;j<n;j++)/第i行到第n-1行
{
for (int k=0;k<n;k++)//每一行的相同列
dp[k]+=matrix[j][k];
max=maxsub(n,dp);
if (max>sum)
sum=max;
}
}
printf("%d\n",sum);
}
return 0;
}