POJ3311 Hie with the Pie ACM解题报告(Floyd+状压dp)

题意是有N个城市(1~N)和一个PIZZA店(0),要求一条回路,从0出发,又回到0,而且距离最短  

也就是TSP(旅行商)问题,首先不难想到用FLOYD先求出任意2点的距离dis[i][j]  

接着枚举所有状态,用11位二进制表示10个城市和pizza店,1表示经过,0表示没有经过  

定义状态DP(S,i)表示在S状态下,到达城市I的最优值  

接着状态转移方程:DP(S,i) = min{DP(S^(1<<i-1),k) + dis[k][j],DP(S,i)},器重S^(1<<i-1)表示未到达城市i的所有状态,1<=k<=n  

对于全1的状态,即S = (1<<n)-1则表示经过所有城市的状态,最终还需要回到PIZZA店0  

那么最终答案就是min{DP(S,i) + dis[i][0]}  

!!!一定要注意这题的情况是dis[i][j]!=dis[j][i],(题目里有说 the time to go directly from location  i  to  j  may not be the same as the time to go directly from location  j  to  i)贡献了8个WA,看了好久题解都没发现,这题是有向图,注意注意!!!

#include<iostream>
#include<cstdio>
#include<cctype>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<string>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<sstream>
#include<stack>
using namespace std;
#define MAX 105
typedef long long LL;
const double pi=3.141592653589793;
const int INF=1e9;
const double inf=1e20;
const double eps=1e-6;
int d[15][2500];
int main(){
   int n;
   int dis[15][15];
   while(cin>>n&&n){
	   for(int i=0;i<=n;i++) for(int j=0;j<=n;j++) cin>>dis[i][j];
	   for(int k=0;k<=n;k++){
		   for(int i=0;i<=n;i++){
			   for(int j=0;j<=n;j++){
				   if(i!=k&&k!=j&&i!=j) dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);//Floyd求出每个点之间的最短路,因为这题可以重复经过某个点
				                                                                     //然后就可以看成不能重复经过的TSP问题了
			   }
		   }
	   }
           for(int s=0;s<(1<<n);s++){//s枚举已经经过的城市的情况
		   for(int i=1;i<=n;i++){
			   if(s&(1<<(i-1))){
				   if(s==(1<<(i-1))) d[i][s]=dis[0][i];//经过情况中只有i城市,所以是从0出发的。
				   else {
					   d[i][s]=INF;
					   for(int j=1;j<=n;j++){
						   if(j!=i&&s&(1<<(j-1))){//枚举已经经过j,并且j不是i
							   d[i][s]=min(d[i][s],dis[j][i]+d[j][s^(1<<(i-1))]);//从j-i的情况,当时经过城市里没有i
                                                                              //dis[j][i]!=dis[i][j],有向图,注意方向,不然这辈子看不出了
						   }
					   }
				   }
			   }
		   }
	   }
	   int ans=INF;
	   for(int i=1;i<=n;i++){
		   ans=min(ans,d[i][(1<<n)-1]+dis[i][0]);
	   }
	   printf("%d\n",ans);
	   
   }
    return 0;
}


你可能感兴趣的:(动态规划,ACM,poj)