当前位置:   article > 正文

算法我也不知道有没有下一个---一个题目的开端(索引堆与图)

ab5,bc4,cd8 火车前端

病痛了一周,折磨来折磨去,终于还是平静了下来,现在能把上周末"贯穿"学到的最后一个基础数据结构的知识给沉淀沉淀了。也是即将再单位分享的东西:图论。这东西,想当年大二,学校的时候,只是理解或者画图,就已经很费劲了,当时我们平时作业包括后来的期末考试也只是到理解原理的层面,会画图,就行,实现什么的,根本别想。我自己当年也是认真看过,死磕过几个算法的实现:最小生成树、最短路径等,然而,是看不懂的。这么多年过去了,自己编程能力和见识多了,回过头来,发现,其实并没那么难,只不过比平时的业务逻辑稍微复杂些,但是相较于框架的源码复杂度来说,还是远远够不到的。这一次,我先以一道题开始,逐步将相关的知识点讲下来,采取这种方式,最后,会拓展的讲讲其他的实现,恩,开干

零、一些话

本次的题目,其实算是一道功能实现题,但是能很直接的看出使用的数据结构和大二数据结构课本中介绍的这种数据结构的基本几个之一的算法,所以非常考察当年本科计算机基础的夯实程度与编程实践能力。还是那句话,不要整天咔咔咔的写个curd就觉得会编程,很牛逼。计算机能解决的问题,远远不止curd,编程能力也远远不止curd。同样,这道题,也是一道面试题(具体出处,我就不说了),成败与否,工资高与低,日子苦逼与否,可能往往取决于一道题,并非危言耸听。有人觉得算法没用,平时用不到,都是api调用,我耗时耗力学,有啥回报呢,我能用来解决什么问题呢?下面我来谈谈我的想法。

算法,我觉得类似于我现在追的一部翻拍武侠剧《倚天屠龙记》中的两部武学经典《九阳神功》与《乾坤大挪移》。张无忌整个过程中似乎并没有学什么具体的武功招式,与外家功夫,只学了这两部内功圣经,就已经干翻六大派高手,干翻赵敏手下的一大堆的援军与玄冥二老,称霸武林。可是这两部武学内功并没有交张无忌一招一式啊!我觉得算法如此。我们现在编程生活中,几乎不用我们去实现链表,更不可能让我们实现一个红黑树,更别提让我们写一个图的深度遍历,因为太多太多的框架都帮助我们实现了。然而,试想一种情况:如果一个开发者,能白板写大顶堆小顶堆,能白板手写红黑树,能立马手动编码实现一道基于数据结构的中等难度的LeetCode上面的习题,难道他不会调用spring的api?难道他不会写个sql去查询数据库?我相信即使他不会,也比一般开发者学的快!更进一步,如果你发现你要根据具体的业务场景实现自己的业务负载均衡模型、在电信的业务模型中对具体的一个表达式进行解析(恩,我弄过)、在游戏开发领域快速判断两点之间的连通性等等这些,都是必须要牢靠掌握数据结构与算法能力的。哪怕是,我们想要更佳平滑的弥补HashMap的多线程坑,不懂数据结构,那是绝对不可能的。

也许百度百度,可能大家都能知道怎么弄。可是自我思考得来的解决与百度来的解决,那是截然不同的两种能力等级!解决问题的速度,也是不一样的。更佳的,当面临系统优化、制定规则引擎、写一些系统工具,那算法与数据结构就是无论如何都无法绕开的了。所以为什么很多大厂面试,为啥愿意面算法,让你手写算法,原因在此。本来就是更上了一个level啊!是否录取一个应聘者,如何开工资,能力等级是多少,高下立判!

一、算法原题

Kiwiland市的铁路系统服务于多个城镇。由于资金问题,所有线路均为“单向线路”。例如,由Kaitaia向Invercargill的方向有一条铁路线,但反方向却没有路线。事实上,即使这两个方向都有铁路线,他们也可能是不同的线路轨迹,里程数也不一定相同。

这道题就是要帮助铁路系统为乘客提供线路信息。具体来讲,你需要计算出某条线路的里程数,两个城镇之间的线路数量以及其中的最短线路。

**Input:**一个有向图(directed graph),其中的节点表示一个城镇,线表示两个城镇间的线路,线的权重表示距离。一条线路不得出现两次或以上,且起点城镇和终点城镇不得相同。

**Output:**测试1至5,如果该线路不存在,则输出 'NO SUCH ROUTE',否则就按照给定的线路行进,注意不要添加任何其他站点! 例如,在下面第1条线路中,必须从A城开始,直接到达B城(距离为5),然后直接到达C城(距离为4)。

  1. 线路A-B-C的距离
  2. 线路A-D的距离
  3. 线路A-D-C的距离
  4. 线路A-E-B-C-D的距离
  5. 线路A-E-D的距离
  6. 从C点开始,在C点结束,沿途最多可以有3个站点,符合该要求的线路有几条? 在下面给出的示例数据中,共有2条符合线路:C-D-C (2 站),C-E-B-C (3 站)
  7. 从A点开始,在C点结束,要求线路上必须有4个站点,符合该要求的线路有几条? 在下面给出的示例数据中,共有3条符合线路:A 到 C (通过 B,C,D); A 到 C (经过 D,C,D); 以及 A 到 C (经过 D,E,B).
  8. A 到 C的最短路线长度
  9. B 到 B的最短路线长度
  10. 从C点开始,在C点结束,要求距离小于30,符合该要求的线路有几条? 在下面给出的示例数据中,符合条件的线路有:CDC, CEBC, CEBCDC, CDCEBC, CDEBC, CEBCEBC, CEBCEBCEBC.

