//最短路徑(狄克斯特拉算法)PshortP1.cpp
//從一個頂點到其余各頂點的最短路徑
#include<iostream.h>
#include<iomanip.h>
#include<stdlib.h>
#include "graph.cpp"
//網G從下標v0到其他頂點的最短距離dist和最短路徑下標path
void PShortPath(AdjMatrix &G,int v0,int dist[],int path[])
{int n=G.NumV();
 int *s=new int[n];
 int mindis,i,j,u;
 for(i=0;i<n;i++)
 {dist[i]=G.GetWeight(v0,i);
  s[i]=0;
  if(i!=v0&&dist[i]<MaxValue) path[i]=v0;
  else path[i]=-1;
 }
 s[v0]=1;//標記頂點v0已從集合T加入到集合S中
 //在當前還未找到最短路徑的頂點集中選取具有最短距離的頂點u
 for(i=1;i<n;i++)
 {mindis=MaxValue;
  for(j=0;j<n;j++)
   if(s[j]==0&&dist[j]<mindis)
    {u=j;
     mindis=dist[j];
    }
  //當已不再存在路徑時算法結束;此語句對非連通圖是必需的
  if(mindis==MaxValue) return;
  s[u]=1;//標記頂點u已從集合T加入到集合S中
  //修改從v0到其他頂點的最短距離和最短路徑
  for(j=0;j<n;j++)
   if(s[j]==0&&G.GetWeight(u,j)<MaxValue&&
    dist[u]+G.GetWeight(u,j)<dist[j])
   {//頂點v0經頂點u到其他頂點的最短距離和最短路徑
    dist[j]=dist[u]+G.GetWeight(u,j);
    path[j]=u;
   }
 }
}
//算法測試
void main()
{cout<<"PShortP1.cpp運行結果:\n";
 int n=6,k1=1,k2=1;
 AdjMatrix g(n,k2);
 g.CreateMatrix(n,k1,k2);
 cout<<"\n輸出鄰接矩陣相應圖的每個頂點:\n";
 g.Creatgraph(n,k2);
 int m=g.NumV();
 int *dist=new int[m];
 int *path=new int[m];
 int v0=0;
 PShortPath(g,v0,dist,path);
 cout<<"從頂點"<<g.GetValue(v0)
     <<"到其他各頂點的最短距離為:\n";
 for(int i=0;i<m;i++)
  cout<<"到頂點"<<g.GetValue(i)
      <<"的最短距離為:"<<dist[i]<<endl;
 cout<<"從頂點"<<g.GetValue(v0)
     <<"到其他各頂點的最短路徑的前一頂點為:\n";
 for(i=0;i<m;i++)
  if(path[i]!=-1)
   cout<<"到頂點"<<g.GetValue(i)<<"的前一頂點為:"
       <<g.GetValue(path[i])<<endl;
 cin.get();cin.get();}