Hihocoder #1081 最短路径一 dijkstra

#1081 : 最短路径·一

时间限制:10000ms

单点时限:1000ms

内存限制:256MB

描述

万圣节的早上，小Hi和小Ho在经历了一个小时的争论后，终于决定了如何度过这样有意义的一天——他们决定去闯鬼屋！

在鬼屋门口排上了若干小时的队伍之后，刚刚进入鬼屋的小Hi和小Ho都颇饥饿，于是他们决定利用进门前领到的地图，找到一条通往终点的最短路径。

鬼屋中一共有N个地点，分别编号为1..N，这N个地点之间互相有一些道路连通，两个地点之间可能有多条道路连通，但是并不存在一条两端都是同一个地点的道路。那么小Hi和小Ho至少要走多少路程才能够走出鬼屋去吃东西呢？

/*

小Ho想了想说道：“唔……我觉得动态规划可以做，但是我找不到计算的顺序，如果我用f[i]表示从S到达编号为i的节点的最短距离的话，我并不能够知道f[1]..f[N]的计算顺序。”

“所以这个问题不需要那么复杂的算法啦，我就稍微讲讲你就知道了！”小Hi道：“路的长度不可能为负数对不对？”

“那是自然，毕竟人类还没有发明时光机器……”小Ho点点头。

于是小Hi问道：“那么如果就看与S相邻的所有节点中与S最近的那一个S'，并且从S到S'的距离为L，那么有可能存在另外的道路使得从S到S'的距离小于L么？”

“不能，因为S'是与S相邻的所有节点中与S最近的节点，那么从S到其他相邻点的距离一定是不小于L的，也就是说无论接下来怎么走，回到L点时总距离一定大于L。”小Ho思考了一会，道。

“也就是说你已经知道了从S到S'的最短路径了是么？”小Hi继续问道。

“是的，这条最短路径的长度是L。”小Ho答道。

小Hi继续道：“那么现在，我们不妨将S同S'看做一个新的节点？称作S1，然后我就计算与S相邻或者与S'相邻的所有节点中，与S最近的哪一个节 点S''。注意，在这个过程中，与S相邻的节点与S的距离在上一步就已经求出来了，那么我要求的只有与S'相邻的那些节点与S的距离——这个距离等于S与 S'的距离加上S'与这些结点的距离，对于其中重复的节点——同时与S和S'相邻的节点，取两条路径中的较小值。”

小Ho点了点头：“那么同之前一样，与S1（即S与S'节点）相邻的节点中与S'距离最近的节点如果是S''的话，并且这个距离是L2，那么我们可以知道S到S''的最短路径的长度便是L2，因为不可能存在另外的道路比这个更短了。”

于是小Hi总结道：“接下来的问题不就很简单了么，只需要以此类推，每次将与当前集合相邻（即与当前集合中任意一个元素）的所有节点中离S最近的节 点（这些距离可以通过上一次的计算结果推导而出）选出来添加到当前集合中，我就能够保证在每一个节点被添加到集合中时所计算的离S的距离是它与S之间的最 短路径！”

“原来是这样！但是我的肚子更饿了呢！”言罢，小Ho的肚子咕咕叫了起来。

*/

输入

每个测试点（输入文件）有且仅有一组测试数据。

在一组测试数据中：

第1行为4个整数N、M、S、T，分别表示鬼屋中地点的个数和道路的条数，入口（也是一个地点）的编号，出口（同样也是一个地点）的编号。

接下来的M行，每行描述一条道路：其中的第i行为三个整数u_i, v_i, length_i，表明在编号为u_i的地点和编号为v_i的地点之间有一条长度为length_i的道路。

对于100%的数据，满足N<=10^3，M<=10^4, 1 <= length_i <= 10^3, 1 <= S, T <= N, 且S不等于T。

对于100%的数据，满足小Hi和小Ho总是有办法从入口通过地图上标注出来的道路到达出口。

输出

对于每组测试数据，输出一个整数Ans，表示那么小Hi和小Ho为了走出鬼屋至少要走的路程。

样例输入

5 23 5 4
1 2 708
2 3 112
3 4 721
4 5 339
5 4 960
1 5 849
2 5 98
1 4 99
2 4 25
2 1 200
3 1 146
3 2 106
1 4 860
4 1 795
5 4 479
5 4 280
3 4 341
1 4 622
4 2 362
2 3 415
4 1 904
2 1 716
2 5 575

样例输出

123

//Dijstra
//hiho23
#include<algorithm>
#include<cassert>
#include<cstdio>
#include<cmath>
#include<ctime>
#include<cstring>
using namespace std;

int S;
int T;
int N,M,i,j,k,m[1024][1024];
bool v[1024];

int main(){
    scanf("%u%u%u%u",&N,&M,&S,&T);
    {
        memset(m,66,sizeof(m));
        while(M--){
            scanf("%u%u%u",&i,&k,&j);
            if(j<m[k][i]){
                m[k][i]=m[i][k]=j;
            }
        }
        memset(&v[1],0,N);
        v[S]=true;
        for(;;){
            k=0;
            for(i=1;i<=N;i++){
                if(v[i]==0&&m[S][i]<m[S][k]){
                    k=i;
                }
            }
            if(k==T){
                printf("%u\n",m[S][T]);
                break;
            }
            v[k]=1;
            for(i=1;i<=N;i++){
                if(v[i]==0&&m[S][i]>(m[k][i]+m[S][k])){
                    m[S][i]=(m[k][i]+m[S][k]);
                }
            }
        }
    }
    return 0;
}