测试输入:

在测试输入中,城镇名字由A、B、C、D代表。例如AB5表示从A到B,线路距离为5.

图表: AB5, BC4, CD8, DC8, DE6, AD5, CE2, EB3, AE7

要求输出:

Output #1: 9

Output #2: 5

Output #3: 13

Output #4: 22

Output #5: NO SUCH ROUTE

Output #6: 2

Output #7: 3

Output #8: 9

Output #9: 9

Output #10: 7

从题目可以看出,其实就是要实现一个最简单的图,并最终实现一个最短路径的算法。当然这个过程中,会涉及一些数据结构与算法:堆(索引最小堆)、邻接矩阵表示图、图的深度优先遍历、有环图的遍历、Dijkastra最短路径算法(正权边)。下面会一个个来说。

二、最小堆

首先我们围绕这道题,从最最基础的结构 — 堆,说起。以往,各种教学,都从大顶堆讲,然后自己私下去拓展学小顶堆,这回我们直接从小顶堆,来入手。首先让我们看看堆的定义(来源于维基百科):

(英语:Heap)是计算机科学中的一种特别的树状数据结构。若是满足以下特性,即可称为堆:“给定堆中任意节点 P 和 C,若 P 是 C 的母节点,那么 P 的值会小于等于(或大于等于) C 的值”。若母节点的值恒小于等于子节点的值,此堆称为最小堆(min heap);反之,若母节点的值恒大于等于子节点的值,此堆称为最大堆(max heap)。在堆中最顶端的那一个节点,称作根节点(root node),根节点本身没有母节点(parent node)。

总结下:

  • 树状结构
  • 完全二叉(多叉)树(除了叶子节点,其他节点必须所有孩子都有值)
  • 任意父亲节点大于(小于)或者等于孩子节点
  • 大于:最大堆;小于:最小堆

下面是堆的基础模型图:

最小堆的示意图

相应的,我们使用的存储堆的数据结构,更多的使用一个数组,如下:

使用数组具体的母子节点的换算公式是:

  • 获取当前节点的父亲节点索引:parentNodeIndex = (currentNodeIndex-1)/2
  • 获取当前节点的孩子索引:
    • 左孩子:leftNodeIndex = currentNodeIndex*2 + 1
    • 右孩子:rightNodeIndex = currentNodeIndex*2 + 2

下面是我总结的大体上堆这种数据结构的优点:

  • 快速的获取最大(最小)值:O(1)复杂度
  • 简单、快速的底层存储逻辑
  • 堆排序(O(nlogn))

三、最小堆的核心实现(基于Java)

先给出基础的骨架代码

  1. public class MinPriorityQueue<K extends Comparable<K>> {
  2. // 堆的基础存储数据结构
  3. private Object[] elements;
  4. // 当前堆中元素个数
  5. private int size;
  6. public MinPriorityQueue() {
  7. this.elements = new Object[Common.INIT_SIZE];
  8. this.size = 0;
  9. }
  10. private void checkLength(int index) {
  11. if (index < 0 || index >= this.elements.length) {
  12. throw new IllegalArgumentException("超出限制");
  13. }
  14. }
  15. public int lchild(int index) {
  16. checkLength(index);
  17. return (index * 2) + 1;
  18. }
  19. public int rchild(int index) {
  20. checkLength(index);
  21. return (index * 2) + 2;
  22. }
  23. public int parent(int index) {
  24. checkLength(index);
  25. return (index - 1) / 2;
  26. }
  27. public void swap(int i, int j) {
  28. checkLength(i);
  29. checkLength(j);
  30. K temp = (K) elements[i];
  31. elements[i] = elements[j];
  32. elements[j] = temp;
  33. }
  34. @Override
  35. public String toString() {
  36. return "MinPriorityQueue{" +
  37. "elements=" + Arrays.toString(elements) +
  38. ", size=" + size +
  39. '}';
  40. }
  41. }

有几点:

  • 内部存储由于擦除原因,使用了Object
  • 每个泛型参数类型都要是可比较类型(Comparable),因为要比较父子节点大小
  • size表示当前堆中元素大小
  • 当然,这里可以直接使用Comparable类型的基础数组存储

1、上浮操作与添加元素

上浮是保证整个堆保证原先数据规则的一种手段,主要用于添加元素的时候,每次添加元素,我们会直接放到内部基础数组最后一个索引位置,然后做上浮操作,下面我们随便添加一个数字(9),然后看看整体上浮过程是如何变化的:

整个过程比较简单,主要是当前节点和父亲加点比较:

  1. /**
  2. * 上浮
  3. */
  4. public void swam(int index) {
  5. // 防止索引溢出
  6. checkLength(index);
  7. int p = index;
  8. while (p > 0) {
  9. K currentEle = (K) elements[p];
  10. int parentIndex = parent(p);
  11. // 核心,用于比对当前节点与父亲节点的大小
  12. if (currentEle.compareTo((K) elements[parentIndex]) < 0) {
  13. // 交换当前索引与父亲节点索引中的值
  14. swap(p, parentIndex);
  15. p = parentIndex;
  16. } else {
  17. // 一旦发现当前值大于父亲节点,就停止循环
  18. break;
  19. }
  20. }
  21. }

