HDU 5067 简单tsp

原文链接: HDU 5067 简单tsp

上一篇: Python TSP旅行商问题 暴力求解

下一篇: Python 遗传算法 求解tsp问题

hdu5067

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

参考

https://www.cnblogs.com/youmuchen/p/6879579.html

dp[v][s] 表示由s出发，经过v中的点回到0的最小值

对于不在v中的点k

dp[v | (1<<k) ][k] = min(

dp[v | (1<<k) ][k] ,

dis[k][s]+dp[v][s],        由k先到s，然后经过v到达0

)

#include <iostream>
#include <cstdio>
#include <cmath>
#include <queue>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#define LL __int64
#define INF 0x3f3f3f3f
const int maxn = 50+10;
using namespace std;
struct node
{
    int x, y;
} p[maxn];

int a[maxn][maxn], c[maxn][maxn], d[(1<<11)+10][15];

int main()
{
    int i, j, k, n, m, cnt;
    while(~scanf("%d%d", &n, &m))
    {
        memset(d, INF, sizeof(d));
        cnt = 0;
        for(i = 0; i < n; i++)
            for(j = 0; j < m; j++)
            {
                scanf("%d", &a[i][j]);
                if(i == 0 && j == 0)
                {
                    p[cnt].x = i;
                    p[cnt++].y = j;
                }
                else if(a[i][j]>0)
                {
                    p[cnt].x = i;
                    p[cnt++].y = j;
                }
            }
        for(i = 0; i < cnt; i++)
        {
            for(j = i+1; j < cnt; j++)
            {
                c[i][j] = c[j][i] = abs(p[i].x-p[j].x)+abs(p[i].y-p[j].y);
            }
            c[i][i] = 0;
        }
        d[0][0] = 0;
        for(i = 0; i < (1<<cnt); i++)
        {
            for(j = 0; j < cnt; j++)
                if(d[i][j]!=INF)  //注意保证该状态存在
                    for(k = 0; k < cnt; k++)
                    {
                        if(i&(1<<k)) continue;  //注意保证该点还没有加入
                        //if(c[j][k]==INF) continue;  //两个点之间可达

                        d[i|(1<<k)][k] = min(d[i|(1<<k)][k], d[i][j]+c[j][k]);
                    }
        }
        int ans = INF;
        for(i = 0; i < cnt; i++)
            ans = min(d[(1<<cnt)-1][i]+c[0][i], ans);

        cout<<ans<<endl;


        // for(int i=0;i<cnt;i++)
          // for(int j=0;j<cnt;j++)
            // printf("%3d  %c",c[j][i],j==(i<<cnt)-1? '\n':' ');


    }
    return 0;
}