ZCMU 1398: 工程【最短路径】

ZCMU 1398: 工程

Time Limit: 1 Sec  Memory Limit: 128 MB

Description

某省自从实行了很多年的畅通工程计划后，终于修建了很多路。不过路多了也不好，每次要从一个城镇到另一个城镇时，都有许多种道路方案可以选择，而某些方案要比另一些方案行走的距离要短很多。这让行人很困扰。

现在，已知起点和终点，请你计算出要从起点到终点，最短需要行走多少距离。

Input

本题目包含多组数据，请处理到文件结束。
每组数据第一行包含两个正整数N和M(0<N<200,0<M<1000)，分别代表现有城镇的数目和已修建的道路的数目。城镇分别以0～N-1编号。
接下来是M行道路信息。每一行有三个整数A,B,X(0<=A,B<N,A!=B,0<X<10000),表示城镇A和城镇B之间有一条长度为X的双向道路。
再接下一行有两个整数S,T(0<=S,T<N)，分别代表起点和终点。

Output

对于每组数据，请在一行里输出最短需要行走的距离。如果不存在从S到T的路线，就输出-1.

Sample Input

3 3

0 1 1

0 2 3

1 2 1

0 2

3 1

0 1 1

1 2

Sample Output

2

-1

思路

普通的求点间最短路径问题，没有负权边，负权回路之类的骚操作，可随意选用Floyd,Dijkstra,Bellman-Ford多种算法A之，这里我用了Floyd算法。

方法：插点法，任意两个城市之间的距离要想缩短，只能通过引入第三点，通过这个点进行中转。注意此处不能用两点间线段最短来思考。

        for(int k=1;k<=N;k++)//在当前最短路径允许通过k点中转的前提下
        {
            for(i=1;i<=N;i++)//从i点到j点
            {
                for(int j=1;j<=N;j++)
                {
                    //如果当前两点间最短路径小于通过k点中转的路径，就更新当前的最短路径
                    if (map[i][j] > map[i][k] + map[k][j])
                        map[i][j] = map[i][k] + map[k][j];
                }
            }
        }

基本思想：逐轮更新从i号顶点到j号顶点允许经过k号顶点中转的最短路程（对从i号顶点到j号顶点的边进行松弛）。

注意每一轮更新后最短路径会变为允许中转后的最短路径，也就是说经过k轮更新后（对所有的边进行k轮松弛后），所得的最短路径就是允许通过任意一点中转后得到的两点间最短路径。接下来就看题目想要从哪去哪了，读入start和end，输出map[start][end]即可，由于本题任何道路都是双向的，没有市际单行道这种骚操作，因此输出map[end][start]也是可以的。

AC代码

#include <iostream>
#include <cstdio>
#include <vector>
#include <cstring>
#include <cmath>
#include <cstring>
#include <string>
#include <list>
#include <algorithm>
#include <set>
#include <list>
#define INF 0x3f3f3f3f
 
using namespace std;
 
int map[200][200];
 
int main()
{
    int N, M;
    while(scanf("%d%d",&N,&M)!=EOF)
    {
        int i, a, b;
        //初始化，将由城市本身到本身的路程设为0，其余均设为无穷大INF
        memset(map, INF, sizeof(int) * 40000);
        for (i = 0; i < N; i++)map[i][i] = 0;
 
 
        for(i=0;i<M;i++)
        {
            scanf("%d%d", &a, &b);
            scanf("%d", &map[a][b]);
            //双向道路，距离相同
            map[b][a] = map[a][b];
        }


        //Floyd-Warshall算法核心代码
        for(int k=0;k<N;k++)
        {
            for(i=0;i<N;i++)
            {
                for(int j=0;j<N;j++)
                {
                    if (map[i][j] > map[i][k] + map[k][j])
                        map[i][j] = map[i][k] + map[k][j];
                }
            }
        }


        int start, end;
        scanf("%d%d", &start, &end);
        if (map[start][end] != INF)printf("%d\n", map[start][end]);
        else printf("-1\n");
    }
    //system("pause");
    return 0;
}