有了上浮操作,那么向堆中添加元素的代码也是比较简单的了:

  1. /**
  2. * 添加元素
  3. */
  4. public void addElement(K element) {
  5. if (size == elements.length) {
  6. // 添加堆的容量
  7. resize();
  8. }
  9. elements[size] = element;
  10. swam(size);
  11. // 维护当前堆数量大小
  12. size++;
  13. }
  14. /**
  15. * 扩容
  16. */
  17. public void resize() {
  18. int length = elements.length;
  19. Object[] tempArr = new Object[length + Common.GROW_SIZE];
  20. System.arraycopy(elements, 0, tempArr, 0, length);
  21. this.elements = tempArr;
  22. }

2、下沉操作与删除元素

相对比较来说,下沉会更加复杂一些,复杂点在于,每次都要进行当前节点与几个孩子节点的比较操作,可是当前节点是否有孩子节点,还是要判断的。总的来说,要花心思去保护索引溢出的情况。首先,让我们看看删除一个堆内元素是如何进行的:

  • 保存堆顶元素的值
  • 交换堆顶和最后的结点
  • 删除最后一个节点的值,指向null
  • 对堆顶(索引为0)做下沉操作

下面是下沉操作的基础代码:

  1. /**
  2. * 下沉
  3. */
  4. public void sink(int index) {
  5. checkLength(index);
  6. int p = index;
  7. /**
  8. * 如果左孩子大于或者等于当前堆中元素个数的话,
  9. * 表明当前节点已经是叶子节点,可以不用继续遍历了
  10. */
  11. while (p < size && lchild(p) < size) {
  12. K currentEle = (K) elements[p];
  13. // 获取左孩子索引
  14. int lchild = lchild(p);
  15. // 获取右孩子索引
  16. int rchild = rchild(p);
  17. // 获取左孩子的值
  18. K lelement = (K) elements[lchild];
  19. // 获取右孩子的值
  20. K relement = (K) elements[rchild];
  21. // 核心,获取左右孩子最小值的那个节点索引,注意,这里右孩子有可能为空
  22. int minChildIndex = Objects.isNull(relement) ? lchild :
  23. (lelement.compareTo(relement) > 0 ? rchild : lchild);
  24. // 使用当前节点值与最小值比较,如果当前值还要小,那就交换
  25. if (currentEle.compareTo((K) elements[minChildIndex]) > 0) {
  26. swap(p, minChildIndex);
  27. p = minChildIndex;
  28. } else {
  29. break;
  30. }
  31. }
  32. }

在这里,主要对节点索引进行了详尽的判断与保护,接下来,就是删除元素的方法代码:

  1. /**
  2. * 删除堆顶元素
  3. */
  4. public K delElemet() {
  5. if(size == 0){
  6. throw new IllegalArgumentException("当前堆已经空了");
  7. }
  8. K result = (K) elements[0];
  9. if (size > 1) {
  10. swap(0, size - 1);
  11. elements[--size] = null;
  12. sink(0);
  13. } else {
  14. size--;
  15. }
  16. return result;
  17. }

3、调整数组成为一个最小堆

通常情况下,我们接收到的是一整个数组的值,然后需要我们整理,才能变成一个最小堆,这个过程我们称它为:heapify。就是重新调整。让我们来看下一组数组Integer[] arr = {90,34,99,57,11,67,55,23,76,33,45}如何进行从新调整的:

做一些介绍:

  • 每次从第一个非叶子节点开始进行每个节点的下沉操作
  • 第一个非叶子节点的索引 = (最后一个节点索引-1)/2
  • 整体时间复杂度是:O(n)

下面就是基础的代码实现:

  1. // 每次传入数组的大小
  2. public MinPriorityQueue(Comparable[] comps) {
  3. if (Objects.isNull(comps) || comps.length == 0) {
  4. throw new IllegalArgumentException("数组不能为空");
  5. }
  6. this.elements = new Object[comps.length];
  7. this.size = 0;
  8. for (int i = 0; i < comps.length && comps[i] != null; i++) {
  9. this.elements[i] = comps[i];
  10. this.size++;
  11. }
  12. heapify();
  13. }
  14. private void heapify() {
  15. if (this.size == 0) {
  16. return;
  17. }
  18. int index = (this.size - 2) / 2;
  19. for (int i = index; i >= 0; i--) {
  20. // 下沉
  21. sink(i);
  22. }
  23. }

到此,我们对于一个最小堆的全部实现,就完成了。下面开始介绍基于堆的一种拓展数据结构 — 索引最小堆

四、索引最小堆

经历了最小堆的"洗礼",发现如下几个问题:

  • 每次交换都是使用原始元素进行交换
  • 没提供修改固定索引元素的方法
  • 如果每个节点元素是一个复杂且庞大的值,那么交换过程会导致很多问题(慢、内存溢出)

就此,索引堆概念出现了:

  • 每个实际元素的数组中的索引位置不变
  • 另申请一个与堆同等大小的数组,存储每个元素在堆中实际位置
  • 再来一个同等大小的数组,反向索引具体堆中的位置(数组索引 —> 堆上的索引)

(最后一点最不好理解,后面会有介绍)如此一来,每次元素添加到堆中,具体的元素就不会随着上浮下沉而变换位置,真正变换的,是索引值(就是一个整型)。

1、索引最小堆演示

下面我根据上面最小堆,进一步扩展成索引最小堆,一步步演示,我们使用的原始数组是:

Integer[] arr = {90,34,99,57,11,67,55,23,76,33,45}

(1)不加反向索引

可以看到全程没有进行存储元素的移动,全部是元素所对应的索引值进行移动,这样一来,就能很好的解决上面元素复杂,移动缓慢的问题。

