Dijkstra算法详解

Dijkstra算法详解:https://blog.csdn.net/longshengguoji/article/details/10756003

在解决单源点最短路径的问题时,常常用到经典的Dijkstra算法,其算法的本质思想是:按路径长度递增依次产生最短路径

下面给出算法的大致流程:

1.初始化所有结点并将起始点设为标记,进入以下循环

2.在到达某点的最短路径中找最小且未标记的点(可以用一维数组表示)

  如:

数组下标:0  1   2   3   4   5

   Len      :-   0   5  10  2   -

这个数组表示 1号节点为初始节点,1号节点到达2号节点的最短路径为5,到3号为10,无法到达5号(具体可以以较大的数表示其路径)。

从中找到一个未标记且Len最短的一个,未标记用另一数组记录

如:

数组下标:0   1   2   3   4   5

   标记     :0   1   0   0   1   0

此数组表示从初始节点到达4号节点的最短路径已找到

从以上两个数组中可以得出:此次循环找到的点为2号节点,进入下一步

3.标记找到的点,以此标记点为中间点重新计算所有未标记点的最短路径(更新最短路径表)

4.循环1.2步至n-1次(n为顶点数,若最后还有未被标记的,说明无法到达此点)

下面是核心代码:

[cpp] view plain copy

while(count

      {  

          tempmin=INFINITE;  

for(i=1;i<=n;i++)  

          {  

if(in[i]==0&&Len[i]

              {  

                  tempmin=Len[i];  

                  si=i;  

              }  

          }  

          in[si]=1;  

          count++;  

for(i=1;i<=n;i++) //updata the length  

          {  

if(in[i]==0&&(tempmin+mGraph.matrix[si][i])

              {  

                  Len[i]=tempmin+mGraph.matrix[si][i];  

              }  

          }      

      }  

核心部分是上面的两个for循环,第一个for循环遍历所有结点,找到未标记并且长度最短的路径;第二个for循环也是遍历所有结点以上面找到的最短路径结点为中间点更新最短路径表;注意在第一个for循环之后要标记找到的点(说明到达此点的最短路径已找到)

下面用一个实例来具体说明:

若有这样一张有向图(此图为《数据结构》严蔚敏 p188页,书中没有详细讲解,现以此为实例)

V0-V5 共6个节点,节点间路径已标出,现要求从V0到其余各节点的最短路径;

有上面的算法流程可知,在使用Dijkstra算法时需要几个结构来保存我们想要的信息:

1.保存这张图的结构体

2.记录V0到其余各节点最短路径的数组(这里设为Len[n+1])

3.记录某节点是否已找到最短路径的数组(这里设为in[n+1])

接下来就是算法实现部分:

1.标记V0 -- in[1]=1   初始化 Len[]={INFINITE , 0 , INFINITE , 10, INFINITE , 30 , 100}  这里数组首元素未用到,数组下标从1开始表示V0以此类推

2.第一次循环与V0相邻的有V2、V4、V5,其中V2距离最短为10,标记V2,并且以V2为中间点(路径为V0->V2->Vx)更新最短路表,此时V3被更新为10+50=60,。

3.进入第二次循环,此时未被标记的有V1、V3、V4、V5,,其中从V0到这些的临时最短路分别为INFINITE、60、30、100,从中找到最小的即V4,将V4标记为1,以V4为中间点(路径为V0->V4->Vx)更新最短路表,此时V3被更新为50,V5被更新为90。

4.进入第三次循环,此时未被标记的有V1、V3、V5,其中临时最短路分别为INFINITE、50、90,从中找到最小的即V3,将V3标记为1,以V3为中间点(路径为V0->V4->V3->Vx)更新最短路表,此时V5被更新成60。

5.进入第四次循环,此时未被标记的有V1、V5,其中临时最短路分别为INFINITE、60,找到V5,标记为1,以V5为中间点更新最短路表,此时没有元素被更新

6.进入第五次循环,这次循环没找到任何东西

7.退出循环,Len表中即为V0到其余各个节点的最短路径。

以上就是Dijkstra算法最基本的思想,当然,在找最短路时我们常常会需要求出到达某点的最短路径,那么下面,我们就来看看要怎样记录下到达某点的最短路径:

在Dijkstra算法的前提下加入查询最短路径其实很简单,只要在每次更新最短路时保存在该顶点的父节点序号即可,最后输出时回退中间节点然后用堆栈输出即可

最后直接给出完整代码(写的很一般,还望指教):

[cpp] view plain copy

#include   

#include   

#define MAX_LEN 100  

#define INFINITE 1000  

typedef struct graph  

{  

int nodenum;  

int edgenum;  

int matrix[MAX_LEN][MAX_LEN];  

}Graph;  


typedef struct stack  

{  

int  bottom;  

int  top;  

int  printout[MAX_LEN];  

}mstack;  


int in[MAX_LEN];  

int Len[MAX_LEN];  

int path[MAX_LEN];  


void InitStack(mstack *s)  

{  

    s->bottom=0;  

    s->top=0;  

memset(s->printout,0,sizeof(int)*MAX_LEN);  

}  


void push(mstack *s,int m)  

{  

    s->printout[s->top++]=m;  

}  


int pop(mstack *s)  

{  

return s->printout[--s->top];  

}  


void InitGraph(Graph *g,int n)  

{  

int i,j;  

for(i=1;i<=n;i++)  

for(j=1;j<=n;j++)  

          {  

if(i==j)g->matrix[i][j]=0;  

else g->matrix[i][j]=INFINITE;       

          }  


for(i=1;i<=n;i++)  

     {  

         in[i]=0;  

         Len[i]=INFINITE;  

         path[i]=0;  

     }  

}  


int main()  

{  


int n,m,i,A,B,templen,count,min,tempmin,si,temp;  

while(scanf("%d %d",&n,&m))  

    {  

        count=0;  

        Graph mGraph;  

        mGraph.edgenum=m;  

        mGraph.nodenum=n;  

        InitGraph(&mGraph,n);  

for(i=0;i

        {  

scanf("%d %d %d",&A,&B,&templen);  

            mGraph.matrix[A][B]=templen;  

        }  



        in[1]=1;  

path[1]=1;//sava path  

        Len[1]=0;  


for(i=2;i<=n;i++)  

        {  

Len[i]=mGraph.matrix[1][i];//Init the len  

if(Len[i]!=INFINITE)path[i]=1;  

        }  


        min=0;  

        si=1;  


while(count

        {  

            tempmin=INFINITE;  

for(i=1;i<=n;i++)  

            {  

if(in[i]==0&&Len[i]

                {  

                    tempmin=Len[i];  

                    si=i;  

                }  

            }  

            in[si]=1;  


for(i=1;i<=n;i++) //updata the length  

            {  

if(in[i]==0&&(tempmin+mGraph.matrix[si][i])

                {  

                    Len[i]=tempmin+mGraph.matrix[si][i];  

                    path[i]=si;  

                }  

            }  

            count++;   

        }  


        mstack s;  

for(i=1;i<=n;i++)  

        {  

            temp=i;  

            InitStack(&s);  

if(path[temp]==0)  

            {   

printf("no path\n");  

continue;  

            }  

while(path[temp]!=1)  

            {  

                push(&s,path[temp]);  

                temp=path[temp];  

            }  

printf("1-->");  

while(s.bottom!=s.top)  

            {  

printf("%d-->",pop(&s));  

            }  

printf("%d  min length is %d\n",i,Len[i]);  


        }  


    }  

return 0;  

}  

附上上面那张图的结果(V0节点为1,V1节点为2.。。。以此类推)

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,547评论 6 477
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,399评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,428评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,599评论 1 274
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,612评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,577评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,941评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,603评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,852评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,605评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,693评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,375评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,955评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,936评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,172评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 43,970评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,414评论 2 342

推荐阅读更多精彩内容