状态压缩位运算之最短Hamilton路径

最短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≤n≤20,0≤a[i,j]≤10^7

分析都在代码里了。

代码：

//二进制0和1表示存在不存在，状态压缩
#include<iostream>
#include<cstring>
#include<algorithm>

using namespace std;

const int N=20,M=1<<20;

int n;
int f[M][N],weight[N][N];

int main()
{
    cin>>n;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            cin>>weight[i][j];
    
    memset(f,0x3f,sizeof f);//(初始化为正无穷)地址f，大小sizeof f,初始化为0x3f
    
    f[1][0]=0;//已经走过的点有0，现在在0位置
    
    for(int i=0;i<1<<n;i++)
        for(int j=0;j<n;j++)
            if(i>>j&1)
                for(int k=0;k<n;k++)
                    if(i-(1<<j)>>k&1)
                        f[i][j]=min(f[i][j],f[i-(1<<j)][k]+weight[k][j]);//f[i][j]=min(f[i][j],f[这里第一次写错了][k]+weight[k][j])
    
    cout<<f[(1<<n)-1][n-1]<<endl;//for中定义的i在这里不能用所以写成1<<n，最后还忘了加括号，减号优先级高于移位运算符
    
    return 0;
}