给定一张 n 个点的带权无向图,点从 0∼n−1 标号,求起点 0 到终点 n−1的最短 Hamilton 路径。
Hamilton 路径的定义是从 0到 n−1不重不漏地经过每个点恰好一次。
输入格式
第一行输入整数 n。接下来 n行每行 n 个整数,其中第 i 行第 j 个整数表示点 i 到 j 的距离(记为 a[i,j])。
对于任意的 x,y,z,数据保证 a[x,x]=0,a[x,y]=a[y,x] 并且 a[x,y] + a[y,z] ≥ a[x,z]。
输出格式
输出一个整数,表示最短 Hamilton 路径的长度。
数据范围
1≤n≤20
0≤a[i,j]≤107
输入样例:
5
0 2 4 5 1
2 0 6 5 3
4 6 0 8 3
5 5 8 0 5
1 3 3 5 0
输出样例:
18
本题直接输入的是一个邻接矩阵,直接用数组存储就行,接下来考虑如何表达遍历n个点的的状态,我们最先考虑的是用一个bool数组来储存对应点的状态,但是题目要求我们要将所有的点都遍历一遍,那我们还需要判断bool数组是否都是true,在这里我们直接使用状态dp,将每个点的遍历状态转换成对应的二进制代码,如走了 0 1 3 4这4个点,对应的状态就是11011,ok,这样我们的遍历状态就确定了,使用一个二进制码来表示遍历状态。
这里借用闫式dp分析法(具体内容可以去acwing上查看)
(图片来源AcWing 91. 最短Hamilton路径 - AcWing)
i表示的是遍历状态,j表示遍历到哪个点,f[i][j]的定义如图。由邻接矩阵知,图中的每一个点之间都有一条边,我们可以从一个点走到任意一个点,状态计算一般是由前一个状态来计算现在的状态,如现在我们遍历到j点,我们可以从1到n这些点到j点,我们只需要找出这些点中f[state][k]最小的就行了。
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N = 20;
int f[1 << N][N];
int g[N][N];
int n;
int main()
{
cin >> n;
for(int i = 0; i < n; i ++)
for(int j = 0; j < n; j ++)
scanf("%d", &g[i][j]);
//遍历第一个点,状态为0...001,也就是1;因为要找最小值,所以将所有f[i][j]都初始化为无穷大
memset(f, 0x3f, sizeof f);
f[1][0] = 0;
for(int state = 0; state < 1 << n; state ++)//遍历每一种状态
for(int j = 0; j < n; j ++)//表示走到j点
if(state >> j & 1)//如果这个点被走过
for(int k = 0; k < n; k ++)//表示走到j点之前,k点为路径的最终点
f[state][j] = min(f[state][j], f[state - (1 << j)][k] + g[k][j]);//state - (1 << j)表示将j点从路径中移除
cout << f[(1 << n) - 1][n - 1];
return 0;
}