[bzoj3245]最快路线

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define N 205
 4 struct ji{
 5     int nex,to,len,v;
 6 }edge[N*N];
 7 struct ji2{
 8     int k,v;
 9 }from[N][505];
10 queue<ji2>q;
11 int E,n,m,k,x,y,v,z,head[N],ans[N],vis[N][505];
12 double d[N][505];
13 void add(int x,int y,int z,int v){
14     edge[E].nex=head[x];
15     edge[E].to=y;
16     edge[E].len=z;
17     edge[E].v=v;
18     head[x]=E++;
19 }
20 void spfa(){
21     for(int i=1;i<=n;i++)
22         for(int j=1;j<=500;j++)d[i][j]=0x3f3f3f3f;
23     d[0][70]=0;
24     q.push(ji2{0,70});
25     while (!q.empty()){
26         ji2 k=q.front();
27         q.pop();
28         vis[k.k][k.v]=0;
29         for(int i=head[k.k];i!=-1;i=edge[i].nex){
30             int v=edge[i].to,flag=0;
31             if (!edge[i].v){
32                 flag=1;
33                 edge[i].v=k.v;
34             }
35             double x=1.0*edge[i].len/edge[i].v;
36             if (d[k.k][k.v]+x<d[v][edge[i].v]){
37                 from[v][edge[i].v]=k;
38                 d[v][edge[i].v]=d[k.k][k.v]+x;
39                 if (!vis[v][edge[i].v]){
40                     vis[v][edge[i].v]=1;
41                     q.push(ji2{v,edge[i].v});
42                 }
43             }
44             if (flag){
45                 k.v=edge[i].v;
46                 edge[i].v=0;
47             }
48         }
49     }
50 }
51 int main(){
52     scanf("%d%d%d",&n,&m,&k);
53     memset(head,-1,sizeof(head));
54     for(int i=1;i<=m;i++){
55         scanf("%d%d%d%d",&x,&y,&v,&z);
56         add(x,y,z,v);
57     }
58     spfa();
59     double s=d[k][1];
60     for(int i=2;i<=500;i++)s=min(s,d[k][i]);
61     for(int i=1;i<=500;i++)
62         if (d[k][i]==s){
63             ans[++ans[0]]=k;
64             ji2 t=from[k][i];
65             while (t.k){
66                 ans[++ans[0]]=t.k;
67                 t=from[t.k][t.v];
68             }
69         }
70     printf("0");
71     for(int i=ans[0];i;i--)printf(" %d",ans[i]);
72 }