旅行商问题(TSP)的 C++ 动态规划解法教学攻略

一、问题描述

旅行商问题(TSP)是一个经典的组合优化问题。给定一个无向图,图中的顶点表示城市,边表示两个城市之间的路径,边的权重表示路径的距离。一个售货员需要从驻地出发,经过所有城市后回到驻地,要求总的路程最短。

二、输入输出形式

输入形式

输入的第一行包含两个整数 n 和 m,分别表示顶点个数和边数。接下来的 m 行中,每行包含三个整数 u、v 和 w,表示顶点 u 和顶点 v 之间有一条边,边的权重为 w。

输出形式

输出一个整数,表示最短路程的总长度。

三、样例输入输出

样例输入

复制代码
5 10
1 2 3
1 3 1
1 4 5
1 5 8
2 3 6
2 4 7
2 5 9
3 4 4
3 5 2
4 5 3

样例输出

16

四、算法分析与设计

1. 动态规划算法

旅行商问题可以通过动态规划(DP)来解决。其核心思想是利用状态压缩来记录已经访问过的城市集合,并使用动态规划表来存储到达每个城市的最短路径。

2. 状态表示

定义 dp[S][i] 表示已经访问了集合 S 中的城市,最后到达城市 i 的最短路径。其中,S 是一个位掩码,表示已经访问过的城市集合。例如,S = (1 << n) - 1 表示所有城市都被访问过。

3. 状态转移

对于每个状态 S,遍历所有可能的城市 i,如果 i 在集合 S 中,则尝试将未访问的城市 j 加入集合 S,更新新的状态 S | (1 << j) 的最短路径。

4. 初始状态和结果提取

初始状态为 dp[full][i] = dist[i][0],其中 full 表示所有城市都被访问过,dist[i][0] 是城市 i 回到起点 0 的距离。最终结果从 dp[1][0] 中提取,表示从起点 0 出发,访问所有城市后回到起点的最短路径。

五、代码实现

复制代码
#include <iostream>
#include <vector>
#include <algorithm>
#include <climits>
using namespace std;

const int INF = 0x3f3f3f3f;

int main() {
    int n, m;
    cin >> n >> m;

    // 初始化邻接矩阵
    vector<vector<int>> dist(n, vector<int>(n, INF));
    for (int i = 0; i < n; i++) {
        dist[i][i] = 0;
    }

    // 输入边信息并构建邻接矩阵
    for (int i = 0; i < m; i++) {
        int u, v, w;
        cin >> u >> v >> w;
        u--; // 转换为从 0 开始的索引
        v--;
        dist[u][v] = w;
        dist[v][u] = w;
    }

    // 动态规划表
    int full = (1 << n) - 1;
    vector<vector<int>> dp(1 << n, vector<int>(n, INF));

    // 初始化动态规划表,full 表示所有城市都被访问过
    for (int i = 0; i < n; i++) {
        dp[full][i] = dist[i][0];
    }

    // 填充动态规划表
    for (int S = full; S >= 0; S--) {
        if ((S & 1) == 0) continue; // 必须包含起点(城市 0)

        for (int i = 0; i < n; i++) {
            if (!((S >> i) & 1)) continue; // 城市 i 不在集合 S 中时跳过

            if (S == full) {
                continue; // full 已经初始化
            }

            for (int j = 0; j < n; j++) {
                if ((S >> j) & 1) continue; // 城市 j 已经在集合 S 中时跳过

                dp[S][i] = min(dp[S][i], dp[S | (1 << j)][j] + dist[i][j]);
            }
        }
    }

    // 输出结果
    cout << dp[1][0] << endl;

    return 0;
}

六、代码解析

1. 邻接矩阵初始化

首先,我们创建一个大小为 n x n 的邻接矩阵 dist,初始化所有元素为 INF(表示无穷大)。然后将对角线元素设置为 0,因为每个城市到自身的距离为 0。接下来,读取边信息并更新邻接矩阵。

2. 动态规划表初始化

定义动态规划表 dp,其中 dp[S][i] 表示已经访问了集合 S 中的城市,最后到达城市 i 的最短路径。初始状态为 dp[full][i] = dist[i][0],其中 full 表示所有城市都被访问过,此时需要回到起点 0。

3. 状态转移填充

从 full 开始,逐步减少集合的大小,填充动态规划表。对于每个状态 S 和城市 i,如果 i 在集合 S 中,则尝试将未访问的城市 j 加入集合 S,更新新的状态 S | (1 << j) 的最短路径。

4. 输出结果

最终结果从 dp[1][0] 中提取,表示从起点 0 出发,访问所有城市后回到起点的最短路径。

七、复杂度分析

1. 时间复杂度

动态规划的状态数为 O(2^n * n),每个状态需要遍历所有可能的城市 n,因此总时间复杂度为 O(2^n * n^2)。

2. 空间复杂度

动态规划表的空间复杂度为 O(2^n * n),用于存储每个状态和城市的最短路径。

八、总结

旅行商问题是一个典型的组合优化问题,通过动态规划算法可以有效地解决。该算法利用状态压缩和动态规划表,能够在合理的时间内找到最短路径。对于较小规模的问题,这种方法是可行的,但对于大规模问题,可能需要更高效的算法或启发式方法。

希望本攻略能够帮助你理解旅行商问题的动态规划解法,并能够在实际问题中应用。如果有任何疑问或需要进一步的解释,请随时提问!

相关推荐
weixin_424183617 小时前
开放—封闭原则实战(下):基于策略的流程编排与编译期多态
c++·分布式·系统架构
Logic1017 小时前
C语言/数据结构位运算题解:异或XOR找出学习计划中的“独特学习时间“——只出现一次的数字
c语言·数据结构·数组·位运算·时间复杂度·算法题·异或性质
wuminyu8 小时前
C++协程实现接收端的零拷贝Buffer管理原理剖析
java·linux·c语言·jvm·c++
小蒋学算法8 小时前
算法-交替删除操作后最后剩下的整数-跟着灵茶山艾府学算法
数据结构·算法
m0_380743878 小时前
C++中的freopen的用法实例详解
开发语言·c++
荆棘鸟智能8 小时前
林业遥感长势评估怎么做?从NDVI时序分析到LiDAR蓄积量回归
人工智能·算法·智慧林业
jimy19 小时前
polymorph.cpp里面的“析构函数、智能指针”问题
开发语言·c++
hhb_6189 小时前
AIRAGDebug:一键定位RAG链路异常
人工智能·python·算法
朝朝辞暮i9 小时前
C++ 第 6 课:switch —— 多状态选择
开发语言·c++·算法
朝朝辞暮i10 小时前
C++ 第 4 课:if / else —— 让程序自己做决定
开发语言·c++