加载中...

地址发布 老王说明书 宣传中心
此板块只作为纯讨论

正经话题,不搞色情!贤者时间必备
查看: 4772|回复: 29
收起左侧

[学习讨论] RPG游戏的寻路算法——从绯月讲起

[复制链接]
 楼主| 发表于 2023-12-25 18:31:47 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有帐号?免费注册

x
本帖最后由 navebayes 于 2023-12-26 12:02 编辑 2 _7 p( i( G! t- V  I. \  u(欢迎访问老王论坛:laowang.vip)
: ?/ n& {9 F' k: x, o! V. W(欢迎访问老王论坛:laowang.vip)
最近正在玩个游戏,也算是国产之光绯月仙行录,这个游戏哪里都好就是bug太多,并且作者过于摆烂,以至于有很多玩家都认为这个游戏就是故意拖着吃赞助的(bushi)
3 j# p' J  ?* V2 p3 V' i言归正传,在游玩这个游戏的过程中,我在一个评论区里看到这样一段话:自从玩了绯月之后,对于其他RPG游戏都看不上眼了,因为这个游戏独创了自动寻路的功能,可以说是RPG类游戏的里程碑式壮举。) e4 H8 j; K& X# n% v(欢迎访问老王论坛:laowang.vip)
我对此感到好奇,因为从前从来没有游玩RPG类游戏的经验,但我学过一点点算法,于是我打算用一些浅显易懂的方式说说自动寻路这样一个功能的实现。5 z  ^& J8 _! ?7 T(欢迎访问老王论坛:laowang.vip)
+ ]6 h; x4 |( q8 j" D# r6 I8 ?% u(欢迎访问老王论坛:laowang.vip)
主流的寻路算法:深度优先,广度优先,Dijkstra,A* 等,我这里主要讲讲后两种。
2 D9 f! X$ B: @4 ?
9 Y& C+ ~" ?) l' I" T" SDijkstra算法:这个算法是目前很多地图软件都在使用的算法,采用OPEN,CLOSE表的方式实现寻路功能。
' l2 F2 Q( K' R  I# P    创建两个表,OPEN, CLOSE。5 H" V3 y/ E9 v: `(欢迎访问老王论坛:laowang.vip)
OPEN表保存所有已生成而未考察的节点,CLOSED表中记录已访问过的节点。
5 h" \2 o3 T- M; e# P  j  x1. 访问路网中里起始点最近且没有被检查过的点,把这个点放入OPEN组中等待检查。
+ {; n# A: h9 S1 G& ~1 l2. 从OPEN表中找出距起始点最近的点,找出这个点的所有子节点,把这个点放到CLOSE表中。
7 S( f" ]/ q5 u2 k& b+ p3. 遍历考察这个点的子节点。求出这些子节点距起始点的距离值,放子节点到OPEN表中。
' J$ e; S/ b# Q) h4. 重复2,3,步。直到OPEN表为空,或找到目标点。
- S# p3 q" V6 T- ]1 o9 r* h* o8 S1 i) y2 ^' y! q(欢迎访问老王论坛:laowang.vip)
实际写代码还是比较占行数的,直接给出链接如下。
1 U5 W2 ~* v6 X$ C参考:https://blog.csdn.net/YiYeZhiNian/article/details/1222174504 Q6 D5 ]" u4 ~( V' e(欢迎访问老王论坛:laowang.vip)

) a* L( ]/ _$ o; {: m3 G; {6 ^3 g* C用这个算法,我写过一个课程作业,具体就是对于各个城市的地铁最短路径规划,大致还是比较成功的。先说说个人对于Dijkstra算法设计地铁线路规划:) C8 V: f6 q- P+ w. t# a* S(欢迎访问老王论坛:laowang.vip)
1.首先用爬虫爬到该城市的地铁网络,包含站点名称,该站点经纬度和线路图,并生成excel表格。用该excel表格生成pickle文件,方便直接调用。
" U' K- F8 M- Y- p2.注册高德地图开发者账号(该功能需要实名),得到一个key用于调用api(每日上限100次,超过付费,我调试到后面不让我调试了...)
" F1 j/ c6 M' V$ `4 i3.输入始发地和目的地,并通过api返回两个位置的经纬度坐标。
0 C/ [  M0 ?5 ?6 s- X4.比较两地理位置之间的最近地铁站,并根据Dijkstra算法实现路径规划。
' o3 Q  Q' i! h- D5.Dijkstra算法的本质就是不断选择新顶点并更新已处理表,并将他和邻居节点进行比对,当所有节点都被处理后即为最短路径
$ e0 G& f' ~( }" E至此,已初步完成具体工作流程。! J5 n/ N/ X# Y7 [- a5 x(欢迎访问老王论坛:laowang.vip)
  1. def get_nearest_subway(data,longitude1,latitude1):: l$ \- Y' K9 Z9 l: |6 L5 b(欢迎访问老王论坛:laowang.vip)
  2.   #找最近的地铁站6 t3 E0 a3 E8 x5 I$ M0 \/ ](欢迎访问老王论坛:laowang.vip)
  3.   longitude1=float(longitude1); t# j2 Y# U  {; m& J2 y(欢迎访问老王论坛:laowang.vip)
  4.   latitude1=float(latitude1)
    $ P2 ?& ?9 ?! b
  5.   distance=float('inf')
    ' J/ }: |) c7 M0 t1 f$ X% w: v
  6.   nearest_subway=None- ~" K; R6 t) H- L  m( u( F4 a(欢迎访问老王论坛:laowang.vip)
  7.   for i in range(data.shape[0]):
    - p" M+ z7 o3 N6 E* {7 x' D
  8.     site1=data.iloc[i]['name']   
    2 e7 C) q* |3 q6 y2 q& s/ P8 M
  9.     longitude=float(data.iloc[i]['longitude']); L2 e' v1 [; q) V(欢迎访问老王论坛:laowang.vip)
  10.     latitude=float(data.iloc[i]['latitude'])    #分别将经纬度代入,计算与目标之间的欧氏距离8 D8 X7 Z) ^, F+ f7 G" ~6 ^(欢迎访问老王论坛:laowang.vip)
  11.     temp=geodesic((latitude1,longitude1), (latitude,longitude)).m   #temp对其遍历即可,这里对比各个地铁站的欧氏距离
    ( g  ]1 s$ U3 B
  12.     if temp<distance:* V: B) L& ?$ u0 l% t3 T* k(欢迎访问老王论坛:laowang.vip)
  13.      distance=temp9 X# p) H# u: q! x$ d8 q0 ?+ [(欢迎访问老王论坛:laowang.vip)
  14.      nearest_subway=site15 f# r1 }3 P( ~" M& Q(欢迎访问老王论坛:laowang.vip)
  15.      return nearest_subway  
复制代码
  1. def subway_line(start,end):    #创建点之间的距离
    $ ]7 G/ e- z5 u$ V3 X+ ~
  2.   file=open('graph.pkl','rb')
    * @! ?8 U0 g. q/ k/ T; k1 S
  3.   graph=pickle.load(file)   #现在我们有了各个地铁站之间的距离存储在graph
    , [& Z$ Q1 H. f
  4.   costs={}    #创建节点的开销表,cost是指从start到该节点的距离
    & C8 U8 b! L/ t% S) D0 `3 q
  5.   parents={}
    * C- J8 v) k% P! ]) F5 w
  6.   parents[end]=None9 t" l: i& I. @1 p& p8 s+ V(欢迎访问老王论坛:laowang.vip)
  7.   for node in graph[start].keys():( r7 c* n; e  k" Z(欢迎访问老王论坛:laowang.vip)
  8.     costs[node]=float(graph[start][node])
    1 g( X  t6 S! D; F9 I! x, Z( Y6 M
  9.     parents[node]=start
    & n3 l9 P  H: u/ |% {
  10.     costs[end]=float('inf')    #终点到起始点距离为无穷大6 M) C: ^( F% e0 J, |3 F" v+ j/ C; y. u(欢迎访问老王论坛:laowang.vip)
  11.     processed=[]      #记录处理过的节点list
    # y3 R# Q$ p/ Y! f. J
  12.     shortest_path=dijkstra(start,end,graph,costs,processed,parents)) Y' t, u; o5 }- W$ Z(欢迎访问老王论坛:laowang.vip)
  13.     return shortest_path
复制代码
  1. #计算图中从start到end的最短路径
    9 u6 s6 `+ P( H2 r. n  M0 ?
  2. def dijkstra(start,end,graph,costs,processed,parents):
      i7 ?/ A# |/ d, a3 d
  3.     #查询到目前开销最小的节点
    9 V  ~  Z. K& [( ^% W4 ]
  4.     node=find_lowest_cost_node(costs,processed)
      _, T% t9 ]. b6 e$ ]
  5.     #使用找到的开销最小节点,计算它的邻居是否可以通过它进行更新: U8 X" M4 ?' L(欢迎访问老王论坛:laowang.vip)
  6.     #如果所有的节点都在processed里面 就结束* E2 M# R8 O5 c- A7 @3 b8 h- \: k(欢迎访问老王论坛:laowang.vip)
  7.     while node is not None:
    4 _2 Y5 _* S/ F2 G( q1 v( L
  8.     #获取节点的cost
    - w2 E$ h" ?- U: M0 Q; K
  9.       cost=costs[node]  #cost 是从node 到start的距离
    2 J4 M3 b4 Y3 f8 p+ q! e
  10.       #获取节点的邻居
    ) K5 Q/ t4 G3 O2 y' [. `
  11.       neighbors=graph[node]
    9 O/ T5 o$ A+ W3 F4 i
  12.       #遍历所有的邻居,看是否可以通过它进行更新0 d" K; d. I3 A/ O% _9 [$ f  J& ^(欢迎访问老王论坛:laowang.vip)
  13.       for neighbor in neighbors.keys():% A" ]# P; p, E8 s(欢迎访问老王论坛:laowang.vip)
  14.       #计算邻居到当前节点+当前节点的开销5 ^; |5 {+ P& j4 B3 ]1 D(欢迎访问老王论坛:laowang.vip)
  15.         new_cost=cost+float(neighbors[neighbor])
    ; u0 S% x: [1 F
  16.         if neighbor not in costs or new_cost<costs[neighbor]:( w/ w; H: F3 n: l- V(欢迎访问老王论坛:laowang.vip)
  17.         costs[neighbor]=new_cost" p6 {" A( ?5 O/ e3 d(欢迎访问老王论坛:laowang.vip)
  18.         #经过node到邻居的节点,cost最少6 o4 @8 P; K6 I8 j+ M(欢迎访问老王论坛:laowang.vip)
  19.         parents[neighbor]=node
    8 R6 }/ F/ e; {4 f0 B
  20.         #将当前节点标记为已处理! N. {4 c8 L" ^" o% v(欢迎访问老王论坛:laowang.vip)
  21.   processed.append(node)
    3 V1 p6 n+ V) H) `1 ]9 p, k3 n4 q% @) ^
  22.   #下一步继续找U中最短距离的节点  costs=U,processed=S9 C; w$ g4 S, Q' A0 \6 I(欢迎访问老王论坛:laowang.vip)
  23.   node=find_lowest_cost_node(costs,processed)
复制代码
  1. def find_lowest_cost_node(costs,processed):7 l' d3 _. G* l/ G" S! F& U(欢迎访问老王论坛:laowang.vip)
  2.     #初始化数据0 q" x! l- I, V; q(欢迎访问老王论坛:laowang.vip)
  3.   lowest_cost=float('inf') #初始化最小值为无穷大  f2 \) }6 K4 @(欢迎访问老王论坛:laowang.vip)
  4.   lowest_cost_node=None# K; ^2 @# H1 e* S+ z+ Q( i  s2 n6 o(欢迎访问老王论坛:laowang.vip)
  5.     #遍历所有节点
    3 P; v1 c( b1 j/ J, d8 W1 Q8 |
  6.   for node in costs:
    9 O+ O: h+ j; {& H6 F
  7.     #如果该节点没有被处理
    % y& n6 ^0 q4 S' E
  8.     if not node in processed:: `& U. A& L$ ^1 H6 Q% c(欢迎访问老王论坛:laowang.vip)
  9.       #如果当前的节点的开销比已经存在的开销小,那么更新该节点为最小开销的节点+ f$ o: j* X9 z(欢迎访问老王论坛:laowang.vip)
  10.       if costs[node]<lowest_cost:
    5 e" N6 Z. Q) V4 w
  11.       lowest_cost=costs[node]' t+ r( N. _3 A2 k(欢迎访问老王论坛:laowang.vip)
  12.       lowest_cost_node=node
    4 z# }4 n/ \. O) a
  13.     return lowest_cost_node
复制代码
上面这段基本上搬运的,主要是他代码已经写的很好了,注释也写的不错,但我写的时候爬虫调不出来(反爬虫技术可以的),最后是我手动去地图里找经纬度得到的结果。. E* d% Z0 E7 s( r8 y(欢迎访问老王论坛:laowang.vip)
引用:https://blog.csdn.net/fengdu78/article/details/111570695
& O/ Q3 ^( K  Y- K9 I  ~. S& I8 i

6 B: v; D* G# ?% W) O
2 d. y, ]+ r6 Q/ `& U: {3 h. R/ Q
. r* q2 X2 y% I( k7 ^6 P* k( ](欢迎访问老王论坛:laowang.vip)
A*算法:这个算法也是非常著名的算法,与Dijkstra算法相比,增加了启发式函数 ---- 启发函数的好坏直接决定了算法的效率和结果。由此衍生的D*(动态A算法)算法也被广泛运用于各类游戏中。D*的搜索效率高的原因就在于,当计划路径上某点被阻碍,因为是反向指针,可以定位到被堵节点的上一节点。也就是说只需重新搜索很小的一部分,其余部分仍然使用初始规划出的路径,大大提高了重规划的效率。0 |; l# d% B% [1 E(欢迎访问老王论坛:laowang.vip)

8 E  E9 C5 Z/ v0 w2 z
3 h0 j8 _" n, t额外补充-dfs&bfs逻辑
# t# I( b, o$ g( e  X6 M9 A深度搜索(dfs)和广度搜索(bfs)的算法逻辑可以说是最具有代表性的,基本任何有关于寻路的算法都没法绕开他俩。我担心可能没了解这2个算法直接去看Djkta和A*会有些懵逼
' |8 D+ _7 d* e
% j( z, o* z0 [4 o2 i5 q深度搜索(dfs) ! R' p+ `3 V( s7 c(欢迎访问老王论坛:laowang.vip)

1 i! P8 e8 T- M2 K8 Odfs就和它字面意思一样,是往更深的地方找(deepfound)
7 W+ N2 b/ P- C9 x7 }2 H它的核心思想很简单:
# V+ n2 j$ A" D5 \一直往前走,走不通就回头" X0 ]: V. T7 D. H1 @(欢迎访问老王论坛:laowang.vip)
dfs.png
: M/ p0 n8 o* q/ q) l" |. f顺序?当然是长幼啦,有长立长无长立幼 (1会先找2而非6,在2时会先找4而非5 ,直到4发现3 发现无路可走再back回去)9 H1 [/ C! G  Q(欢迎访问老王论坛:laowang.vip)
大致伪代码如下! c5 G8 R% b/ y; L/ x- }(欢迎访问老王论坛:laowang.vip)
  1. input 地图/ i4 L2 e; b- r8 i: E; o(欢迎访问老王论坛:laowang.vip)
  2. create 已经过点% u$ A3 h% a) w5 y) J- |! `; ?(欢迎访问老王论坛:laowang.vip)
  3. create 结果存储
    2 m5 \9 w4 @( Q3 s) _

  4. $ a# r' @% j( k6 m
  5. type node{
    7 c* Y5 E" t% l, O7 O' n
  6. node nextNode[];//下一节点们' L; D3 n2 U- Y# [; b5 _: H" W- X7 _(欢迎访问老王论坛:laowang.vip)
  7. nodeval;//节点标记物
    / M$ u; N  ]  a8 K: _
  8. }* {, z9 N) S0 C0 U, B1 |(欢迎访问老王论坛:laowang.vip)
  9. / A2 Q' B, |0 R& Z% t- U( z(欢迎访问老王论坛:laowang.vip)
  10. 3 d5 L7 h0 T/ z& f(欢迎访问老王论坛:laowang.vip)
  11. //这里开始是函数
    ' Q7 T( L% c$ L
  12. fun dfs(node* nowPoint)0 w2 l# U4 e. Z" T4 E(欢迎访问老王论坛:laowang.vip)
  13. {+ l; C1 ^" ]. j. h* V(欢迎访问老王论坛:laowang.vip)
  14. define u 为 nextNode[] size& b" s) j" h- }- b(欢迎访问老王论坛:laowang.vip)
  15. int  key;
    * `  s* s9 R! m) K8 r; Q
  16. % p( J, Y( B$ h- `8 C0 a(欢迎访问老王论坛:laowang.vip)
  17. for i in u {( B: w4 ]/ ^6 p5 r3 @(欢迎访问老王论坛:laowang.vip)
  18.       if (nowPoint.noteval == target)9 y/ E* ?- ^! F3 r0 ~  E6 o+ H7 m(欢迎访问老王论坛:laowang.vip)
  19.         {4 \6 `2 Q7 V7 F(欢迎访问老王论坛:laowang.vip)
  20.             结果存储[0] = nowPoint.noteval;
    , ?3 g" L2 F8 Y4 Z. J1 b6 S
  21.             return 0 ;/ r# A2 W! v  I& ]- F; A; ]8 x' G(欢迎访问老王论坛:laowang.vip)
  22.         }  
    ' v) w9 f$ D. G) l& W4 E/ _" i
  23.       else if( key = (dfs(nowPoint->nextNode[i])) != -1 ) //如果dfs这次没有返回负1(即 找到终点了)
    ; p2 v$ }& Z+ G: x
  24.             {
    6 Y1 F/ y5 L0 N+ B' F4 ]! F
  25.                 结果存储[key] = nowPoint.noteval;
    - |7 y1 g5 `0 D5 [( }8 w! d- ?+ h
  26.                 return key+1;         
    ( R  E7 \, h2 m  t1 h
  27.             }   , w5 X( k! D; y) ~& Y(欢迎访问老王论坛:laowang.vip)
  28.         else
    3 e  M$ t; a9 I, t/ E* g
  29.             {
    " O9 J# r% y4 s" d
  30.                 /********************/
    2 g" n5 s8 p! t: c6 m9 y
  31.                 /**        nothing          **/0 l+ r  u" T+ g3 j2 M(欢迎访问老王论坛:laowang.vip)
  32.                 /********************/0 Q" M# ~2 @1 S(欢迎访问老王论坛:laowang.vip)
  33.             }    : i: b# ]4 n$ V/ B3 B0 r(欢迎访问老王论坛:laowang.vip)

  34. - Q/ G8 o! T6 I/ J8 W
  35.         / V1 x& B# B% r" O) l2 L0 j(欢迎访问老王论坛:laowang.vip)
  36.    }   
    % r* F8 h# A- Q2 H5 d" X$ \( L
  37.     return -1;
    / Y) {2 x3 ?3 M6 U; S
  38. }* O+ [% ?" I6 n5 E9 z2 B(欢迎访问老王论坛:laowang.vip)
复制代码
就那么短,你只需要确定是不是就行了
4 [: A2 w. k7 Z- o是不是很简单:p 但是这就是一个比较原始的寻路算法的模样-顺其自然流
; Y) j2 y( {( ^' ]4 l' D6 k2 [
# X6 |2 |" S9 M$ [' ]- q* D在dfs算法中,你需要做的就是 询问+上抛
' m) A7 S' T8 R当然,dfs算法唯一能保证的就是‘找得到路’,这也是为什么纯粹的dfs算法不常用于生产环境中9 t; o. _1 M1 N0 t, s: n" |(欢迎访问老王论坛:laowang.vip)
6 K5 _5 Q  _5 ?, b% [(欢迎访问老王论坛:laowang.vip)

; x. {0 U+ k4 y, n7 Z广度搜索(bfs)
# [5 q- k) V4 B2 y知道深度,广度就很容易联想了 先找周围嘛 无论如何,先找周围0 C3 E! _5 A5 i( Z& Q0 }(欢迎访问老王论坛:laowang.vip)
bfs.png 9 M+ x+ I- e( F4 f$ S(欢迎访问老王论坛:laowang.vip)
这里不进行代码补充了,只简单地说一下逻辑
5 \$ c% I# ^0 e
! K3 Z# A7 Y" ]# k7 G6 u! z1 z1 l6 ?这个算法分以下几步:
0 R6 R7 }* j( |" C) C0,准备一个队列(就是那个queue),然后丢入我们的起点 命运的齿轮开始拨动
$ F  g) P( Y! w. A( I" p5 J
# q, c* p: t( S6 y3 D# f4 b1,询问当前节点是否终点?
* V, {5 E" B* b: B/ A' X2,将自己的nextNode 塞入队列中(除了访问过的)
/ m: z5 K) r8 z8 b7 C* P; X3,从队列里output一个节点,作为下一个‘当前节点’# S% ]- X* w# F; Q" ^(欢迎访问老王论坛:laowang.vip)
* {2 \8 e) C! _8 X3 Z(欢迎访问老王论坛:laowang.vip)
然后就是循环1~3
# z% J( z0 X7 L0 r! o
. p, H5 @' }2 q5 B4 b9 i# }是不是很简单?
6 P4 S- z" q& K9 R% }2 B& |
7 x: V: l. ~3 S+ i' D0 W这2个算法都属于随性流,一个适合终点远一个适合终点近。但无论如何这俩都暂时没有比较最优的功能
8 }* E  @4 Y- m- n. `6 W0 v因为他们刚~满~十~八~岁~的审敛条件就只是找得到路  ^5 j8 H9 U" |* O+ a(欢迎访问老王论坛:laowang.vip)

5 t9 \# u- G! H! g但你可以发现哦?如果将dfs的‘长幼判断’换成‘最优判断’,将呆值传递换为矩阵存储 就是dj算法了诶(a*也是类似的): d' S) L3 Z. T3 b(欢迎访问老王论坛:laowang.vip)
而bfs作为‘扩散式搜索’显然地在进行多点寻路时效用更加明显5 ?# P( y/ _. y+ J(欢迎访问老王论坛:laowang.vip)
如果觉得寻路算法很难的话,不妨先从dfs&bfs开始了解
0 R" s- r3 t0 C4 f
4 ~& F) n4 }2 G; K& b: D/ n
+ T% e+ r7 O4 \
3 e$ r, C) [6 F( Q: J7 N; V
. x0 F4 O* \0 V2 t( z; H

A*寻路算法

A*寻路算法

评分

参与人数 1软妹币 +235 收起 理由
navebayes + 235 cheese!!

查看全部评分

本帖被以下淘专辑推荐:

回复

使用道具 举报

发表于 2023-12-25 20:31:04 手机版 | 显示全部楼层
Applcu 发表于 2023-12-25 20:19
1 s/ a* Q  u/ L0 I我在想怎么算说明运作逻辑.....真贴上去么?这个编辑器到底怎么用啊/ f! ]+ H8 F: |  d(欢迎访问老王论坛:laowang.vip)
...
! H2 _- z3 S! I% B+ r(欢迎访问老王论坛:laowang.vip)
有一个代码框功能,可以用那个。主要编辑的话建议用正八经的md编辑器(我用的是国产的typora)
回复 支持 3 反对 0

使用道具 举报

发表于 2023-12-25 20:25:06 | 显示全部楼层
navebayes 发表于 2023-12-25 20:11
& M0 j+ V( U# E) u申请失败,原因:对算法本身描述过少,请更具体地说明算法的运作逻辑和方式,如果可以最好给出部分低代码
8 Y* e. [; v, |8 ]0 d3 P3 R3 _3 |$ b  g3 Q. b, Q, Y' d( h! F) b3 E  d(欢迎访问老王论坛:laowang.vip)
...

5 s9 L2 Q; N: [8 v- R4 T啊这。不如解释一下单源最短路径问题中对于负权值和负权值回路为什么不能使用上述算法4 ]! v8 h2 V9 F1 V. p6 i(欢迎访问老王论坛:laowang.vip)
+ k6 V: Q& N) t4 A6 ]. r(欢迎访问老王论坛:laowang.vip)
这玩意考试考得头昏脑涨,我tm 深度生成树画错了,越想越气 怎么会这么蠢
9 v* I$ M- X4 v  o3 t; C) P
回复 支持 2 反对 0

