状压dp解决旅行商问题TSP_dp 商旅问题 stp-程序员宅基地

技术标签: 算法  c++  动规  动态规划  

毕业旅行问题

小明目前在做一份毕业旅行的规划。打算从北京出发,分别去若干个城市,然后再回到北京,每个城市之间均乘坐高铁,且每个城市只去一次。由于经费有限,希望能够通过合理的路线安排尽可能的省一些路上的花销。给定一组城市和每对城市之间的火车票的价钱,找到每个城市只访问一次并返回起点的最小车费花销。

时间限制:C/C++ 1秒,其他语言2秒
空间限制:C/C++ 32M,其他语言64M

输入描述:

城市个数n(1<n≤20,包括北京)

城市间的车票价钱 n行n列的矩阵 m[n][n]

输出描述:

最小车费花销 s

示例1

输入例子:

4
0 2 6 5
2 0 4 4
6 4 0 2
5 4 2 0

输出例子:

13

例子说明:

共 4 个城市,城市 1 和城市 1 的车费为0,城市 1 和城市 2 之间的车费为 2,城市 1 和城市 3 之间的车费为 6,城市 1 和城市 4 之间的车费为 5,依次类推。假设任意两个城市之间均有单程票可购买,且票价在1000元以内,无需考虑极端情况。

题解

我们最终的目的是为找一个包含所有点的最小环。因此可以任意选一个节点(城市)作为出发节点(城市),不会影响最终的效果。这里我们挑选第一个节点(城市:北京)作为出发点节。

首先,我们考虑暴力结题的思路。从一个点出发,任意其他没走过的城市都可以作为备选,从一开始有N-1个备选直到最后只有一个备选。那么时间复杂度就是 O ( n ! ) O(n!) O(n!),这个时间复杂度是非常高的。在规定的时间内无法完成最高20个城市的旅行商问题。

因此我们考虑使用状态压缩的动态规划来解决这个问题,每个城市是否去过用一个bit表示,那么就需要 2 20 = 1048576 2^{20} = 1048576 220=1048576大小的数组就可以表示所有去过哪些城市的状态。我们还需要另一个状态表示当前处于哪一个节点(城市)。

d i s [ i ] [ j ] dis[i][j] dis[i][j] 表示 i i i 节点(城市)和 j j j 节点(城市)的距离。
d p [ v ] [ S ] dp[v][S] dp[v][S] 表示从 1 1 1 节点出发到过 S S S 这些城市且最后停留在 v v v 这个城市所需要的最短路程。接下来我们结合具体数值解释一下 d p [ v ] [ S ] dp[v][S] dp[v][S] 索引的意思,比如 S = 6 S=6 S=6,j换成二进制就是0000 0000 0000 0000 0000 0110,表示 20 20 20个城市目前从 1 1 1 节点出发到过 2 2 2 节点(城市)和 3 3 3 节点(城市); v = 3 v=3 v=3 表示目前停留在 3 3 3 节点(城市)。

