2017年西安邮电大学第十二届数学建模竞赛B题论文

       笔者参加了2017年西安邮电大学第十二届数学建模竞赛,想到自己参赛前想了解一下学校往年的赛题,但能找到的资料很少,在这里我把我们参加B题论文放上来,希望能给大家一些帮助.
       2017年的校赛B题同时也是同一时间举行的2017年中兴捧月杯算法精英挑战赛赛题。

       当时我们三个大一学生,用整整五一的三天假期,完成了这篇论文。当时并没有掌握多少知识,可以说是初生牛犊不怕虎了,就大胆地做,乱搞。
       很幸运的是居然拿了校一等奖,当时我们大概是B题排名前20名的唯一一个大一组,我深知自己当时的水平,能拿一等奖应该说是思路比较正确,正好和老师对上了。总之,lucky,还是lucky。

赛题:

B题

        最强大脑中的收官蜂巢迷宫变态级挑战,相信大家都叹为观止!最强大脑收官战打响后,收视率节节攀升,就连蚁后也不时出题难为一下她的子民们。在动物世界中,称得上活地图的,除了蜜蜂,蚂蚁当仁不让。在复杂多变的蚁巢中, 蚂蚁总是能以最快、最高效的方式游历在各个储藏间(存储食物)。今天,她看完最新一期节目,又发布了一项新任务:

小蚁同学,我需要玉米库的玉米,再要配点水果,去帮我找来吧。小蚁正准备出发,蚁后又说:哎呀,回来,我还没说完呢,还有若干要求如下:

1.小蚁同学,你需要尽可能以最少的花费拿到食物(附件图中路线上的数值表示每两个储物间的花费);

2.小蚁同学,你最多只能经过9个储藏间拿到食物(包含起止两个节点,多次通过同一节点按重复次数计算);

3.小蚁同学,你必须经过玉米间,水果间(附件图中标绿色节点);

4.别忘了,食蚁兽也在路上活动呢,一旦与食蚁兽相遇,性命危矣!不过小蚁微信群公告已经公布了敌人信息(附件图中标红色路段)

5.最后,千万别忘了,还有两段路是必须经过的,那里有我准备的神秘礼物等着你呢(附件图中标绿色路段)

这下小蚁犯难了,这和它们平时找食物的集体活动规则不一样嘛,看来这次需要单独行动了。要怎么选路呢?小蚁经过一番苦思冥想,稿纸堆了一摞,啊,终于找到了!亲爱的同学们,你们能否也设计一种通用的路径搜索算法,来应对各种搜索限制条件,找到一条最优路径,顺利完成蚁后布置的任务呢?

注:

1、蚁巢,有若干个储藏间(附件图中圆圈表示),储藏间之间有诸多路可以到达(各储藏间拓扑图见附件);

2、节点本身通行无花费;

3、该图为无向图,可以正反两方向通行,两方向都会计费,并且花费相同;

4、起止节点分别为附件图中S点和E点。

5、最优路径:即满足限制条件的路径。

论文

                                    蚂蚁的最优路径分析

        针对空间给定两点,在限定条件下,本文对蚂蚁移动过程进行分析,根据图上两点间的对应关系建立模型——带权无向图模型,并讨论模型的精度和稳定性,最后又根据建立的模型改进得到了确定的路线图。

       基于带权无向图数学模型,运用了Dijkstra算法,建立了最佳路线模型,通过C++编程,利用每个条件对结果的影响程度大小对条件进行重要性排序,使结果优先满足重要性高的条件。并逐步增加条件来修改Dijkstra算法,将原问题分解成几个易于用Dijkstra算法求解的子问题,先对各个子问题逐一求局部最优解,再在此基础上求全局最优解,得出一个满足尽量多的条件下的结果

        最后针对计算结果中的误差,验证估计结果是否正确。

        求解得出总路线最短、总花费最少且相对均衡的最优组路线,路线如下:

                                             0→2→4→5→6→7→8→14→13→12→16→17

       具体路线如图所示:

                                

关键字:带权无向图模型  两点之间最短路径  条件重要性排序     Dijkstra算法  花费最少

 

                              蚂蚁的最优路径分析

                                             一、问题重述

        最强大脑中的收官蜂巢迷宫变态级挑战,相信大家都叹为观止!最强大脑收官战打响后,收视率节节攀升。在动物世界中,称得上活地图的,除了蜜蜂,蚂蚁当仁不让。在复杂多变的蚁巢中,蚂蚁总是能以最快、最高效的方式游历在各个储藏间(存储食物)。小蚁同学现需要通过蚁巢贮藏间取得玉米库的玉米,水果库的水果,还有要满足若干要求,如下:

1.最少的花费拿到食物(图中路线上的数值表示每两个储物间的花费);

2.最多只能经过9个储藏间拿到食物(包含起止两个节点,多次通过同一节点按重复次数计算);

3.必须经过玉米间,水果间(附件图中标绿色节点);

4.食蚁兽也在路上活动(附件图中标红色路段),一旦与食蚁兽相遇,会遭遇生命危险;

5.有两段路有蚁后准备的神秘礼物(附件图中标绿色路段),是必须经过的。

     请设计一种通用的路径搜索算法,来应对各种搜索限制条件,找到一条最优路径,顺利完成蚁后布置的任务。

注:

1、蚁巢,有若干个储藏间(图中圆圈表示),储藏间之间有诸多路可以到达;

2、节点本身通行无花费;

3、该图为无向图,可以正反两方向通行,两方向都会计费,并且花费相同;

4、起止节点分别为附件图中S点和E点。

5、最优路径:即满足限制条件的路径是能以最快、最高效的方式游历在各个储藏间(存储食物)。蚁巢有若干个储藏间(下图中圆圈表示),储藏间之间有诸多路可以到达(各储藏间及路径如图所 示)。

                                      

                                                                      图1.1  蚁巢中储藏间及路线图

 

                                   二、模型假设与符号说明

2.1 模型假设

  1. 假定不同普通储物间(非玉米间、水果间)除位置外其他条件都是相同的;
  2. 假定每两个储物间之间的花费为图中两点之间边的权值;

     3.假定第一个条件求最短路径设定为解决此问题的大前提;

     4.根据题意,假定不同的条件对最终结果的影响程度是不同的,是可以按照对结果的影响大小(条件的重要性)对条件进行分级排列的;

     5.经分析论证,本题的所有条件不可能全部实现(分析论证见下文的条件分析),假定将本题的最优路径含义理解为是最接近于结果的一条路径(即题中所说 蚂蚁总是能以最快、最高效的方式游历在各个储藏间)。

2.2 符号说明

  1. 将点N1、N2...N16设定为点1、2...16,将S点和E点设定为点0和点17,如下表所示:

符号

意义

符号

意义

符号

意义

0

点S

6

点N6

12

点N12

1

点N1

7

点N7

13

点N13

2

点N2

8

点N8

14

点N14

3

点N3

9

点N9

15

点N15

4

点N4

10

点N10

16

点N16

5

点N5

11

点N11

17

点E

 

    2.  I表示条件对问题结果的影响程度,即条件的重要性。按数值大小重要性递增。如条件2的重要性为I_{2}

    3.d(j)表示j点到起点的距离

    4.a_{j,k}表示j,k两点之间的距离

    5.INF代表无穷大的值

    6.ans代表结果

                                                三、问题分析

3.1条件分析

1.以最少的花费拿到食物(附件图中路线上的数值表示每两个储物间的花费);

*分析:在本题的分析、建模、求解过程中,最少的花费是这个问题的核心约束条件。我们将其作为解决这个问题的大前提,不与其他4个条件并列分析。

2.最多只能经过9个储藏间拿到食物(包含起止两个节点);

*分析:在本题中,若满足题目要求同时过玉米仓、水果仓、两条绿色路径、起点、终点,则已经需要经过8个固定的点,而小蚁同学最多只能经过9个储藏间拿到食物(包含起止两个节点,多次通过同一节点按重复次数计算),显然不能成立。所以本条件是已经确定不能成立的了,所以将其重要性设定为I_{2}=1

3.必须经过玉米间,水果间(附件图中标绿色节点);

*分析:蚁后分配的任务就是要求小蚁通过蚁巢贮藏间取得玉米间的玉米、水果间的水果。显然,与蚂蚁路径的长短相比,完成任务是更重要的。所以,我们将其重要性设定为I_{3}=3

4. 食蚁兽也在路上活动,一旦与食蚁兽相遇,会遭遇生命危险;(附件图中标红色路段);

*分析:红色路段为不可控因素,如果走红色路段,虽然不确定蚂蚁一定能和食蚁兽相遇,但一旦相遇,蚂蚁就会有生命危险,这样蚂蚁任务的其他的条件全都无法实现,在确定最优计划的实际问题中,我们应该尽可能的避免不可控因素以及其所造成的影响。所以,我们将其重要性设定为I_{4}=4

