动态规划——求最大子序列
1.先来分析一维的:
#include<stdio.h> int main() { int t,m,n,sum,maxSum; scanf("%d",&m); while(m--) { scanf("%d%d",&n,&sum); maxSum=sum; while(--n) { scanf("%d",&t); if(sum<0) sum=t; else sum+=t; if(maxSum<sum) maxSum=sum; } printf("%d\n",maxSum); } return 0; }
这里需要两个结论。首先,对于array[1...n],如果array[i...j]就是满足和最大的子串,那么对于任何k(i<=k<=j),我们有array[i...k]的和大于0。因为如果存在k使得array[i...k]的和小于0,那么我们就有array[k+1...j]的和大于array[i...j],这与我们假设的array[i...j]就是array中和最大子串矛盾。
其次,我们可以将数组从左到右分割为若干子串,使得除了最后一个子串之外,其余子串的各元素之和小于0,且对于所有子串array[i...j]和任意k(i<=k<j),有array[i...k]的和大于0。此时我们要说明的是,满足条件的和最大子串,只能是上述某个子串的前缀,而不可能跨越多个子串。我们假设array[p...q],是array的和最大子串,且array[p...q],跨越了array[i...j],array[j+1...k]。根据我们的分组方式,存在i<=m<j使得array[i...m]的和是array[i...j]中的最大值,存在j+1<=n<k使得array[j+1...n]的和是array[j+1...k]的最大值。由于array[m+1...j]使得array[i...j]的和小于0。此时我们可以比较array[i...m]和array[j+1...n],如果array[i...m]的和大于array[j+1...n]则array[i...m]>array[p...q],否array[j+1...n]>array[p...q],无论谁大,我们都可以找到比array[p...q]和更大的子串,这与我们的假设矛盾,所以满足条件的array[p...q]不可能跨越两个子串。对于跨越更多子串的情况,由于各子串的和均为负值,所以同样可以证明存在和更大的非跨越子串的存在。对于单元素和最大的特例,该结论也适用。
根据上述结论,我们就得到了Kadane算法的执行流程,从头到尾遍历目标数组,将数组分割为满足上述条件的子串,同时得到各子串的最大前缀和,然后比较各子串的最大前缀和,得到最终答案。我们以array={−2, 1, −3, 4, −1, 2, 1, −5, 4}为例,来简单说明一下算法步骤。通过遍历,可以将数组分割为如下3个子串(-2),(1,-3),(4,-1,2,1,-5,4),这里对于(-2)这样的情况,单独分为一组。各子串的最大前缀和为-2,1,6,所以目标串的最大子串和为6。
南阳理工44就是一个一维的字串和:
子串和
描述
- 给定一整型数列{a1,a2...,an},找出连续非空子串{ax,ax+1,...,ay},使得该子序列的和最大,其中,1<=x<=y<=n。
输入
- 第一行是一个整数N(N<=10)表示测试数据的组数)
每组测试数据的第一行是一个整数n表示序列中共有n个整数,随后的一行里有n个整数I(-100=<I<=100),表示数列中的所有元素。(0<n<=1000000) - 输出
- 对于每组测试数据输出和最大的连续子串的和。
- 样例输入
-
1 5 1 2 -1 3 -2
- 样例输出
-
5
#include<stdio.h> int main() { int i,j,n,m,sum,maxsum; scanf("%d",&n); while(n--) { scanf("%d%d",&m,&j); sum=maxsum=j; for(i=1;i<m;i++) { scanf("%d",&j); if(sum<0) sum=j; else sum+=j; if(maxsum<sum) maxsum=sum; } printf("%d\n",maxsum); } return 0; }
最大和
时间限制:1000 ms | 内存限制:65535 KB难度:5- 描述
-
给定一个由整数组成二维矩阵(r*c),现在需要找出它的一个子矩阵,使得这个子矩阵内的所有元素之和最大,并把这个子矩阵称为最大子矩阵。
例子:
0 -2 -7 0
9 2 -6 2
-4 1 -4 1
-1 8 0 -2
其最大子矩阵为:9 2
-4 1
-1 8
其元素总和为15。
- 输入
- 第一行输入一个整数n(0<n<=100),表示有n组测试数据;
每组测试数据:
第一行有两个的整数r,c(0<r,c<=100),r、c分别代表矩阵的行和列;
随后有r行,每行有c个整数; - 输出
- 输出矩阵的最大子矩阵的元素之和。
- 样例输入
-
1 4 4 0 -2 -7 0 9 2 -6 2 -4 1 -4 1 -1 8 0 -2
- 样例输出
-
15
-
//我们已经知道一维数组的最值求法,现在我们可以把二位数组转化为一维数组 //我们可以把二维数组按照行数进行压缩 #include<stdio.h> #include<string.h> int max(int n,int *a)//这是求一维数组的最大值 { int i,sum,maxsum; sum=maxsum=-9999999; for(i=0;i<n;i++) { if(sum<0) sum=a[i]; else sum+=a[i]; maxsum=maxsum>sum?maxsum:sum; } return maxsum; } int main() { int i,j,n,m,sum,maxsum,t,a[110][110],b[110],p; scanf("%d",&p); while(p--) { maxsum=-9999999; scanf("%d%d",&n,&m); for(i=0;i<n;i++) for(j=0;j<m;j++) scanf("%d",&a[i][j]); for(i=0;i<n;i++)//行号从0到n { memset(b,0,sizeof(b)); for(j=i;j<n;j++)//把从i到n行压缩成一行 { for(t=0;t<m;t++) b[t]+=a[j][t]; sum=max(m,b);//求每行的最值 maxsum=maxsum>sum?maxsum:sum;//总得最值 } } printf("%d\n",maxsum); } return 0; }