Project Euler 83

发表于 更新于 分类于 阅读次数： 23 本文字数： 1.1k 阅读时长 ≈ 1 分钟

Project Euler 83

题目

Path sum: four ways

NOTE: This problem is a significantly more challenging version of Problem 81.

In the 5 by 5 matrix below, the minimal path sum from the top left to the bottom right, by moving left, right, up, and down, is indicated in bold red and is equal to 2297.

(131673234103182019634296515063080374642211153769949712195680573252437331)

Find the minimal path sum from the top left to the bottom right by moving left, right, up, and down in matrix.txt (right click and "Save Link/Target As..."), a 31K text file containing an 80 by 80 matrix.

解决方案

本题不再适用于动态规划的做法，因为不能够很好地分阶段。

注意到，本题求的是最小值。因此，可以将整个模型转化为图论的单源最短路去完成。

令矩阵大小为 n。把矩阵中的每一个元素视为一个节点，那么就是求节点 (1,1) 到节点 (n,n) 的最短路径。注意到这里使用的是有向图，边集如下构造： 对于矩阵中所有元素，都有 4 条出边指向其四个方向的相邻元素，边权为指向的元素的值。

最终，计算出 (1,1) 到 (n,n) 的最短路径长度再加上 (1,1) 元素本身的值，就是最终答案。

单源最短路以 Dijkstra，Bellman-Ford 为代表。多源最短路则以 Floyd 为代表。

本代码实现主要使用 networkx 库利用 Dijkstra 算法导出结果。

代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
import networkx as nx

ls = open('p081_matrix.txt', 'r').readlines()
n = len(ls)
a = []
for i in range(n):
    b = [int(x) for x in ls[i][:-1].split(',')]
    a.append(b)
g = nx.DiGraph()
d = [[-1, 0], [1, 0], [0, -1], [0, 1]]

for i in range(n):
    for j in range(n):
        for dx, dy in d:
            x, y = i + dx, j + dy
            if 0 <= x < n and 0 <= y < n:
                g.add_edge((i, j), (x, y), weight=a[x][y])
ans = nx.dijkstra_path_length(g, (0, 0), (n - 1, n - 1)) + a[0][0]
print(ans)