5.有两段路是必须经过的,那里有神秘礼物(附件图中标绿色路段)。

*分析:在实现小蚁以最少花费拿到玉米和水果的过程中,神秘礼物对该目的的实现造成的影响不可预测,应尽量降低它对解决问题影响程度,所以将其重要性设定为I_{5}=2 

3.2条件权重排序

通过对条件的逐条分析,确定权重排序:

    1.食蚁兽也在路上活动,一旦与食蚁兽相遇,会遭遇生命危险;(附件图中标红色路段);I_{4}

    2.必须经过玉米间,水果间(附件图中标绿色节点);I_{3}

    3.有两段路是必须经过的,那里有神秘礼物(附件图中标绿色路段)。I_{5}

    4.最多只能经过9个储藏间拿到食物(包含起止两个节点);I_{2}

                                                        I_{4}=4>I_{3}=3>I_{5}=2>I_{2}=1

                           


                                                                     图3.2.1 条件权重排序图

       将图中各个储藏间及之间的路径抽象为带权无向图模型,将问题转化为求带权无向图中经过指定中间节点集的最短路径,按照3.2中条件重要性顺序设计算法求解满足题意的路径,确定最优解。

                           

                                                            图3.2.2 建模思路流程图

四、模型建立与求解

4.1模型准备

       将蚂蚁巢简化为一个带权无向图模型,每两个储物间之间的花费简化为边的权值。于是,问题转化成了寻找S点和E点之间的最短路径问题经过指定中间节点集和边集的最短路径问题。

      首先,我们满足重要性最大的条件即条件4:不与食蚁兽相遇。将题目附带的图整理简化为下图,直接将红色的边删除。

          

         并将图转换为数据,以便程序读入存入邻接矩阵

                              

 

4.2模型的求解

(1) Dijkstra算法求带权无向图的单源最短路

       在重要性最高的第四个条件满足后,我们用深度优先搜索的方法来寻找所有可能的路径,共找到了148895条路径,因为数量太多,无法一一列举,搜索代码见附录。

      之后我们运用Dijkstra算法来寻求最短路径。

      Dijkstra算法的基本思路是:先将与起点有边直接相连的节点到起点的距离记为对应的边的权重值,将与起点无边直接相连的节点到起点的距离记为无穷大。然后以起点为中心向外层层扩展,计算所有节点到起点的最短距离。每次新扩展到一个距离最短的点后,更新与它有边直接相邻的节点到起点的最短距离。当所有点都扩展进来后,所有节点到起点的最短距离将不会再被改变,因而保证了算法的正确性

[1] Dijkstra算法求解最短路径问题的基本步骤如下:

      ①设立 Y 和 N 两个集合, Y 用于保存所有等待访问的节点,N记录所有已经访问过的节点。

      ②访问网络节点中距离起始节点最近且没有被访问过的节点,把这个节点放入 Y 中等待访问。

      ③从 Y 中找出距离起点最近的节点,放入 N 中,更新与这个节点有边直接相连的相邻节点到起始节点的最短距离,同时把这些相邻节点加入 Y 中。

                                                             d(j)=min(d(j),a_{j,k}+d(k))

      ④重复步骤(2)和(3),直到 Y 集合为空, N 集合为网络中所有节点为止。

 

[2]由于题中图节点较多,下面用该图来演示一下Dijkstra算法

                                      

                                 图4.2 Dijkstra算法演示样图

开始点(源点):start

D[i]:顶点i到start的最短距离。

初始:

     D[start]=0

     D[i]=a[start,i]  (无边设为INF)

 

       集合1:已求点                                              集合2:未求点

           

       

算法步骤:

1、在集合2中找一个到start距离最近的顶点k ,距离=d[k]

2、把顶点k加到集合1中,同时修改集合2 中的剩余顶点j的d[j]是否经过k后变短。如果变短修改d[j]

If  ( d[k]+a[k,j]<d[j])  -> d[j]=d[k]+a[k,j])

3、重复1,直至集合2空为止。

                   

                  

                  

 

[3]至此,我们用可以用典型的Dijkstra算法可以找出满足重要性最大的条件4的最短路径(实现代码见附录)

               

                                                                       路径①

              

                                                                         路径②

            路径①,②皆经过7个点,花费为6

