【bzoj5047】空间传送装置

Description

太空中一共有n座星球，它们之间可以通过空间传送装置进行转移。空间传送装置分为m种，第i种装置可以用4个参
数a_i,b_i,c_i,d_i来描述。因为时空抖动的问题，在非整数时刻禁止使用空间传送装置。如果在整数s时刻使用装
置，那么需要花费((a_i*s+b_i) mod c_i)+d_i单位时间才能完成传送。现在是s时刻，小Q位于1号星球，请写一个
程序计算从1号星球到每个星球最少需要的时间。
Input

第一行包含4个正整数n,m,s,e(2<=n<=100000,1<=m<=50,1<=s<=2000,1<=e<=200000)
分别表示星球的个数、空间传送装置的种类数、当前的时间以及空间传送装置的个数。
接下来m行，每行4个正整数a_i,b_i,c_i,d_i(1<=a_i,b_i,c_i,d_i<=2000)，依次描述每种装置的参数。
接下来e行，每行3个正整数u_i,v_i,w_i(1<=u_i,v_i<=n,u_i!=v_i,1<=w_i<=m)
表示从星球u_i可以使用第w_i种装置单向传送到星球v_i。
Output

输出n-1行，每行一个整数，第i行表示从1到i+1的最少所需时间，若无解输出-1。
Sample Input

3 2 1 3

1 1 5 1

2 2 7 1

1 2 1

2 3 2

3 1 1
Sample Output

3

6

HINT

1到3：在时刻1使用第一种装置从1传送到2，花费时间3，再等待2单位时间，于时刻6使用第二种装置到达3，花费时间1。

题解
记录f[i][j]表示装置i，在j%c_i到达时刻所花费的最小时间，然后跑Dijkstra

代码

#include<bits/stdc++.h>
#define mod 1000000007
#define inf 1000000000
#define N 100005
#define pa pair<int,int>
typedef long long ll;
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,s,e,tot,a[N],b[N],c[N],d[N];
int Head[N],ret[2*N],Next[2*N],opt[2*N];
int f[55][2005],dis[N];
int t[2005],l[2005],r[2005];
bool vis[N];
priority_queue<pa,vector<pa>,greater<pa> >q;
inline void ins(int u,int v,int w)
{
    ret[++tot]=v;opt[tot]=w;
    Next[tot]=Head[u];Head[u]=tot;
}
void dijkstra()
{
    for (int i=1;i<=n;i++) dis[i]=inf;
    q.push(make_pair(s,1));dis[1]=s;
    while (!q.empty())
    {
        int now=q.top().second;q.pop();
        if (vis[now]) continue; vis[now]=1;
        for (int i=Head[now];i;i=Next[i])
            if (dis[now]+f[opt[i]][dis[now]%c[opt[i]]]<dis[ret[i]])
            {
                dis[ret[i]]=dis[now]+f[opt[i]][dis[now]%c[opt[i]]];
                q.push(make_pair(dis[ret[i]],ret[i]));
            }
    }
}
int main()
{
    n=read();m=read();s=read();e=read();
    for (int i=1;i<=m;i++)
    {
        a[i]=read(),b[i]=read(),c[i]=read(),d[i]=read();
        for (int j=0;j<c[i];j++) t[j]=(a[i]*j+b[i])%c[i]+d[i];
        l[0]=t[0];r[c[i]-1]=t[c[i]-1]+c[i]-1;
        for (int j=1;j<c[i];j++) l[j]=min(l[j-1],t[j]+j);
        for (int j=c[i]-2;j>=0;j--) r[j]=min(r[j+1],t[j]+j);
        for (int j=0;j<c[i];j++) f[i][j]=min(l[j]+c[i]-j,r[j]-j);
    }
    for (int i=1;i<=e;i++)
    {
        int u=read(),v=read(),w=read();
        ins(u,v,w);
    }
    dijkstra();
    for (int i=2;i<=n;i++)
        printf("%d\n",dis[i]==inf?-1:dis[i]-s);
    return 0;
}