Dijkstra的matlab算法.docVIP

//Dijkstra算法 ? 时间复杂度 ? O（NV^2） //D[i]为v到i的最短路径长度 //P[i]为此最短路径中i的前趋 //注意： //1。i，j不相邻时E[i][j]应为infinity //2。不使用AV //3。边权不能为负值 void ? Graph::Dijkstra(int ? v,vector int ? D,vector int ? P)const { vector bool ? M ? = ? VM; for(int ? i=0;i NV;i++) { ? ? ? D[i] ? = ? E[v][i]; ? ? P[i] ? = ? v; ? ? } M[v] ? = ? false; ? D[v] ? = ? 0; for(;;) { int ? min,w; min ? = ? infinity; for(int ? j=0;j NV;j++) { if(!M[j])continue; if(D[j] ? min){ ? ? ? min ? = ? D[j]; ? ? w ? = ? j; ? } } if(min ? == ? infinity)break; ? ? ? M[w] ? = ? false; for(int ? j ? = ? 0;j NV;j++) { if(!M[j])continue; if(min ? + ? E[w][j] D[j]){ ? ? ? D[j] ? = ? min ? + ? E[w][j]; ? ? P[j] ? = ? w; ? ? } } } } //Floyd算法 ? 时间复杂度O（NV*|E| ? + ? NV^2） //D[i][j]为i到j的最短距离 //P[i][j]为求出D[i][j]时j的前趋 //注意： //1。i，j不相邻时E[i][j]应为infinity //2。边权不能为负值 //3。不使用AV void ? Graph::Floyd_Warshall(VVI ? D,VVI ? P)const { int ? i,j,k; for ? (i=0;i NV;i++) for ? (j=0;j NV;j++) { ? ? ? D[i][j]=E[i][j]; P[i][j]=i; ? ? ? ? ? ? ? } for ? (i=0;i NV;i++) ? { ? ? ? D[i][i]=0;P[i][i]=0; ? } for ? (k=0;k NV;k++) { if(!VM[k])continue; for ? (i=0;i NV;i++) { if(!VM[i])continue; if(D[i][k] ? == ? infinity)continue; for ? (j=0;j NV;j++) { if(!VM[j])continue; if ? (D[i][j] D[i][k]+D[k][j]) { ? ? ? D[i][j]=D[i][k]+D[k][j]; ? ? P[i][j]=P[k][j]; ? ? } } } } } 求一个点到其他各点的最短路径 function [d,DD]=dijkstra(D,s)%Dijkstra最短路算法Matlab程序用于求从起始点s到其它各点的最短路%D为赋权邻接矩阵%d为s到其它各点最短路径的长度%DD记载了最短路径生成树[m,n]=size(D);d=inf.*ones(1,m);d(1,s)=0;Vdd=zeros(1,m);dd(1,s)=1;y=s;DD=zeros(m,m);DD(y,y)=1;counter=1;while length(find(dd==1))m? ? for i=1:m? ?? ? if dd(i)==0? ?? ?? ? d(i)=min(d(i),d(y)+D(y,i));? ?? ? end? ? end? ? ddd=inf;? ? for i=1:m? ? if dd(i)==0d(i)ddd? ?? ?ddd=d(i);? ?? ?end? ? end? ? yy=find(d==ddd);? ? counter=counter+1;? ? DD(y,yy(1,1))=counter;? ? DD(yy(1,1),y)=counter;? ? y=yy(1,1);? ? dd(1,y)=1;end 匈牙利算法的matlab实现 function [y,fval,flag]=Hungary(C)%y返回的是最优值的行列号%******