(2)基于Dijkstra算法的带权无向图经过指定边及节点集的单源最短路径模型

       在重要性最高的两个条件4、3都满足的情况下(蚂蚁既不与食蚁兽相遇,又经过水果间和玉米间),我们用深度优先搜索的方法来寻找所有可能的路径,共找到了57242条路径,因为数量太多,无法一一列举,搜索代码见附录。

      之后我们改进典型的Dijkstra算法来寻求最短路径

      上文展示了Dijkstra算法的主要思想,但考虑本题中的限制条件后,典型Dijkstra算法就不再适用了,需要对其进行改进。我们再考虑满足重要性第二大的条件3,即寻找包含7、12点的最短路径。

1.经过指定节点集的改进

        贪心算法,是一种通过分级处理某些最优解问题的方法。贪心算法,在求解过程的每一步中都采取在当前状态下看来是最好或最优的选择,从而希望最终的结果是最好或最优的。如果一个复杂的问题能够分解成几个子问题来解决,并且子问题的最优解能递推到最终问题的最优解,简言之,如果局部最优解能最终推导出全局最优解,那么贪心算法在解决这类问题中将会非常有效。

       贪心算法求解问题的基本步骤:

       ①把原问题分解成若干个易于求解的简化的子问题。

       ②对每一子问题分别求解,得到所有子问题的局部最优解。

       ③通过对全部子问题的某个局部最优解进行分析、聚合等处理,生成原问题的一个解。

       ④从上一步所求出的原问题的解集或某个初始解出发,通过筛选、迭代等手段求出原问题的解(全局最优解)。

       

       我们算法的改进就是基于贪心算法思想的,将原问题分解成几个易于Dijkstra算法求解的子问题,先对各个子问题逐一求局部最优解,再在此基础上求全局最优解。具体说来,就是先将预处理后的网络拓扑图中的所有节点分为三个子集合,分别为包含起点的起点集,包含全部的需要经过的中间节点的中间节点集,包含终点的终点集。通过Dijkstra算法依次求起点集到中间节点集之间的局部最短路径,连通中间节点集中所有节点的局部最短路径,中间节点集到终点集之间的局部最短路径,以此求出一条从起点出发经过指定的所有中间节点后到达终点的待选全局最短路径。通过对有限的待选全局最短路径进行筛选,选出其中距离最短的路径,即为满足要求的全局最短路径。

对本题来说

                                

                                                                                       或

                               

       二者结果是一样的。

       这样,我们可以找出满足条件4、3的最短路径(实现代码见附录)

                                  

 

2.经过指定节点集和边集的改进

       在经过指定节点集的Dijkstra算法的基础上,将指定边的两节点加入到指定节点的集合中,并且该边的权值记录下来,图中的权值赋值为一个极大值,即不能重复走这条边。

       这样,指定节点集中就有了6个点,对这6个点的排列顺序进行全排列,对每一种组合进行判断,如果两条特殊边的两顶点不靠在一起,就舍弃这种排列。

      假设E中以i,j为起点终点的路径最短

                                         

      这样,我们就只剩最后一个条件没有满足了。因为这个条件在本题中其他条件限制下是不能满足的,所以将其最多经过的储藏间个数扩大为10个、11个、12个······通过深度优先搜索算法来寻找所有可能的路径。

          ①扩大为10个时仍然是无解的

          ②扩大为11个点时有1条路径

            0->2->4->5->12->6->7->8->14->13->17      14(指总花费)

         ③扩大为12个点时共有63条路径

路径

总花费

0->1->2->4->5->12->6->7->8->14->13->17

17

0->1->4->2->3->7->6->12->13->14->15->17

19

0->1->4->2->3->7->6->14->13->12->16->17

17

0->1->4->2->3->7->8->14->13->12->16->17

15

0->1->4->2->5->12->6->7->8->14->13->17

17

0->2->4->5->3->7->6->12->13->14->15->17

18

0->2->4->5->3->7->6->14->13->12->16->17

16

0->2->4->5->3->7->8->14->13->12->16->17

14

0->2->4->5->6->7->8->14->13->12->16->17

13

0->2->4->5->10->12->6->7->8->14->13->17

14

0->2->4->5->12->6->3->7->8->14->13->17

16

0->2->4->5->12->6->7->8->14->13->15->17

19

0->2->4->5->12->6->7->8->14->13->16->17

16

0->2->4->5->12->6->7->8->14->13->17

14

0->2->4->5->12->6->7->8->15->14->13->17

17