(2)加上反向索引

上面的索引堆还存在一个问题,就是,我更改堆中对应一个具体索引的元素值,之后,我并不知道这个值对应的具体在堆中的位置,例如:

  1. 改变索引4中元素11的值:data[4] = 12
  2. 我们想要确定11这个元素堆中的位置(其实是在索引数组的第一个位置上,就是0上面)
  3. 然后根据这个位置,对这个位置前面的进行上浮操作,这个位置后面的进行下沉操作,这样就能从新调整堆

这个时候,反向索引就有用了!我下面展现下最终加上反向索引的最终形式:

如果i表示具体索引值,recs表示反向索引数组,indexes表示索引数组,这几个的关系是:

  • recs[indexes[i]] = i
  • indexes[recs[i]] = i

反向数组不太好理解,可以多看看,多看看。其实就是一个存储当前索引对应的数据,在堆中的位置,这个位置其实也是一个索引,只不过这个索引指向的是索引数组

2、索引最小堆的基础实现

首先是基础框架的代码:

  1. public class IndexMinPriorityQueue<K extends Comparable<K>> {
  2. // 索引数组
  3. private int[] indexes;
  4. // 原始数据的数组
  5. private Object[] elements;
  6. // 反向索引数组
  7. private int[] recs;
  8. // 当前堆元素的大小
  9. private int size;
  10. public IndexMinPriorityQueue() {
  11. this.indexes = new int[Common.INIT_SIZE];
  12. this.elements = new Object[Common.INIT_SIZE];
  13. this.recs = new int[Common.INIT_SIZE];
  14. this.size = 0;
  15. for (int i = 0; i < Common.INIT_SIZE; i++) {
  16. this.indexes[i] = -1;
  17. this.recs[i] = -1;
  18. }
  19. }
  20. public IndexMinPriorityQueue(int count) {
  21. this.indexes = new int[count];
  22. this.elements = new Object[count];
  23. this.recs = new int[count];
  24. this.size = 0;
  25. for (int i = 0; i < count; i++) {
  26. // 默认初始化时候,堆内没有元素时候两个索引数组的每个字段都初始化成-1
  27. this.indexes[i] = -1;
  28. this.recs[i] = -1;
  29. }
  30. }
  31. }

加入了两个辅助性质的数组,我们来看看其他的一些修改点,首先是交换数据的方法:

  1. public void swap(int i, int j) {
  2. if (i < 0 || i >= size) {
  3. throw new IllegalArgumentException("超出限制");
  4. }
  5. if (j < 0 || j >= size) {
  6. throw new IllegalArgumentException("超出限制");
  7. }
  8. // 可以看到,每次只交换索引数组中的值,真实数据数组不变
  9. int temp = indexes[i];
  10. indexes[i] = indexes[j];
  11. indexes[j] = temp;
  12. // 反向索引数组内部值的确认
  13. recs[indexes[i]] = i;
  14. recs[indexes[j]] = j;
  15. }

相应的,上浮与下沉的方法也是要变动的:

  1. /**
  2. * 下沉
  3. */
  4. public void sink(int index) {
  5. if (index < 0 || index > size) {
  6. throw new IllegalArgumentException("超出限制");
  7. }
  8. int p = index;
  9. while (p < size && lchild(p) < size) {
  10. // 可以看到,每次取元素,都是要据具体的索引数组的值来定位真实数据数组的位置
  11. K currentEle = (K) elements[indexes[p]];
  12. int lchild = lchild(p);
  13. int rchild = rchild(p);
  14. K lelement = (K) elements[indexes[lchild]];
  15. // 索引数组的值为-1,表示当前元素被删除或者就是没存入
  16. int minChildIndex = indexes[rchild] == -1 ? lchild :
  17. (lelement.compareTo((K) elements[indexes[rchild]]) > 0
  18. ? rchild : lchild);
  19. if (currentEle.compareTo((K) elements[indexes[minChildIndex]]) > 0) {
  20. swap(p, minChildIndex);
  21. p = minChildIndex;
  22. } else {
  23. break;
  24. }
  25. }
  26. }
  27. /**
  28. * 上浮
  29. */
  30. public void swam(int index) {
  31. if (index < 0 || index >= size) {
  32. throw new IllegalArgumentException("超出限制");
  33. }
  34. int p = index;
  35. while (p > 0) {
  36. K currentEle = (K) elements[indexes[p]];
  37. int parentIndex = parent(p);
  38. K pelement = (K) elements[indexes[parentIndex]];
  39. if (Objects.nonNull(pelement) && currentEle.compareTo(pelement) < 0) {
  40. swap(p, parentIndex);
  41. p = parentIndex;
  42. } else {
  43. break;
  44. }
  45. }
  46. }

添加是对添加蒜素的修改:

  1. /**
  2. * 添加元素
  3. */
  4. public void addElement(K element) {
  5. if (size == elements.length) {
  6. // 扩容
  7. resize();
  8. }
  9. // 索引数组的末尾添加这个元素的索引值
  10. indexes[size] = size;
  11. // 反向索引也是最后一位添加
  12. recs[size] = size;
  13. // 其实这也是最后一位添加,因为此时indexes[size] = size
  14. elements[indexes[size]] = element;
  15. // 先对最后一位的索引进行上浮操作,然后再将size加一
  16. swam(size++);
  17. }
  18. /**
  19. * 扩容
  20. */
  21. public void resize() {
  22. int length = elements.length;
  23. Object[] tempArr = new Object[length + Common.GROW_SIZE];
  24. int[] tempIndexes = new int[length + Common.GROW_SIZE];
  25. int[] tempRecs = new int[length + Common.GROW_SIZE];
  26. System.arraycopy(elements, 0, tempArr, 0, length);
  27. System.arraycopy(indexes, 0, tempIndexes, 0, length);
  28. System.arraycopy(recs, 0, tempRecs, 0, length);
  29. for (int i = length; i < tempIndexes.length; i++) {
  30. tempIndexes[i] = -1;
  31. tempRecs[i] = -1;
  32. }
  33. this.elements = tempArr;
  34. this.indexes = tempIndexes;
  35. this.recs = tempRecs;
  36. }

