首篇先了解一下什么是动态规划,动态规划(Dynamic Programming)俗称dp,是一种解决多阶段决策问题的优化方法。它通过将问题分解为多个重叠子问题,并利用子问题的解来构建原问题的解。动态规划算法通常用于具有重复子问题的优化问题。
那什么是数字三角形模型,其实就是二维的类三角形的图,通过根据已推出的量继续向后推,直到求出我们需要的路线或者花费之类的,我们所需要考虑的就是当前点能由哪些点推来即可写出状态转移方程,最后按照dp的老套路,一般分为下面几个步骤进行即可:
- 确定dp数组(dp table)以及下标的含义
- 确定递推公式
- 初始化dp数组
- 确定遍历顺序
- 举例推导dp数组
Hello Kitty想摘点花生送给她喜欢的米老鼠。
她来到一片有网格状道路的矩形花生地(如下图),从西北角进去,东南角出来。
地里每个道路的交叉点上都有种着一株花生苗,上面有若干颗花生,经过一株花生苗就能摘走该它上面所有的花生。
Hello Kitty只能向东或向南走,不能向西或向北走。
问Hello Kitty最多能够摘到多少颗花生。
输入格式
第一行是一个整数T,代表一共有多少组数据。
接下来是T组数据。
每组数据的第一行是两个整数,分别代表花生苗的行数R和列数 C。
每组数据的接下来R行数据,从北向南依次描述每行花生苗的情况。每行数据有C个整数,按从西向东的顺序描述了该行每株花生苗上的花生数目M。
输出格式
对每组输入数据,输出一行,内容为Hello Kitty能摘到得最多的花生颗数。
数据范围
1≤T≤100,
1≤R,C≤100,
0≤M≤1000
输入样例:
2
2 2
1 1
3 4
2 3
2 3 4
1 6 5
输出样例:
8
16
思路:题目要求我们从左上角走到右下角,并且只能向右和向下,求出路上经过点的和的最大值。由于行进方向被确定了,因此对于任意一个除起点外的点(i, j),它只能通过(i - 1, j),或者(i, j - 1),推来,因此可以得出我们的状态转移方程dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]) ,然后再加上这个点的花生数即可,即dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]) + a[i][j]。数组默认值为0,所以初始化直接从(1, 1)开始,最后dp[n][m]即为我们所要的结果
#include <iostream> #include <cstring> #include <algorithm> using namespace std; const int N = 110; int T; int n, m; int dp[N][N]; int main() { scanf("%d", &T); while (T -- ) { scanf("%d%d", &n, &m); int a[n + 1][m + 1]; for(int i = 1; i <= n; i ++ ) for(int j = 1; j <= m; j ++ ) cin >> a[i][j]; for(int i = 1; i <= n; i ++ ) for(int j = 1; j <= m; j ++ ) dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]) + a[i][j]; printf("%d\n", dp[n][m]); } return 0; }
一个商人穿过一个 N×N 的正方形的网格,去参加一个非常重要的商务活动。
他要从网格的左上角进,右下角出。
每穿越中间 11 个小方格,都要花费 11 个单位时间。
商人必须在 (2N−1) 个单位时间穿越出去。
而在经过中间的每个小方格时,都需要缴纳一定的费用。
这个商人期望在规定时间内用最少费用穿越出去。
请问至少需要多少费用?
注意:不能对角穿越各个小方格(即,只能向上下左右四个方向移动且不能离开网格)。
输入格式
第一行是一个整数,表示正方形的宽度 N。
后面 N 行,每行 N 个不大于 100100 的正整数,为网格上每个小方格的费用。
输出格式
输出一个整数,表示至少需要的费用。
数据范围
1≤N≤100
输入样例:
5
1 4 6 8 10
2 5 7 15 17
6 8 9 18 20
10 11 12 19 21
20 23 25 29 33
输出样例:
109
样例解释
样例中,最小值为 109=1+2+5+7+9+12+19+21+33109=1+2+5+7+9+12+19+21+33。
思路:这题和上题思路类似,只不过这题要求的是最小费用,初始化和遍历顺序方面有所不同。题目给的方向说是四个方向都可以,但是又给了时间限制,要求在2 * N - 1内走到右下角,可以发现我们只能向右和向下,因为从左上角到右下角至少需要2 * N - 2步,所以如果向左或者向上走了就一定会超时。
#include <iostream> #include <cstring> #include <algorithm> using namespace std; const int N = 110; int n; int a[N][N]; int dp[N][N]; int main() { scanf("%d", &n); for(int i = 1; i <= n; i ++ ) for(int j = 1; j <= n; j ++ ) scanf("%d", &a[i][j]); memset(dp, 0x3f, sizeof dp); for(int i = 1; i <= n; i ++ ) for(int j = 1; j <= n; j ++ ) { if(i + j == 2) dp[i][j] = a[i][j]; else if(i == 1) dp[i][j] = dp[i][j - 1] + a[i][j]; else if(j == 1) dp[i][j] = dp[i - 1][j] + a[i][j]; else dp[i][j] = min(dp[i - 1][j], dp[i][j - 1]) + a[i][j]; } printf("%d\n", dp[n][n]); return 0; }