专题题解:基本算法

 

问题 C: 最短Hamilton路径

时间限制: 2 Sec  内存限制: 128 MB
[提交] [状态]

题目描述

给定一张 n(n≤20) 个点的带权无向图,点从 0~n-1 标号,求起点 0 到终点 n-1 的最短Hamilton路径。 Hamilton路径的定义是从 0 到 n-1 不重不漏地经过每个点恰好一次。

输入

第一行一个整数n。
接下来n行每行n个整数,其中第i行第j个整数表示点i到j的距离(一个不超过10^7的正整数,记为a[i,j])。
对于任意的x,y,z,数据保证 a[x,x]=0,a[x,y]=a[y,x] 并且 a[x,y]+a[y,z]>=a[x,z]。

输出

一个整数,表示最短Hamilton路径的长度。

样例输入 Copy

4
0 2 1 3
2 0 2 1
1 2 0 1
3 1 1 0

样例输出 Copy

4

提示

从0到3的Hamilton路径有两条,0-1-2-3和0-2-1-3。前者的长度为2+2+1=5,后者的长度为1+2+1=4

[提交][状态]

题解

用二进制上的数代表一个点的状态,取(1)或不取(0)。题目让求从点1到n的最短汉密顿路径,即经过每个点一次,这时的状态用二进制表示就是 (1<<n)-1 (n个1)。用dp[i][j]表示在状态 i 下,从1到 j 的最短汉密顿路径。

dp[i][j]可由上一个状态(上一状态就是把 j从当前状态中去掉)dp[i^(1<<(j-1))][k]得到,其中保证k是中存在的点,即 (i>>k)&1。

表示 i 的第 k 位是1,即经过点 k。  注意是 i>>k 不是 i<<k

则状态转移方程为:dp[i][j]=min{dp[i^(1<<j)][k]+Map[k][j]}(k=1~n); 其含义就是枚举到达点j之前的前一个点k,取其最短;
 

 

哈密顿图: 经过所有顶点且一次。

最开始想的用普通dp,但是只能保证两点之间的最短路程,并不能保证经过所有顶点。

第二,此题不能写ll,详情见第一个代码,解析见博客链接https://blog.csdn.net/yu121380/article/details/81041476。所以

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a[25][25],dis[(1<<20)+5][25];***
int main()
{
    int n;
    scanf("%d",&n);
    for(int i=0; i<n; i++)
    {
        for(int j=0; j<n; j++)
        {
            scanf("%lld",&a[i][j]);
        }
    }
    memset(dis,0x3f,sizeof(dis));
    dis[1][0]=0;
    for(int i=1; i<=(1<<n)-1; i++)
    {
        for(int j=0; j<n; j++)
        {
            if((i>>j)&1)//等价于i除以2的j次方,然后与1做与运算,也就是如果i/(2^j)得到的数装换为2进制后,如果最后一位为1,就输出1,最后一位为0就输出0。 
            {
                for(int k=0; k<n; k++)
                {
                    if((i>>k)&1)
                    {
                        dis[i][j]=min(dis[i][j],dis[i^(1<<j)][k]+a[k][j]);
                    }
                }
            }
        }
    }
 
    printf("%lld",dis[(1<<n)-1][n-1]);
 
}
 
/**************************************************************
    Problem: 7780
    User: 2019UPC110
    Language: C++
    Result: 运行错误
****************************************************************/
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int a[25][25],dis[(1<<20)+5][25];***避免溢出
int main()
{
    int n;
    scanf("%d",&n);
    for(int i=0; i<n; i++)
    {
        for(int j=0; j<n; j++)
        {
            scanf("%d",&a[i][j]);
        }
    }
    memset(dis,0x3f,sizeof(dis));
    dis[1][0]=0;
    for(int i=1; i<=(1<<n)-1; i++)
    {
        for(int j=0; j<n; j++)
        {
            if((i>>j)&1)
            {
                for(int k=0; k<n; k++)
                {
                    if((i>>k)&1)
                    {
                        dis[i][j]=min(dis[i][j],dis[i^(1<<j)][k]+a[k][j]);
                    }
                }
            }
        }
    }
 
    printf("%d",dis[(1<<n)-1][n-1]);
 
}
 
/**************************************************************
    Problem: 7780
    User: 2019UPC110
    Language: C++
    Result: 正确
    Time:1835 ms
    Memory:104424 kb
****************************************************************/

 

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值