bzoj3245 最快路线(二维？spfa)

只记到每个点的最小时间显然是不行的。我们还需要知道到这个点的速度是多少。因此我们用d[i][j]表示到i，速度为j的最小时间。最后因为要输出路径，所以我们再记录一个from数组，表示d[i][j]取到最小值的前一条边是谁。最后输出时，递归输出即可。

#include <bits/stdc++.h>
using namespace std;
#define pa pair<int,int>
#define ll long long
#define inf 0x3f3f3f3f
#define N 155
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,t,h[N],num=0;
bool inq[N][505];
double d[N][505];
pa from[N][505];
struct edge{
    int to,next,v,l;
}data[N*N];
void spfa(){
    queue<pa>q;memset(d,127,sizeof(d));
    q.push(make_pair(0,70));inq[0][70]=1;d[0][70]=0;
    while(!q.empty()){
        int x=q.front().first,vx=q.front().second;q.pop();inq[x][vx]=0;
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to,vy=data[i].v,l=data[i].l;
            if(vy&&d[x][vx]+(double)l/vy<d[y][vy]){
                d[y][vy]=d[x][vx]+(double)l/vy;
                from[y][vy]=make_pair(x,vx);
                if(!inq[y][vy]) inq[y][vy]=1,q.push(make_pair(y,vy));
            }
            if(!vy&&d[x][vx]+(double)l/vx<d[y][vx]){
                d[y][vx]=d[x][vx]+(double)l/vx;
                from[y][vx]=make_pair(x,vx);
                if(!inq[y][vx]) inq[y][vx]=1,q.push(make_pair(y,vx));
            }
        }
    }
}
void print(int x,int vx){
    if(x!=0) print(from[x][vx].first,from[x][vx].second);
    printf("%d",x);
    if(x!=t) printf(" ");
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();t=read();
    while(m--){
        int x=read(),y=read(),v=read(),l=read();
        data[++num].to=y;data[num].next=h[x];h[x]=num;
        data[num].v=v;data[num].l=l;
    }
    spfa();
    int ans=0;
    for(int i=1;i<=500;++i) if(d[t][i]<d[t][ans]) ans=i;
    print(t,ans);
    return 0;
}