考虑一下最短路,但是有的边边权是不确定的,我们在spfa的时候把这些边推到另一个队列里,spfa完事之后取出边权不定的边的队列里的一个边,找到能使他边权最小的前一个边后暂时定下这条边的边权,然后从这条边的终点再跑一次spfa(每次跑spfa都把当前不在边权不定队列里的边权不定边推进去)(即使不能更新终点的最短路也要跑,因为可能更新终点的连出去的不限速边的数量),除非这个点已经在这一部分入队超过50次了
记录一下每一条边的前一条边是哪一个(记录点不行,因为可能走的不是到这个点的最短路的边),然后倒叙输出即可
这个算动点spfa吧-_-还是spfa套spfa-_-并不清楚
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<bitset>
using namespace std;
#define MAXN 160
#define MAXM 22510
#define ll long long
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
struct vec{
int to;
int fro;
int v;
double l;
int num;
vec(){
}
vec(int _to,int _fro,int _v,double _l,int _num){
to=_to;
fro=_fro;
v=_v;
l=_l;
num=_num;
}
};
vec mp[MAXM];
int tai[MAXN],cnt;
int n,m,T;
int v1[MAXM],v2[MAXM],v[MAXM],l[MAXM];
vector<vec>r[MAXN];
queue<vec>qm;
double dis[MAXN];
int q[MAXN],hd,tl;
bool iqm[MAXM];
int rd[MAXM];
int Rd[MAXN];
bool iq[MAXN];
int ans[MAXN];
int mxv[MAXN];
int tcn[MAXN];
inline void be(int x,int y,int v,double l,int num){
mp[++cnt].to=y;
mp[cnt].fro=tai[x];
tai[x]=cnt;
mp[cnt].v=v;
mp[cnt].l=l;
mp[cnt].num=num;
}
void spfa(){
int i,x,y;
while(hd!=tl){
iq[x=q[hd++]]=0;
hd%=MAXN;
for(i=tai[x];i;i=mp[i].fro){
y=mp[i].to;
if(!mp[i].v){
if(!iqm[mp[i].num]){
qm.push(vec(mp[i].to,x,0,mp[i].l,mp[i].num));
iqm[mp[i].num]=1;
}
}else{
rd[mp[i].num]=Rd[x];
if(dis[y]>dis[x]+1.0*mp[i].l/mp[i].v){
dis[y]=dis[x]+1.0*mp[i].l/mp[i].v;
Rd[y]=mp[i].num;
if(!iq[y]){
iq[q[tl++]=y]=1;
tl%=MAXN;
}
}
}
}
}
}
int main(){
int i,x,xx,y;
scanf("%d%d%d",&n,&m,&T);
T++;
for(i=1;i<=m;i++){
scanf("%d%d%d%d",&v1[i],&v2[i],&v[i],&l[i]);
v1[i]++;
v2[i]++;
be(v1[i],v2[i],v[i],l[i],i);
if(v[i]){
r[v2[i]].push_back(vec(v2[i],v1[i],v[i],l[i],i));
}
}
r[1].push_back(vec(1,0,70,0,0));
for(i=2;i<=n;i++){
dis[i]=1e60;
}
iq[q[tl++]=1]=1;
spfa();
while(!qm.empty()){
vec t=qm.front();
iqm[t.num]=0;
qm.pop();
x=t.fro;
y=t.to;
double d=1e60;
int v;
for(i=0;i<r[x].size();i++){
if(dis[r[x][i].fro]+1.0*(r[x][i].l+t.l)/r[x][i].v<d){
d=dis[r[x][i].fro]+1.0*(r[x][i].l+t.l)/r[x][i].v;
v=r[x][i].v;
xx=r[x][i].num;
}
}
bool flag=0;
r[y].push_back(vec(y,0,v,v*d,t.num));
if(dis[y]>d){
dis[y]=d;
Rd[y]=t.num;
flag=1;
}
rd[t.num]=xx;
if(tcn[y]<=50){
tcn[y]++;
iq[q[tl++]=y]=1;
tl%=MAXN;
spfa();
}
}
int tot=0;
ans[++tot]=T;
for(i=Rd[T];i;i=rd[i]){
ans[++tot]=v1[i];
}
for(i=tot;i!=1;i--){
printf("%d ",ans[i]-1);
}
printf("%d\n",ans[1]-1);
return 0;
}
/*
4 6 3
0 1 0 5
1 2 0 4
2 3 0 6
3 2 1 0
3 1 0 1
0 2 2 0
*/