洛谷P1266——速度限制

速度限制

好恶心啊……调的累死了

写法比较鬼畜，不过应该是跑的最快的做法吧。。。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]);
}