AcWing 91 最短Hamilton路径

哈密顿路径，恰好经过所有顶点一次。比如n=4，从0-3，起点终点固定，有0-1-2-3和0-2-1-3两种路径选择，我们选择其中长度最短的一条就可以了，n最大为20，从0-19，状态数有18！这么多，暴力求解不能承受。0-n-1的中间状态一共2^n种，采用状态压缩，即用一串n位的二进制整数表示中间状态，为了确定dp数组的每一种状态，需要两个参数，其一是哪些顶点被选择了，其二是当前处在哪个顶点。状态转移方程为f[state][j] = f[state_k][k] + v[k][j]. f[state][j]表示当前已选择顶点的情况为state且目前处在顶点j的状态，f[state_k][k]表示当前在k顶点且从原来的顶点集中除去顶点k的情况。

代码设计的注意点：其一，memset(f,0x3f,sizeof f);sizeof运算符优先级挺高，不用加括号也可，还有一般用memset给每个字节填充0x3f也就是0011 1111，表示INF，相加不会溢出；其二，状态转移方程化为代码。一共2^n乘以n种状态，头两层循环表示这么多种状态，我们要做的就是从f[1][0]转移到f[2^n-1][n-1]。需要注意的是，上面这么多种状态有些状态是不用扩散的，比如f[state][j] j不在state中就不用考虑。还有就是状态转移方程：f[i][j] = min(f[i][j],f[i-(1<<j)][k] + v[k][j]) 从依次从i中遍历顶点作为k，但是得先去掉当前顶点j，因为j=k时相当于什么都没去掉。

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 20,M = 1 << 20;
int n;
int f[M][N],v[N][N];
int main(){
	cin>>n;
	for(int i = 0;i < n;i++){
		for(int j = 0;j < n;j++){
			cin>>v[i][j];
		}
	}
	memset(f,0x3f,sizeof f);
	f[1][0] = 0;
	for(int i = 1;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] + v[k][j]);
					}
				}
			}
		}
	}
	cout<<f[(1<<n)-1][n-1];
	return 0;
}