状态压缩动态规划：最短Hamilton路

题目描述

给定一张 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\le n\le 20$
$0\le a[i,j]\le 10^7$

输入样例

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

算法思想（状压DP）

状态表示

f[s][j]表示从0到j点所有经过点为集合s的二进制形式的最短Hamilton路径的长度。数s的二进制形式表示点0 ~ n - 1每个点是否被走过
例如 s = 10，二进制表示为$(1010{)}_2$，表示1、3号点走过，0、2号点没有走过。

状态计算

考虑要走到j点之前，要经过的前一个点（设其编号为i），将集合分为n份，取最小值。那么
f [ s ] [ j ] = m i n { f [ s − { j } ] [ i ] } + a [ i ] [ j ] , i ∈ [ 0 , n − 1 ] f[s][j] = min\{f[s-\{j\}][i]\} + a[i][j], i ∈ [0,n-1] f[s][j]=min{f[s−{j}][i]}+a[i][j],i∈[0,n−1]

初始状态

因为求最小值，所以首先将f数组初始化为正无穷。其次，求从起点0出发终点 n-1 的最短Hamilton路径，所以经过编号为0的点并且到达0点，此时的最短距离为0，即f[1][0] = 0。

时间复杂度

状态数$2^n\times n$，转移过程要循环$n$次，所以时间复杂度为$O(2^n\times n^2)$

代码实现

#include <iostream>
#include <cstring>
using namespace std;
const int N = 20, M = 1 << 20, INF = 0x3f3f3f3f;
int g[N][N], f[M][N];
int main()
{
    int n;
    cin >> n;
    for(int i = 0; i < n; i ++)
        for(int j = 0; j < n; j ++)
            cin >> g[i][j];
    
    //求最小值，所以将f数组初始化无穷大      
    memset(f, 0x3f, sizeof f);
    //从起点0出发到达0时，最短Hamilton路径的长度为0
    f[1][0] = 0;
    
    //枚举所有子集状态
    for(int s = 0; s < (1 << n); s ++)
    {
        //枚举当前到达的点i
        for(int i = 0; i < n; i ++)
        {
            //如果i点包含在集合中
            if((s >> i) & 1)
            {
            	//枚举上一次到达的点j
                for(int j = 0; j < n; j ++)
                {
                    //i和j都在集合中，并且i和j不相同
                    if((s >> j & 1) && i != j)
                    {
                        f[s][i] = min(f[s][i], f[s - (1 << i)][j] + g[i][j]);
                    }
                }
            }
        }
    }    
    cout << f[(1 << n) - 1][n - 1] << endl;    
    return 0;
}