圖解Dijkstra演算法+程式碼實現

Carol淋發表於2022-05-20

簡介

Dijkstra(迪傑斯特拉)演算法是典型的單源最短路徑演算法,用於計算一個節點到其他所有節點的最短路徑。主要特點是以起始點為中心向外層層擴充套件,直到擴充套件到終點為止。Dijkstra演算法是很有代表性的最短路徑演算法,在很多專業課程中都作為基本內容有詳細的介紹,如資料結構,圖論,運籌學等等。注意該演算法要求圖中不存在負權邊。
對應問題:在無向圖G=(V,E)中,假設每條邊E(i)的長度W(i),求由頂點V0到各節點的最短路徑。
image-20220520143004325

工作過程

Dijkstra演算法將頂點集合分為兩組,一組記錄已經求得最短路徑的頂點記為finallyNodes,一組正在求解中的頂點記為processNodes,
step1:finallyNodes中頂點最開始只有源節點,最短路徑長度為0,而processNodes中包含除源節點以外的節點,並初始化路徑長度,與源節點直接相連的記路徑長度為權重,不相連的記為♾️。
step2:從process中選擇路徑長度最小的頂點,加入finallyNodes,並且更新processNodes,將與當前頂點相連的頂點路徑長度更新為min(當前權重,當前頂點最短路徑長度+當前頂點與頂點相連邊權重)。
step3:重複step2,直至processNodes陣列為空。
image-20220520125520468

總體思路

這次我想先描述一下自己的大概思路,下面再寫具體實現。
首先為了方便,我採用的是鄰接表儲存圖結構,鄰接表是一個二維陣列,值儲存權重。根據上面工作過程中描述的內容,我們會有兩個中間集合記錄,finallyNodes記錄的是最終結果,我們只需要將計算的結果往裡面塞即可。但是processNodes卻是一個不斷變化更新的集合,其中的操作包括刪除節點,更改節點值,查詢節點值,同時我們每次需要拿出processNodes中記錄的距離最小的值,所以ProcessNodes準備用最小堆來做,那再刪除節點,更改節點值之後都需要調整堆為最小堆,java自帶的優先佇列沒有提供更改節點值的操作,因此我們這裡需要自己實現一個小根堆,支援以上操作。
然後就中規中矩實現dijkstra演算法即可。

實現

小根堆

如果對堆不太熟悉的可以先看看這篇文章:堆(優先佇列),這裡就不過多解釋了,直接貼程式碼。
這裡堆中存的資料格式為int二維陣列,儲存節點下標位置和對應距離,排序按儲存的距離進行排序。

public class MinHeap {
        List<int[][]> heap ;
        /**
         * 獲取並移除堆頂元素,並調整堆
         * @return
         */
        public int[][] pop() {
            int[][] top = heap.get(0);
            heap.set(0, heap.get(heap.size() - 1));
            heap.remove(heap.size() - 1);
            //調整堆
            this.adjust(0, heap.size() - 1);
            return top;
        }

        /**
         * 判斷是否為空
         * @return true/false
         */
        public boolean isEmpty() {
            if (null == this.heap) {
                return true;
            }
            if (this.heap.size() == 0) {
                return true;
            }
            return false;
        }
        /**
         * 修改index位置節點的value值,並調整最小堆(Java priorityQueue未提供)
         * @param index 修改節點位置
         * @param value 修改值
         */
        public void changeValue(int index, int value) {
            int src = heap.get(index)[0][1];
            heap.get(index)[0][1] = value;
            //直接比較當前值是變大還是變小,然後考慮是向上調整還是向下調整
            //則當前值可能往上移動
            if (src > value) {
                this.upAdjust(index);
                return;
            }
            this.adjust(index, heap.size() - 1);
        }

        public void upAdjust(int index) {
            //依次與雙親節點進行比較,小於雙親節點就直接交換。一直到根節點
            while (index > 0) {
                int parent = index >> 1;
                //雙親節點本來小於當前節點不需要進行調整
                if (heap.get(parent)[0][1] <= heap.get(index)[0][1]) {
                    break;
                }
                swap(index, parent);
                index = parent;
            }
        }
        
        /**
         * 初始化一個最小堆
         * @param nums
         */
        public void init(int[][] nums) {
            heap = new ArrayList<>(nums.length);
            for (int i = 0 ; i < nums.length; i ++) {
                int[][] temp = new int[1][2];
                temp[0][0] = nums[i][0];
                temp[0][1] = nums[i][1];
                heap.add(temp);
            }
            //從最後一個雙親節點開始將堆進行調整
            for (int i = nums.length / 2 ; i >= 0 ; -- i) {
                this.adjust(i, nums.length - 1);
            }
        }
        /**
         * 從當前index開始調節為最小堆
         * @param index 當前節點下標
         * @param end 最後一個節點下標
         */
        private void adjust(int index, int end) {
            //找到當前節點的孩子節點,將較小的節點與當前節點交換,一直往下,直至end
            while (index <= end) {
                //左孩子節點
                int left = index << 1;
                if (left + 1 <= end && heap.get(left + 1)[0][1] < heap.get(left)[0][1] ) {
                    //找到當前較小的節點
                    ++ left;
                }
                //沒有孩子節點,或者當前的孩子節點均已大於當前節點,已符合最小堆,不需要進行調整
                if (left > end || heap.get(index)[0][1] <= heap.get(left)[0][1]) {
                    break;
                }
                swap(index, left);
                index = left;
            }
        }
        private void swap(int i, int j) {
            int[][] temp = heap.get(i);
            heap.set(i, heap.get(j));
            heap.set(j, temp);
        }
    }

