第十三周 项目2:Kruskal算法的验证

来源:互联网 发布:韩国漫画下载软件 编辑:程序博客网 时间:2024/06/05 04:41
  1. * Copyright (c)2016,烟台大学计算机与控制工程学院      
  2. * All rights reserved.      
  3. * 文件名称:jiang.cpp      
  4. * 作    者:姜孝龙     
  5. * 完成日期:2016年12月8日      
  6. * 版 本 号:v1.0       
  7. *问题描述:验证Kruskal算法 
  8. *输入描述:无      
  9. *程序输出:测试数据      
  10. */
程序中graph.h是图存储结构的“算法库”中的头文件

  1. #include <stdio.h>  
  2. #include <malloc.h>  
  3. #include "graph.h"  
  4. #define MaxSize 100  
  5. typedef struct  
  6. {  
  7.     int u;     //边的起始顶点  
  8.     int v;     //边的终止顶点  
  9.     int w;     //边的权值  
  10. } Edge;  
  11.   
  12. void InsertSort(Edge E[],int n) //对E[0..n-1]按递增有序进行直接插入排序  
  13. {  
  14.     int i,j;  
  15.     Edge temp;  
  16.     for (i=1; i<n; i++)  
  17.     {  
  18.         temp=E[i];  
  19.         j=i-1;              //从右向左在有序区E[0..i-1]中找E[i]的插入位置  
  20.         while (j>=0 && temp.w<E[j].w)  
  21.         {  
  22.             E[j+1]=E[j];    //将关键字大于E[i].w的记录后移  
  23.             j--;  
  24.         }  
  25.         E[j+1]=temp;        //在j+1处插入E[i]  
  26.     }  
  27. }  
  28.   
  29. void Kruskal(MGraph g)  
  30. {  
  31.     int i,j,u1,v1,sn1,sn2,k;  
  32.     int vset[MAXV];  
  33.     Edge E[MaxSize];    //存放所有边  
  34.     k=0;                //E数组的下标从0开始计  
  35.     for (i=0; i<g.n; i++)   //由g产生的边集E  
  36.         for (j=0; j<g.n; j++)  
  37.             if (g.edges[i][j]!=0 && g.edges[i][j]!=INF)  
  38.             {  
  39.                 E[k].u=i;  
  40.                 E[k].v=j;  
  41.                 E[k].w=g.edges[i][j];  
  42.                 k++;  
  43.             }  
  44.     InsertSort(E,g.e);      //采用直接插入排序对E数组按权值递增排序  
  45.     for (i=0; i<g.n; i++)   //初始化辅助数组  
  46.         vset[i]=i;  
  47.     k=1;    //k表示当前构造生成树的第几条边,初值为1  
  48.     j=0;    //E中边的下标,初值为0  
  49.     while (k<g.n)       //生成的边数小于n时循环  
  50.     {  
  51.         u1=E[j].u;  
  52.         v1=E[j].v;      //取一条边的头尾顶点  
  53.         sn1=vset[u1];  
  54.         sn2=vset[v1];   //分别得到两个顶点所属的集合编号  
  55.         if (sn1!=sn2)   //两顶点属于不同的集合  
  56.         {  
  57.             printf("  (%d,%d):%d\n",u1,v1,E[j].w);  
  58.             k++;                     //生成边数增1  
  59.             for (i=0; i<g.n; i++)   //两个集合统一编号  
  60.                 if (vset[i]==sn2)   //集合编号为sn2的改为sn1  
  61.                     vset[i]=sn1;  
  62.         }  
  63.         j++;               //扫描下一条边  
  64.     }  
  65. }  
  66.   
  67. int main()  
  68. {  
  69.     MGraph g;  
  70.     int A[6][6]=  
  71.     {  
  72.         {0,6,1,5,INF,INF},  
  73.         {6,0,5,INF,3,INF},  
  74.         {1,5,0,5,6,4},  
  75.         {5,INF,5,0,INF,2},  
  76.         {INF,3,6,INF,0,6},  
  77.         {INF,INF,4,2,6,0}  
  78.     };  
  79.     ArrayToMat(A[0], 6, g);  
  80.     printf("最小生成树构成:\n");  
  81.     Kruskal(g);  
  82.     return 0;  
  83. }  


运行结果:

0 0
原创粉丝点击