第十三周项目二（Kruskal算法的验证）

1. *Copyright(c)2016,烟台大学计算机与控制工程学院    
2. *All right reserved.    
3. *文件名称：利用遍历思想求解图问题.cpp    
4. *作者：宋雨静 
5. *时间：11月25日    
6. *版本号；v1.0    
7. *问题描述：    
8.             
9.    Prim算法的验证和Kruskal算法的验证    
10. *输入描述：带权图的邻接矩阵    
11. *程序输出：最小生成树各边以及权值。    
12. /    
13.   
14. [cpp] view plain copy  
15. #include <stdio.h>    
16. #include <malloc.h>    
17. #define MAXV 100                //最大顶点个数    
18. #define INF 32767       //INF表示∞    
19. #define MaxSize 100    
20. typedef int InfoType;    
21.     
22. //以下定义邻接矩阵类型    
23. typedef struct    
24. {    
25.     int no;                     //顶点编号    
26.     InfoType info;              //顶点其他信息，在此存放带权图权值    
27. } VertexType;                   //顶点类型    
28.     
29. typedef struct                  //图的定义    
30. {    
31.     int edges[MAXV][MAXV];      //邻接矩阵    
32.     int n,e;                    //顶点数，弧数    
33.     VertexType vexs[MAXV];      //存放顶点信息    
34. } MGraph;                       //图的邻接矩阵类型    
35.     
36. //以下定义邻接表类型    
37. typedef struct ANode            //弧的结点结构类型    
38. {    
39.     int adjvex;                 //该弧的终点位置    
40.     struct ANode *nextarc;      //指向下一条弧的指针    
41.     InfoType info;              //该弧的相关信息,这里用于存放权值    
42. } ArcNode;    
43.     
44. typedef int Vertex;    
45.     
46. typedef struct Vnode            //邻接表头结点的类型    
47. {    
48.     Vertex data;                //顶点信息    
49.     int count;                  //存放顶点入度,只在拓扑排序中用    
50.     ArcNode *firstarc;          //指向第一条弧    
51. } VNode;    
52.     
53. typedef VNode AdjList[MAXV];    //AdjList是邻接表类型    
54.     
55. typedef struct    
56. {    
57.     AdjList adjlist;            //邻接表    
58.     int n,e;                    //图中顶点数n和边数e    
59. } ALGraph;                      //图的邻接表类型    
60. typedef struct    
61. {    
62.     int u;     //边的起始顶点    
63.     int v;     //边的终止顶点    
64.     int w;     //边的权值    
65. } Edge;    
66. void ArrayToMat(int *Arr, int n, MGraph &g);    
67. void Prim(MGraph g,int v);    
68. void InsertSort(Edge E[],int n);    
69. void Kruskal(MGraph g);    
70. void ArrayToMat(int *Arr, int n, MGraph &g)    
71. {    
72.     int i,j,count=0;  //count用于统计边数，即矩阵中非0元素个数    
73.     g.n=n;    
74.     for (i=0; i<g.n; i++)    
75.         for (j=0; j<g.n; j++)    
76.         {    
77.             g.edges[i][j]=Arr[i*n+j]; //将Arr看作n×n的二维数组，Arr[i*n+j]即是Arr[i][j]，计算存储位置的功夫在此应用    
78.             if(g.edges[i][j]!=0)    
79.                 count++;    
80.         }    
81.     g.e=count;    
82. }    
83.     
84.     
85.     
86. void InsertSort(Edge E[],int n) //对E[0..n-1]按递增有序进行直接插入排序    
87. {    
88.     int i,j;    
89.     Edge temp;    
90.     for (i=1; i<n; i++)    
91.     {    
92.         temp=E[i];    
93.         j=i-1;              //从右向左在有序区E[0..i-1]中找E[i]的插入位置    
94.         while (j>=0 && temp.w<E[j].w)    
95.         {    
96.             E[j+1]=E[j];    //将关键字大于E[i].w的记录后移    
97.             j--;    
98.         }    
99.         E[j+1]=temp;        //在j+1处插入E[i]    
100.     }    
101. }    
102.     
103. void Kruskal(MGraph g)    
104. {    
105.     int i,j,u1,v1,sn1,sn2,k;    
106.     int vset[MAXV];    
107.     Edge E[MaxSize];    //存放所有边    
108.     k=0;                //E数组的下标从0开始计    
109.     for (i=0; i<g.n; i++)   //由g产生的边集E    
110.         for (j=0; j<g.n; j++)    
111.             if (g.edges[i][j]!=0 && g.edges[i][j]!=INF)    
112.             {    
113.                 E[k].u=i;    
114.                 E[k].v=j;    
115.                 E[k].w=g.edges[i][j];    
116.                 k++;    
117.             }    
118.     InsertSort(E,g.e);      //采用直接插入排序对E数组按权值递增排序    
119.     for (i=0; i<g.n; i++)   //初始化辅助数组    
120.         vset[i]=i;    
121.     k=1;    //k表示当前构造生成树的第几条边,初值为1    
122.     j=0;    //E中边的下标,初值为0    
123.     while (k<g.n)       //生成的边数小于n时循环    
124.     {    
125.         u1=E[j].u;    
126.         v1=E[j].v;      //取一条边的头尾顶点    
127.         sn1=vset[u1];    
128.         sn2=vset[v1];   //分别得到两个顶点所属的集合编号    
129.         if (sn1!=sn2)   //两顶点属于不同的集合    
130.         {    
131.             printf("  (%d,%d):%d\n",u1,v1,E[j].w);    
132.             k++;                     //生成边数增1    
133.             for (i=0; i<g.n; i++)   //两个集合统一编号    
134.                 if (vset[i]==sn2)   //集合编号为sn2的改为sn1    
135.                     vset[i]=sn1;    
136.         }    
137.         j++;               //扫描下一条边    
138.     }    
139. }    
140.     
141. int main()    
142. {    
143.     MGraph g;    
144.     int A[6][6]=    
145.     {    
146.         {0,6,1,5,INF,INF},    
147.         {6,0,5,INF,3,INF},    
148.         {1,5,0,5,6,4},    
149.         {5,INF,5,0,INF,2},    
150.         {INF,3,6,INF,0,6},    
151.         {INF,INF,4,2,6,0}    
152.     };    
153.     ArrayToMat(A[0], 6, g);    
154.     printf("最小生成树构成:\n");    
155.     Kruskal(g);    
156.     return 0;    
157. }    

运行结果：