0->2->4->5->12->13->14->6->7->8->15->17

22

0->2->4->9->5->12->6->7->8->14->13->17

17

0->2->4->9->10->12->6->7->8->14->13->17

14

0->3->2->4->5->12->6->7->8->14->13->17

15

0->2->4->5->3->7->6->12->13->14->15->17

18

0->2->4->5->3->7->6->14->13->12->16->17

16

0->2->4->5->3->7->8->14->13->12->16->17

14

0->2->4->5->6->7->8->14->13->12->16->17

13

0->2->4->5->10->12->6->7->8->14->13->17   

14

0->2->4->5->12->6->3->7->8->14->13->17

16

0->2->4->5->12->6->7->8->14->13->15->17

19

0->2->4->5->12->6->7->8->14->13->16->17

16

0->2->4->5->12->6->7->8->14->13->17

14

0->2->4->5->12->6->7->8->15->14->13->17

17

0->2->4->5->12->13->14->6->7->8->15->17

22

0->2->4->9->5->12->6->7->8->14->13->17

17

0->2->4->9->10->12->6->7->8->14->13->17

14

0->3->2->4->5->12->6->7->8->14->13->17

15

0->1->2->4->5->12->6->7->8->14->13->17

17

0->1->4->2->3->7->6->12->13->14->15->17

19

0->1->4->2->3->7->6->14->13->12->16->17

17

0->1->4->2->3->7->8->14->13->12->16->17

15

0->1->4->2->5->12->6->7->8->14->13->17

17

0->2->4->5->3->7->6->12->13->14->15->17

18

0->2->4->5->3->7->6->14->13->12->16->17

16

0->2->4->5->3->7->8->14->13->12->16->17

14

0->2->4->5->6->7->8->14->13->12->16->17

13

0->2->4->5->10->12->6->7->8->14->13->17

14

0->2->4->5->12->6->3->7->8->14->13->17

16

0->2->4->5->12->6->7->8->14->13->15->17

19

0->2->4->5->12->6->7->8->14->13->16->17

16

0->2->4->5->12->6->7->8->14->13->17

14

0->2->4->5->12->6->7->8->15->14->13->17

17

0->2->4->5->12->13->14->6->7->8->15->17

22

0->2->4->9->5->12->6->7->8->14->13->17

17

0->2->4->9->10->12->6->7->8->14->13->17

14

0->1->2->4->5->12->6->7->8->14->13->17

17

0->1->4->2->5->12->6->7->8->14->13->17

17

0->2->4->5->6->7->8->14->13->12->16->17

13

0->2->4->5->10->12->6->7->8->14->13->17

14

0->2->4->5->12->6->7->8->14->13->15->17

19

0->2->4->5->12->6->7->8->14->13->16->17

16

0->2->4->5->12->6->7->8->14->13->17

14

0->2->4->5->12->6->7->8->15->14->13->17

17

0->2->4->5->12->13->14->6->7->8->15->17

22

0->2->4->9->5->12->6->7->8->14->13->17

17

0->2->4->9->10->12->6->7->8->14->13->17

14

0->3->2->4->5->12->6->7->8->14->13->17

15

           ④扩大到13个点时有331条路径,因为路径太多,无法一一列举

          ⑤······

 

        用改进后的Dijkstra算法求得最优路线是

                              

          共经过12个点,总花费是13

        结合上述数据验证,这条路径的确是花费最少的路径。

       虽然经过的点不是最少的,但实现了其他条件下花费最少,又因为经过储藏间的个数是对结果影响最小的一个条件,所以,此道路是最优解。

所以,本题的最优解为

                                      

五、模型评价

5.1模型优缺点

5.1.1模型优点

     (1)逐步求解方法可以保证模型中最后结果的精确性,使误差变得尽可能的小。

   (2)抽象出的模型具有可转移性,当研究问题中的限制条件发生改变时,可更改所建模型中相关代码实现方法的通用应用。

     (3)运用 C++语言编写程序,方便高效地完成求解带权无向模型。

5.1.2模型缺点

     (1)缺点:条件筛选所用的条件增加逐步求解方法,过程略微繁琐。

   (2)在运行代码搜索路径的过程中,若有多种花费相同的符合条件可选择路径,只能输出其中一条。只能通过后续分析筛选出最优路径。

5.2 模型改进

(1)通过lingo的敏感分析(如分析已知量的变化,分析最优函数的反应)来进行更深层次的数据分析。

(2)可以用链表和结构体对算法进行优化。

