第十三周项目（3）——Dijkstra算法的验证

1. /*      
2. * Copyright (c)2016,烟台大学计算机与控制工程学院      
3. * All rights reserved.      
4. * 文件名称：wu.cpp      
5. * 作    者：武昊      
6. * 完成日期：2016年12月8日      
7. * 版 本 号：v1.0       
8. *问题描述：验证Dijkstra算法 
9. *输入描述：无      
10. *程序输出：测试数据      
11. */

程序中graph.h是图存储结构的“ 算法 库”中的头文件

1. #include <stdio.h>  
2. #include <malloc.h>  
3. #include "graph.h"  
4. #define MaxSize 100  
5. void Ppath(int path[],int i,int v)  //前向递归查找路径上的顶点  
6. {  
7.     int k;  
8.     k=path[i];  
9.     if (k==v)  return;          //找到了起点则返回  
10.     Ppath(path,k,v);            //找顶点k的前一个顶点  
11.     printf("%d,",k);            //输出顶点k  
12. }  
13. void Dispath(int dist[],int path[],int s[],int n,int v)  
14. {  
15.     int i;  
16.     for (i=0; i<n; i++)  
17.         if (s[i]==1)  
18.         {  
19.             printf("  从%d到%d的最短路径长度为:%d\t路径为:",v,i,dist[i]);  
20.             printf("%d,",v);    //输出路径上的起点  
21.             Ppath(path,i,v);    //输出路径上的中间点  
22.             printf("%d\n",i);   //输出路径上的终点  
23.         }  
24.         else  printf("从%d到%d不存在路径\n",v,i);  
25. }  
26. void Dijkstra(MGraph g,int v)  
27. {  
28.     int dist[MAXV],path[MAXV];  
29.     int s[MAXV];  
30.     int mindis,i,j,u;  
31.     for (i=0; i<g.n; i++)  
32.     {  
33.         dist[i]=g.edges[v][i];      //距离初始化  
34.         s[i]=0;                     //s[]置空  
35.         if (g.edges[v][i]<INF)      //路径初始化  
36.             path[i]=v;  
37.         else  
38.             path[i]=-1;  
39.     }  
40.     s[v]=1;  
41.     path[v]=0;              //源点编号v放入s中  
42.     for (i=0; i<g.n; i++)               //循环直到所有顶点的最短路径都求出  
43.     {  
44.         mindis=INF;                 //mindis置最小长度初值  
45.         for (j=0; j<g.n; j++)       //选取不在s中且具有最小距离的顶点u  
46.             if (s[j]==0 && dist[j]<mindis)  
47.             {  
48.                 u=j;  
49.                 mindis=dist[j];  
50.             }  
51.         s[u]=1;                     //顶点u加入s中  
52.         for (j=0; j<g.n; j++)       //修改不在s中的顶点的距离  
53.             if (s[j]==0)  
54.                 if (g.edges[u][j]<INF && dist[u]+g.edges[u][j]<dist[j])  
55.                 {  
56.                     dist[j]=dist[u]+g.edges[u][j];  
57.                     path[j]=u;  
58.                 }  
59.     }  
60.     Dispath(dist,path,s,g.n,v);     //输出最短路径  
61. }  
62.   
63. int main()  
64. {  
65.     MGraph g;  
66.     int A[7][7]=  
67.     {  
68.         {0,4,6,6,INF,INF,INF},  
69.         {INF,0,1,INF,7,INF,INF},  
70.         {INF,INF,0,INF,6,4,INF},  
71.         {INF,INF,2,0,INF,5,INF},  
72.         {INF,INF,INF,INF,0,INF,6},  
73.         {INF,INF,INF,INF,1,0,8},  
74.         {INF,INF,INF,INF,INF,INF,0}  
75.     };  
76.     ArrayToMat(A[0], 7, g);  
77.     Dijkstra(g,0);  
78.     return 0;  
79. }