好惡心啊……調的累死了
寫法比較鬼畜,不過應該是跑的最快的做法吧。。。lc233大佬老早就A掉了我特麼還調的要死
裂點,dist[x][y]表示從起點出發,上一步是y,走到x的最短路徑
這樣記錄有什麼用呢?
因爲我們更新x->t只需要用到y->x這條邊時的限速,所以y之前的怎麼走怎麼走其實是不需要知道的,只需要知道一個最小值就行了。
然後直接跑spfa記錄方案
記錄二維狀態,隊列中存儲sped表示這一個位置上的狀態的速度
莫名其妙跑的飛快
#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <queue>
#include <string>
#include <map>
#include <cstring>
#include <ctime>
#include <vector>
#define maxn 1000001
#define inf 1e9
#define ll long long
#define For(i,j,k) for(int i=j;i<=k;i++)
#define Dow(i,j,k) for(int i=k;i>=j;i--)
using namespace std;
int poi[maxn],nxt[maxn],f[maxn];
double dist[501][501],v[maxn],sp[maxn],sped[maxn*4];
int inq[501][501],ans[501][501],Ans[501],tot;
int tx[maxn],ty[maxn],cnt,n,m,d,x,y,z,t;
inline void add(int x,int y,int z,int t)
{
poi[++cnt]=y;nxt[cnt]=f[x];v[cnt]=z;sp[cnt]=t;f[x]=cnt;
}
inline void dfs(int x,int y)
{
Ans[++tot]=x;if(!x) return;dfs(y,ans[x][y]);
}
int main()
{
scanf("%d%d%d",&n,&m,&d);
For(i,0,n) For(j,0,n) dist[i][j]=1e9;
For(i,1,m)
scanf("%d%d%d%d",&x,&y,&t,&z),add(x,y,z,t);
tx[1]=0;ty[1]=n+1;dist[0][0]=0;sped[1]=70;
int l=1,r=1;
while(l<=r)
{
int x=tx[l],y=ty[l];
inq[x][y]=0;
for(int i=f[x];i;i=nxt[i])
{
if(sp[i]==0)
{
if(dist[poi[i]][x]>dist[x][y]+(double)v[i]/sped[l])
{
if(!inq[poi[i]][x])
{
tx[++r]=poi[i];ty[r]=x;
sped[r]=sped[l];
inq[poi[i]][x]=1;
}
dist[poi[i]][x]=dist[x][y]+(double)v[i]/sped[l];
ans[poi[i]][x]=y;
}
}
else
{
if(dist[poi[i]][x]>dist[x][y]+(double)v[i]/sp[i])
{
if(!inq[poi[i]][x])
{
tx[++r]=poi[i];ty[r]=x;
sped[r]=sp[i];
inq[poi[i]][x]=1;
}
dist[poi[i]][x]=dist[x][y]+(double)v[i]/sp[i];
ans[poi[i]][x]=y;
}
}
}
l++;
}
int x=0;
For(i,1,n-1) if(dist[d][i]<dist[d][x]) x=i;
dfs(d,x);
Dow(i,1,tot) printf("%d ",Ans[i]);
}