洛谷P1004 [NOIP2000 提高组] 方格取数题解

本文介绍了如何使用四维动态规划解决一个NOIP2000提高组的竞赛题目,该题要求找到两条从方格图左上角到右下角的路径,使途中取到的数字和最大化。状态转移方程为dp[i][j][k][l]=max(dp[i-1][j][k-1][l],dp[i][j-1][k][l-1],dp[i-1][j][k][l-1],dp[i][j-1][k-1][l])+a[i][j]+a[k][l],其中a[i][j]表示点(i,j)的数值。通过四重循环计算,可以找到最佳路径组合,例如样例中给出的最优解为67。" 109480555,10261843,C# 实践:字符提取与整数整除方法,"['C#编程', '整型运算', '字符串处理']

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

四维动态规划,属于是模版题

题目描述

设有 N × N 的方格图 ( N ≤ 9 ) ,我们将其中的某些方格中填入正整数,而其他的方格中则放入数字0。如下图所示(见样例):

A
 0  0  0  0  0  0  0  0
 0  0 13  0  0  6  0  0
 0  0  0  0  7  0  0  0
 0  0  0 14  0  0  0  0
 0 21  0  0  0  4  0  0
 0  0 15  0  0  0  0  0
 0 14  0  0  0  0  0  0
 0  0  0  0  0  0  0  0
                         B

某人从图的左上角的A点出发,可以向下行走,也可以向右走,直到到达右下角的B点。在走过的路上,他可以取走方格中的数(取走后的方格中将变为数字0)。

此人从A点到B点共走两次,试找出2条这样的路径,使得取得的数之和为最大。

输入格式

输入的第一行为一个整数 N(表示 N × N 的方格图),接下来的每行有三个整数,前两个表示位置,第三个数为该位置上所放的。一行单独的0表示输入结束。

样例输入:
8
2 3 13
2 6  6
3 5  7
4 4 14
5 2 21
5 6  4
6 3 15
7 2 14
0 0  0
输出格式

只需输出一个整数,表示2条路径上取得的最大的和。

推状态转移方程

这是道很明显的多维DP题目。因为是2条路径,所以同时需要dp两个点,使用四维数组。

定义一个四维数组,dp[ i ][ j ][ k ][ l ]的值表示第一条路径到点( i , j ),第二条路径到点( k , l )时所取的最大数之和,
因为只能向下或向右行走,于是就能得到状态转移方程:

dp[ i ][ j ][ k ][ l ] = max(dp[ i-1 ][ j ][ k-1 ][ l ],dp[ i ][ j-1 ][ k ][ l-1 ],dp[ i-1 ][ j ][ k ][ l-1 ],dp[ i ][ j-1 ][ k-1 ][ l ]) + a[ i ][ j ]+a[ k ][ l ]

PS:a[ i ][ j ]代表点( i , j )上的数字大小

这里注意,第一条路径经过后会将格子中的数拿走,当i=k且j=l时,因为已经被第一条路径取走了,不需要再增加a[ k ][ l ]的值

蒟个栗子

还是看样例:

A 
 0  0  0  0  0  0  0  0
 0  0 13  0  0  6  0  0
 0  0  0  0  7  0  0  0
 0  0  0 14  0  0  0  0
 0 21  0  0  0  4  0  0
 0  0 15  0  0  0  0  0
 0 14  0  0  0  0  0  0
 0  0  0  0  0  0  0  0
                         B

通过四重循环,我们可以得到两条数字和最大的路径分别为,穿过13-14-4和穿过21-15,总和为67。

代码实现看这里

#include 
### NOIP 2000 提高 方格问题分析 该问题属于经典的动态规划类题目,目标是在给定的一个二维矩阵中选若干个不相邻的字使得总和最大。此问题可以通过状态转移的方式解决。 #### 动态规划的核心思路 定义 `dp[i][j]` 表示到达第 `i` 行第 `j` 列时能够得的最大值之和[^1]。由于每次移动仅能向右或者向下,因此可以得出如下状态转移方程: 对于任意位置 `(i, j)` 的状态更新方式为: \[ dp[i][j] = \text{matrix}[i][j] + \max(dp[i-1][j], dp[i][j-1]) \] 其中需要注意边界条件处理以及初始值设定。当处于第一行或第一列时,路径的选择会受到限制[^2]。 #### 边界情况考虑 如果当前单元位于网格的第一行,则只能由左侧进入;同理,若当前位置处在首列,则唯一可能来自上方。这些特殊情况需单独初始化以确保计算准确性[^3]。 ```python def max_sum(matrix): m, n = len(matrix), len(matrix[0]) # 初始化 DP dp = [[0]*n for _ in range(m)] # 设置起点 dp[0][0] = matrix[0][0] # 填充第一行 for j in range(1, n): dp[0][j] = dp[0][j-1] + matrix[0][j] # 填充第一列 for i in range(1, m): dp[i][0] = dp[i-1][0] + matrix[i][0] # 完成剩余部分填充 for i in range(1,m): for j in range(1,n): dp[i][j] = max(dp[i-1][j], dp[i][j-1]) + matrix[i][j] return dp[-1][-1] ``` 以上代码片段展示了如何通过构建辅助来记录每一步的最佳决策过程,并最终返回全局最优解。 #### 复杂度评估 时间复杂度主要决于遍历整个输入矩阵所需的操作次,即 O(M*N),这里 M 和 N 分别代表矩阵的高度与宽度。空间复杂度同样为 O(M*N),因为我们需要额外的空间存储中间结果以便后续访问[^4]。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值