1726: [Usaco2006 Nov]Roadblocks第二短路

1726: [Usaco2006 Nov]Roadblocks第二短路

Time Limit: 5 Sec   Memory Limit: 64 MB
Submit: 1178   Solved: 564
[ Submit][ Status][ Discuss]

Description

贝茜把家搬到了一个小农场，但她常常回到FJ的农场去拜访她的朋友。贝茜很喜欢路边的风景，不想那么快地结束她的旅途，于是她每次回农场，都会选择第二短的路径，而不象我们所习惯的那样，选择最短路。 贝茜所在的乡村有R(1<=R<=100,000)条双向道路，每条路都联结了所有的N(1<=N<=5000)个农场中的某两个。贝茜居住在农场1，她的朋友们居住在农场N（即贝茜每次旅行的目的地）。 贝茜选择的第二短的路径中，可以包含任何一条在最短路中出现的道路，并且，一条路可以重复走多次。当然咯，第二短路的长度必须严格大于最短路（可能有多条）的长度，但它的长度必须不大于所有除最短路外的路径的长度。

Input

* 第1行: 两个整数，N和R，用空格隔开

* 第2..R+1行: 每行包含三个用空格隔开的整数A、B和D，表示存在一条长度为 D(1 <= D <= 5000)的路连接农场A和农场B

Output

* 第1行: 输出一个整数，即从农场1到农场N的第二短路的长度

Sample Input

4 4
1 2 100
2 4 200
2 3 250
3 4 100

Sample Output

450

输出说明:

最短路：1 -> 2 -> 4 (长度为100+200=300)
第二短路：1 -> 2 -> 3 -> 4 (长度为100+250+100=450)

HINT

Source

Gold

[ Submit][ Status][ Discuss]
先用SPFA算出从起点到每个点的最短路，记为f[i] 然后是终点到每个点的最短路，记为g[i] 最后暴力枚举一下每条边走几次，用f,g更新下答案就行了 显然即使不是简单路径，次短路上一条边也不会走很多次的

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<vector>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<bitset>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;
 
const int maxn = 5050;
const int INF = ~0U>>1;
 
struct E{
    int to,w; E(){}
    E(int to,int w): to(to),w(w){}
};
 
int n,m,Ans = INF,f[maxn],g[maxn];
bool vis[maxn];
 
vector <E> v[maxn];
queue <int> Q;
 
void SPFA(int S,int *dis)
{
    for (int i = 1; i <= n; i++) dis[i] = INF;
    Q.push(S); vis[S] = 1; dis[S] = 0;
    while (!Q.empty())
    {
        int k = Q.front(); Q.pop(); vis[k] = 0;
        for (int i = 0; i < v[k].size(); i++)
        {
            E e = v[k][i];
            if (dis[e.to] > dis[k] + e.w)
            {
                dis[e.to] = dis[k] + e.w;
                if (!vis[e.to]) vis[e.to] = 1,Q.push(e.to);
            }
        }
    }
}
 
int main()
{
    #ifdef DMC
        freopen("DMC.txt","r",stdin);
    #endif
     
    cin >> n >> m;
    while (m--)
    {
        int x,y,w;
        scanf("%d%d%d",&x,&y,&w);
        v[x].push_back(E(y,w));
        v[y].push_back(E(x,w));
    }
    SPFA(1,f); SPFA(n,g);
    for (int i = 1; i <= n; i++)
        for (int j = 0; j < v[i].size(); j++)
        {
            E e = v[i][j];
            int now = f[i] + e.w + g[e.to];
            if (now > g[1]) Ans = min(Ans,now);
            now = f[i] + 2 * e.w + g[i];
            if (now > g[1]) Ans = min(Ans,now);
            now = f[i] + 3 * e.w + g[e.to];
            if (now > g[1]) Ans = min(Ans,now);
        }
    cout << Ans << endl;
    return 0;
}