推广:本文的求解方法可以适用于交通路线的选择等实际问题。

参考文献

[1] 黄书力,胡大裟,蒋玉明. 经过指定的中间节点集的最短路径算法. [计算机工程与应用]  41-45  2015,51(11)

[2] Thomas H.CormenCharles E.Leiserson 《算法导论》 北京 机械工业出版社 2014

[3] 刘汝佳 《算法竞赛入门经典》 北京 清华大学出版社 2015

 

附录:

  • Dijstra算法(满足不经过食蚁兽红边)

/*

程序会输出最优解和最优路径

*/

#include <iostream>

#include<stdio.h>

#include<cstring>

#include<stack>

#define maxn 536870912

#define up(x,y) for(int i=x;i<=y;i++)

#define mem(x) memset(x,0,sizeof(x))

using namespace std;

 int mapp[20][20];

    bool b[20];

    int d[20],fa[20];

stack<int>q;

int main()

{

    freopen("input.txt","r",stdin);

    int t,n,x,y,e,v,minn;

    n=17;

    t=42;

    up(0,n)

      for(int j=0;j<=n;j++)

        if(i==j)mapp[i][j]=0;else mapp[i][j]=maxn;

    while(t--)

    {

        cin>>x>>y>>e;

        mapp[x][y]=mapp[y][x]=e;

 

    }

    mem(b);mem(fa);

    up(0,n)

      {

          d[i]=mapp[0][i];

      }

    mapp[11][12]=mapp[12][11]=maxn;//不经过食蚁兽红边

    mapp[2][4]=mapp[4][2]=0;

    mapp[13][14]=mapp[14][13];

    up(1,n-1)//从每次找最小的边开始,用这条边松弛别的边

        {

        minn=maxn;

        for(int j=0;j<=n;j++)

          if(b[j]==0&&d[j]<minn)

          {

              minn=d[j];

              v=j;

          }

        b[v]=1;

        for(int j=0;j<=n;j++)

            if(b[j]==0&&d[j]>mapp[j][v]+d[v])//读入图的边,存入邻接矩阵

              {

                  d[j]=mapp[j][v]+d[v];

                  fa[j]=v;

              }

    }

    cout<<d[n]<<"\n";

    while(fa[n]!=0)

    {

        q.push(fa[n]);

        n=fa[n];

    }

    q.push(0);

    while(q.empty()==0)//输出最短路径

    {

        cout<<q.top()<<"->";

        q.pop();

    }

    cout<<"17\n";

   return 0;

}

 

  • Dijstra算法经过指定节点集的改进

#include <iostream>

#include<stdio.h>

#include<cstring>

#define maxn 536870912

#define n 17

#define t 42

using namespace std;

int mapp[20][20];

bool b[20];

int d[20],fa[20];

int x,y,e,v,minn;

int dijkstra(int s,int e)

{

    memset(b,0,sizeof(b));

    memset(fa,0,sizeof(fa));

    for(int i=0;i<=n;i++)

      {

          d[i]=mapp[s][i];

      }

    for(int i=1;i<=n-1;i++)//从每次找最小的边开始,用这条边松弛别的边

        {

        minn=maxn;        for(int j=0;j<=n;j++)

          if(b[j]==0&&d[j]<minn)

          {

              minn=d[j];

              v=j;

          }

        b[v]=1;

        for(int j=0;j<=n;j++)

            if(b[j]==0&&d[j]>mapp[j][v]+d[v])

              {

                  d[j]=mapp[j][v]+d[v];

                  fa[j]=v;

              }

    }

   return d[e];

 

}

int main()

{

    freopen("input.txt","r",stdin);

    for(int i=0;i<=n;i++)

      for(int j=0;j<=n;j++)

        if(i==j)mapp[i][j]=0;else mapp[i][j]=maxn;//读入图的边,存入邻接矩阵

    for(int i=1;i<=t;i++)

    {

        cin>>x>>y>>e;

        mapp[x][y]=mapp[y][x]=e;

    }

    mapp[11][12]=mapp[12][11]=maxn;;//去掉食蚁兽红边

    cout<<dijkstra(0,7)+dijkstra(7,12)+dijkstra(12,17)<<endl;

   return 0;

}

 

 

 

 

  • Dijkstra算法经过指定节点集和边集的改进

/*

程序最终会输出最优解,通过最优解的值可以从上面输出的路径中找出

最优解的路径的中间点集的顺序

输入数据即上文4.1模型准备中将图转换为的数据

*/