下面看看删除一个堆顶元素的修改:

  1. /**
  2. * 删除堆顶元素
  3. *
  4. * @return 堆顶的具体元素值
  5. */
  6. public K delElemet() {
  7. K result = (K) elements[indexes[0]];
  8. delEl();
  9. return result;
  10. }
  11. /**
  12. * 删除堆顶元素
  13. *
  14. * @return 堆顶的具体元素索引值
  15. */
  16. public int delElemetIndex() {
  17. int result = indexes[0];
  18. delEl();
  19. return result;
  20. }
  21. private void delEl() {
  22. if (size > 1) {
  23. // 这里其实是交换索引数组的第一位和最后一位的值
  24. swap(0, size - 1);
  25. }
  26. // 此时要把末尾的索引值对应的元素置空,代表删除原始数据
  27. elements[indexes[--size]] = null;
  28. // 当然,反向索引数组的值也要删除,置位-1
  29. recs[indexes[size]] = -1;
  30. // 当然,最后要把索引数组值删除,其实就是最后一位
  31. indexes[size] = -1;
  32. // 对索引数组,就是第一位开始,做下沉
  33. sink(0);
  34. }

最后的最后,我们就可以看看如何修改一个堆中的元素了

  1. public void change(int index, K element) {
  2. if (index < 0 || index > size) {
  3. throw new IllegalArgumentException("超出限制");
  4. }
  5. this.elements[index] = element;
  6. // 这时候,反向索引数组就显示作用了:获取这个修改值对应的堆中的索引值
  7. int currentHeapIndex = this.recs[index];
  8. swam(currentHeapIndex);
  9. sink(currentHeapIndex);
  10. }

五、图的基础

好吧,到了图这里,已经深入到计算机科学领域了,我作为一个工程狗,真心无法,也没有能力一箭到底!所以在此,我先上定义,然后我们截取最简单,切合具体的实际问题(这里就是第一章那道题)来说相关的图论的基础结构与算法,其他的,有关图论的东西太多太多,有兴趣可以自己深入研究。下面是wiki上面的一个定义:

图有多种变体,包括简单图、多重图、有向图、无向图等,但大体上有以下两种定义方式。

二元组的定义

一张 G 是一个二元组(V,E),其中V称为顶点集,E称为边集。它们亦可写成V(G)E(G)E的元素是一个二元组数对,用(x,y)表示,其中x,y \in V

三元组的定义

一张 G 是一个三元组(V,E,I),其中V称为顶集(Vertices set),E称为边集(Edges set),EV不相交;I称为关联函数,IE中的每一个元素映射到V\times V。如果I(e)=(u,v) (e\in E, u,v \in V)那么称边e连接顶点u,v,而u,v则称作e的端点,u,v此时关于e相邻。同时,若两条边有一个公共顶点,则称关于相邻。

恩,说实话,我也看不懂,具体对于我们这次要解决的实际问题,相关的数据结构的知识点,提取出以下几点:

  • 我们要实现一个有向图
  • 而且是有向加权图
  • 会有环
  • 边用E(edge)来表示
  • 节点用V(Vertices)来表示
  • 我们这里固定使用一种表示方法来进行图的存储:邻接表

下面我们就开搞!

1、有向图

我们来看看题目中,要我们实现的有向图的基本输入是:

AB5, BC4, CD8, DC8, DE6, AD5, CE2, EB3, AE7

根据此文字表述,我们画出如下的示意图:

有话说:不要在意边的长短与权值大小的对应比例,只是一种展示、展示、展示!

接下来我们就用代码来实现相应的V(定点)与E(边)

  1. // 定点的抽象数据定义
  2. public interface IVertices<V> {
  3. // 获取这个节点的值
  4. V getData();
  5. // 获取当前节点的索引值
  6. int getIndex();
  7. }
  8. // 边的抽象数据定义
  9. public interface IEdge<V> {
  10. // 获取一条边起始顶点的对象
  11. V getV1();
  12. // 获取一条边结束顶点的对象
  13. V getV2();
  14. // 获取一条边的权值
  15. int getWeight();
  16. }

2、邻接表存储表示法

最常用的图的表示法有两种:邻接矩阵、邻接表。响应的适用场景如下:

  • 邻接矩阵:适用于稠密图
  • 邻接表:适用于稀疏图

一般来说,我们解决实际问题,都是稀疏图,所以常用邻居表,本次我们只看邻接表的表示法,来存储一个图。下面是基本的邻接表示意图:

根据上面这个图的展示,整个图使用邻接表存储,会有下面的几点:

  • 一个图对象保存一个map
  • key值是图中的顶点
  • value值是list,list中的每个值是以key值顶点为起始的边
  • 顶点对象要实现hashCode与equals方法,因为要当key