Dijsktra

資料結構

圖節點僅儲存節點值,一個Node陣列nodes,儲存圖中所有節點,一個二維陣列adjacencyMatrix,儲存圖中節點之間邊的權重,行和列下標與nodes陣列下標對應。

 //節點
 Node[] nodes;

 //鄰接矩陣
 int[][] adjacencyMatrix;
public class Node {
        private char value;
        Node(char value) {
            this.value = value;
        }
    }

初始化

初始化圖
values標誌的圖中所有節點值,edges標誌圖中邊,資料格式為(node1的下標,node2的下標,邊權重)

private void initGraph(char[] values, String[] edges) {
        nodes = new Node[values.length];
  			//初始化node節點
        for (int i = 0 ; i < values.length ; i ++) {
            nodes[i] = new Node(values[i]);
        }
        adjacencyMatrix = new int[values.length][values.length];
  			//初始化鄰接表,同一個節點權重記為0,不相鄰節點權重記為Integer.MAX_VALUE
        for (int i = 0 ; i < values.length ; i++) {
            for (int j = 0 ; j < values.length ; j ++) {
                if (i == j) {
                    adjacencyMatrix[i][j] = 0;
                    continue;
                }
                adjacencyMatrix[i][j] = Integer.MAX_VALUE;
                adjacencyMatrix[j][i] = Integer.MAX_VALUE;
            }
        }
  			//根據edges更新相鄰節點權重值
        for (String edge : edges) {
            String[] node = edge.split(",");
            int i = Integer.valueOf(node[0]);
            int j = Integer.valueOf(node[1]);
            int weight = Integer.valueOf(node[2]);
            adjacencyMatrix[i][j] = weight;
            adjacencyMatrix[j][i] = weight;
        }
        visited = new boolean[nodes.length];

    }

初始化dijsktra演算法必要的finallyNodes和processNodes

    
		/**
		* 標誌對應下標節點是否已經處理,避免二次處理
		*/
		boolean[] visited;
    /**
     * 記錄已經求得的最短路徑 finallyNodes[0][0]記錄node下標,finallyNodes[0][1]記錄最短路徑長度
     */
    List<int[][]> finallyNodes;
    /**
     * 記錄求解過程目前的路徑長度,因為每次取當前已知最短,所以最小堆進行記錄
     * 但是java優先佇列沒有實現改變值,這裡需要自己實現
     * 首先每次取出堆頂元素之後,堆頂元素加入finallyNodes,此時需要更新與當前元素相鄰節點的路徑長度
     * 然後重新調整小根堆
     * 首先:只會更新變小的資料,所以從變小元素開始往上進行調整,或者直接呼叫調整方法,從堆頂往下進行調整
     */
    MinHeap processNodes;
	/**
     * 初始化,主要初始化finallyNodes和processNodes,finallyNodes加入源節點,processNodes加入其他節點
     * @param nodeIndex
     */
    private void initDijkstra(int nodeIndex) {
        finallyNodes = new ArrayList<>(nodes.length);
        processNodes = new MinHeap();
        int[][] node = new int[1][2];
        node[0][0] = nodeIndex;
        node[0][1] = adjacencyMatrix[nodeIndex][nodeIndex];
        finallyNodes.add(node);
        visited[nodeIndex] = true;
        int[][] process = new int[nodes.length - 1][2];
        int j = 0;
        for (int i = 0 ; i < nodes.length ; i++) {
            if (i == nodeIndex) {
                continue;
            }
            process[j][0] = i;
            process[j][1] = adjacencyMatrix[nodeIndex][i];
            ++ j;
        }
        //初始化最小堆
        processNodes.init(process);
    }

dijsktra演算法實現

public void dijkstra() {
        //1。堆頂取出最小元素,加入finallyNodes
  			//2。將與堆頂元素相連節點距離更新,
        while (!processNodes.isEmpty()) {
            int[][] head = processNodes.pop();
            finallyNodes.add(head);
            int nodeIndex = head[0][0];
            visited[nodeIndex] = true;
            //跟堆頂元素相鄰的元素
            for (int j = 0 ; j < nodes.length ; j ++) {
                //找到相鄰節點
                if (visited[j] || Integer.MAX_VALUE == adjacencyMatrix[nodeIndex][j]) {
                    continue;
                }
                for (int i = 0 ; i < processNodes.heap.size() ; i++) {
                    int[][] node = processNodes.heap.get(i);
                    //找到節點並且值變小,需要調整
                    if (node[0][0] == j && node[0][1] > head[0][1] + adjacencyMatrix[nodeIndex][j]) {
                        processNodes.changeValue(i, head[0][1] + adjacencyMatrix[nodeIndex][j]);
                        break;
                    }
                }
            }

        } 
    }

測試

public static void main(String[] args) {
        char[] values = new char[]{'A','B','C','D','E','F','G','H'};
        String[] edges = new String[]{"0,1,2","0,2,3","0,3,4","1,4,6","2,4,3","3,4,1","4,5,1","4,6,4","5,7,2","6,7,2"};
        Dijkstra dijkstra = new Dijkstra();
        dijkstra.initGraph(values, edges);
        int startNodeIndex = 0;
        dijkstra.initDijkstra(startNodeIndex);
        dijkstra.dijkstra();
        for (int[][] node : dijkstra.finallyNodes) {
            System.out.println(dijkstra.nodes[node[0][0]].value + "距離" + dijkstra.nodes[startNodeIndex].value + "最短路徑為:" + node[0][1]);
        }
    }

image-20220520153626856

相關文章