最小生成树问题---Prim算法与Kruskal算法实现(MATLAB语言实现)

  2015-12-17晚,复习,甚是无聊,阅《复杂网络算法与应用》一书,得知最小生成树问题(Minimum spanning tree)问题。记之。

  何为:连通且不含圈的图称为树。

  图T=(V,E),|V|=n,|E|=m,下列关于树的说法等价:

  1. T是一个树。
  2. T无圈,且m=n-1。
  3. T连通,且m=n-1。
  4. T无圈,但每加一新边记得到唯一一个圈。
  5. T连通,但任舍去一边就不连通。
  6. T中任意两点,有唯一道路相连。

  何为生成树:若图G=(V,E)的生成子图是一棵树,则称该树为图G的生成树,也称支撑树,简称为图G的数。图G中属于生成树的边称为数枝(Branch)。

  何为最小生成树:连通图G=(V,E),每条边上有非负权L(e)。一棵树上所有树枝权的总和,称为这个生成树的权。具有最小权的生成树称为最小生成树,也就是说最小支撑树,简称最小树。

  私以为,两种算法其实都是贪心,所以需要严格的证明。由于最近时间零散、数学久置未学、对算法领域没有系统了解。所以不进行深入探讨(也就是说证明),仅以一个简单实例做一个入门级的了解。

  Prim算法:

  给定连通赋权图G=(V,E,W),其中W为邻接矩阵,设两个集合P和Q,其中P用于存放G的最小生成树中的节点,集合Q存放G的最小G的最小生成树中的边。另集合P的初值为P={v1}(假设构造最小生成树时从v1出发),集合Q的初值为P={空集}。

  (1)P = {v1},Q = {空集};

  (2)while P ~= Q

    找到最小边pv,其中p∈P,v∈V-P;

    P = P + {v};

    Q = Q + {pv};

    end

  Kruskal算法

  (1)选e1∈E(G),使得w(e1) = min(选e1的权值最小)。

  (2)e1,e2,...,ei已选好,则从E(G)-{e1,e2,...,ei}中选取ei+1,使得G[{e1,e2,...,ei,ei+1}]中无圈,且,w(ei+1) = min。

  (3)直到选得en-1为止。

  以下是问题:

一个乡有9个自然村,其间道路及各道路长度如图所示,各边上数字代表距离。问如何架设电线最短。

  最小生成树问题---Prim算法与Kruskal算法实现(MATLAB语言实现)

预输入:

 A = zeros();
A(,:) = [ ];
A(,[, ]) = [ ];
A(, ) = ;
A(, ) = ;
A(, ) = ;
A(, ) = ;
A(, ) = ;
A(, ) = ;

Prim算法实现:

 A = A+A';
A(A==) = Inf;
P = zeros(, );
P(,) = ;
V = :;
V_P = V - P;
link = zeros(,);
k=;
while k<
p = P(P~=);
v = V_P(V_P~=);
pv = min(min(A(p,v)));
[x, y] = find(A==pv);
for i=:length(x)
if any(P==x(i)) && any(V_P==y(i))集合判断,关键!
P(,y(i)) = y(i);
V_P = V - P;
link(k, :) = [x(i), y(i)];
k = k+;
break;
end
end
end

输出:

>> link
link =

Kruskal算法实现:

 A = sparse(A');
A(A==) = Inf; B = sparse(, );
link = graphallshortestpaths(B, 'directed', false);
while sum(sum(link)) == Inf%如果不连通,则和无穷大
ind = find(A==min(min(A)));
[x, y] = ind2sub(size(A), ind);%寻找最短边
for i=:length(x)
if link(x(i), y(i))==Inf%判断是否连通,关键!!
B(x(i), y(i)) = A(x(i), y(i));
A(x(i), y(i))=Inf;%取差集
end
end
link = graphallshortestpaths(B, 'directed', false);
end

输出:

B =
(,)
(,)
(,)
(,)
(,)
(,)
(,)
(,)

所以结果是:

  最小生成树问题---Prim算法与Kruskal算法实现(MATLAB语言实现)

(算法实现写于2015-12-18)

  下一步工作,用其他语言实现,复杂度分析,深入论证算法正确性

  另外,代码中使用了graphallshortestpaths函数,这是最短路径的问题。关于最短路径将在下一篇讲述。其实MATLAB有关于求求最小生成树的函数graphminspantree,关于这一套函数,将结合官网资料,以后做梳理。

  A = sparse(A');
>> graphminspantree(A)
ans =
(,)
(,)
(,)
(,)
(,)
(,)
(,)
(,)

  如此,可直接求得最小生成树。

上一篇:APNS .p12文件转换为 .pem文件


下一篇:最小生成树问题---Prim算法学习