算法第二章上机实践报告
题目:最大子列和问题
问题描述:
给定K个整数组成的序列{ N1, N2, ..., NK },“连续子列”被定义为{ Ni, Ni+1, ..., Nj },其中 1≤i≤j≤K。“最大子列和”则被定义为所有连续子列元素的和中最大者。例如给定序列{ -2, 11, -4, 13, -5, -2 },其连续子列{ 11, -4, 13 }有最大的和20。现要求你编写程序,计算给定整数序列的最大子列和。
本题旨在测试各种不同的算法在各种数据情况下的表现。各组测试数据特点如下:
- 数据1:与样例等价,测试基本正确性;
- 数据2:102个随机整数;
- 数据3:103个随机整数;
- 数据4:104个随机整数;
- 数据5:105个随机整数;
输入格式:
输入第1行给出正整数K (≤100000);第2行给出K个整数,其间以空格分隔。
输出格式:
在一行中输出最大子列和。如果序列中所有整数皆为负数,则输出0。
输入样例:
6
-2 11 -4 13 -5 -2
输出样例:
20
算法描述:
首先,我们可以把整个序列平均分成左右两部分,答案则会在以下三种情况中:
1、所求序列完全包含在左半部分的序列中。
2、所求序列完全包含在右半部分的序列中。
3、所求序列刚好横跨分割点,即左右序列各占一部分。
我们只要计算出:以分割点为起点向左的最大连续序列和、以分割点为起点向右的最大连续序列和,这两个结果的和就是第三种情况的答案。因为已知起点,所以这两个结果都能在O(N)的时间复杂度能算出来。
递归不断减小问题的规模,直到序列长度为1的时候,那答案就是序列中那个数字。
#include<iostream> using namespace std; int MaxSubSum(int* a, int left,int right ) { int sum=0; if(left==right) { sum=a[left]>0?a[left]:0; } else { int mid=(left+right)/2; int leftsum=MaxSubSum(a,left,mid); int rightsum=MaxSubSum(a,mid+1,right); int sum1=0; int lefts=0; for(int i=mid;i>=left;i--) { lefts +=a[i]; if(lefts>sum1) sum1=lefts; } int sum2=0; int rights=0; for(int i=mid+1;i<=right;i++) { rights +=a[i]; if(rights>sum2) sum2=rights; } sum=sum1+sum2; if(sum<leftsum) sum=leftsum; if(sum<rightsum) sum=rightsum; } return sum; } int main() { int n,a[100000]; cin>>n; for(int i=0;i<n;i++) cin>>a[i]; cout<<MaxSubSum(a,0,n-1); return 0; }
算法时间及空间复杂度分析:
设总的时间复杂度为T(N),分成两份每边都是T(N/2),还要加上算跨越分界线部分O(N)(表示每一个元素都扫描了一遍,即N的常数倍,每一次求跨越边界的子列都要扫描左右两边的元素,每个元素多会扫 描一遍以上),一直分到N/=1,也就是递归的终止,左右两边相等。T(1)=O(1),k =logN,
=N,ckN=cNlogN=O(NlogN)
当两个复杂度在一起时我们取大的那个O(NlogN)。