100道动态规划10—— SCU 2842 The TSP Problem 记忆化搜索 状态压缩

这就是传说中的状态压缩的题了吧！虽然说在暑假集训的时候做过一点。。但是当时完全没有感觉。。现在再来看一看，还是有点感觉了的。

TSP问题，恩，旅行商问题，算法导论上也提到过，这个问题目前还没有多项式时间的解，我之前说过一道双调巡游的题目就是在TSP问题上加以限制，才能利用动态规划在多项式时间内解出来。假若不加以限制的话，那么利用动态规划就是指数级的时间复杂度了。

分析一下题目吧

因为我们的线路构成一个环，因此我们以任意一个点为起点都是可以的，因此无妨以0号城市为起点

定义状态dp[i][s]，其中s是一个二进制数，代表城市的集合，表示当前在第i个城市，还需要访问s里面的城市各一次之后的最小值，因此状态转移方程就是

dp[i][s]=min(dis(i,j)+dp[j][s-{j}])for all j ∈ s

Tips：本题用的是记忆化搜索，感觉在这种状态定义下，自底向上感觉有点不好写。。

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

int board[25][25],n,dp[25][1<<21],t,_find(int now,int s),limi;

int main(){
	scanf("%d",&t);
	while(t--){
		memset(dp,-1,sizeof dp);
		scanf("%d",&n);

		for(int i=0;i<n;limi|=(1<<i++))
		for(int j=0;j<n;++j)
			scanf("%d",&board[i][j]);

		--limi;
		for(int i=0;i<n;++i)
			dp[i][0]=board[i][0];
		printf("%d\n",_find(0,limi));
		limi=0;
	}
	return 0;
}

int _find(int now,int s){
	if(dp[now][s]!=-1)
		return dp[now][s];

	dp[now][s]=0x3f3f3f3f;
	for(int i=0;i<n;++i)
	if(s&(1<<i))
		dp[now][s]=min(dp[now][s],board[now][i]+_find(i,s^(1<<i)));

	return dp[now][s];
}