void ShortestPath_DIJ(MGraph G,int v0,PathMatrix &P,ShortPathTable &D)
{ // 算法7.15
  // 用Dijkstra算法求有向網G的v0頂點到其余頂點v的最短路徑P[v]
  // 及其帶權長度D[v]。
  // 若P[v][w]為TRUE，則w是從v0到v當前求得最短路徑上的頂點。
  // final[v]為TRUE當且僅當v∈S,即已經求得從v0到v的最短路徑。
  int i=0,j, v,w,min;
  bool final[MAX_VERTEX_NUM];
  for (v=0; v<G.vexnum; ++v) {
    final[v] = FALSE;  
    D[v] = G.arcs[v0][v].adj;
    for (w=0; w<G.vexnum; ++w)  P[v][w] = FALSE;  // 設空路徑
    if (D[v] < INFINITY) { P[v][v0] = TRUE;  P[v][v] = TRUE; }
  }
  D[v0] = 0;  final[v0] = TRUE;        // 初始化，v0頂點屬于S集
  //--- 開始主循環，每次求得v0到某個v頂點的最短路徑，并加v到S集 ---
  for (i=1; i<G.vexnum; ++i) {         // 其余G.vexnum-1個頂點
    min = INFINITY;                    // 當前所知離v0頂點的最近距離
    for (w=0; w<G.vexnum; ++w)
      if (!final[w])                           // w頂點在V-S中
        if (D[w]<min) { v = w;  min = D[w]; }  // w頂點離v0頂點更近
    final[v] = TRUE;                       // 離v0頂點最近的v加入S集
    for (w=0; w<G.vexnum; ++w)             // 更新當前最短路徑及距離
      if (!final[w] && (min+G.arcs[v][w].adj<D[w])) { 
        // 修改D[w]和P[w], w∈V-S
        D[w] = min + G.arcs[v][w].adj;
        for(j=0;j<G.vexnum;j++) P[w][j] = P[v][j]; //第v行賦值于第w行
        P[w][w] = TRUE;   // P[w] = P[v]+[w]
      }//if
  }//for
} // ShortestPath_DIJ