#include <iostream>

#include<stdio.h>

#include<cstring>

#include<algorithm>

#include<stack>

#define maxn 536870912

#define n 17

#define t 42

using namespace std;

int mapp[20][20];

int a[10]={2,4,7,13,14,12};

bool b[20];

int d[20],fa[20];

int x,y,e,v,minn;

stack<int>q;

stack<int>o;

int dijkstra(int s,int e)

{

    memset(b,0,sizeof(b));

    memset(fa,0,sizeof(fa));

    for(int i=0;i<=n;i++)

      {

          d[i]=mapp[s][i];

      }

    for(int i=1;i<=n-1;i++)//从每次找最小的边开始,用这条边松弛别的边

        {

        minn=maxn;//先吧min给赋初值

        for(int j=0;j<=n;j++)

          if(b[j]==0&&d[j]<minn)

          {

              minn=d[j];

              v=j;

          }

        b[v]=1;

        for(int j=0;j<=n;j++)

            if(b[j]==0&&d[j]>mapp[j][v]+d[v])

              {

                  d[j]=mapp[j][v]+d[v];

                  fa[j]=v;

              }

    }

   return d[e];

 

}

int main()

{

    //freopen("input.txt","r",stdin);

    for(int i=0;i<=n;i++)

      for(int j=0;j<=n;j++)

        if(i==j)mapp[i][j]=0;else mapp[i][j]=maxn;//读入图的边,存入邻接矩阵

    for(int i=1;i<=t;i++)

    {

        cin>>x>>y>>e;

        mapp[x][y]=mapp[y][x]=e;

    }

    mapp[11][12]=mapp[12][11]=maxn;//去掉食蚁兽红边

    mapp[2][4]=mapp[4][2]=maxn;//绿色特殊便

    mapp[13][14]=mapp[14][13]=maxn;//绿色特殊边

    int sf,sum,cnt(100000);

    while(next_permutation(a,a+5))//中间节点集的6个点的全排列

    {

        sf=0;

        for(int i=0;i<=5;i++)

        {

            if(a[i]==4)//判断绿色特殊边

              if(a[i+1]==2||a[i-1]==2)

                sf++;

            if(a[i]==2)

              if(a[i+1]==4||a[i-1]==4)

                sf++;

            if(a[i]==14)

              if(a[i+1]==13||a[i-1]==13)

                sf++;

            if(a[i]==13)

              if(a[i+1]==14||a[i-1]==14)

                sf++;

        }

        if(sf==4)//如果两条特殊边都经过

        {

            cout<<0<<' ';

            for(int i=0;i<=5;i++)cout<<a[i]<<' ';cout<<17<<endl;

            sum=0;

            for(int i=0;i<=4;i++)

            {

                if((a[i]==4&&a[i+1]==2)||(a[i]==2&&a[i+1]==4))

                {

                    sum+=2;

                }

                else

                if((a[i]==13&&a[i+1]==14)||(a[i]==14&&a[i+1]==13))

                {

                    sum+=1;

                }

                else

                {

                sum+=dijkstra(a[i],a[i+1]);

                }

            }

        sum+=dijkstra(0,a[0]);

        sum+=dijkstra(a[5],17);

        cout<<sum<<"----------------"<<endl;

        if(sum<cnt)

            cnt=sum;

    }

    }

    cout<<"\n最优解:"<<cnt<<endl;

   return 0;

}

  • 满足不经过食蚁兽红线的路径搜索代码

#include <iostream>

#include<stdio.h>

#include<cstring>

#include<stack>

#define maxn 536870912

using namespace std;

int mapp[20][20];

bool b[20];

int d[20],fa[20];

int t(42),n(17),x,y,e,v,minn;

void print(int x)

{

    stack<int>q;

    int flag(0);

    while(x!=0)

    {

        q.push(x);

        x=fa[x];

    }

    q.push(0);

    while(q.empty()==0)

    {

        cout<<q.top()<<"->";

        q.pop();

    }

    cout<<"17\n";

}

void dfs(int x)

