bzoj:1726: [Usaco2006 Nov]Roadblocks第二短路

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

 

 

 

SPFA，记录最大和次大，然后按更新位置搜就好了

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;

struct na{
    int y,z,ne;
    na(){
        ne=0;
    }
}b[200001];
struct ma{
    int x,t;
};
int x,y,z,p,n,m,num=0,l[5001],r[5001],dis[5001][2];
bool bo[5001][2];
char ch;
inline int read(){
    p=0;ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') p=p*10+ch-48, ch=getchar();
    return p;
}
inline void in(int x,int y,int z){
    num++;
    if (l[x]==0) l[x]=num;else b[r[x]].ne=num;
    b[num].y=y;b[num].z=z;r[x]=num;
}
inline void ins(int x,int y,int z){
    in(x,y,z);in(y,x,z);
}
const int INF=1e8;
queue <ma> q;
inline void pu(ma x){
    if (!bo[x.x][x.t]){
        bo[x.x][x.t]=1;
        q.push(x);
    }
}
int main(){
    register int i;
    n=read();m=read();
    for (i=1;i<=m;i++){
        x=read();y=read();z=read();
        ins(x,y,z);
    }
    for (i=1;i<=n;i++) dis[i][0]=dis[i][1]=INF;
    ma cmp;cmp.x=1;cmp.t=0;
    q.push(cmp);bo[1][0]=1;dis[1][0]=0;
    while(!q.empty()){
        ma k=q.front(),cmp;bo[k.x][k.t]=0;
        q.pop();
        if (k.x==n) continue;
        for (i=l[k.x];i;i=b[i].ne){
            if (!k.t){
                if (dis[b[i].y][0]>dis[k.x][0]+b[i].z) dis[b[i].y][1]=dis[b[i].y][0],dis[b[i].y][0]=dis[k.x][0]+b[i].z,cmp.x=b[i].y,cmp.t=0,pu(cmp);else
                if (dis[b[i].y][1]>dis[k.x][0]+b[i].z&&((b[i].y!=n)||(dis[b[i].y][0]!=dis[k.x][0]+b[i].z)))
                dis[b[i].y][1]=dis[k.x][0]+b[i].z,cmp.x=b[i].y,cmp.t=1,pu(cmp);
            }else
                if (dis[b[i].y][1]>dis[k.x][1]+b[i].z&&((b[i].y!=n)||(dis[b[i].y][0]!=dis[k.x][1]+b[i].z)))
                dis[b[i].y][1]=dis[k.x][1]+b[i].z,cmp.x=b[i].y,cmp.t=1,pu(cmp);
        }
    }
    printf("%d\n",dis[n][1]);
}