最短路径算法 2.Dijkstra算法

Dijkstra 算法解决的是带权重的有向图上单源最短路径问题，该算法要求所有边的权重都为非负值。该算法的时间复杂度是O(N²)，相比于处理无负权的图时，比Bellmad-Ford算法效率更高。

算法描述：

首先引用《算法导论》中的一段比较官方的话，如果可以看懂，那下一部分就可以跳过了：

“Dijkstra算法在运行过程中维持的关键信息是一组结点集合S。从源结点s到该集合中每个结点之间的最短路径已经被找到。算法重复从结点集 V - S 中算则最短路径估计的最小的结点 u ,将 u 加入到集合S，然后对所有从 u 出发的边进行松弛。” 所谓松弛操作，简单的说就是更新两点间的最短距离。

不是很好理解对吧，那么下面的描述是更容易理解的一种描述：

设起始点为s，dis[v]表示s点到v点的最短路径，pre[v]是v的前驱结点，用来输出路径。

1、初始化：dis[v]=∞(v≠s) dis[s]=0,pre[s]=0;

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

(1)在没有被访问过的点中，即上述的V - S集合，找到一个点 u 使得dis[u]是最小的。

(2)标记 u 为已确定的最短路径。

(3)for(每个与 u 相连且没有确定过最短距离的点 v)

if(dis[u]+m[u][v]<dis[v]){

dis[v]=dis[u]+m[u][v];

pre[v]=u;

}

3、结束：结果dis[v]就是s到v的最短距离。

算法理解：

其中自以为有几点理解需要说明：

1、为什么用到中间点？

2、取s到中间点的距离时采用什么策略？

第一个问题，从起点到另一个点的最短路径至少会经历一个中间点，所以我们要求出经过这个中间的到另一个点的路径，就要先求出起点到中间点的最短路径。

第二个问题，其实这里采用的是一中贪心的策略。当然这个策略可以被严格证明是正确的，但是我也一知半解，只知道是可以被证明的，在这里也就不浪费时间了。(详细可以参考《算法导论》)

最后，解释一下为什么有负权边的时候不可以：

连接矩阵如下(图可以自己在旁边画一下)：

1 2 3

1 \ 2 1

2 2 \ -4

3 1 -4 \

那么第一次标记的点就为3并且把dis[3]记为1，但实际上dis[3]应该时-2，因此就会出现错误。

最后附上一段不那么标准的代码：

 #include<stdio.h>
 #include<stdlib.h>
 int m[][],e,dist[],n,b[],pre[],dist[];
 void dij(int s){
      b[s]=;
      int i,j;
      for(i=;i<=n;i++)
       dist[i]=m[s][i];
      dist[s]=;
      pre[s]=;

      for(i=;i<=n;i++){
       int min=,k=;
       for(j=;j<=n;j++)
        if(b[j]!= && dist[j]<min)
         {min=dist[j];k=j;}
        b[k]=;
        for(j=;j<=n;j++)
         if(min+m[k][j]<dist[j]&&b[j]!=)
          {
           dist[j]=min+m[k][j];
           pre[j]=i;
           }
       }
       for(i=;i<=n;i++)
        if(i!=s)
         printf("%d ",dist[i]);
      }
 int main(){
     int i,j;
     scanf("%d%d",&n,&e);
     memset(b,,sizeof(b));
     memset(m,,sizeof(m));
     for(i=;i<=e;i++){
      int x,y;
      scanf("%d%d",&x,&y);
      scanf("%d",&m[x][y]);
      }
     int w;
     scanf("%d",&w);
     dij(w);
     system("pause");
     return ;
     }