Dijkstra算法.doc

Dijkstra算法Dijkstra算法是典型最短路算法，用于計(jì)算一個(gè)節(jié)點(diǎn)到其他所有節(jié)點(diǎn)的最短路徑。主要特點(diǎn)是以起始點(diǎn)為中心向外層層擴(kuò)展，直到擴(kuò)展到終點(diǎn)為止。Dijkstra算法能得出最短路徑的最優(yōu)解，但由于它遍歷計(jì)算的節(jié)點(diǎn)很多，所以效率低。Dijkstra算法是很有代表性的最短路算法，在很多專業(yè)課程中都作為基本內(nèi)容有詳細(xì)的介紹，如數(shù)據(jù)結(jié)構(gòu)，圖論，運(yùn)籌學(xué)等等。Dijkstra一般的表述通常有兩種方式，一種用永久和臨時(shí)標(biāo)號(hào)方式，一種是用OPEN,CLOSE表方式，Drew為了和下面要介紹的A*算法和D*算法表述一致，這里均采用OPEN,CLOSE表的方式。其采用的是貪心法的算法策略大概過程：創(chuàng)建兩個(gè)表，OPEN,CLOSE。OPEN表保存所有已生成而未考察的節(jié)點(diǎn)，CLOSED表中記錄已訪問過的節(jié)點(diǎn)。1訪問路網(wǎng)中距離起始點(diǎn)最近且沒有被檢查過的點(diǎn)，把這個(gè)點(diǎn)放入OPEN組中等待檢查。2從OPEN表中找出距起始點(diǎn)最近的點(diǎn)，找出這個(gè)點(diǎn)的所有子節(jié)點(diǎn)，把這個(gè)點(diǎn)放到CLOSE表中。3遍歷考察這個(gè)點(diǎn)的子節(jié)點(diǎn)。求出這些子節(jié)點(diǎn)距起始點(diǎn)的距離值，放子節(jié)點(diǎn)到OPEN表中。4重復(fù)第2和第3步,直到OPEN表為空，或找到目標(biāo)點(diǎn)。編輯本段算法實(shí)現(xiàn)#include#defineMaxNum765432100usingnamespacestd；ifstreamfin(Dijkstra.in)；ofstreamfout(Dijkstra.out)；intMap501501；boolis_arrived501；intDist501,From501,Stack501；intp,q,k,Path,Source,Vertex,Temp,SetCard；intFindMin()intp,Temp=0,Minm=MaxNum；for(p=1；p=Vertex；p+)if(DistpSourceVertex；for(p=1；p=Vertex；p+)for(q=1；qMappq；if(Mappq=0)Mappq=MaxNum；for(p=1；p=Vertex；p+)Distp=MapSourcep；if(Distp!=MaxNum)Fromp=Source；elseFromp=p；is_arrivedSource=true；SetCard=1；doTemp=FindMin()；if(Temp!=0)SetCard=SetCard+1；is_arrivedTemp=true；for(p=1；pDistTemp+MapTempp)&(!is_arrivedp)Distp=DistTemp+MapTempp；Fromp=Temp；elsebreak；while(SetCard!=Vertex)；for(p=1；p=Vertex；p+)if(p!=Source)fout=n；foutSource:SourcenTarget:pn；if(Distp=MaxNum)foutDistance:Infinityn；foutPath:NoWay!；elsefoutDistance:Distpn；k=1；Path=p；while(FromPath!=Path)Stackk=Path；Path=FromPath；k=k+1；foutPath:=1；q-)foutStackq；fout1=Source:2Target:3Distance:25Path:2-3=Source:2Target:4Distance:50Path:2-1-4=Source:2Target:5Distance:50Path:2-3-5=Source:2Target:6Distance:60Path:2-3-5-6=Source:2Target:7Distance:InfinityPath:NoWay!=示例程序及相關(guān)子程序：voidDijkstra(intn,intDistance,intiPath)intMinDis,u；inti,j；/從鄰接矩陣復(fù)制第n個(gè)頂點(diǎn)可以走出的路線，就是復(fù)制第n行到Distancefor(i=0；iVerNum；i+)Distance=Arcn,i；Visited=0；/第n個(gè)頂點(diǎn)被訪問，因?yàn)榈趎個(gè)頂點(diǎn)是開始點(diǎn)Visitedn=1；/找到該頂點(diǎn)能到其他頂點(diǎn)的路線、并且不是開始的頂點(diǎn)n、以前也沒走過。/相當(dāng)于尋找u點(diǎn)，這個(gè)點(diǎn)不是開始點(diǎn)nfor(i=0；iVerNum；i+)u=0；MinDis=No；for(j=0；jVerNum；j+)if(Visitedj=0&(DistancejMinDis)MinDis=Distancej；u=j；/如范例P1871圖G6，Distance=No,No,10,No,30,100，第一次找就是V2，所以u(píng)=2/找完了，MinDis等于不連接，則返回。這種情況類似V5。if(MinDis=No)return；/確立第u個(gè)頂點(diǎn)將被使用，相當(dāng)于Arcv,u+Arcu,w中的第u頂點(diǎn)。Visitedu=1；/尋找第u個(gè)頂點(diǎn)到其他所有頂點(diǎn)的最小路，實(shí)際就是找Arcu,j、j取值在0，VerNum。/如果有Arci,u+Arcu,jArci,j，則Arci,j=Arci,u+Arcu,jArci,j/實(shí)際中，因?yàn)镈istance是要的結(jié)果，對(duì)于起始點(diǎn)確定的情況下，就是：/如果(Distanceu+Arcu,j)=Distancej則：/Distancej=Distanceu+Arcu,j；/而iPath保存了u點(diǎn)的編號(hào)；/同理：對(duì)新找出的路線，要設(shè)置Visitedj=0，以后再找其他路，這個(gè)路可能別利用到。例如V3for(j=0；jVerNum；j+)if(Visitedj=0&Arcu,jNo&u!=j)if(Distanceu+Arcu,j)=Distancej)Distancej=Distanceu+Arcu,j；Visitedj=0；iPathj=u；/輔助函數(shù)voidPrim()inti,m,n=0；for(i=0；i0)if(m=MinAdjNode(n)!=-1)Tn.Nodes.Add(Tm)；n=m；Visitedn+；elsen=MinNode(0)；if(n0)TMin2.Nodes.Add(TMin1)；Visitedn+；listBox1.Items.Add(Vn)；treeView1.Nodes.Add(T0)；voidTopoSort()inti,n；listBox1.Items.Clear()；StackS=newStack()；for(i=0；i=0；i-)if(InDegree(i)=0)S.Push(i)；Visited+；while(S.Count!=0)n=(int)S.Pop()；listBox1.Items.Add(Vn)；ClearLink(n)；for(i=VerNum-1；i=0；i-)if(Visited=0&InDegree(i)=0)S.Push(i)；Visited+；voidAOETrave(intn,TreeNodeTR,intw)inti,w0；if(OutDegree(n)=0)return；for(i=0；iVerNum；i+)if(w0=Arcn,i)!=0)listBox1.Items.Add(V+t+(w+w0).ToString()+t+i.ToString()+t+n.ToString()；TreeNodeT1=newTreeNode()；T1.Text=V+W=+(w+w0).ToString()+；TR.Nodes.Add(T1)；AOETrave(i,T1,w+w0)；voidAOE()inti,w=0,m=1；TreeNodeT1=newTreeNode()；for(i=0；iVerNum；i+)Visited=0；T1.Text=V0；listBox1.Items.Add(雙親表示法顯示這個(gè)生成樹：)；listBox1.Items.Add(VtWtIDtPID)；for(i=0；iVerNum；i+)if(w=Arc0,i)!=0)listBox1.Items.Add(V+t+w.ToString()+t+i.ToString()+t0)；TreeNodeT2=newTreeNode()；T2.Text=V+W=+w.ToString()+；AOETrave(i,T2,w)；T1.Nodes.Add(T2)；listBox1.Items.Add(tt樹+m.ToString()；m+；treeView1.Nodes.Clear()；treeView1.Nodes.Add(T1)；intIsZero()inti；for(i=0；i=0)returni；return-1；in