(无向图)邮路.doc

（无向图）中国邮路 一、相关概念 二、实际意义 一笔画问题：求欧拉回路。 邮路问题：经过每条边至少一次，且重复的边的总长最小。 三、算法及步骤 附录 参考资料 [1]甘应爱,田丰等. 运筹学（第三版）[M]. 北京: 清华大学出版社, 2005年6月第三版: 276-280 [2]戴一奇,胡冠章,陈卫. 图论与代数结构[M]. 北京: 清华大学出版社, 1995年8月第1版: 32-33 function [P,LEN]=PostRoad(C) %========================================================================== %相关概念： % 1、连通图G=（V,E）中的一条经过所有边的简单回路（道路）称为G的欧拉回路（道 %路）。 % 2、无向连通图G存在欧拉回路的充要条件是G中各点的度都是偶数。 % 3、若无向连通图G只有两个度为奇数的结点，则G存在欧拉道路。 % 4、若有向连通图G中各结点的正度都等于该结点的负度，则G存在欧拉回路。 % 5、L是无向连通图G的最佳邮路的充要条件是： % （1）G的每条边最多重复一次； % （2）在G的任意一个回路上，重复边的长度之和不超过该回路长度的一半。 %========================================================================== %实际意义： % 一笔画问题：求欧拉回路。 % 邮路问题：经过每条边至少一次，且重复的边的总长最小。 %========================================================================== %算法及步骤： % Edmonds最小权匹配算法。 % step1：求出图G中度为奇数的结点，构成V0。 % step2：求V0中任意两点在G中的最短路径P(i,j)及其长度len(i,j)。 % step3：对V0中的结点进行最小权匹配，即选出|V0|/2个len(i,j)，保证V0中的每个结 %点只出现1次且总和最小。 % step4：在最小权匹配中选中的各个路径的边在G中重复一次。 %========================================================================== %输入数据：正权无向图G。 %输出结果：需要重复的路径及其总长度。 %========================================================================== [x,y]=size(C); N=max(x,y); V0=[]; %存储度为奇数的结点编号 for i=1:N C(i,i)=0; %消除自环 l=length(find(C(i,:)~=0)); %统计结点i的度，即第i行非零元素个数 if floor(l/2)~=l/2 %度为奇数 V0(length(V0)+1)=i; end; end; L=length(V0); %统计度为奇数的结点的个数 len=zeros(L); %存储V0中任意两个结点之间的最短路径长度 road=zeros(L*(L-1)/2,N); %记录V0中任意两个结点之间的最短路径 m=0; for i=1:L-1 len(i,i)=sum(sum(abs(C)))*N*N; %将结点到自身的距离设置为无穷以便于最小权匹配的进行 for j=i+1:L m=m+1; [len(i,j),r]=ShortestPath(C,V0(i),V0(j)); %求最短路径 road(m,1:length(r))=r; len(j,i)=len(i,j); end; end; len(L,L)=sum(sum(abs(C)))*N*N; [LEN,match]=PerfectMatch(len,min); %最小权匹配 m=0; no=0; P=[]; %存储被选出的重复路径 for i=1:L-1 for j=i+1:L m=m+1; if match(i,j)==1 no=no+1; P(no,:)=road(m,:); end; end; end; LEN=LEN/2; 《图论的MAT