求图的连通,8.docVIP

“离散数学”实验报告 （求图的连通性） 专 业 网络工程 班 级 1202班 学 号 姓 名 张敏慧 2013.12.14 目录 TOC \o 1-3 \h \z \u HYPERLINK \l _Toc276474920 一.实验目的 3 HYPERLINK \l _Toc276474921 二.实验原理 3 HYPERLINK \l _Toc276474922 1. Warshell算法 3 HYPERLINK \l _Toc276474923 2. 矩阵幂算法 3 HYPERLINK \l _Toc276474924 三.实验环境 3 HYPERLINK \l _Toc276474925 四. 实验流程图 3 HYPERLINK \l _Toc276474926 1.实验内容 3 HYPERLINK \l _Toc276474927 2.实验流程图 4 HYPERLINK \l _Toc276474928 五.实验代码 5 HYPERLINK \l _Toc276474929 1.MATLAB语言. 5 HYPERLINK \l _Toc276474930 2.C语言 7 HYPERLINK \l _Toc276474931 六. 实验结果 8 HYPERLINK \l _Toc276474934 七. 实验总结 9 一、实验目的 用计算机语言编写图的连通性判断算法，可输入图的邻接矩阵，判断图是否连通以及确定连通分支的个数，掌握Warshell算法或矩阵幂算法的实现方法。 二、实验原理 1、Warshell算法 Warshell算法可解决图是否连通的问题, 而且效率很高。在该算法中，矩阵是判断矩阵，表示从到连通，表示从到不连通。 (1)置新矩阵 P:= C; (2)置 = 1; (3)对所有的,若, 则对k=1,2,…,n, 有; (4) ; (5) 如转向步骤(3), 否则停止。 2、矩阵幂算法 由于邻接阵包含了图的所有信息，和关联阵一样，是图的等价表示。可以通过对邻接阵C做一些计算，得到图G的一些性质。例如考虑中的的元素，如果它不为零，由于，则至少存在一组或一个长度为3的链使端和端相连。从而，通过计算C的各阶幂次可得到关于图是否连通的信息。 三、实验环境 使用visual C++6.0为编程软件，采用C语言为编程语言实现。 使用MATLAB等语言实现图的连通性判断算法。 四、实验流程图 实验内容： 利用MATLAB等语言实现图的连通性判断算法，可对输入的邻接阵进行连通性以及连通分支数的判断。 实验流程图 矩阵幂算法： YesNoC1=C^kP=P+C1 Yes No C1=C^k P=P+C1 k=k+1 得到C的阶数n k=1 开始 输入邻接矩阵C k=n? 连通分支数S=n-连通矩阵P的秩 结束 Warshell算法： i++若p(j,i)=1 i++ 若p(j,i)=1 p(i,j)=p(i,k)+*p(k,j) j++ i=1,j=1 输入邻接矩阵C 得到C的阶数n 边数m k=1 开始 i=n? j=i? Yes No Yes No 结束 连通分支数S 五、实验代码 MatLab 源代码： clear,clc; %输入邻接矩阵 disp(图的连通性以及连通分支数的判断); C = input(请输入图的邻接矩阵(格式如：[1 1 0;1 1 1;0 1 1]) C=); %矩阵幂算法 n=size(C,1);%邻接矩阵阶数 P=zeros(n,n);%构造连通矩阵P k=1; for k=1:n %计算矩阵幂的和 C1=C^k; P = P + C1; end S=n-rank(P);%连通分支数为0特征值个数 %Warshell算法 S1=0;a=1; G=zeros(n,1); for i=1:n for j=(i+1):n if C(i,j)==1%若两端之间有边连通 if G(i)==G(j)%若两端之间有连通链，说明二者在同一连通分支 if G(i)==0 G(i)=a;G(j)=a; a=a+1; S1=S1+1; end else if G(i)==0 G(i)=G(j);%若与i不连通，则与j在同一连通分支 elseif G(j)==0