資料結構：點對之間最短距離--Floyd演算法

標籤：資料結構   動態規劃   floyd   圖   演算法   

                           Floyd演算法

Floyd演算法

    Dijkstra演算法是用於解決單源最短路徑問題的，Floyd演算法則是解決點對之間最短路徑問題的。Floyd演算法的設計策略是動態規劃，而Dijkstra採取的是貪心策略。當然，貪心演算法就是動態規劃的特例。

演算法思想

    點對之間的最短路徑只會有兩種情況：

1. 兩點之間有邊相連，weight(Vi,Vj)即是最小的。
2. 通過另一點：中介點，兩點相連，使weight(Vi,Vk)+weight(Vk,Vj)最小。

Min_Distance(Vi,Vj)=min{weight(Vi,Vj),weight(Vi,Vk)+weight(Vk,Vj)}。正是基於這種背後的邏輯，再加上動態規劃的思想，構成了Floyd演算法。故當Vk取完所有頂點後，Distance(Vi,Vj)即可達到最小。
題外話：代碼本身不重要，演算法思想才是精髓。思想極難得到，而有了思想，稍加經驗即可寫出代碼。向思想的開創者致敬！
思想很難，代碼卻比較簡單，直接上代碼代碼類定義

#include<iostream>  #include<iomanip>#include<stack>using namespace std;#define MAXWEIGHT 100#undef INFINITY#define INFINITY 1000class Graph{private://頂點數  int numV;//邊數  int numE;//鄰接矩陣  int **matrix;public:Graph(int numV);//建圖  void createGraph(int numE);//析構方法  ~Graph();//Floyd演算法void Floyd();//列印鄰接矩陣  void printAdjacentMatrix();//檢查輸入  bool check(int, int, int);};

類實現

//建構函式，指定頂點數目Graph::Graph(int numV){//對輸入的頂點數進行檢測while (numV <= 0){cout << "頂點數有誤！重新輸入 ";cin >> numV;}this->numV = numV;//構建鄰接矩陣，並初始化matrix = new int*[numV];int i, j;for (i = 0; i < numV; i++)matrix[i] = new int[numV];for (i = 0; i < numV; i++)for (j = 0; j < numV; j++){if (i == j)matrix[i][i] = 0;elsematrix[i][j] = INFINITY;}}void Graph::createGraph(int numE){/*對輸入的邊數做檢測一個numV個頂點的有向圖，最多有numV*(numV - 1)條邊*/while (numE < 0 || numE > numV*(numV - 1)){cout << "邊數有問題！重新輸入 ";cin >> numE;}this->numE = numE;int tail, head, weight, i;i = 0;cout << "輸入每條邊的起點(弧尾)、終點(弧頭)和權值" << endl;while (i < numE){cin >> tail >> head >> weight;while (!check(tail, head, weight)){cout << "輸入的邊不正確！請重新輸入 " << endl;cin >> tail >> head >> weight;}matrix[tail][head] = weight;i++;}}Graph::~Graph(){int i;for (i = 0; i < numV; i++)delete[] matrix[i];delete[]matrix;}/*弗洛伊德演算法求各頂點對之間的最短距離及其路徑*/void Graph::Floyd(){//為了不修改鄰接矩陣，多用一個二維數組int **Distance = new int*[numV];int i, j;for (i = 0; i < numV; i++)Distance[i] = new int[numV];//初始化for (i = 0; i < numV; i++)for (j = 0; j < numV; j++)Distance[i][j] = matrix[i][j];//prev數組int **prev = new int*[numV];for (i = 0; i < numV; i++)prev[i] = new int[numV];//初始化prevfor (i = 0; i < numV; i++)for (j = 0; j < numV; j++){if (matrix[i][j] == INFINITY)prev[i][j] = -1;elseprev[i][j] = i;}int d, v;for (v = 0; v < numV; v++)for (i = 0; i < numV; i++)for (j = 0; j < numV; j++){d = Distance[i][v] + Distance[v][j];if (d < Distance[i][j]){Distance[i][j] = d;prev[i][j] = v;}}//列印Distance和prev數組cout << "Distance..." << endl;for (i = 0; i < numV; i++){for (j = 0; j < numV; j++)cout << setw(3) << Distance[i][j];cout << endl;}cout << endl << "prev..." << endl;for (i = 0; i < numV; i++){for (j = 0; j < numV; j++)cout << setw(3) << prev[i][j];cout << endl;}cout << endl;//列印頂點對最短路徑stack<int> s;for (i = 0; i < numV; i++){for (j = 0; j < numV; j++){if (Distance[i][j] == 0);else if (Distance[i][j] == INFINITY)cout << "頂點 " << i << " 到頂點 " << j << " 無路徑！" << endl;else{s.push(j);v = j;do{v = prev[i][v];s.push(v);} while (v != i);//列印路徑cout << "頂點 " << i << " 到頂點 " << j << " 的最短路徑長度是 "<< Distance[i][j] << " ，其路徑序列是...";while (!s.empty()){cout << setw(3) << s.top();s.pop();}cout << endl;}}cout << endl;}//釋放空間for (i = 0; i < numV; i++){delete[] Distance[i];delete[] prev[i];}delete[]Distance;delete[]prev;}//列印鄰接矩陣  void Graph::printAdjacentMatrix(){int i, j;cout.setf(ios::left);cout << setw(7) << " ";for (i = 0; i < numV; i++)cout << setw(7) << i;cout << endl;for (i = 0; i < numV; i++){cout << setw(7) << i;for (j = 0; j < numV; j++)cout << setw(7) << matrix[i][j];cout << endl;}}bool Graph::check(int tail, int head, int weight){if (tail < 0 || tail >= numV || head < 0 || head >= numV|| weight <= 0 || weight >= MAXWEIGHT)return false;return true;}

主函數

int main(){cout << "******Floyd***by David***" << endl;int numV, numE;cout << "建圖..." << endl;cout << "輸入頂點數 ";cin >> numV;Graph graph(numV);cout << "輸入邊數 ";cin >> numE;graph.createGraph(numE);cout << endl << "Floyd..." << endl;graph.Floyd();system("pause");return 0;}

運行


小結Floyd演算法代碼看似很長，其實並不難。代碼中很多都是用於準備工作和輸出。
完整代碼下載：Floyd演算法

轉載請註明出處，本文地址：http://blog.csdn.net/zhangxiangdavaid/article/details/38366923
若有所協助，頂一個哦！
專欄目錄：

• 資料結構與演算法目錄
• c指標