传送门biu~
dis[i][j]代表到第i个点速度为j时需要的最少时间,spfa转移即可。
#include<bits/stdc++.h>
#define pa pair<int,int>
#define N 155
#define M 505
using namespace std;
int n,m,d;pa from[N][M];queue<pa>q;
int head[N],to[N*N],nex[N*N],v[N*N],l[N*N],tp;
inline void add(int x,int y,int p,int w){
nex[++tp]=head[x];
head[x]=tp;
to[tp]=y;
v[tp]=p;
l[tp]=w;
}
double dis[N][M];bool in[N][M];
inline void spfa(){
memset(dis,127,sizeof dis);
dis[0][70]=0; q.push(pa(0,70));
while(!q.empty()){
int x=q.front().first,w=q.front().second;
q.pop(); in[x][w]=false;
for(int i=head[x];i;i=nex[i]){
if(!v[i] && dis[x][w]+1.0*l[i]/w<dis[to[i]][w]){
from[to[i]][w]=pa(x,w);
dis[to[i]][w]=dis[x][w]+1.0*l[i]/w;
if(!in[to[i]][w]) q.push(pa(to[i],w)),in[to[i]][w]=true;
}
if(v[i] && dis[x][w]+1.0*l[i]/v[i]<dis[to[i]][v[i]]){
from[to[i]][v[i]]=pa(x,w);
dis[to[i]][v[i]]=dis[x][w]+1.0*l[i]/v[i];
if(!in[to[i]][v[i]]) q.push(pa(to[i],v[i])),in[to[i]][v[i]]=true;
}
}
}
}
void print(int x,int w){
if(x) print(from[x][w].first,from[x][w].second);
printf("%d ",x);
}
int main(){
scanf("%d%d%d",&n,&m,&d);
for(int i=1;i<=m;++i){
int u,v,p,w;
scanf("%d%d%d%d",&u,&v,&p,&w);
add(u,v,p,w);
}
spfa();
int ans=0;
for(int i=1;i<=500;++i) if(dis[d][i]<dis[d][ans]) ans=i;
print(d,ans);
return 0;
}