AcWing 91. 最短Hamilton路径（状压dp+位运算）

AcWing 91. 最短Hamilton路径（状压dp+位运算）

给定一张 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⁷

输入样例：
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-2）！种方案，复杂度是很大的，会超时，所以这里我们用状压dp来解决这个问题

首先假设有5个点0，1，2，3，4，我们要计算0 -> 4(每个点都遍历到了)的最短距离，那我们就要比较0 -> 1（除4都遍历到了，后面同理），0 -> 2,0 -> 3这几个距离那个最短，要计算这几个距离最短，以计算0 -> 1怎样走最短为例，那我们就需要计算0 -> 2(除1，4都遍历到了)，0 -> 3，那个最短，以此类推，就是我们计算最短距离的方法

由此得出，要计算0 -> j 的最短距离，要通过 0-> k 的最短距离+ a[k][j] （k到j的距离）得出，此时就必须保证0 -> j 必须经过k

定义dp[i][j]表示从0到当前点j经过i（i为二进制表示的一个集合，1表示经过这个点）的最短距离, 那状态转移方程就为

	dp[i][j]=min(dp[i][j],dp[i-(1<<j)][k]+a[k][j]);
	//1<<j表示2^j ，i-(1<<j)表示从i中删除j，即不经过j

AC代码：

#include<bits/stdc++.h>
using namespace std;
const int N=1<<20,M=30;
int dp[N][M],a[M][M];
//dp[i][j]表示从0走到j经过了i的最短路径
//i是二进制表示的一个集合，1表示经过这个点
//如1001表示，经过0，3这两个点 
int main()
{
	int i,j,k,n;
	cin>>n;
	for(i=0;i<n;i++)
	for(j=0;j<n;j++)
	cin>>a[i][j];
	memset(dp,0x3f3f3f3f,sizeof(dp));//求最小值，初始最大化
	dp[1][0]=0;//初始状态从0出发，当前在0这个点，且只经过0这个点，则i为00..1，为2^0=1 
	for(i=0;i<(1<<n);i++)//枚举2^n种情况
	{
		for(j=0;j<n;j++)//当前在那个点 
		{
			if(i>>j&1)//i的第j位为1，即j在i中 
			{
				for(k=0;k<n;k++)//k表示在走到j之前，以k为终点的最短距离
				{
					if(i>>k&1)//i的第k位为1，即k在i中
					dp[i][j]=min(dp[i][j],dp[i-(1<<j)][k]+a[k][j]);//后者表示从0到k的最短距离+从k到j的距离 
				 } 
			 } 
		}
	 } 
	 cout<<dp[(1<<n)-1][n-1]<<"\n";//(1<<n)-1表示所有点都走过了，即11..1,（n-1个1），为2^n-1; 
	return 0;
}