Kruskal 最小生成树的总结

写了3，4天的最小生成树，里边写了报告和没写的加起来也有十题左右了。发现这些简单的算法里头有不少的pin。特别是2421让人抓狂，bug一个接一个，wa了不下十数次后，在波波的指导下成功的ac了，后来自己又找到了原来程序的几个bug。

原先写这篇报告的时候以为Kruskal比Prim优秀许多。后来波波跟我说Prim也是很强的,于是我分别用Kruskal和Prim写了Poj1789,发现原来Prim也是很strong的,用的时间几乎只有Kruskal的20%,囧。于是现在把这篇报告改为Kruskal最小生成树总结,顺便把标程贴出来以后好找~~

 

下面的是Kruskal：

1. //这段代码用Kruskal来解决
3. #include <stdio.h>
4. #include <stdlib.h>
5. #include <memory.h>
6. #include <algorithm>
7. using namespace std;
8. // 边类型
9. typedef struct Edge{
10.     int x;
11.     int y;
12.     int weigh;
13.     bool operator() ( Edge a,Edge b) {
14.         return (a.weigh < b.weigh);
15.     }
16. }Edge;
17. //最多可能的边
18. Edge edge[251001];
19. int parent[501];
20. //用于qsort()的比较函数
21. int compare(const void*a,const void *b)
22. {
23.     return ((Edge*)a)->weigh - ((Edge*)b)->weigh;
24. }
25. // 找父节点并压缩路径
26. int Find(int x)
27. {
28.     if( parent[x] >= 0)
29.         parent[x] = Find(parent[x]);
30.     else
31.         return x;
32.     return parent[x];
33. }
34. // 合并集合
35. void Union(int x,int y)
36. {
37.     x = Find(x);
38.     y = Find(y);
39.     // 注意:这里的重边判断很重要
40.     if( x == y)
41.         return ;
42.     if(parent[x] < parent[y]) 
43.         parent[y] = x;
44.     else
45.         if(parent[x] > parent[y])
46.             parent[x] = y;
47.         else {
48.             parent[x]--;
49.             parent[y] = x;
50.         }
51. }
52. int main()
53. {
54.     int m,n;
55.     int count = 0;
56.     scanf("%d",&n);
57.     for( int i = 0; i < n; i++) {
58.         scanf("%d",&m);
59.         memset(parent,-1,m * sizeof(int));
60.         for( int j = 0; j < m; j++)
61.             for( int k = 0; k < m; k++) {
62.                 scanf("%d",&edge[count].weigh);
63.                 edge[count].x = j;
64.                 edge[count++].y = k;
65.             }
66.         sort(edge,edge + m * m,Edge());
67.         count = 1;
68.         int max;
69.         // 忽略那些为0的节点,有m个,然后从小到大去边
70.         for( int j = m; j < (m * m) && count < m; j++)
71.             if(Find(edge[j].x) != Find(edge[j].y)) {
72.                 count ++;
73.                 Union(edge[j].x,edge[j].y);
74.                 max = edge[j].weigh;
75.             }
76.         printf("%d/n",max);
77.     }
78. }

唉,还是有parent好啊~~