1. /*
    2. prim算法在用邻接矩阵方式存储的图中的实现
    3. */
    4. #define _CRT_SECURE_NO_WARNINGS
    5. #define MAXSIZE 100 //数组最大值
    6. #include <stdio.h>
    7. #include <stdlib.h>
    8. typedef struct Graph {
    9. char Vertex[MAXSIZE];
    10. int Edge[MAXSIZE][MAXSIZE];
    11. int numV, numE;//顶点、边数量
    12. }adjMatrix;
    13. int getSum(adjMatrix *G, int *prims) {//获得最小生成树的权值
    14. int sum = 0;
    15. for (int i = 1; i < G->numV; i++) {
    16. int min = 32767;
    17. for (int j = 0; j < i; j++) {
    18. if (G->Edge[prims[j]][prims[i]] < min) {
    19. min = G->Edge[prims[j]][prims[i]];//疯了,把prims[j]写成了j,害我调试了好久!!!
    20. }
    21. }
    22. sum += min;
    23. }
    24. return sum;
    25. }
    26. void prim(adjMatrix *G, int start) {
    27. int prims[MAXSIZE];//存储最小生成树结果数组
    28. int lowCost[MAXSIZE];//当前已加入最小生成树的邻接顶点的边的权值
    29. int min, k, index = 0;
    30. lowCost[start] = 0;//自己到自己的距离为0
    31. prims[index++] = start;
    32. for (int i = 0; i < G->numV; i++) {
    33. lowCost[i] = G->Edge[start][i];//把当前传入顶点的所有连接顶点的边的权值存入
    34. }
    35. for (int i = 0; i < G->numV; i++) {//进行遍历,遍历一次加入一个顶点
    36. if (start == i)
    37. continue;
    38. min = 32767;
    39. for (int j = 0; j < G->numV; j++) {
    40. if (lowCost[j] != 0 && lowCost[j] < min) {//如果当前顶点未曾加入最小树中且小于目前最小值,更新
    41. min = lowCost[j];
    42. k = j;//记录位置
    43. }
    44. }
    45. //已找到最小值,更新prims数组
    46. prims[index++] = k;
    47. lowCost[k] = 0;//将第k个顶点置为已访问,即代表它已加入最小生成树
    48. //如有顶点未处理,则看情况需更新lowCost数组
    49. for (int j = 0; j < G->numV; j++) {
    50. if (lowCost[j] && G->Edge[k][j] < lowCost[j]) {
    51. lowCost[j] = G->Edge[k][j];
    52. /*
    53. 对于这里的理解,为什么要替换?然后为什么可以这样替换?
    54. 首先回答第一个:因为普里姆算法是每次根据我们加入的点去寻找目前可纳入的最小边,
    55. 所以我们没加入一个顶点都要进行判断,加入的顶点是否让某些顶点的
    56. 可达距离变短
    57. 第二个问题: 我们每次确定一个顶点,也就意味着一条边也确定了下来,所以此时新加入顶点
    58. 到其余顶点的距离若小于之前的距离,则需要进行更换
    59. 比如这里的 1 2 3,当1 3 确定,1 3 边的权值为1,我们的目标是让最后的权值
    60. 最小,我们知道之前1到2权值为6,总权值就为1 + 6 = 7,而此时3 到 2仅为5
    61. 小于 6,所以要进行更新替换。
    62. */
    63. }
    64. }
    65. }
    66. printf("%d ", getSum(G, prims));
    67. for (int i = 0; i < G->numV; i++) {
    68. printf("%c ", G->Vertex[prims[i]]);
    69. }
    70. }
    71. int main() {
    72. void createGraphFromFile(adjMatrix *);
    73. void dispGraph(adjMatrix *G);
    74. adjMatrix *G = (adjMatrix *)malloc(sizeof(adjMatrix));
    75. createGraphFromFile(G);
    76. dispGraph(G);
    77. prim(G, 0);
    78. return 0;
    79. }