#include<cstdio>
#include<cstring>
using namespace std;
int n;
int dist[15][15];//输入的距离矩阵
int min_dist[15][15];//求出的最短路径距离矩阵
int d[15][1<<15];//用于表示DP状态
bool vis[15][1<<15];
int dp(int i,int s){//记忆话搜索,这个dp函数不会计算dp(0,0)的切记
if(vis[i][s])return d[i][s];//这个是结束条件,所以VIS[i][1<<i]与d[i][1<<i]都要预更新
vis[i][s]=true;
int &ans = d[i][s];//本次DP要更新的状态答案
ans = 1e9;
for(int j=0;j<=n;j++)//集合s中的一个节点号j
if(s&(1<<j)&&j!=i)//假设是从点集的节点J走到当前点I
if(ans>dp(j,s^(1<<i))+min_dist[j][i])//当前距离小于子集生成的距离
ans = dp(j,s^(1<<i))+min_dist[j][i];//递归求解,回溯更新答案
return ans;
}
void floyd(){//计算最短路径距离
for(int k=0;k<=n;k++)
for(int i=0;i<=n;i++)
for(int j=0;j<=n;j++)
if(min_dist[i][j]>min_dist[i][k]+min_dist[k][j])
min_dist[i][j] = min_dist[i][k]+min_dist[k][j];
}
int main(){
while( scanf("%d",&n)==1&&n ){
for(int i=0;i<=n;i++)
for(int j=0;j<=n;j++){
scanf("%d",&dist[i][j]);
min_dist[i][j]= dist[i][j];
}
floyd();//计算最短路径距离
memset(vis,0,sizeof(vis));
for(int i=0;i<=n;i++){//注意这里一定要单独初始化,否则dp函数计算出来的结果错误
vis[i][1<<i]=true;
d[i][1<<i]=min_dist[0][i];
}
printf("%d\n",dp(0,(1<<(n+1))-1 ));
}
return 0;
}
旅行商问题
POJ 3311 状态压缩
Sample Input
3
0 1 10 10
1 0 1 2
10 1 0 10
10 2 10 0
0
Sample Output
8
题意:求走过所有点并回到原点的最短路,可以走一个点多次.
思路:
因为可以走一个点多次,所以可以先求出每两个点之间的最短路,然后用经典的旅行商问题的状态压缩DP做法