下面就是一个图的基本代码实现,使用邻接表的方式:

  1. public class Graph<V extends IVertices, E extends IEdge<V>> {
  2. // 当前图中的结点数目
  3. private int vsize;
  4. // 当前图中的边的数目
  5. private int esize;
  6. // 邻接表
  7. private Map<V, List<E>> vectors;
  8. public Graph() {
  9. this.vsize = 0;
  10. this.esize = 0;
  11. this.vectors = new HashMap<>();
  12. }
  13. // 根据所有的结点来初始化
  14. public Graph(V[] datas) {
  15. if (Objects.isNull(datas) || datas.length == 0) {
  16. throw new IllegalArgumentException("初始化失败");
  17. }
  18. this.vsize = datas.length;
  19. this.esize = 0;
  20. this.vectors = new HashMap<>();
  21. for (int i = 0; i < this.vsize; i++) {
  22. this.vectors.put(datas[i], new LinkedList<>());
  23. }
  24. }
  25. // 添加一条边
  26. public void addEdge(E w) {
  27. List<E> ts = this.vectors.get(w.getV1());
  28. if (ts == null) {
  29. throw new IllegalArgumentException("这个节点不存在");
  30. }
  31. if (!ts.contains(w)) {
  32. ts.add(w);
  33. this.esize++;
  34. }
  35. }
  36. // 获取总节点数
  37. public int getVsize() {
  38. return vsize;
  39. }
  40. // 获取总边数
  41. public int getEsize() {
  42. return esize;
  43. }
  44. public Map<V, List<E>> getVectors() {
  45. return vectors;
  46. }
  47. }

六、图的遍历(深度优先)

根据题目,我这里介绍深度优先遍历的整个过程与实现。当然还有广度优先遍历,会在后面的最短路径上深入介绍。

1、题目节点遍历展示

我们随便找出一个小题:

从C点开始,在C点结束,沿途最多可以有3个站点,符合该要求的线路有几条? 在下面给出的示例数据中,共有2条符合线路:C-D-C (2 站),C-E-B-C (3 站)

从上面的图,我们来看下整体深度遍历是如何进行的:

细心的看这个动态图,很好的展示了,整个深度遍历的求解过程。

2、代码实现

下面我就用代码来实现这一过程。相对来说,要实现这个过程并非那么容易,因为要控制类似于:"最多","刚好","最多长度"等这些个区间求解。所以,要有响应的访问计数值,来进行辅助。首先我们来看看深度遍历图的对象的基础:

  1. public class WordDepestPath {
  2. // 要遍历的图
  3. private Graph<WordVector, WordEdge> graph;
  4. // 记录当前边被访问到的次数
  5. private Map<WordEdge, Integer> edgeVisitedCount;
  6. // 记录当前顶点被访问到的次数
  7. private Map<WordVector, Integer> verticesVisitedCount;
  8. public WordDepestPath(Graph<WordVector, WordEdge> graph) {
  9. this.graph = graph;
  10. this.edgeVisitedCount = new HashMap<>();
  11. this.verticesVisitedCount = new HashMap<>();
  12. // 获取图的邻接表
  13. Map<WordVector, List<WordEdge>> vectors = graph.getVectors();
  14. // 遍历所有的边,进行初始化
  15. for (Map.Entry<WordVector, List<WordEdge>> entry : vectors.entrySet()) {
  16. List<WordEdge> value = entry.getValue();
  17. for (WordEdge it : value) {
  18. // 起始状态顶点和边都没有被访问到一次
  19. edgeVisitedCount.put(it, 0);
  20. verticesVisitedCount.put(it.getV1(), 0);
  21. }
  22. }
  23. }
  24. }

这里的WordVector、WordEdge分别是IVertices、IEdge的实现,主要顶点内包装了个字符串,例如:"A"。当然,因为要作为map的key值,必须要实现hashcode与equals方法,如下代码:

  1. public class WordEdge implements IEdge<WordVector>, Comparable<WordEdge> {
  2. private WordVector v1, v2;
  3. private int weight;
  4. public WordEdge() {
  5. }
  6. public WordEdge(WordVector v1, WordVector v2, int weight) {
  7. this.v1 = v1;
  8. this.v2 = v2;
  9. this.weight = weight;
  10. }
  11. @Override
  12. public WordVector getV1() {
  13. return v1;
  14. }
  15. @Override
  16. public WordVector getV2() {
  17. return v2;
  18. }
  19. @Override
  20. public int getWeight() {
  21. return weight;
  22. }
  23. @Override
  24. public boolean equals(Object o) {
  25. if (this == o) return true;
  26. if (o == null || getClass() != o.getClass()) return false;
  27. WordEdge wordEdge = (WordEdge) o;
  28. return weight == wordEdge.weight &&
  29. Objects.equals(v1, wordEdge.v1) &&
  30. Objects.equals(v2, wordEdge.v2);
  31. }
  32. @Override
  33. public int hashCode() {
  34. return Objects.hash(v1, v2, weight);
  35. }
  36. @Override
  37. public String toString() {
  38. return "WordEdge{" +
  39. "v1=" + v1 +
  40. ", v2=" + v2 +
  41. ", weight=" + weight +
  42. '}';
  43. }
  44. @Override
  45. public int compareTo(WordEdge o) {
  46. return this.weight - o.getWeight();
  47. }
  48. }
  49. public class WordVector implements IVertices<String>, Comparable<WordVector>{
  50. private int index;
  51. private String word;
  52. public WordVector() {
  53. }
  54. public WordVector(int index, String word) {
  55. this.index = index;
  56. this.word = word;
  57. }
  58. @Override
  59. public int getIndex() {
  60. return index;
  61. }
  62. @Override
  63. public boolean equals(Object o) {
  64. if (this == o) return true;
  65. if (o == null || getClass() != o.getClass()) return false;
  66. WordVector that = (WordVector) o;
  67. return index == that.index &&
  68. Objects.equals(word, that.word);
  69. }
  70. @Override
  71. public int hashCode() {
  72. return Objects.hash(index, word);
  73. }
  74. @Override
  75. public String toString() {
  76. return word;
  77. }
  78. @Override
  79. public int compareTo(WordVector o) {
  80. return this.getIndex() - o.getIndex();
  81. }
  82. @Override
  83. public String getData() {
  84. return this.word;
  85. }
  86. }

