商旅问题TSP——动态规划(c++ 动态规划)

1.定义

  TSP问题(旅行商问题)是指旅行家要旅行n个城市,要求各个城市经历且仅经历一次然后回到出发城市,并要求所走的路程最短。

本文以2019字节跳动春招笔试第5题,毕业旅行问题为例:

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

输入描述:
城市个数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

输出例子1:
13

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

2 动态规划思想

1 出发城市:0 ;经过{1,2,3}各城市有且仅有一次,再回到出发城市0
2 dp[i]{v}:表示从城市i出发,经过集合v,回到出发点0 的最段路径
3 最后结果dp[0]{1,2,3};

dp[0]{1,2,3}=min(dis[0][1]+dp[1]{2,3},dis[0][2]+dp[2]{1,3},dis[0][3]+dp[3]{1,2})

dp[i]{v}=min(dis[i][k]+dp[k]{v-k});(i不属于v,k属于v)
初始条件:
dp[i]{}=dis[i][0];//从i出发,不经过任何城市,回到城市0

难点:

1 怎么表示集合V?
除出发城市0 外剩下的城市为{1,2,3},则剩下城市所包含的所有子集为
{},{1},{2},{3},{1,2},{1,3},{2,3},{1,2,3}正好2^(n-1)个

设dp[i][j]中j为除出发点之外剩下的城市集合,则j的二进制中第i为正好对应第i个城市是否在j中。如:dp[0][7]=dp[0]{3,2,1} dp[0][5]=dp[0]{3,1}。
则dp表的列数为2^(n-1)列:(可写为1<<(n-1))

2 怎么查看集合j中某一位的状态?(判断城市k是否在集合中)
采用(j>>(k-1))&1==1判断,等于1则在,等于0则不在
如:判断城市1,2是否在集合j=5={3,1}中;
(j>>(1-1))&1=5&1=101&001=1:则城市1在集合j=5中
(j>>(2-1))&1=2&1=10&01=0:则城市2不在集合j=5中

3 怎么将集合j中某个城市置0?(即从集合j中取出某个城市k)
j^(1<<(k-1)) (和0异或不变,和1异或取反)
如从集合j=7={3,2,1}中取出城市2变为j={3,1}=101=5;
7^(1<<(2-1))=7 ^ 2=111 ^ 010=101=5;

代码如下:

#include<iostream>
#include<vector>
using namespace std;
const int MAX = 100000;
int main
  • 2
    点赞
  • 24
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
TSP问题旅行家旅行n个城市要求各个城市经历经历一次,然后回到出发城市,并要求所走路程最短。动态规划是解决TSP问题的一种方法。下面是一个使用C++实现TSP问题动态规划的例子: ```C++ #include <iostream> #include <cstring> using namespace std; const int MAXN = 20; const int INF = 0x3f3f3f3f; int c[MAXN][MAXN]; // 城市间距 int d[1 << MAXN][MAXN]; // TSP最终结果 int col; // TSP最终结果表的列数 int sub[1 << MAXN]; // 子集 void setDistance(int num) { // 创建城市间间距矩阵C for (int i = 0; i < num; i++) { for (int j = 0; j < num; j++) { cin >> c[i][j]; } } } int subContain(int n) { // 判断元素是否存在于子集 int i = 0; while (sub[i] != -1) { if (sub[i] == n) { return 1; } i++; } return 0; } int subGather(int n, int num) { // 将数字n转化为子集形式存放在集合sub[]中 int i = 0; int length = 0; memset(sub, -1, sizeof(sub)); while (n > 0) { if (n & 1) { sub[i] = num - 1; length++; } n >>= 1; i++; } return length; } int deSub(int k, int j) { // 求集合{1,3}剔除元素3后的子集 int t = 1 << j; return (k & (~t)); } int doTSP(int num) { memset(d, INF, sizeof(d)); d[1][0] = 0; for (int i = 1; i < (1 << num); i++) { int length = subGather(i, num); for (int j = 0; j < num; j++) { if (subContain(j)) { for (int k = 0; k < num; k++) { if (subContain(k) && k != j) { int t = deSub(i, j); d[i][j] = min(d[i][j], d[t][k] + c[k][j]); } } } } } return d[(1 << num) - 1][0]; } void showResult(int num) { // 输出结果 int i = (1 << num) - 1; int j = 0; col = 0; int result[20]; memset(result, -1, sizeof(result)); while (i > 0) { result[col++] = j; int t = i; i = deSub(i, j); j = 0; while (t > i) { t >>= 1; j++; } } cout << "最短路径为:" << endl; for (int i = col - 1; i >= 0; i--) { cout << result[i] << " "; } cout << endl; } int main() { int num; cout << "请输入城市数量:" << endl; cin >> num; cout << "请输入城市间距矩阵C:" << endl; setDistance(num); int shortestPath = doTSP(num); cout << "最短路径长度为:" << shortestPath << endl; showResult(num); return 0; } ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值