319
技术社区[云栖]
【算法小总结】最大连续子序列和最大连续子矩阵的关系与实现
最大连续子序列和最大连续子矩阵的关系与实现
求最长连续子序列的优化方法（非DP）
//求最大连续子序列和与对应的开头元素和结束元素
实现代码：
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<limits.h>
#define MAXN 100002
using namespace std;
int main()
{
int i,j,n,a[MAXN],Max,sum,max;
int ps,pe,ts,te;//记录开头元素和结束元素
scanf("%d",&n);
for(i=0;i<n;i++)
scanf("%d",&a[i]);
sum=0;
Max=INT_MIN;
for(i=0;i<n;i++)
{
if(sum<=0)//如果序列和为负数，初始化sum与开始结束下标
{
sum=a[i];
ts=i;
te=i;
}
else//如果序列和为正数，继续累加，并记录新的结尾坐标
{
sum+=a[i];
te=i;
}
if(sum>Max)//一旦发现更大的连续子序列，就更新最大值和开始结束下标
{
Max=sum;
ps=ts;
pe=te;
}
}
printf("%d %d %d\n",Max,a[ps],a[pe]);
return 0;
}
例子：
8
1 -3 -5 2 6 -1 4 9
20 2 9
类似的是，在矩阵中求最大子矩阵的方法可以转化成一维求最大连续子序列和的方法。
将矩阵的行或列压缩。这里我选择列压缩
如：
1 -1 2
4 9 15
-7 1 10
计算每一递进行的最长连续子序列和，什么是递进行？举个例子
全局变量MAX=-999999999
刚刚的矩阵
1 -1 2
4 9 15
-7 1 10
————————————————————————————————————
第一次记算取 1 -1 2
求出最长连续子序列是2，MAX=2
第二次取
1 -1 2
4 9 15
这里两行的每一列元素之和：5 8 17（从这里开始就体现矩阵压缩了）
求出最长连续子序列是30，MAX=30
第三次取
1 -1 2
4 9 15
-7 1 10
这里三行的每一列元素之和：-2 9 27
求出最长连续子序列是36，MAX=36
————————————————————————————————————
第四次取：4 9 15
求出最长连续子序列是28，MAX=36
第五次取
4 9 15
-7 1 10
这里两行的每一列元素之和：-3 10 25
求出最长连续子序列是35，MAX=36
————————————————————————————————————
第六次次取：-7 1 10
求出最长连续子序列是11，MAX=36
————————————————————————————————————
到此为止，所有可能的子矩阵和全都计算完毕，得出最大的子矩阵元素之和为36
即是
-1 2
9 15
1 10
这个子矩阵
实现代码：
#include <stdio.h>
#include <string.h>
#include <limits.h>
#define MAXN 10002
int a[MAXN][MAXN],temp[MAXN];
int n;
int getsum()
{
int i,sum=0,max=0;
for(i=0;i<n;i++)
{
sum+=temp[i];
if(sum>max) max=sum;
if(sum<0) sum=0;
}
return max;
}
int main()
{
int i,j,ans,k;
while(scanf("%d",&n)!=EOF)
{
memset(a,0,sizeof(a));
for(i=0;i<n;i++)
for(j=0;j<n;j++)
scanf("%d",&a[i][j]);
ans=INT_MIN;
for(i=0;i<n;i++)//从第i行开始
{
memset(temp,0,sizeof(temp));
for(j=i;j<n;j++)//从 i行到 n-1行都尝试一次
{
for(k=0;k<n;k++)//把 j至 k行的每一列都加起来，就是矩阵压缩
temp[k]+=a[j][k];
int pre=getsum();//计算压缩矩阵形成的一维数组的最长连续序列
if(ans<pre)//更新最大值
ans=pre;
}
}
printf("%d\n",ans);
}
return 0;
}
刚刚的例子答案：
3
1 -1 2
4 9 15
-7 1 10
36
验证方法无误～
欢迎访问程序员之洞https://blog.csdn.net/acmman/
最后更新：2017-04-03 05:39:54