P9751 [CSP-J 2023] 旅遊巴士

纯粹的發表於2024-08-11

原題連結

分析

很逆天的一道題

設 \(dp[i][j]\) 為到達第 \(i\) 個點的時刻 \(t\) 且滿足 \(t\mod k=j\) 的最小 \(t\)

則有答案為 \(dp[n][0]\)

更新也很簡單，設當前點為 \(u\)，當前時間為 \(t\) 需要遍歷的下一個點 \(v\)，則有 \(dis[v][(t+1)\%k]=dis[u][t\%k]+1\)

如果道路還沒開放，我們可以原地等 \(k\) 的倍數時間（等價於在起點等），因為我們只更新 \([(t+1)\% k]\)，其他的時間自然會有其他時刻來更新

這個設計非常巧妙，避開了 \(a_i\) 的限制

code

#include<bits/stdc++.h>
using namespace std;
/*
#define int long long
#define double long double
#define lowbit(x) ((x)&(-x))
const int inf=1e18;
const int mod=1e9+7;

const int N=4e5;
int qpow(int a,int n)
{
    int res=1;
    while(n)
    {
        if(n&1) res=res*a%mod;
        a=a*a%mod;
        n>>=1;
    }
    return res;
}
int inv(int x)
{
    return qpow(x,mod-2);
}
int fa[2000005];
int finds(int now){return now==fa[now]?now:finds(fa[now]);}

vector<int> G[200005];

int dfn[200005],low[200005];
int cnt=0,num=0;
int in_st[200005]={0};
stack<int> st;
int belong[200005]={0};

void scc(int now,int fa)
{
    dfn[now]=++cnt;
    low[now]=dfn[now];
    in_st[now]=1;
    st.push(now);

    for(auto next:G[now])
    {
        if(next==fa) continue;

        if(!dfn[next])
        {
            scc(next,now);
            low[now]=min(low[now],low[next]);
        }
        else if(in_st[next])
        {
            low[now]=min(low[now],dfn[next]);
        }
    }

    if(low[now]==dfn[now])
    {
        int x;
        num++;
        do
        {
            x=st.top();
            st.pop();
            in_st[x]=0;
            belong[x]=num;
        }while(x!=now);
    }
}
vector<int> prime;
bool mark[200005]={0};
void shai()
{
    for(int i=2;i<=200000;i++)
    {
        if(!mark[i]) prime.push_back(i);

        for(auto it:prime)
        {
            if(it*i>200000) break;

            mark[it*i]=1;
            if(it%i==0) break;
        }
    }
}
*/

struct node
{
    int to,cost;
};
vector<node> G[10004];

struct fresh
{
    int place,mod,t;
    bool operator<(const fresh&c)const
    {
        return c.t<t;
    }
};


int dis[10005][105];
void solve()
{
    memset(dis,0x3f,sizeof dis);
    int n,m,k;
    cin>>n>>m>>k;

    for(int i=1;i<=m;i++)
    {
        int x,y,t;
        cin>>x>>y>>t;

        G[x].push_back({y,t});
    }


    priority_queue<fresh> q;
    q.push({1,0,0});

    while(q.size())
    {
        auto [now,mod,t]=q.top();
        q.pop();
        if(dis[now][mod]<=t) continue;

        dis[now][mod]=t;
        for(auto next:G[now])
        {
            auto [to,wait]=next;

            if(wait<=t)
            {
                q.push({to,(mod+1)%k,t+1});
            }
            else
            {
                int add=(wait-t)/k+((wait-t)%k!=0);
                q.push({to,(mod+1)%k,t+add*k+1});
            }
        }
    }

    if(dis[n][0]>1e9) cout<<-1;
    else cout<<dis[n][0];
}
signed main()
{
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    int TT=1;
    //cin>>TT;
    while(TT--) solve();
    return 0;
}

相關文章