HDU-1874 畅通工程续

畅通工程续

Time Limit: 3000/1000 MS (Java/Others)
Memory Limit: 32768/32768 K (Java/Others)

Problem 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

Reference Code

Dijkstra算法，0MS，1584K：

#include <cstdio>
#include <cstring>
#include <vector>
#include <queue>
const int maxn=1e4+3;
const int INF=0x3f3f3f3f;
struct node{
    int v,w;
    bool operator<(const node &oth)const{
        return oth.w<w;
    }
};
std::vector<node> G[maxn];
int n,m,d[maxn];
bool vis[maxn];
void init(){
    for (int i=0;i<maxn;++i) G[i].clear();
    memset(vis,false,sizeof(vis));
    memset(d,INF,sizeof(d));
}
void addedge(){
    int u,v,w;
    scanf("%d%d%d",&u,&v,&w);
    G[u].push_back((node){v,w});
    G[v].push_back((node){u,w});
}
int dijkstra(int st,int ed){
    std::priority_queue<node> pq;
    pq.push((node){st,0});
    d[st]=0;
    while (!pq.empty()){
        node now=pq.top();pq.pop();
        int v=now.v;
        if (vis[v]) continue;
        vis[v]=true;
        for (int i=0;i<G[v].size();++i){
            int vi=G[v][i].v,w=G[v][i].w;
            if (!vis[vi]&&d[vi]>w+d[v]){
                d[vi]=w+d[v];
                pq.push((node){vi,d[vi]});
            }
        }
    }
    return d[ed];
}
int main(){
    while(~scanf("%d%d",&n,&m)){
        init();
        for(int i=0;i<m;i++) addedge();
        int st,ed;
        scanf("%d%d",&st,&ed);
        int res=dijkstra(st,ed);
        if (INF!=res) printf("%d\n",res);
        else printf("-1\n");
    }
    return 0;
}

Bellman-ford算法，15MS，1236K：

#include <cstdio>
#include <cstring>
#include <algorithm>
const int maxn=1e3+3;
const int INF=0x3f3f3f3f;
struct edge{int u,v,w;}edge[maxn*2];
int n,m,dist[maxn];
void inputedge(int E){
    for (int i=0;i<E;++i){
        scanf("%d%d%d",&edge[i].u,&edge[i].v,&edge[i].w);
        edge[i+E].v=edge[i].u;
        edge[i+E].u=edge[i].v;
        edge[i+E].w=edge[i].w;
    }
}
void bellmanford(int st,int V,int E) {
    memset(dist,INF,sizeof(dist));
    dist[st]=0;
    for (int i=0;i<V;++i)
        for (int j=0;j<2*E;++j)
            dist[edge[j].u]=std::min(dist[edge[j].u],dist[edge[j].v]+edge[j].w);
}

int main() {
    int s, t;
    while (~scanf("%d%d",&n,&m)) {
        inputedge(m);
        scanf("%d%d",&s,&t);
        bellmanford(s,n,m);
        if (dist[t]<INF) printf("%d\n",dist[t]);
        else printf("-1\n");
    }
    return 0;
}

Floyd算法，78MS，5156K：

#include <cstdio>
#include <cstring>
#include <algorithm>
const int maxn=1e3+3;
const int INF=0x3f3f3f3f;
int maze[maxn][maxn];
void init(){
    memset(maze,INF,sizeof(maze));
    for (int i=0;i<maxn;++i) maze[i][i]=0;
}
void inputedge(int E){
    for (int i=0;i<E;++i){
        int u,v,w;
        scanf("%d%d%d",&u,&v,&w);
        if (maze[u][v]>w) maze[u][v]=maze[v][u]=w;
    }
}
void floyd(int V) {
    for (int i=0;i<V;++i)
        for (int j=0;j<V;++j)
            for (int k=0;k<V;++k)
                maze[j][k]=
                std::min(maze[j][k],maze[j][i]+maze[i][k]);
}
int main(){
    int n,m,s,t;
    while (~scanf("%d%d",&n,&m)) {
        init();
        inputedge(m);
        floyd(n);
        scanf("%d%d",&s,&t);
        if (maze[s][t]<INF) printf("%d\n",maze[s][t]);
        else printf("-1\n");
    }
    return 0;
}

SPFA算法，46MS，5220K：

#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
const int maxn=1e3+3;
const int INF=0x3f3f3f3f;
int maze[maxn][maxn],vis[maxn],dist[maxn];
void init(){
    memset(vis,false,sizeof(vis));
    memset(dist,INF,sizeof(dist));
    memset(maze,INF,sizeof(maze));
    for (int i=0;i<maxn;++i) maze[i][i]=0;
}
void inputedge(int E){
    for (int i=0;i<E;++i){
        int u,v,w;
        scanf("%d%d%d",&u,&v,&w);
        if (maze[u][v]>w) maze[u][v]=maze[v][u]=w;
    }
}
void SPFA(int st,int n){
    std::queue<int> q;
    q.push(st);
    dist[st]=0;
    vis[st]=true;
    while (!q.empty()){
        int top=q.front();q.pop();
        vis[top]=false;
        for (int i=0;i<n;++i){
            int buf=dist[top]+maze[top][i];
            if (dist[i]>buf){
                dist[i]=buf;
                if (!vis[i]){
                    q.push(i);
                    vis[i]=true;
                }
            }
        }
    }
}
int main() {
    int n,m,s,t;
    while (~scanf("%d%d",&n,&m)) {
        init();
        inputedge(m);
        scanf("%d%d",&s,&t);
        SPFA(s,n);
        if (dist[t]<INF) printf("%d\n",dist[t]);
        else printf("-1\n");
    }
    return 0;
}

Tips

SSSP的经典裸题，现给出了四种算法，以本题来看，Dijkstra算法用时最短，Bellman-ford算法占内存最少。综合来看还是Dijkstra较优。