hdu 5418 Victor and World（旅行商问题）

题目链接：

http://acm.hdu.edu.cn/showproblem.php?pid=5418

解题思路：

旅行商问题（Traveling Salesman Problem，TSP）又译为旅行推销员问题、货郎担问题，简称为TSP问题，是最基本的路线问题，该问题是在寻求单一旅行者由起点出发，通过所有给定的需求点之后，最后再回到原点的最小路径成本。最早的旅行商问题的数学规划是由Dantzig（1959）等人提出。

TSP问题是NP问题，没有已知的多项式时间的高效算法可以解决这一问题。不过在程序设计竞赛中还是有可能出现这种范围较小的题目的。

所有可能的路线共有（n-1）!种。这是一个非常大的值，即使n（假如为16）已经很小了，仍然无法试遍每一种情况。对于这个问题，我们使用DP来解决。首先让我们试着写出它的递推式。

假设现在已经访问过的顶点的集合（起点0当做还未访问过的顶点）为S，当前所在的顶点为v，用dp[S][v]表示从v出发访问剩余的所有顶点，最终回到顶点0的路径的权重总和的最小值。由于从v出发移动到任意的一个节点u不属于S，因此有如下递推式

dp[V][0] = 0

dp[S][v] = min{dp[S U {u}] + dis(v,u) | u 不属于S}

我们只要按照这个递推式进行计算就可以了。由于在这个递推式中，有一个下标是集合而不是普通的整数，因此需要稍加处理。首先我们试着使用记忆化搜索求解。虽然有一个下标不是整数，但是我们可以把它编码成一个整数，或者给它们定义一个全序关系并用二叉搜索树存储，从而可以使用记忆化搜索来求解。特别的，对于集合我们可以把每一个元素的选取与否对应到一个二进制位里，从而把状态压缩成一个整数，大大方便了计算和维护。对于不是整数的情况，很多时候很难确定一个合适的递推式，因此使用记忆化搜索可以避免这个问题。不过在这个问题中，对于任意两个整数i和j，如果它们对应的集合S(i)属于S(j),就有i<=j，因此可以像下面的写法一样，通过循环求出答案。

int dp[1<<N][N];//记忆化搜索使用的数组

void solve(){
    //用足够大的值初始化数组
    memset(dp,INF,sizeof(dp));
    dp[(1<<n)-1][0] = 0;

    //根据递推式依次计算
    for(int s = (1<<n)-2; s >= 0; s--){
        for(int v = 0; v < n; v++){
            for(int u = 0; u < n; u++){
                if(!(s>>u & 1)){
                    dp[s][v] = min(dp[s][v],dp[s|1<<u][u]+dis[v][u]);
                }
            }
        }
    }
    printf("%d\n",dp[0][0]);
}

官方题解：

我们首先需要预处理出任意两个国家之间的最短距离，因为数据范围很小，所以直接用Floyd就行了。之后，我们用f[S][i]表示访问国家的情况为S，当前最后访问的一个国家是i所需要的最小总油量，其中，S的二进制表示记录了访问国家的情况，S在二进制表示下的第i位（不管是从左往右还是从右往左都可以）如果是1则表示第i个国家被访问过了，否则表示第i个国家没有被访问过，那么f[S|(1<<i)][i]=min(f[S][j]+f[i][j])，其中i和j满足S&(1<<j)=1且S&(1<<i)=0。最开始时，除了f[1][1]是0，其他情况都是无穷大，之后先枚举S，再枚举i（我验题的时候因为这里搞反结果WA了），那么最终的答案就是min(f[(1<<n)-1][i]+f[i][1])，其中i$\in$ [2,n]。总复杂度为O(n^3+n^2*2^n)O(n​3​​+n​2​​∗2​n​​)。

AC代码：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;

const int INF = 0x3f3f3f;
const int N = 18;
//输入
int n,m;
int dis[N][N];
int dp[1<<N][N];//记忆化搜索使用的数组

void floyd(){
    for(int k=0; k<n; k++)
        for(int i=0; i<n; i++)
            for(int j=0; j<n; j++)
                if(dis[i][k]+dis[k][j]<dis[i][j])
                    dis[i][j]=dis[i][k]+dis[k][j];
}

void solve(){
    //用足够大的值初始化数组
    memset(dp,INF,sizeof(dp));
    dp[(1<<n)-1][0] = 0;

    //根据递推式依次计算
    for(int s = (1<<n)-2; s >= 0; s--){
        for(int v = 0; v < n; v++){
            for(int u = 0; u < n; u++){
                if(!(s>>u & 1)){
                    dp[s][v] = min(dp[s][v],dp[s|1<<u][u]+dis[v][u]);
                }
            }
        }
    }
    printf("%d\n",dp[0][0]);
}

int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        int u,v,w;
        memset(dis,INF,sizeof(dis));
        for(int i = 0; i < n; i++)
            dis[i][i] = 0;
        for(int i = 1; i <= m; i++){
            scanf("%d%d%d",&u,&v,&w);
            u--;v--;
            dis[u][v] = dis[v][u] = min(dis[u][v],w);
        }
        floyd();
        solve();
    }
    return 0;
}