{

    if(x==17)

    {

        print(fa[x]);

        return;

    }

    for(int i=0;i<=n;i++)

        if(b[i]==0&&mapp[x][i]<maxn)

        {

            b[i]=1;

            fa[i]=x;

            dfs(i);

            b[i]=0;

        }

}

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);//若采用文件输入输出可以去掉注释

    for(int i=0;i<=n;i++)

      for(int j=0;j<=n;j++)

        if(i==j)mapp[i][j]=0;else mapp[i][j]=maxn;//读入图的边,存入邻接矩阵

    while(t--)

    {

        cin>>x>>y>>e;

        mapp[x][y]=mapp[y][x]=e;

    }

    mapp[11][12]=mapp[12][11]=maxn;

    memset(b,0,sizeof(b));

    memset(fa,0,sizeof(fa));

    dfs(0);

   return 0;

}

  • 满足不经过食蚁兽红线和经过玉米水果间的路径搜索代码

#include <iostream>

#include<stdio.h>

#include<cstring>

#include<stack>

#define maxn 536870912

using namespace std;

int mapp[20][20];

bool b[20];

int d[20],fa[20];

int t(42),n(17),x,y,e,v,minn;

void print(int x)

{

    stack<int>q;

    int flag(0);

    while(x!=0)

    {

        q.push(x);

        if(x==7||x==12)flag++;

        x=fa[x];

    }

    q.push(0);

    if(flag!=2)return ;

    while(q.empty()==0)

    {

        cout<<q.top()<<"->";

        q.pop();

    }

    cout<<"17\n";

}

void dfs(int x)

{

    if(x==17)

    {

        print(fa[x]);

        return;

    }

    for(int i=0;i<=n;i++)

        if(b[i]==0&&mapp[x][i]<maxn)

        {

            b[i]=1;

            fa[i]=x;

            dfs(i);

            b[i]=0;

        }

}

int main()

{

    //freopen("input.txt","r",stdin);

    //freopen("output.txt","w",stdout);//若采用文件输入输出可以去掉注释

    for(int i=0;i<=n;i++)

      for(int j=0;j<=n;j++)

        if(i==j)mapp[i][j]=0;else mapp[i][j]=maxn;//读入图的边,存入邻接矩阵

    while(t--)

    {

        cin>>x>>y>>e;

        mapp[x][y]=mapp[y][x]=e;

    }

    mapp[11][12]=mapp[12][11]=maxn;

    memset(b,0,sizeof(b));

    memset(fa,0,sizeof(fa));

    dfs(0);

   return 0;

}

⑥储物间个数扩大的是搜索路径的代码

#include <iostream>

#include<stdio.h>

#include<cstring>

#include<stack>

#define maxn 536870912

using namespace std;

int mapp[20][20];

bool b[20];

int d[20],fa[20];

int t(42),n(17),x,y,e,v,minn;

void print(int x)

{

    stack<int>q;

    int flag(0);

    while(x!=0)

    {

        q.push(x);

        if(x==7||x==12)flag++;

        if((x==4&&fa[x]==2)||(x==2&&fa[x]==4))flag++;

        if((x==14&&fa[x]==13)||(x==13&&fa[x]==14))flag++;

        x=fa[x];

    }

    q.push(0);

    if(flag!=4||q.size()>11)return ;//要实现寻找10、11、12个点的路径,只需把q.size()>后面的数字修改即可

    while(q.empty()==0)

    {

        cout<<q.top()<<"->";

        q.pop();

    }

    cout<<"17\n";

}

void dfs(int x)

{

    if(x==17)//如果搜索到终点就打印路径并输出

    {

        print(fa[x]);

        return;

    }

    for(int i=0;i<=n;i++)

        if(b[i]==0&&mapp[x][i]<maxn)

        {

            b[i]=1;//选择这个点的方向走

            fa[i]=x;//记录父节点

            dfs(i);

            b[i]=0;//不选这个点的情况

        }

}

int main()

{

    //freopen("input.txt","r",stdin);//若采用文件输入输出可以去掉注释

    //freopen("output.txt","w",stdout);

    for(int i=0;i<=n;i++)

      for(int j=0;j<=n;j++)

        if(i==j)mapp[i][j]=0;else mapp[i][j]=maxn;//读入图的边,存入邻接矩阵

    while(t--)

    {

        cin>>x>>y>>e;

        mapp[x][y]=mapp[y][x]=e;

    }

    mapp[11][12]=mapp[12][11]=maxn;

    memset(b,0,sizeof(b));

    memset(fa,0,sizeof(fa));

    dfs(0);

   return 0;

}

       

        后来了解到更好的做法是用模拟退火算法来寻优,但是对于本题,似乎贪心算法求出的初始解基本接近最优解,所以......嘻嘻O(∩_∩)O

  • 4
    点赞
  • 24
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值