使用道具 举报

 楼主| 发表于 2023-12-25 23:28:47 | 显示全部楼层
昔有佳人公孙氏 发表于 2023-12-25 23:08! ]3 I6 @$ C1 t8 T% ]3 ?3 a. V(欢迎访问老王论坛:laowang.vip)
其实目前,本科生能接触到的寻路算法 无外乎迪杰斯特拉和弗洛伊德 一个不吃负值 一个不吃负环
3 }4 }' s/ F/ e. W# Y9 X但其实寻找最 ...
3 O* V; E* w; A* a+ W3 a! f9 Q(欢迎访问老王论坛:laowang.vip)
黏菌那个早有耳闻,确实是大自然的力量。我之后还想写点加密算法,比如AES,RSA,ECC之类的,其实这个地铁的算法也是我的一个课设,我大概花了两天就搞完了,我的主专业课是自控/单片机/数电/模电/嵌入式之类的,更偏自动化这个方向,算法我研究的不深。这个当毕业论文可能还是有一些勉强的,个人感觉工科是要做出一点实物的东西,当然不同学校可能要求不尽相同。
/ U8 p6 W. i; K
; K1 P( p) U* P" Q& z) x2 z) ?
# Y7 |' \' M- Q  I; D/ H1 L9 s
回复 支持 1 反对 0