下面是深度遍历的入口方法:

  1. public List<List<WordEdge>> depestPath(WordVector src, WordVector dest) {
  2. List<List<WordEdge>> result = new ArrayList<>();
  3. Map<WordVector, List<WordEdge>> vectors = this.graph.getVectors();
  4. List<WordEdge> wordEdgesSrc = vectors.get(src);
  5. List<WordEdge> wordEdgesDest = vectors.get(dest);
  6. if (wordEdgesSrc == null || wordEdgesDest == null) {
  7. throw new IllegalArgumentException("没有此次搜索路径的结点");
  8. }
  9. for (WordEdge edge : wordEdgesSrc) {
  10. Stack<WordEdge> stack = new Stack<>();
  11. stack.push(edge);
  12. edgeVisitedCount.put(edge, edgeVisitedCount.get(edge) + 1);
  13. depestPathLargest3(edge, dest, 1, stack, result);
  14. edgeVisitedCount.put(stack.peek(), 0);
  15. stack.pop();
  16. }
  17. return result;
  18. }

最后是我们核心深度递归遍历图的方法:

  1. /**
  2. * 最多三站地的路径求解方法
  3. * 解决思路:
  4. * 如果整个路径上面最多只允许有三个顶点(除开起始节点)
  5. * 那么就是说,如果是一个有环有向图,那么这个解的边
  6. * 最多只能被访问一次。如果多被访问一次,那么就会超出
  7. * 三个顶点的题目要求,所以整个过程,重点要控制这个边
  8. * 的访问次数
  9. * @param currentEdge 当前遍历到的边
  10. * @param dest 目标顶点
  11. * @param largest 当前遍历到的站点的数量
  12. * @param stack 栈,用于辅助保存遍历节点
  13. * @param result 结果集,保存了所有符合条件的路径
  14. */
  15. private void depestPathLargest3(WordEdge currentEdge, WordVector dest, int largest,
  16. Stack<WordEdge> stack, List<List<WordEdge>> result) {
  17. // 递归终止条件
  18. if (edgeVisitedCount.get(stack.peek()).intValue() >= 2) {
  19. stack.pop();
  20. return;
  21. }
  22. WordVector v2 = currentEdge.getV2();
  23. // 注意这里:如果到了目的顶点,但是不满足最长路径值,同样不能成为结果
  24. if (v2.getData().equals(dest.getData()) && largest <= 3) {
  25. // 此分支,就是表示一条我们求解的路径
  26. ArrayList<WordEdge> rightPath = new ArrayList<>();
  27. result.add(rightPath);
  28. for (WordEdge it : stack) {
  29. rightPath.add(it);
  30. }
  31. }
  32. List<WordEdge> edges = this.graph.getVectors().get(v2);
  33. // 这种情况表示下个顶点没有临边,那就要栈顶的边不可走,要出栈
  34. if (edges.isEmpty()) {
  35. stack.pop();
  36. return;
  37. }
  38. // 对当前层遍历到的边的结束顶点,进行邻接边的遍历
  39. for (WordEdge it : edges) {
  40. // 每每遍历到了一个邻接边,要进行入栈+计数
  41. stack.push(it);
  42. edgeVisitedCount.put(it, edgeVisitedCount.get(it) + 1);
  43. // 开始递归
  44. depestPathLargest3(it, dest, largest + 1, stack, result);
  45. }
  46. /**
  47. * 注意这里,走到了这个地方,表示遍历邻边已经结束,
  48. * 每次弹出一个值的之前,要把这个邻边对应的访问计数
  49. * 清零,为了不影响后面的递归结束判断
  50. */
  51. edgeVisitedCount.put(stack.peek(), 0);
  52. stack.pop();
  53. }

根据这种思路,我们要求解下面两个问题,也是依葫芦画瓢了:

  1. 从A点开始,在C点结束,要求线路上必须有4个站点,符合该要求的线路有几条? 在下面给出的示例数据中,共有3条符合线路:A 到 C (通过 B,C,D); A 到 C (经过 D,C,D); 以及 A 到 C (经过 D,E,B).
  2. 从C点开始,在C点结束,要求距离小于30,符合该要求的线路有几条? 在下面给出的示例数据中,符合条件的线路有:CDC, CEBC, CEBCDC, CDCEBC, CDEBC, CEBCEBC, CEBCEBCEBC.

七、图的最短路径问题

我们这里求最短路径,使用最经典的迪杰克斯拉算法(Dijkstra),这个算法,用来求单源最短路径使用,有一定的限制:

  • 所有边的权值只能是正数,不能有负权边
  • 每次要使用索引堆进行辅助

下面我们一步步来讲解下,如何求上面算法题中的一小题,就是最短路径的问题:

A 到 C的最短路线长度

1、迪杰克斯拉算法展示

下面是动图展示,然后我们一步步说如何求解:

一些前置条件:

  1. 索引堆中,索引值对应每个节点:1->A,2->B,3->C,4->D,5->E
  2. 索引堆中,具体的真实值存储选取的单源节点到当前节点最短的路径值:A->0,B->5,C->9,D->5,E->7

下面是具体的求解过程:

  1. 首先将原点加入到堆中,那原点对应的值是0
  2. 开始循环从堆中拿堆顶的数据,直到堆为空为止
  3. 每次拿到数据之后,遍历拿到的节点的所有邻边
  4. 拿到每个邻边的另一端的结点,判断这个节点是否被确定下来
  5. 如果没有确定下来,那就看当前从此路径过来的总路径,是否比当前节点的路径要短
  6. 如果短,则跟新当前节点的路径值,如果不短,则跳过
  7. 当然,这个结点有可能还没有统计过路径长度,路径为空,这种时候直接将节点,对应的路径值,压入索引堆中
  8. 如果当前节点已经是确定路径的结点了,那就跳过当前节点与邻边的遍历
  9. 注意:每次从堆顶拿到的,是最小的路径值,那这个路径对应另一端的节点就被确认下来,因为我们没有考虑负权边,所以不可能存在另外一个路径,到这个节点还要短,所以每次堆顶拿到的值,必须要被标记成已确定
  10. 如此遍历,直到堆为空的情况结束

这就是全过程,下面我们来看看代码实现

2、迪杰克斯拉算法的代码实现

  1. public class Dijsktra {
  2. // 我们要求解的原图(邻接表存储)
  3. private Graph<WordVector, WordEdge> graph;
  4. // 索引小顶堆
  5. private IndexMinPriorityQueue<Integer> queue;
  6. // 保存每个节点是否被确认下来的映射,默认是false
  7. private Map<WordVector, Boolean> isMarked;
  8. // 保存源节点到每个节点的最小路径值
  9. private Number[] distTo;
  10. // 保存每个节点的最短路径是从那个邻边到达的
  11. private WordEdge[] from;
  12. public Number[] getDistTo() {
  13. return distTo;
  14. }
  15. public Dijsktra(Graph<WordVector, WordEdge> graph, WordVector src) {
  16. this.graph = graph;
  17. this.queue = new IndexMinPriorityQueue<>(graph.getVsize());
  18. this.isMarked = new HashMap<>();
  19. distTo = new Number[graph.getVsize()];
  20. from = new WordEdge[graph.getVsize()];
  21. // 默认将每个节点的确认映射,设置成false,就是都是未确认状态
  22. for (Map.Entry<WordVector, List<WordEdge>> entry : graph.getVectors().entrySet()) {
  23. WordVector key = entry.getKey();
  24. isMarked.put(key, false);
  25. }
  26. // 初始化最短路径保存数组,与最短路径的邻边数组
  27. for (int i = 0; i < graph.getVsize(); i++) {
  28. distTo[i] = 0.0;
  29. from[i] = null;
  30. }
  31. // 第一个将源节点加入到结构中
  32. from[src.getIndex()] = new WordEdge(src, src, 0);
  33. // 注意,加入的索引值与对应的结点,还有值是路径的长度
  34. this.queue.insert(src.getIndex(), from[src.getIndex()].getWeight());
  35. this.isMarked.put(src, true);
  36. // 开始遍历索引堆
  37. while (!this.queue.isEmpty()) {
  38. // 获取堆顶的元素索引
  39. Integer nodeIndex = this.queue.delElemetIndex();
  40. // 通过索引值与邻边数组,获取对应的当前遍历的堆顶定点
  41. WordVector v2 = from[nodeIndex].getV2();
  42. // 当前节点就是最短路径了,所以标记已被确认
  43. this.isMarked.put(v2, true);
  44. // 开始遍历当前节点的所有邻边
  45. List<WordEdge> edges = this.graph.getVectors().get(v2);
  46. for (WordEdge it : edges) {
  47. // 查询邻边另一边的结点,看看路径情况
  48. WordVector nextNode = it.getV2();
  49. if (!this.isMarked.get(nextNode)) {
  50. int nextNodeIndex = nextNode.getIndex();
  51. /**
  52. * 核心:首先的if逻辑判断很关键,
  53. * 看看当前节点索引对应的路径值有没有开始统计,
  54. * 并且如果开始统计有值的话,就计算从当前邻边到
  55. * 达当前节点的路径长度,是否小于已经存在的路径
  56. * 如果小,就要更新,不小就略过。
  57. */
  58. if (from[nextNodeIndex] == null
  59. || distTo[v2.getIndex()].intValue()
  60. + it.getWeight() < distTo[nextNodeIndex].intValue()) {
  61. // 内部,表示要更新当前节点的最短路径,那就要改各个数组与索引堆中的值
  62. distTo[nextNodeIndex] = distTo[v2.getIndex()].intValue()
  63. + it.getWeight();
  64. from[nextNodeIndex] = it;
  65. // 索引堆也有两种情况,有这个值,没这值
  66. if (queue.contain(nextNodeIndex))
  67. queue.change(nextNodeIndex, distTo[nextNodeIndex].intValue());
  68. else
  69. queue.insert(nextNodeIndex, distTo[nextNodeIndex].intValue());
  70. }
  71. }
  72. }
  73. }
  74. }
  75. }

如此,其实已经解决了我们题目中的所有问题

转载于:https://my.oschina.net/UBW/blog/3039949

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/小舞很执着/article/detail/943796
推荐阅读
相关标签
  

闽ICP备14008679号