哈密顿路径 状压DP

对于某类问题， 求解经过所有点的最短路径（哈密顿路径）或者求经过所有点最短路径后再回到开始的地方（哈密顿回路）例如旅行商问题（TSP）。

首先我们假设 每个点与点之间的权值(路的长度)已经求出 记为 w[i][j].
（对于这个部分有些时候需要进行处理，自己去把权值求出来）

我们把每一个点的状态压缩到应该二进制的位 1代表走这个点，0代表没有走这个点。

对于某一个状态，我们可以得到它可以由哪一个状态转移过来。
（设有7个点）例如： 当前我们得到了一个状态 1000001 表示第1个点和第7个点走过，
而它的上一个状态可以为0000001，1000000推出。
0000001 -> 1000001 说明此时是走了 7-> 1 这条边。

我们设dp[i][j] 表示当前第i个状态且正在访问第j个节点的最小花费。
那么上述结果为dp[1000001][1] = min(dp[0000001][7]+w[7][1],dp[1000001][1]);

可以得到主要步骤
1 ： 枚举一个状态从 1-（1<<n）。
2 ： 找到上一个走的点。
3 ： 枚举上一个状态
4 ： 更新当前状态的最小值
5： 最后输出一个dp[(1<<n)-1][…]或者枚举这个状态的所有点输出最小值或者再加上w[…][1]的权值计算哈密顿回路 （看具体情况）。

因为之后右移运算符的原因，我们让dp的第二维的大小减1 即 1号点为 0 也可以不用，看个人喜好。
一个例题： swustoj

题目说经过每个点一次，且权值已经给好，我们直接套用上述步骤即可。

详细解释看代码：

#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#include<bits/stdc++.h>
#define int long long
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;
 
const int INF = 0x3f3f3f3f;
const double eps = 1e-4;
const int mod = 999911659;
const int N = 25;

int dp[1<<20][N],a[N][N];

signed main(){
	IOS;
	#ifdef ddgo
		freopen("C:/Users/asus/Desktop/ddgoin.txt","r",stdin);
	#endif
	
	int n; cin>>n;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			cin>>a[i][j];
	//初始化，好像太大了memset比较慢，所以直接用的 for 循环
	for(int i=0;i<(1<<20);i++)
		for(int j=0;j<25;j++)
			dp[i][j] = INF;
	// 从第一个点开始，所以我们把 dp[1][0] 表示1(0000000000000000...0001)这个状态且当前
	//正在走1(1-1=0)这个节点初始化为0
	dp[1][0] = 0;
	for(int i=0;i< 1<<n;i++)//枚举所有状态
		for(int j=0;j<n;j++)// 找到当前状态的终点，可以理解为当前这个位置是1，但上一个状态这个位置是0
			if(i>>j & 1)
				for(int k=0;k<n;k++)//找上一个状态 
					if((i-(1<<j))>>k & 1)//上一个状态中为1的点都可以为起点。
						dp[i][j] =min(dp[i][j],dp[i-(1<<j)][k]+a[k][j]);
	int mi = INF;
	// 哈密顿回路再更新一下最小值。
	for(int i=0;i<n;i++) mi = min(dp[(1<<n)-1][i]+a[i][0],mi);
	cout<<mi<<endl;
    return 0;
}

题目地址：hdu5418
题目说经过点至少一次，所以用floyd处理每个点与点之间的最短距离

#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#include<bits/stdc++.h>
#define int long long
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;
 
const int INF = 0x3f3f3f3f;
const double eps = 1e-4;
const int mod = 999911659;
const int N = 20;

int dp[1<<16][N],a[N][N];

signed main(){
	IOS;
	#ifdef ddgo
		freopen("C:/Users/asus/Desktop/ddgoin.txt","r",stdin);
	#endif
	int tt; cin>>tt;
	while(tt --){
		memset(a,INF,sizeof(a));
		int n,m; cin>>n>>m;
		for(int i=0;i<m;i++){
			int u,v,w; cin>>u>>v>>w;
			a[u-1][v-1] = min(a[u-1][v-1],w);
			a[v-1][u-1] = min(a[v-1][u-1],w);
		}
		
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				for(int k=0;k<n;k++)
					a[i][j] = min(a[i][j],a[i][k]+a[k][j]);
		
		for(int i=0;i<(1<<16);i++)
			for(int j=0;j<20;j++)
				dp[i][j] = INF;
		dp[1][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)
							dp[i][j] =min(dp[i][j],dp[i-(1<<j)][k]+a[k][j]);
		int mi = INF;
		a[0][0] = 0; // 可以是直接到达起点 
		for(int i=0;i<n;i++) mi = min(dp[(1<<n)-1][i]+a[i][0],mi);
		cout<<mi<<endl;
	}
    return 0;
}