使用道具 举报

 楼主| 发表于 2023-12-25 21:58:28 | 显示全部楼层
先说说个人对于Dijkstra算法设计地铁线路规划:" p8 f/ O' w  v+ R(欢迎访问老王论坛:laowang.vip)
1.首先用爬虫爬到该城市的地铁网络,包含站点名称,该站点经纬度和线路图,并生成excel表格。用该excel表格生成pickle文件,方便直接调用。
  _  B/ ?, v7 p0 W. _9 M. x% P2.注册高德地图开发者账号(该功能需要实名),得到一个key用于调用api(每日上限100次,超过付费,我调试到后面不让我调试了...)
, I; k+ z2 y) e& T' e8 \3.输入始发地和目的地,并通过api返回两个位置的经纬度坐标。
: x/ D) q# r9 Y1 `; T$ L: Q! T' r4.比较两地理位置之间的最近地铁站,并根据Dijkstra算法实现路径规划。1 ~: P; `3 i0 l6 ^(欢迎访问老王论坛:laowang.vip)
5.Dijkstra算法的本质就是不断选择新顶点并更新已处理表,并将他和邻居节点进行比对,当所有节点都被处理后即为最短路径0 C' X, f) d7 e; p* X! q(欢迎访问老王论坛:laowang.vip)
至此,已初步完成具体工作流程。  u$ k& y  s9 g6 i5 f(欢迎访问老王论坛:laowang.vip)
  1. def get_nearest_subway(data,longitude1,latitude1):) C- j, l  @. ](欢迎访问老王论坛:laowang.vip)
  2.     #找最近的地铁站
    : N+ @0 j" m' e; j
  3.     longitude1=float(longitude1)( @" _' o' {+ W0 ^: q2 w/ b(欢迎访问老王论坛:laowang.vip)
  4.     latitude1=float(latitude1)7 k4 _) F/ _4 ](欢迎访问老王论坛:laowang.vip)
  5.     distance=float('inf')
    9 Z, W& H3 ~: I& B0 a$ b+ y
  6.     nearest_subway=None4 ~3 n1 z5 i2 J1 F(欢迎访问老王论坛:laowang.vip)
  7.     for i in range(data.shape[0]):
    ; u0 o. K3 @8 w6 h6 h# l. ^
  8.   site1=data.iloc[i]['name']    ' K" h/ E. k$ o9 ?7 U3 \% E! h- M(欢迎访问老王论坛:laowang.vip)
  9.   longitude=float(data.iloc[i]['longitude'])
    ( k/ X6 {& E$ Q7 ~2 s$ l1 [& r
  10.   latitude=float(data.iloc[i]['latitude'])    #分别将经纬度代入,计算与目标之间的欧氏距离! K' Y3 O6 G. R(欢迎访问老王论坛:laowang.vip)
  11.   temp=geodesic((latitude1,longitude1), (latitude,longitude)).m   #temp对其遍历即可,这里对比各个地铁站的欧氏距离8 U+ b. d; W. b+ B' x& U5 z* A  \(欢迎访问老王论坛:laowang.vip)
  12.   if temp<distance:
    ( y4 O% N. Z4 ^( m8 _
  13.     distance=temp) U% S0 x( ^' Q  j' @(欢迎访问老王论坛:laowang.vip)
  14.     nearest_subway=site18 @% d' F+ B8 K* [(欢迎访问老王论坛:laowang.vip)
  15.     return nearest_subway  
复制代码
  1. def subway_line(start,end):    #创建点之间的距离1 [1 s. |- B) ^, ?" j(欢迎访问老王论坛:laowang.vip)
  2.     file=open('graph.pkl','rb')) t* K! ^. }8 p/ x& L; P(欢迎访问老王论坛:laowang.vip)
  3.     graph=pickle.load(file)   #现在我们有了各个地铁站之间的距离存储在graph9 n3 B- S" l% H- }(欢迎访问老王论坛:laowang.vip)
  4.     costs={}    #创建节点的开销表,cost是指从start到该节点的距离
    6 K6 ~+ Z$ ^; H# S
  5.     parents={}& Y* \$ `' S; X. @  W(欢迎访问老王论坛:laowang.vip)
  6.     parents[end]=None# Y: V/ m# Y, `6 O(欢迎访问老王论坛:laowang.vip)
  7.     for node in graph[start].keys():
    + v9 @/ s4 X7 ^& O" `# ?
  8.   costs[node]=float(graph[start][node])
    3 {* l, T( G5 f) j) k
  9.   parents[node]=start
    , @2 r4 q. K4 j( s2 q$ m1 T
  10.     costs[end]=float('inf')    #终点到起始点距离为无穷大$ m; I8 u0 k& r0 E8 K" m: g(欢迎访问老王论坛:laowang.vip)
  11.     processed=[]      #记录处理过的节点list- ~( j/ H, b+ X/ Z$ I$ K(欢迎访问老王论坛:laowang.vip)
  12.     shortest_path=dijkstra(start,end,graph,costs,processed,parents), t% G8 X' h0 s1 U2 y, Y4 {(欢迎访问老王论坛:laowang.vip)
  13.     return shortest_path
复制代码
  1. #计算图中从start到end的最短路径) M6 f2 h4 M3 g8 g- @" U(欢迎访问老王论坛:laowang.vip)
  2. def dijkstra(start,end,graph,costs,processed,parents):
    2 k, K8 ~: Z  o/ j3 J
  3.     #查询到目前开销最小的节点
    1 L3 z: \% n$ q7 b3 Y( G1 a
  4.     node=find_lowest_cost_node(costs,processed)
    , v( }3 N7 L: u' l0 C, L) D* w! P7 M
  5.     #使用找到的开销最小节点,计算它的邻居是否可以通过它进行更新
    7 W, j9 ~. p% [. e" {
  6.     #如果所有的节点都在processed里面 就结束/ j1 H7 n+ H9 i; U0 a% [: X1 o4 |; M(欢迎访问老王论坛:laowang.vip)
  7.     while node is not None:. M& V" m4 u$ R. ?9 m# ^(欢迎访问老王论坛:laowang.vip)
  8.   #获取节点的cost
    4 ~9 b( A% s$ s* u# \( U# M
  9.   cost=costs[node]  #cost 是从node 到start的距离: T0 N' {" Y. W! w) H: V(欢迎访问老王论坛:laowang.vip)
  10.   #获取节点的邻居! \; ~) ]! a& r5 x" S# k(欢迎访问老王论坛:laowang.vip)
  11.   neighbors=graph[node]
    ! F4 o; d% o& ]9 F, O0 n
  12.   #遍历所有的邻居,看是否可以通过它进行更新' Z  |' S) u/ I2 p- a) o(欢迎访问老王论坛:laowang.vip)
  13.   for neighbor in neighbors.keys():& ~- v; B: t% `8 Z" \1 E(欢迎访问老王论坛:laowang.vip)
  14.     #计算邻居到当前节点+当前节点的开销
    - i% ~+ M6 f/ ^8 I" z7 N! ]. ^
  15.     new_cost=cost+float(neighbors[neighbor]), Y4 ^1 u( z9 `" p. D(欢迎访问老王论坛:laowang.vip)
  16.     if neighbor not in costs or new_cost<costs[neighbor]:
    ) C& ?% f" j9 |3 M: [0 O! q
  17.       costs[neighbor]=new_cost
    7 ^  c( |. t& E+ v0 K$ E
  18.       #经过node到邻居的节点,cost最少
    - o: c4 H8 h8 A* c" }/ T# N
  19.       parents[neighbor]=node
    ' R5 v* `, K9 R
  20.   #将当前节点标记为已处理
    ) h. X: p" |; s6 G; O: N7 _
  21.   processed.append(node)
    4 w% z: a; [0 s) \$ M1 o% l
  22.   #下一步继续找U中最短距离的节点  costs=U,processed=S
    # b) a7 w3 d' Y; T
  23.   node=find_lowest_cost_node(costs,processed)
复制代码
4 S7 k6 E/ S; F3 M$ a3 o(欢迎访问老王论坛:laowang.vip)

6 E# p9 L. J  d3 P. X
  1. def find_lowest_cost_node(costs,processed):, m9 K# ^# M7 H7 e3 ^5 T(欢迎访问老王论坛:laowang.vip)
  2.     #初始化数据; {, k) U; l9 E; U! ?- d/ {0 w(欢迎访问老王论坛:laowang.vip)
  3.     lowest_cost=float('inf') #初始化最小值为无穷大, o# B, w6 ~% Z& N" d(欢迎访问老王论坛:laowang.vip)
  4.     lowest_cost_node=None
    5 ^/ i' e. i( {/ K- \- s: k3 C
  5.     #遍历所有节点
    & ^9 u2 B5 c& \- ~. K
  6.     for node in costs:
    : \+ }0 V6 v. g
  7.   #如果该节点没有被处理
    9 n' L7 s+ w' D2 d
  8.   if not node in processed:
    ! h! w1 t: a# h
  9.     #如果当前的节点的开销比已经存在的开销小,那么更新该节点为最小开销的节点
    % e& r6 f$ A) p7 D+ _/ r/ L) B
  10.     if costs[node]<lowest_cost:2 e$ q0 W. y/ w(欢迎访问老王论坛:laowang.vip)
  11.       lowest_cost=costs[node]
    " x" k5 L9 R4 s' X1 L+ c: J8 ?* ^  L! n
  12.       lowest_cost_node=node# O6 R9 G# w4 c: I(欢迎访问老王论坛:laowang.vip)
  13.     return lowest_cost_node
复制代码
上面这段基本上搬运的,主要是他代码已经写的很好了,注释也写的不错,但我写的时候爬虫调不出来(反爬虫技术可以的),最后是我手动去地图里找经纬度得到的结果。
% a$ H( M% k# x$ |2 N" B引用:https://blog.csdn.net/fengdu78/article/details/111570695% v8 ^& h$ m, c(欢迎访问老王论坛:laowang.vip)
回复 支持 1 反对 0

使用道具 举报

发表于 2023-12-25 20:11:24 手机版 | 显示全部楼层

修改修改

本帖最后由 navebayes 于 2023-12-25 20:14 编辑
$ \) ~7 }8 \; f
4 @0 t& Q' P: S8 ~申请失败,原因:对算法本身描述过少,请更具体地说明算法的运作逻辑和方式,如果可以最好给出部分低代码" V. m) ?  r) n2 n" u. Q) Z3 i(欢迎访问老王论坛:laowang.vip)

- h. b: A- K# m" g- W/ Xps:只要做滴好,王爷有赏啊
回复 支持 1 反对 0

使用道具 举报

 楼主| 发表于 2023-12-25 20:19:00 | 显示全部楼层
navebayes 发表于 2023-12-25 20:113 k$ {, k7 E+ ?) L4 e(欢迎访问老王论坛:laowang.vip)
申请失败,原因:对算法本身描述过少,请更具体地说明算法的运作逻辑和方式,如果可以最好给出部分低代码/ {/ F/ V' |! @2 ]1 S3 T' D" Q(欢迎访问老王论坛:laowang.vip)
2 b+ }! q5 p  _6 ~& _+ u(欢迎访问老王论坛:laowang.vip)
...

+ t8 M+ Z- Y# f% r" F我在想怎么算说明运作逻辑.....真贴上去么?这个编辑器到底怎么用啊
9 e" X: {9 }& d8 a; `5 U
回复 支持 反对

使用道具 举报

发表于 2023-12-25 20:35:28 手机版 | 显示全部楼层
昔有佳人公孙氏 发表于 2023-12-25 20:25( F- d  N( T& F0 r! O; n% @; k(欢迎访问老王论坛:laowang.vip)
啊这。不如解释一下单源最短路径问题中对于负权值和负权值回路为什么不能使用上述算法+ X9 k+ G6 V9 _6 n: ](欢迎访问老王论坛:laowang.vip)

3 c& b. C* e8 E# l0 e* C0 l7 T这玩意考试考得头 ...

8 S- i' ]' d& A- i( A( e看作者本人嘛,我不想干涉思路... t; X5 p4 L- F# F# e(欢迎访问老王论坛:laowang.vip)
(程序:哥们这可不是什么负权圈,这是时空之门啊啊啊啊啊(然后获得了每过一年减一岁的黄金体验))
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-12-25 20:43:49 | 显示全部楼层
昔有佳人公孙氏 发表于 2023-12-25 20:25
+ t% g! i# \' j) S2 ~, H0 L' Q9 V% `6 X* m啊这。不如解释一下单源最短路径问题中对于负权值和负权值回路为什么不能使用上述算法
- e" b+ }2 V1 L; B
  r# n+ i+ W2 t; ]7 s6 v! Q9 @) W这玩意考试考得头 ...
2 Q& j# v+ J  k7 X$ M(欢迎访问老王论坛:laowang.vip)
这应该是考专业课考破防了2 c3 Y9 l( b( m6 h. L1 E(欢迎访问老王论坛:laowang.vip)
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-12-25 20:48:19 | 显示全部楼层
昔有佳人公孙氏 发表于 2023-12-25 20:250 w* J4 G( k, ?- R! J- u3 D* W$ S(欢迎访问老王论坛:laowang.vip)
啊这。不如解释一下单源最短路径问题中对于负权值和负权值回路为什么不能使用上述算法  ?/ e+ R' I, W(欢迎访问老王论坛:laowang.vip)

. f! f5 p1 W3 H, u' J+ D% N这玩意考试考得头 ...
0 l. U8 C; \. E: Q+ O- l9 }. Y- w0 C(欢迎访问老王论坛:laowang.vip)
负权不就卡bug了么,越来越低就成环了" Z9 w; @4 |8 Q' q  @(欢迎访问老王论坛:laowang.vip)
回复 支持 反对

使用道具 举报

发表于 2023-12-25 20:56:03 手机版 | 显示全部楼层
Applcu 发表于 2023-12-25 20:48' [! Y2 [# d( H1 I' _(欢迎访问老王论坛:laowang.vip)
负权不就卡bug了么,越来越低就成环了
) i5 e5 Z6 @+ h+ L/ E. z5 I(欢迎访问老王论坛:laowang.vip)
宝宝,把剩下的内容补上嘛
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-12-25 21:15:16 | 显示全部楼层
navebayes 发表于 2023-12-25 20:56
9 E( Z' H: a9 @# L6 v  X宝宝,把剩下的内容补上嘛

/ r- T5 [6 R' S3 d2 K; a! R; n真是男童啊.....
1 G& v) u8 k$ S* B2 s
回复 支持 反对

使用道具 举报

发表于 2023-12-25 21:37:16 手机版 | 显示全部楼层
Applcu 发表于 2023-12-25 21:15
# I$ f# o1 L) I真是男童啊.....
, B+ u% @% a. e! d) \(欢迎访问老王论坛:laowang.vip)
诽谤管理封七天
回复 支持 反对

使用道具 举报

发表于 2023-12-25 22:34:33 | 显示全部楼层
Applcu 发表于 2023-12-25 21:58
9 h* p9 i6 W  {/ i先说说个人对于Dijkstra算法设计地铁线路规划:' F4 T% S' j2 p) ]5 h+ Y(欢迎访问老王论坛:laowang.vip)
1.首先用爬虫爬到该城市的地铁网络,包含站点名称,该站点 ...
+ N/ {* t4 b. Q5 W* l(欢迎访问老王论坛:laowang.vip)
可以加入正文喵
( T7 C( Z: q# N
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-12-25 22:42:13 | 显示全部楼层
昔有佳人公孙氏 发表于 2023-12-25 20:25! r6 |- O  j6 `8 L(欢迎访问老王论坛:laowang.vip)
啊这。不如解释一下单源最短路径问题中对于负权值和负权值回路为什么不能使用上述算法
& \) h& n9 n0 p- v7 U, Q. d+ T. W' d6 }/ r2 X0 ?& d(欢迎访问老王论坛:laowang.vip)
这玩意考试考得头 ...
5 k! |: c: g4 L- h(欢迎访问老王论坛:laowang.vip)
我也头大了.....但本质还是因为Dijkstra算法认为,如果有一条最短路线可以直达目的地,那么从旁边走另外的路将会直接不考虑,但在负权边这里,反而会减小cost,所以Dijkstra算法失效了,因为无法发现这条最短的路径" L: Q6 m9 q& l2 b(欢迎访问老王论坛:laowang.vip)
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-12-25 22:51:56 | 显示全部楼层
navebayes 发表于 2023-12-25 22:34: X4 E. f0 U$ l6 v6 Y+ r(欢迎访问老王论坛:laowang.vip)
可以加入正文喵

* l" @- I9 x( [+ D  X* F加了加了
$ O, S/ ?2 e! l. C3 A9 _
回复 支持 反对

使用道具 举报

发表于 2023-12-25 23:02:44 | 显示全部楼层
Applcu 发表于 2023-12-25 22:51
4 r' A3 R. i$ m# h) f8 F2 Q6 c加了加了
$ `. S9 y  r* N$ Y& c4 T# ]% {(欢迎访问老王论坛:laowang.vip)
介意让我改一下吗https://pic.imgdb.cn/item/6589996fc458853aef070f38.gif
8 ^( I+ R2 ~" a8 c' w" L
回复 支持 反对

使用道具 举报

发表于 2023-12-25 23:08:15 | 显示全部楼层
其实目前,本科生能接触到的寻路算法 无外乎迪杰斯特拉和弗洛伊德 一个不吃负值 一个不吃负环, a+ a' s  |! F( m, t/ m' l& b' P(欢迎访问老王论坛:laowang.vip)
但其实寻找最短路径的算法里还有很知名的黏菌算法。+ T7 O( t' ?5 T(欢迎访问老王论坛:laowang.vip)
有兴趣的可以去B站搜搜黏菌走迷宫,看看在大自然顶级的算法编辑下的单源最短路径问题的解法
3 E2 P9 J' n0 N+ ^6 m1 c- C不是说将奖励放置位置对应于日本的各大城市,两天内的黏菌路线就基本和目前日本的城市规划大致相似了吗,当然不太严谨是真的
  @6 l( @) W" B8 k) U就好像楼主说的做了地铁最短路径规划, 这玩意我舍友和你类似,她做的事全国铁路路径规划,涉及到价格,时间,反正做到最后是你输入起始地址和目的地,就能得出花费最少,时间最少,转站最少的路线,C++课设 我记得很清楚,我记得老师说的评价是,如果再加上运筹学的内容,说不定可以当毕业论文了,不出意外她拿了满分。
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-12-25 23:30:40 | 显示全部楼层
navebayes 发表于 2023-12-25 23:02& B, V+ m. E* u% D1 ^4 ^0 U(欢迎访问老王论坛:laowang.vip)
介意让我改一下吗
3 V6 V# m5 o) n/ k* Y- v(欢迎访问老王论坛:laowang.vip)
可以可以) R9 F1 @8 E  Q+ r/ ^- P$ ^. j) \(欢迎访问老王论坛:laowang.vip)
回复 支持 反对

使用道具 举报

发表于 2023-12-26 12:03:53 | 显示全部楼层
Applcu 发表于 2023-12-25 23:30
% E6 g/ H3 m" n6 ?" x/ Y8 ]可以可以
% J' p3 }- @: y7 k, s(欢迎访问老王论坛:laowang.vip)
排版改了下喵,东西加好了喵# p' E0 U: V1 r( d7 z( ](欢迎访问老王论坛:laowang.vip)
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 免费注册
点击进行验证

本版积分规则

我们不生产资源,只做资源的搬运工。

tags标签-春满四合院-AvGood-Archiver-小黑屋- |网站地图