python用动态规划求删除路径_Python | 动态规划求解TSP
解题思路主要有两部分:
i为当前节点(城市),S为还没有遍历的节点(城市集合),表示从第i个节点起,经历S集合中所有的点,到达终点的最短路径长度。
回溯找到最优的路径,需要将S集合一一对应一个数字(类似于编码,一般用二进制),然后比如从节点i等于0开始,未经历集合为S={1,2,3},而下一步最优的节点 j 等于2,那么M[i][s]=j,回溯时只用从M[0][S]向后推即可。
import numpy as np
import itertools
import random
import matplotlib.pyplot as plt
# 设置中文识别
plt.rcParams['font.sans-serif'] = ['SimHei']
# tsp问题
class Solution:
def __init__(self,X,start_node):
self.X = X #距离矩阵
self.start_node = start_node #开始的节点
self.array = [[0]*(2**(len(self.X)-1)) for i in range(len(self.X))] #记录处于x节点,未经历M个节点时,矩阵储存x的下一步是M中哪个节点
def transfer(self, sets):
su = 0
for s in sets:
su = su + 2**(s-1) # 二进制转换
return su
# tsp总接口
def tsp(self):
s = self.start_node
num = len(self.X)
cities = list(range(num)) #形成节点的集合
# past_sets = [s] #已遍历节点集合
cities.pop(cities.index(s)) #构建未经历节点的集合
node = s #初始节点
return self.solve(node, cities) #求解函数
def solve(self, node, future_sets):
# 迭代终止条件,表示没有了未遍历节点,直接连接当前节点和起点即可
if len(future_sets) == 0:
return self.X[node][self.start_node]
d = 99999
# node如果经过future_sets中节点,最后回到原点的距离
distance = []
# 遍历未经历的节点
for i in range(len(future_sets)):
s_i = future_sets[i]
copy = future_sets[:]
copy.pop(i) # 删除第i个节点,认为已经完成对其的访问
distance.append(self.X[node][s_i] + self.solve(s_i,copy))
# 动态规划递推方程,利用递归
d = min(distance)
# node需要连接的下一个节点
next_one = future_sets[distance.index(d)]
# 未遍历节点集合
c = self.transfer(future_sets)
# 回溯矩阵,(当前节点,未遍历节点集合)——>下一个节点
self.array[node][c] = next_one
return d
# 计算两点间的欧式距离
def distance(vector1,vector2):
d=0;
for a,b in zip(vector1,vector2):
d+=(a-b)**2;
return d**0.5;
首先在(10,10)的坐标上,随机生成10个城市:
# 随机生成10个坐标点
n = 10
random_list = list(itertools.product(range(1, n), range(1, n)))
cities = random.sample(random_list, n)
x = []
y = []
for city in cities:
x.append(city[0])
y.append(city[1])
fig = plt.figure()
plt.scatter(x,y,label='城市位置',s=30)
plt.xlabel('x')
plt.ylabel('y')
plt.title('TSP问题 随机初始城市')
plt.legend()
plt.show()
使用欧氏距离计算两两城市之间的距离:
distence_matrix = np.zeros([n,n])
for i in range(0, n):
for j in range(n):
distence = distance(cities[i],cities[j])
distence_matrix[i][j] = distence
使用动态规划求解TSP问题:
S = Solution(distence_matrix,0)
print("最短距离:" + str(S.tsp()))
# 开始回溯
M = S.array
lists = list(range(len(S.X)))
start = S.start_node
city_order = []
while len(lists) > 0:
lists.pop(lists.index(start))
m = S.transfer(lists)
next_node = S.array[start][m]
print(start,"--->" ,next_node)
city_order.append(cities[start])
start = next_node
得到最终的访问路线:
x1 = []
y1 = []
for city in city_order:
x1.append(city[0])
y1.append(city[1])
x2 = []
y2 = []
x2.append(city_order[-1][0])
x2.append(city_order[0][0])
y2.append(city_order[-1][1])
y2.append(city_order[0][1])
plt.plot(x1,y1,label='路线',linewidth=2,marker='o',markersize=8)
plt.plot(x2,y2,label='路线',linewidth=2,color='r',marker='o',markersize=8)
plt.xlabel('x')
plt.ylabel('y')
plt.title('TSP问题 路线图')
plt.legend()
plt.show()
python用动态规划求删除路径_Python | 动态规划求解TSP相关推荐
- python迷宫问题求最短路径_用栈求解迷宫问题的所有路径及最短路径程序
目的:能将栈运用的更为熟练 实验内容:求解迷宫问题程序,要求输出如图所示的迷宫的路径,并求出第一条最短路径的长度以及最短路径. 设计的算法功能: mgpath(int xi,int yi,int xe ...
- python获取当前路径和上一层路径_python获取当前路径和上一级路径
围观人数: 3 标签:getc 当前目录 路径 print bsp 获取 目录 dir dirname """获取当前目录"&q ...
- 启动python项目找不到路径_Python os.startfile找不到FI
我对Python编程有问题.在我的程序中,我使用os.startfile()来启动程序,但它不允许我这样做. 当我的程序执行启动文件的命令时,我在Windows脚本主机中收到以下错误:Script: ...
- Python代码与Linux删除路径文件命令脚本
在使用linux操作进行深度学习的模型训练过程中,我们经常需要进行各种删除以及复制文件的操作 linux系统下进行文件夹的删除命令: -r 就是向下递归,不管有多少级目录,一并删除 -f 就是直接强行 ...
- python根据文件名获取文件路径_python 查看文件名和文件路径
以下是把sourceDir目录下的以.JPG结尾的文件所有拷贝到targetDir目录下: >>>import os >>> import os.path > ...
- python 获取向上两级路径_Python学习第171课--相对路径和绝对路径
[每天几分钟,从零入门python编程的世界!] 这节我们补充2个概念:相对路径和绝对路径. ●绝对路径 就像一棵大树一样,从它的根开始,往上会有大的枝干,在大的枝干上面又会有小一点的树枝,小树枝上面 ...
- python f.write 保存图片到路径_Python实现生成图片路径和对应标签的方式
学计算机的同学都知道图片的路径更改特别麻烦,尤其要对应到相应的标签,这个让人炸花了眼,接下来就会通过python语言直接生成图片对应标签,大家也知道现在的python特别的火爆,原因就是它能够处理很多 ...
- python sys.path.append()添加路径_Python调用CST进行天线建模仿真:环境搭建指南
CST微波工作室支持使用Matlab对其进行调用,但是一直以来官方没有提供python调用接口,但在最新的CST2020更新中添加了对Python3.6的调用支持,在安装路径下就已经附带了python ...
- 用python二重循环求成绩表_python的循环
python的循环 编写程序时经常有代码需要重复运行,python提供了while和for进行循环操作. 一.while循环 1.while循环可以根据条件进行判断,决定是否要循环执行语句块,语法如下 ...
最新文章
- 干货 | 为你解读34篇ACL论文
- android 图片加载 软引用_Android 解决图片大量下载:软引用必须懂4点
- Black Hat 2021上的七大网络威胁趋势
- 《机器学习》 周志华学习笔记第五章 神经网络(课后习题) python实现
- Android ImageButton示例代码
- writing avocado tests(写avocado测试用例)
- CV之Harris特征点检测器-兴趣点检测(详解)
- WinForm 中自定义文件与自己的应用程序相关联
- elasticsearch6 php,elasticsearch 6.x php-client
- Struts2工作流程
- Oracle Instance
- 反编译apk修改v7包_APK反编译
- JSK-A1144 代码填空:阶梯三角形【水题】
- win10连Android机,Windows10电脑可以和安卓手机无缝连接使用了,实用性干翻Mac-win10手机...
- 如何将 CHM 文件翻译成中文
- android 讯飞语音 提示组件未安装,迅雷看看提示“未安装组件”原因及其解决方法...
- 本地唯一985,要去省会了!
- 扫盲:回囙囜囝回囟因囡団囤囥囦囧囨囩囱囫…(认识吗?)
- linux启动jar后运行其他命令,Linux 部署jar 后台持续运行
- Linux一些基本操作(一)