最短Hamilton路径(状压dp)

链接：https://ac.nowcoder.com/acm/problem/50909
来源：牛客网

时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 262144K，其他语言524288K
64bit IO Format: %lld
题目描述
给定一张 n(n \leq 20)(n≤20) 个点的带权无向图，点从0 \sim n-10∼n−1标号，求起点 0 到终点 n-1 的最短Hamilton路径。 Hamilton路径的定义是从 0 到 n-1 不重不漏地经过每个点恰好一次。
输入描述:
第一行一个整数n。
接下来n行每行n个整数，其中第i行第j个整数表示点i到j的距离（一个不超过10^710
7
的正整数，记为a[i,j]）。
对于任意的x,y,z，数据保证 a[x,x]=0，a[x,y]=a[y,x] 并且a[x,y]+a[y,z] \geq a[x,z]a[x,y]+a[y,z]≥a[x,z]。
输出描述:
一个整数，表示最短Hamilton路径的长度。
示例1
输入
复制

4
0 2 1 3
2 0 2 1
1 2 0 1
3 1 1 0

输出
复制

4

说明
从0到3的Hamilton路径有两条，0-1-2-3和0-2-1-3。前者的长度为2+2+1=5，后者的长度为1+2+1=4
/*
状压dp:
用 dp[status][k] 表示 经过status这些点，终点是k的最短路。
dp[status][j] = min(dp[status][j],dp[status-(1<<j)][k] + a[k][j])
dp[status-(1<<j)][k] + a[k][j]:
通过k到达j，所以之前j应该为0，
即status减去j状态（status^(1<<j)也是相同的意思）。

*/

#include <bits/stdc++.h>
using namespace std;
int a[21][21];
int dp[1<<20][21];
int main()
{
    int n;
    cin>>n;
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < n; j++)
        {
            cin>>a[i][j];
        }
    }
    int inf = 0x3f3f3f3f;
    memset(dp,inf,sizeof(dp));
    dp[1][0] = 0; //边界：状态1，终点为0，路径长度为0
    for(int i = 1; i <(1<<n); i++)
    {
        for(int j = 0; j < n; j++)
        {
            if((i>>j)&1)//status中终点为j，此时想再通过其他点到j
                for(int k = 0; k < n; k++)
                {
                    dp[i][j] = min(dp[i][j],dp[i^(1<<j)][k] + a[k][j]);
                }
        }
    }
    cout<<dp[(1<<n)-1][n-1]<<endl;
    return 0;
}