最短路问题

一.五行算法----Floyd-Warshall算法

    假设现在有下面这样一副图：

很容易得到这个图矩阵为

0 2 6 4

* 0 3 * 

7 * 0 1

5 * 12 0（*为不可到达）

不过从图中不难看出如果从1到3 不是直接到达 而是通过2再到达3 那么花费的路程为5 5<6 故此时可以更新从1到3的最小值 以此类推 看从一点到达另外一点借助其他点到达是否会更小一点 找到这个最小值并更新这个最小值。

这样得到的最小矩阵就是

0 2 5 4

9 0 3 4

6 8 0 1

5 7 10 0

接下来是代码实现这个矩阵的转换

for(k=1;k<=n;k++)//k来寻找中间可以借助的点
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
            if(e[i][j]>e[i][k]+e[k][j])//如果原来存在的i与j之间的距离借助k点可以更小 则更新这个距离
                e[i][j]=e[i][k]+e[k][j];

简短的五行算法 可以求出所有点之间的最短距离 注意循环的顺序不能颠倒。

二.单元最短路----Dijkstra算法

假设现在有下面这样一副图：

用e数组储存这个图的权边为

0 1 12 * * *

* 0 9 3 * * 

* * 0 * 5 *

* * 4 0 13 15

* * * * 0 4

* * * * * 0（*为无法到达 记为无穷大）

找单源最短路 我们以从1到其他路为例 先设一个数组dis为1到其他点的权值

得到dis为 0 1 12 * * * 

从1到1就不用说了为0

先对其他个点进行一次遍历 看从1到其他各点的最小值是多少 记录是到哪个点 结果是到2点最小为1 标记2这个点 看其他点当从1经过2到达时的值 如果比原始值小的话 更新dis数组

接下来继续进行一次遍历 遍历找最小值（已经找到过的最小值点不用再记录）这次找到的最小权值是4 是从1到4点的权值 故这次找其他各点从1经过4到达时的值 如果比原始值小的话更新dis数组

继续进行这样的循环 直到所有的点都已找完 

此时的dis数组的值为 0 1 8 4 13 17

这就是从1到其他各点的最小值 即单源最短路的算法

接下来是代码实现：

# include <iostream>
# include <cstdio>
# include <cstdio>
# include <cstring>
# define inf 0x3f3f3f
using namespace std;
int main()
{
    int dis[110];//储存单源权值
    int e[110][110];//储存各点之间的权值
    int n,m;//n为多少个点 m为有多少种关系
    int i,j,k;
    int x,y,z;
    int b[110];//记录点是否被使用过
    memset(b,0,sizeof(b));
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=n;j++)
        {
            if(i==j) e[i][j]=0;
            else e[i][j]=inf;
        }
    }
    for(i=1;i<=m;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        e[x][y]=z;//有向图
    }
    for(i=1;i<=n;i++)
        dis[i]=e[1][i];
    b[1]=1;//从1到1不用进行记录
    for(i=1;i<=n-1;i++)//进行n-1次因为1到1不用进行判断 只有n-1个点
    {
        int min=inf;
        for(j=1;j<=n;j++)
        {
            if(b[j]==0&&dis[j]<min)//寻找没有找过的点 并判断是否是最小值
            {
                min=dis[j];//更新最小值
                k=j;//记录最小值点
            }
        }
        b[k]=1;//表示该点已经找过 下一轮不用再找
        for(j=1;j<=n;j++)
        {
            if(dis[j]>dis[k]+e[k][j])//查看通过最小的k点到j点能否使权值变小
            {
                dis[j]=dis[k]+e[k][j];
            }
        }
    }
    for(i=1;i<=n;i++)
        printf("%d ",dis[i]);
    printf("\n");
}
/*
 6 9
 1 2 1
 1 3 12
 2 3 9
 2 4 3
 3 5 5
 4 3 4
 4 5 13
 4 6 15
 5 6 4
 
 0 1 8 4 13 17
 */