HDU3001 Travelling (三进制状压DP)

该博客讨论了一个经典的图论问题——旅行商问题。通过使用动态规划和三进制状态表示,找到访问N个城市且每个城市不超过两次的最小成本路径。文章详细介绍了算法实现,包括状态转移方程和预处理步骤,并提供了C++代码示例。
摘要由CSDN通过智能技术生成

题面如下
在这里插入图片描述
题意:
N个城市,每个城市访问不超过两次,旅行城市之间需要花费边权cost [u][v]求访问N个城市的最小花费。
分析:因为每个城市最多可以访问两次,我们用三进制来表示它的状态,进制位的数字表示为当前城市已经访问过的次数。
我们记dp[K][j] 表示在K状态下,以城市J作为结束城市所需要的最小花费预处理dp[3x][x] =0 即从这个点出发,花费为0,你所在的城市就是访问了1次,源点无花费。

枚举每个状态K,枚举城市J,枚举城市P,check由J城市能不能走到P城市,若可以,那么更新dp[K+3_P][P] =min(dp[K+3_P][P],dp[i][j]+cost[J][P])
如果枚举J的时候没有任何一个城市没被访问过,那么枚举当前状态的所有城市,更新答案。

ll G[14][14];
ll vis[14][14];
ll mi[14];
ll dp[177147][14];
ll dp_[177147][14];
signed main()
{
    ll n,m;
    mi[0]=1;
    for(int i=1; i<=12; i++)
    {
        mi[i]=mi[i-1]*3;
    }
    while(scanf("%lld%lld",&n,&m)!=EOF)
    {
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<n; j++)
            {
                G[i][j]=inf;
                vis[i][j]=0;
            }
        }
        for(int i=0; i<mi[n]; i++)
        {
            ll ok=i;
            for(int j=0; j<n; j++)
            {
                dp[i][j]=inf;
                dp_[i][j]=ok%3;
                ok=ok/3;
            }
        }
        for(int i=0; i<n; i++)
        {
            dp[mi[i]][i]=0;
        }
        for(int i=1; i<=m; i++)
        {
            ll u,v,w;
            read(u);
            read(v);
            read(w);
            u--;
            v--;
            vis[u][v]=vis[v][u]=1;
            G[u][v]=G[v][u]=min(G[u][v],w);
        }
        ll ans=inf;
        for(int i=0; i<mi[n]; i++)
        {
            ll lcj=1;
            for(int j=0; j<n; j++)
            {
                // 以j 结束
                if(dp_[i][j]==0)
                    lcj=0;
                for(int k=0; k<n; k++)
                {
                    if(j==k)
                        continue;
                    if(dp_[i][k]!=2&&vis[j][k]!=0)
                    {
                        dp[i+mi[k]][k]=min(dp[i][j]+G[j][k],dp[i+mi[k]][k]);
                    }
                }
            }
            if(lcj)
            {
                for(int p=0; p<n; p++)
                {
                    ans=min(ans,dp[i][p]);
                }
            }
        }
        if(ans>=inf)
        {
            printf("-1\n");
            continue;
        }
        printf("%lld\n",ans);
    }
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值