数塔问题及其变形

数塔（一）

如下图所示的一个数字三角形宝塔。数字三角形中的数字为不超过 100 的正整数。现规定从最顶层走到最底层，每一步可沿左斜线向下或右斜线向下走。
假设三角形行数 n≤100，编程求解从最顶层走到最底层的一条路径，使得沿着该路径所经过的数字的总和最大，输出最大值。

基本上每个人学习动态规划的时候都做过这题。
转移方程也比较简单。（滚动数组优化）

#include <cstdio>
#include <cstring>
#define Max(a, b) ((a>b) ? a : b)
using namespace std;
const int maxn=1010;
int opt[2][maxn], n;

void work() {
   int i, j, t, pre=1, now=0;
   for (i=1; i<=n; i++) {
      pre=now; now=1-pre;
      for (j=1; j<=i; j++) {
         scanf("%d", &t);
         opt[now][j]=Max(opt[pre][j], opt[pre][j-1])+t;
      }
   }

   int ans=0;
   for (i=1; i<=n; i++) ans=Max(ans, opt[now][i]);
   printf("%d\n", ans);
}

int main() {
   scanf("%d", &n);
   memset(opt, 0, sizeof(opt));
   work();
   return 0;
}

数塔（2）

要求走到最后mod m最大

注意这里是走到最后 mod m 最大，因此不能简单的求出最大值 mod m。这样并不满足最优子问题。因此这里需要在原来的转移方程上加上一维，$dp[i][j][m]$表示在i, j位置的时候余数 m 是否存在。这是存在性DP问题

#include <iostream>
#include <algorithm>
#include <string.h>
#include <stdio.h>
using namespace std;
#define ll long long
const int maxn = 205;
int a[maxn][maxn];
bool dp[maxn][maxn][maxn];
int main(int argc, char const *argv[])
{
    int n, k;
    scanf("%d %d", &n, &k);
    memset(dp, false, sizeof dp);

    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= i; j++)
            scanf("%d", &a[i][j]);
    dp[1][1][a[1][1]%k] = true;
    int maxa = 0;
    for (int i = 1; i < n; i++) {
        for (int j = 1; j <= i; j++)
            for (int m = 0; m < k; m++) {
                if (dp[i][j][m]) {
                    dp[i+1][j][(m+a[i+1][j])%k] = true;
                    dp[i+1][j+1][(m+a[i+1][j+1])%k] = true;
                }
            }
    }

    for (int i = 1; i <= n; i++)
        for (int m = 0; m < k; m++) {
            if (dp[n][i][m])
                maxa = max(maxa, m);
        }
    cout << maxa << endl;

    return 0;
}

数塔（三）

要求必须经过其中的某一个点，使之走的路程和最大。

这里要求我们必须经过某一个点，只需要将那一个点加上无穷大∞，那么必定会经过那个点，最后输出的时候减去无穷大。

#include <iostream>
#include <algorithm>
#include <string.h>
#include <stdio.h>
using namespace std;
#define ll long long
const int maxn = 1005;
const int inf = 0x3f3f3f3f;
int a[maxn][maxn], dp[maxn][maxn];
int main(int argc, char const *argv[])
{
    int n;
    scanf("%d", &n);
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= i; j++) {
            scanf("%d", &a[i][j]);
            dp[i][j] = a[i][j];
        }
    int x, y;
    scanf("%d %d", &x, &y);
    a[x][y] += inf;
    dp[x][y] = a[x][y];

    for (int i = n-1; i >= 1; i--)
        for (int j = 1; j <= i; j++)
            dp[i][j] = max(dp[i+1][j], dp[i+1][j+1])+a[i][j];
    int maxa = 0;
    
    printf("%d\n", dp[1][1]-inf);


    return 0;
}