由此我们可知只要求得 d p [ 1 ] [ 2 20 − 1 ] dp[1][2^{20}-1] dp[1][2201] 的值就是最终的答案了。这个值我们可以有由最初的状态不断地推导出来,且最初的状态结果是显而易见的。
d p [ v ] [ S ] = { d i s [ 1 ] [ v ] S ⊕ 2 v − 1 = = 0 min ⁡ ( d p [ u ] [ S ⊕ 2 v − 1 ] + d i s [ u ] [ v ] ) S ∣ 2 u − 1 = = 1 \begin{aligned} dp[v][S]=\left\{ \begin{array}{clc} dis[1][v] & & S \oplus 2^{v-1} == 0\\ \min(dp[u][S \oplus 2^{v-1}]+dis[u][v]) & & S | 2^{u-1} == 1 \end{array} \right. \end{aligned} dp[v][S]={ dis[1][v]min(dp[u][S2v1]+dis[u][v])S2v1==0S2u1==1

C++代码


#include <bits/stdc++.h>
using namespace std;

int N;  // 节点总数
int dis[21][21];  // 个节点之间的距离
int dp[21][1<<20];  // 第一个状态是最后一站,第二个状态描述目前旅行的地方

int main() {
    
    cin>>N;
    for(int i=1;i<=N;i++)
    for(int j=1;j<=N;j++){
    
        scanf("%d", &dis[i][j]);
    }
    memset(dp, 0x3f, sizeof(dp));
    for(int i=1;i<=N;i++){
     // 递推公式的第一个条件
        dp[i][(1<<(i-1))] = dis[1][i];
    }
    for(int i=1;i<(1<<N);i++){
    
        for(int v=1;v<=N;v++){
      // 下一个要去的节点
            if(i&(1<<(v-1))) continue;  // 已经去过v节点了
            for(int u=1;u<=N;u++){
      // 当前所在节点
                if(!(i&(1<<(u-1)))) continue; // 还没去去过u
                int nxt = i|(1<<(v-1));  // 下一个状态
                if(dp[v][nxt] > dp[u][i]+dis[u][v]){
      // 发现更近的路程
                    dp[v][nxt] = dp[u][i]+dis[u][v];
                }
            }
        }
    }
    cout<<dp[1][(1<<N)-1];  // 最后到达了1节点,并且经历过了所有节点。

    return 0;
}
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/zfq17796515982/article/details/130526768

智能推荐

从零开始搭建Hadoop_创建一个hadoop项目-程序员宅基地

文章浏览阅读331次。第一部分:准备工作1 安装虚拟机2 安装centos73 安装JDK以上三步是准备工作,至此已经完成一台已安装JDK的主机第二部分:准备3台虚拟机以下所有工作最好都在root权限下操作1 克隆上面已经有一台虚拟机了,现在对master进行克隆,克隆出另外2台子机;1.1 进行克隆21.2 下一步1.3 下一步1.4 下一步1.5 根据子机需要,命名和安装路径1.6 ..._创建一个hadoop项目

心脏滴血漏洞HeartBleed CVE-2014-0160深入代码层面的分析_heartbleed代码分析-程序员宅基地

文章浏览阅读1.7k次。心脏滴血漏洞HeartBleed CVE-2014-0160 是由heartbeat功能引入的,本文从深入码层面的分析该漏洞产生的原因_heartbleed代码分析

java读取ofd文档内容_ofd电子文档内容分析工具(分析文档、签章和证书)-程序员宅基地

文章浏览阅读1.4k次。前言ofd是国家文档标准,其对标的文档格式是pdf。ofd文档是容器格式文件,ofd其实就是压缩包。将ofd文件后缀改为.zip,解压后可看到文件包含的内容。ofd文件分析工具下载:点我下载。ofd文件解压后,可以看到如下内容: 对于xml文件,可以用文本工具查看。但是对于印章文件(Seal.esl)、签名文件(SignedValue.dat)就无法查看其内容了。本人开发一款ofd内容查看器,..._signedvalue.dat

基于FPGA的数据采集系统(一)_基于fpga的信息采集-程序员宅基地

文章浏览阅读1.8w次,点赞29次,收藏313次。整体系统设计本设计主要是对ADC和DAC的使用,主要实现功能流程为:首先通过串口向FPGA发送控制信号,控制DAC芯片tlv5618进行DA装换,转换的数据存在ROM中,转换开始时读取ROM中数据进行读取转换。其次用按键控制adc128s052进行模数转换100次,模数转换数据存储到FIFO中,再从FIFO中读取数据通过串口输出显示在pc上。其整体系统框图如下:图1:FPGA数据采集系统框图从图中可以看出,该系统主要包括9个模块:串口接收模块、按键消抖模块、按键控制模块、ROM模块、D.._基于fpga的信息采集

微服务 spring cloud zuul com.netflix.zuul.exception.ZuulException GENERAL-程序员宅基地

文章浏览阅读2.5w次。1.背景错误信息:-- [http-nio-9904-exec-5] o.s.c.n.z.filters.post.SendErrorFilter : Error during filteringcom.netflix.zuul.exception.ZuulException: Forwarding error at org.springframework.cloud..._com.netflix.zuul.exception.zuulexception

邻接矩阵-建立图-程序员宅基地

文章浏览阅读358次。1.介绍图的相关概念  图是由顶点的有穷非空集和一个描述顶点之间关系-边(或者弧)的集合组成。通常,图中的数据元素被称为顶点,顶点间的关系用边表示,图通常用字母G表示,图的顶点通常用字母V表示,所以图可以定义为:  G=(V,E)其中,V(G)是图中顶点的有穷非空集合,E(G)是V(G)中顶点的边的有穷集合1.1 无向图:图中任意两个顶点构成的边是没有方向的1.2 有向图:图中..._给定一个邻接矩阵未必能够造出一个图

随便推点

MDT2012部署系列之11 WDS安装与配置-程序员宅基地

文章浏览阅读321次。(十二)、WDS服务器安装通过前面的测试我们会发现,每次安装的时候需要加域光盘映像,这是一个比较麻烦的事情,试想一个上万个的公司,你天天带着一个光盘与光驱去给别人装系统,这将是一个多么痛苦的事情啊,有什么方法可以解决这个问题了?答案是肯定的,下面我们就来简单说一下。WDS服务器,它是Windows自带的一个免费的基于系统本身角色的一个功能,它主要提供一种简单、安全的通过网络快速、远程将Window..._doc server2012上通过wds+mdt无人值守部署win11系统.doc

python--xlrd/xlwt/xlutils_xlutils模块可以读xlsx吗-程序员宅基地

文章浏览阅读219次。python–xlrd/xlwt/xlutilsxlrd只能读取,不能改,支持 xlsx和xls 格式xlwt只能改,不能读xlwt只能保存为.xls格式xlutils能将xlrd.Book转为xlwt.Workbook,从而得以在现有xls的基础上修改数据,并创建一个新的xls,实现修改xlrd打开文件import xlrdexcel=xlrd.open_workbook('E:/test.xlsx') 返回值为xlrd.book.Book对象,不能修改获取sheett_xlutils模块可以读xlsx吗

关于新版本selenium定位元素报错:‘WebDriver‘ object has no attribute ‘find_element_by_id‘等问题_unresolved attribute reference 'find_element_by_id-程序员宅基地

文章浏览阅读8.2w次,点赞267次,收藏656次。运行Selenium出现'WebDriver' object has no attribute 'find_element_by_id'或AttributeError: 'WebDriver' object has no attribute 'find_element_by_xpath'等定位元素代码错误,是因为selenium更新到了新的版本,以前的一些语法经过改动。..............._unresolved attribute reference 'find_element_by_id' for class 'webdriver

DOM对象转换成jQuery对象转换与子页面获取父页面DOM对象-程序员宅基地

文章浏览阅读198次。一:模态窗口//父页面JSwindow.showModalDialog(ifrmehref, window, 'dialogWidth:550px;dialogHeight:150px;help:no;resizable:no;status:no');//子页面获取父页面DOM对象//window.showModalDialog的DOM对象var v=parentWin..._jquery获取父window下的dom对象

什么是算法?-程序员宅基地

文章浏览阅读1.7w次,点赞15次,收藏129次。算法(algorithm)是解决一系列问题的清晰指令,也就是,能对一定规范的输入,在有限的时间内获得所要求的输出。 简单来说,算法就是解决一个问题的具体方法和步骤。算法是程序的灵 魂。二、算法的特征1.可行性 算法中执行的任何计算步骤都可以分解为基本可执行的操作步,即每个计算步都可以在有限时间里完成(也称之为有效性) 算法的每一步都要有确切的意义,不能有二义性。例如“增加x的值”,并没有说增加多少,计算机就无法执行明确的运算。 _算法

【网络安全】网络安全的标准和规范_网络安全标准规范-程序员宅基地

文章浏览阅读1.5k次,点赞18次,收藏26次。网络安全的标准和规范是网络安全领域的重要组成部分。它们为网络安全提供了技术依据,规定了网络安全的技术要求和操作方式,帮助我们构建安全的网络环境。下面,我们将详细介绍一些主要的网络安全标准和规范,以及它们在实际操作中的应用。_网络安全标准规范