poj 3311 dp求解TSP

题意：一个开披萨店的要把披萨所购所有购买者。一直披萨店在位置0，n个购买者分别在位置1-n。题目给出这n+1个点之间的直线距离，问送货员从披萨店出发，送给所有的购买者且最后回到披萨店所需要的最小距离。

思路：现用floyd处理一下给定的图，剩下的就是TSP问题了。题目给的数据范围比较小（n<=10），可以考虑动态规划方法。目前为止都送给了哪些购买者用一个状态（一个数）表示，相应位为1表示经过过，为0表示未经过。dp[ i ][ j ]表示状态i^(1<<j) (也就是从状态 i 中除去 j 城市的那个状态)到城市 j 的最短距离。状态转移方程也就是dp[ i ][ j ] = min(dp[ i ] [ j ] , dp[ pre ][ k ]+dis[ k ][ j ])。（其中pre选择为从状态 i 中除去 j 城市的那个状态）

#include <stdio.h>
#include <string.h>
#define INF 0x3fffffff
#define min(a,b) ((a)<(b)?(a):(b))
#define N 15
int n;
int g[N][N],dp[1<<11][N];
void floyd(){
	int i,j,k;
	for(k = 0;k<=n;k++)
		for(i = 0;i<=n;i++)
			for(j = 0;j<=n;j++)
				g[i][j]=min(g[i][j],g[i][k]+g[k][j]);
}
void solve(){
	int i,j,k,m,res = INF;
	m = 1<<n;
	for(i = 0;i<m;i++)
		for(j = 1;j<=n;j++){
			if(i == (1<<(j-1))){
				dp[i][j] = g[0][j];
				continue;
			}
			if((1<<(j-1)) & i){
				dp[i][j] = INF;
				for(k = 1;k<=n;k++)
					if( ((1<<(k-1))&i) && k!=j)
						dp[i][j]=min(dp[i][j],dp[i^(1<<(j-1))][k]+g[k][j]);
			}
		}
	for(i = 1;i<=n;i++)
		res = min(res,dp[m-1][i]+g[i][0]);
	printf("%d\n",res);
}
int main(){
	freopen("a.txt","r",stdin);
	while(scanf("%d",&n) && n){
		int i,j;
		for(i = 0;i<n+1;i++)
			for(j = 0;j<n+1;j++)
				scanf("%d",&g[i][j]);
		floyd();
		solve();
	}
	return 0;
}