快捷導(dǎo)航

python實現(xiàn)狄克斯特拉算法

更新時間：2021年04月07日 11:39:28 作者：果然還是我比較shuai

這篇文章主要介紹了python實現(xiàn)狄克斯特拉算法。想了解數(shù)據(jù)結(jié)構(gòu)和算法朋友可以參考下

數(shù)據(jù)結(jié)構(gòu)

1、路由信息

dictRoute = {}
dictRoute[nodeId] = {}
dictRoute[nodeId][nebrId] = distance
操作：
①根據(jù)nodeId找到該node的路由信息
②根據(jù)nebrId找到某一條路由的距離

2、節(jié)點信息

dictNode = {}
dictNode[nodeId] = [shortDis, fatherId, bIsCheck]
操作：
①找到nodes中最短距離的節(jié)點
②查找節(jié)點的shortDis，根據(jù)情況更新shortDis、fatherId
③檢查過的節(jié)點，更新bIsCheck

功能實現(xiàn)

/* 找到最短距離節(jié)點的Id，已經(jīng)檢查的不計算在內(nèi) */
def FindShortNodeId(dictNode):
return shortNodeId

/* dikstra算法流程 */
1、找到最短距離節(jié)點Id，并標記已檢查過（如果節(jié)點Id不存在，表示查找完成）
2、得到最短距離節(jié)點的距離
3、輪詢最短距離節(jié)點的鄰居節(jié)點
4、計算鄰居節(jié)點的新距離、得到原最短距離，進行比較
5、如果新距離 < 原距離，則更新鄰居節(jié)點最短距離
概括為兩步：步驟1 （1）- 找到當(dāng)前最短距離節(jié)點
步驟2（2~5） - 更新最短距離節(jié)點鄰居節(jié)點信息

代碼實現(xiàn)

import os
import sys

'''
信息輸入：
1、節(jié)點數(shù)目、路由數(shù)目
2、路由信息 
3、開始節(jié)點、結(jié)束節(jié)點
'''
nodeNum = 0 # 節(jié)點數(shù)目
routeNum = 0 # 路由數(shù)目
listRoute = [] # 臨時存儲輸入的路由信息
listNodeId = []# 臨時存儲節(jié)點id 

nodeIdStart = ''
nodeIdEnd = ''
dictRoute = {} # 解析后的路由信息
dictNode = {} # 節(jié)點信息
# 輸入節(jié)點數(shù)目、路由數(shù)目
strInput = input()
list0 = strInput.split(' ')
nodeNum = int(list0[0])
routeNum = int(list0[1])

# 輸入路由信息
for index in range(routeNum):
 strInput = input()
 listRoute.append(strInput)
 
# 輸入開始節(jié)點、結(jié)束節(jié)點
strInput = input()
list0 = strInput.split(' ')
nodeIdStart = list0[0]
nodeIdEnd = list0[1]

# 解析得到節(jié)點Id
listNodeId.append(nodeIdStart)
listNodeId.append(nodeIdEnd)
for index in listRoute:
 list0 = index.split(' ')
 nodeIdA = list0[0]
 nodeIdB = list0[1]
 if nodeIdA not in listNodeId:
  listNodeId.append(nodeIdA) 
 if nodeIdB not in listNodeId:
  listNodeId.append(nodeIdB) 

# 初始化路由信息字典、節(jié)點信息字典
for nodeId in listNodeId:
 # 節(jié)點字典信息
 dictNode[nodeId] = [10000, '', False] # 最短距離、父節(jié)點、是否檢查過
 # 每個路由字典創(chuàng)建
 dictRoute[nodeId] = {}
dictNode[nodeIdStart][0] = 0

# 初始化路由信息
for index in listRoute:
 list0 = index.split(' ')
 nodeIdA = list0[0]
 nodeIdB = list0[1]
 dictRoute[nodeIdA][nodeIdB] = int(list0[2])
 dictRoute[nodeIdB][nodeIdA] = int(list0[2])
 
# 打印輸入信息
def PrintInputInfo():
 print('nodeNum routeNum:')
 print(str(nodeNum) + ' ' + str(routeNum))
 print('nodeStart nodeEnd')
 print(nodeIdStart+' '+nodeIdEnd)
 print('route info:')
 for nodeId in dictRoute.keys():
  for nebrId in dictRoute[nodeId].keys():
   print(nodeId+'->'+nebrId+' = '+str(dictRoute[nodeId][nebrId]))
 print('node info:')
 for nodeId in dictNode.keys():
  print(nodeId+':'+str(dictNode[nodeId][0])+' '+dictNode[nodeId][1]+' '+str(dictNode[nodeId][2]))

#PrintInputInfo()

'''
狄克斯特拉實現(xiàn)
'''
# 找到最短距離節(jié)點id
def FindShortNodeId(dictNode):
 shortNodeId = ''
 shortDis = 10000
 for nodeId in dictNode.keys():
  if dictNode[nodeId][0] < shortDis and dictNode[nodeId][2] == False:
   shortNodeId = nodeId
   shortDis = dictNode[nodeId][0]
 return shortNodeId
 
# 狄克斯特拉算法
shortNodeId = FindShortNodeId(dictNode)
while shortNodeId:
 if shortNodeId == nodeIdEnd:
  break;
 dictNode[shortNodeId][2] = True
 shortDis = dictNode[shortNodeId][0]
 for nebrId in dictRoute[shortNodeId].keys():
  newDis = dictRoute[shortNodeId][nebrId] + shortDis
  if newDis < dictNode[nebrId][0]:
   dictNode[nebrId][0] = newDis
   dictNode[nebrId][1] = shortNodeId
 shortNodeId = FindShortNodeId(dictNode)
 
# 打印結(jié)果
listRst = []
nodeId = nodeIdEnd
while nodeId:
 listRst.append(nodeId)
 nodeId = dictNode[nodeId][1]
listRst.reverse()

strRst = ''
for nodeId in listRst:
 if nodeId == listRst[-1]:
  strRst += nodeId
 else:
  strRst += nodeId + '->'

if dictNode[nodeIdEnd][1] == '':
 print('cant reach '+nodeIdEnd)
else:
 print(strRst)
 print(dictNode[nodeIdEnd][0])