数字三角形(动态规划)--算法学习

问题

7
3 8
8 1 0
2 7 4 4
4 5 2 6 5

在上面的数字三角形中寻找一条从顶部到底边的路径，使得
路径上所经过的数字之和最大。路径上的每一步都只能往左下或
右下走。只需要求出这个最大和即可，不必给出具体路径。

三角形的行数大于1小于等于100，数字为 0 - 99

输入格式：

5 //三角形行数。下面是三角形
7
3 8
8 1 0
2 7 4 4
4 5 2 6 5

输出

要求输出最大和

递归
//也是我能想到的
思路
来自北大
郭 炜

用二维数组存放数字三角形。
D( r, j) : 第r行第 j 个数字(r,j从1开始算)
MaxSum(r, j) : 从D(r,j)到底边的各条路径中，
最佳路径的数字之和。
问题：求 MaxSum(1,1)
典型的递归问题。
D(r, j)出发，下一步只能走D(r+1,j)或者D(r+1, j+1)。故对于N行的三角形：
if ( r == N)
MaxSum(r,j) = D(r,j)
else
MaxSum( r, j) = Max{ MaxSum(r＋1,j), MaxSum(r+1,j+1) } 加
D(r,j)

递归程序:

#include <iostream> 
#include <algorithm>
#define MAX 101 
using namespace std;
int D[MAX][MAX]; 
int n; 
int MaxSum(int i, int j){ 
if(i==n)
return D[i][j]; 
int x = MaxSum(i+1,j); 
int y = MaxSum(i+1,j+1); 
return max(x,y)+D[i][j]; 

#include <iostream> 
#include <algorithm>
#define MAX 101 
using namespace std;
int D[MAX][MAX]; 
int n; 
int MaxSum(int i, int j){ 
if(i==n)
return D[i][j]; 
int x = MaxSum(i+1,j); 
int y = MaxSum(i+1,j+1); 
return max(x,y)+D[i][j]; 
}
}

显然易见，超时。
原因：

重复计算，如果采用递规的方法，深度遍历每条路径，存在大
量重复计算。则时间复杂度为 2n,对于 n = 100 行，
肯定超时

改进

如果每算出一个MaxSum(r,j)就保存起来，下次用
到其值的时候直接取用，则可免去重复计算。那么
可以用O(n2)时间完成计算。因为三角形的数字总
数是 n(n+1)/2

记忆递归型动归

#include <algorithm>
using namespace std;
#define MAX 101 
int D[MAX][MAX]; int n; 
int maxSum[MAX][MAX];
int MaxSum(int i, int j){ 
if( maxSum[i][j] != -1 )
return maxSum[i][j]; 
if(i==n) 
maxSum[i][j] = D[i][j];
else {
int x = MaxSum(i+1,j); 
int y = MaxSum(i+1,j+1); 
maxSum[i][j] = max(x,y)+ 
D[i][j];
}
return maxSum[i][j];

int main(){ 
int i,j; 
cin >> n; 
for(i=1;i<=n;i++) 
for(j=1;j<=i;j++) {
cin >> D[i][j];
maxSum[i][j] = -1;
}
cout << MaxSum(1,1) << endl; 
}

递归转成递推

例如：
7
3 8
8 1 0
2 7 4 4
4 5 2 6 5

先把最后一行写进。
然后判断4上面的位置也就是原本是2的那个数，判断2+4大还是2+大 ，显然2+5大，所以4上面的位置填7.

继续判断，7+5大还是7+2大，显然7+5大，所以12的填到原本7的位置

以此类推到第一行
代码：

#include<iostream>
#include<algorithm>
using namespace std;
#define MAX 101;
int D[101][101]; int n;
int maxSum[101][101];
int main() {
    int i, j;
    cin >> n;
    for (i = 1; i <= n; i++)
        for (j = 1; j <= i; j++)
            cin >> D[i][j];
    for (int i = 1; i <= n; ++i)
        maxSum[n][i] = D[n][i];
    for (int i = n - 1; i >= 1; --i)
        for (int j = 1; j <= i; ++j)
            maxSum[i][j] =
            max(maxSum[i + 1][j], maxSum[i + 1][j + 1]) + D[i][j];
            cout << maxSum[1][1] << endl;
}

优化

没必要用二维maxSum数组存储每一个MaxSum(r,j),只要从底层一行行向上
递推，那么只要一维数组maxSum[100]即可,即只要存储一行的MaxSum值就
可以。
进一步考虑，连maxSum数组都可以不要，直接用D的 第n行替代maxSum即可。
节省空间，时间复杂度不变

代码

#include <iostream> 
#include <algorithm>
using namespace std;
#define MAX 101 
int D[MAX][MAX]; 
int n; int * maxSum;
int main(){ 
int i,j; 
cin >> n; 
for(i=1;i<=n;i++) 
for(j=1;j<=i;j++) 
cin >> D[i][j];
maxSum = D[n]; //maxSum指向第n行
for( int i = n-1; i>= 1; --i ) 
for( int j = 1; j <= i; ++j ) 
maxSum[j] = max(maxSum[j],maxSum[j+1]) + D[i][j];
cout << maxSum[1] << endl; 
}