2020-08-30

                                                                   动态规划 数塔问题
数塔问题，又称数字三角形、数字金字塔问题。数塔问题是多维动态规划问题中一类常见且重要的题型，
一、数塔问题原型
1.1 问题描述
        7
      3   8
    8   1   0
  2   7   4   4
4   5   2   6   5
有一个r行的数塔，数塔上有若干数字。问从数塔的最高点到底部，在所有的路径中，经过的数字的和最大为多少？
如上图，是一个5行的数塔，其中7—3—8—7—5的路径经过数字和最大，为30。

如果使用贪心算法，则会选择相邻数组最大的数字；即7+8+1+7+5=28.但是该数塔的经过最大数字和为30.
如是使用深度搜索，数据很大时，时间复杂度非常高。
可以发现，要想经过一个数字，只能从左上角或右上角的数字往下到达。

所以显然，经过任一数字A时，路径所经过的数字最大和——是这数字A左上方的数字B以及右上方的数字C两个数字中，所能达到的数字最大和中较大的那一个，再加上该数字A。可以发现，要想经过一个数字，只能从左上角或右上角的数字往下到达。所以显然，经过任一数字A时，路径所经过的数字最大和——是这个数字A左上方的数字B以及右上方的数字C两个数字中，所能达到的数字最大和中较大的那一个，再加上数字A。
所以 定义数组 num和dp .num数组是数塔本塔。Dp数组是经过路线的 最大路线数字的和。所以满足方程
Dp[i][j]=max(dp[i-1][j],dp[i-1][j-1])+num[i][j];
其中设 DP[I-1][J]是左上角 dp[i-1][j-1]是右上角。
在经过第三行的数字1时，她3左上角的数字3和右上角的数字8其能达到的最大和。3显然只有7—3一条路径，故最大和是10，8显然也只有7—8一条路径，其最大和是15；两者中较大的是15，故经过1所能达到的最大和是15+1=16。
一行一行的想下遍历。经过每个位置的最大值都储存在dp数组中。寻找最大值输出即可。最后一行的Dp数组的最大值即为答案。
代码：
#include <stdio.h>
#include <algorithm>
using namespace std;
 int num[1005][1005];
 int dp[1005][1005];
int main()
{
    int r;
    scanf("%d",&r);
    for(int i=1;i<=r;i++)
        for(int j=1;j<=i;j++)
            scanf("%d",&num[i][j]);
        for(int i=1;i<=r;i++)
        for(int j=1;j<=i;j++)
            dp[i][j]=max(dp[i-1][j],dp[i-1][j-1])+num[i][j];
    //经过该数字的最大和，为左上角和右上角中的max，再加上该数字
    int ans=0;
    for(int i=1;i<=r;i++)
        ans=max(ans,dp[r][i]);
    printf("%d\n",ans);
    return 0;
}

2.1矩形数塔
 问题：
有一m行，n列（n为奇数）的数字矩阵（数字中存在部分负数）。
以最后一行的正中间下方为出发点，每次移动可以选择向前方、左前方、右前方移动，问从出发点一直到矩阵的另一侧，所经过的最大数字和为多少。
6 7
16 4 3 12 6 0 3
4 -5 6 7 0 0 2
6  0 -1 -2 3 6 8
5 3 4 0 0 -2 7
-1 7 4 0 7 -5 6
0 -1 3 4 12 4 2
       S、
S是出发点，第一次移动可以移动到3 4 12处。
从起始点移动到居住你的另一测，经过的最大数字是41.
相比于三角形的数塔，矩形数塔比三角形数塔多了一个移动的方向，所以dp的方程变为
Dp[i][j]=max(dp[i+1][j],dp[i+1][j-1],dp[[i+1][j+1]]+num[i][j];
从上倒下遍历。
因为
矩阵中存在负数，所以必须将dp初始化为非常小的负数。否则在使用MAX函数时，结果可能产生错误。

代码：
#include <bits/stdc++.h>
using namespace std;
int num[205][205];
int dp[205][205];
int main()
{
    int m,n;
    while(scanf("%d%d",&m,&n)!=EOF)
    {
        memset(dp,-9999,sizeof(dp));//初始化dp数组为绝对值较大的负数
        for(int i=1;i<=m;i++)
            for(int j=1;j<=n;j++)
            scanf("%d",&num[i][j])
        dp[m][n/2+1]=num[m][n/2+1];
        dp[m][n/2]=num[m][n/2];
        dp[m][n/2+2]=num[m][n/2+2];
        for(int i=m-1;i>=1;i--)
            for(int j=1;j<=n;j++)
            dp[i][j]=max(max(dp[i+1][j+1],dp[i+1][j-1]),dp[i+1][j])+num[i][j];
        int ans=0;
        for(int i=1;i<=n;i++)
            if(dp[1][i]>ans) ans=dp[1][i];        printf("%d\n",ans);
    }
    return 0;
}

电脑应该明天就能修好，下次一定。