Submission #709777

#TimeUsernameProblemLanguageResultExecution timeMemory
709777Jarif_RahmanEscape Route (JOI21_escape_route)C++17
70 / 100
9072 ms296720 KiB
#include "escape_route.h"
#include <bits/stdc++.h>
#define pb push_back
#define f first
#define sc second
using namespace std;
typedef long long int ll;
typedef string str;

const ll inf = 5e17;
const int N = 90;
vector<tuple<int, ll, ll>> v[N];
bitset<N> bl;
priority_queue<pair<ll, int>> pq;

ll dis2[N][N];

vector<ll> dis1[N][N];
vector<ll> sth[N];
ll change[N];


vector<ll> calculate_necessary_time(int n, int m, ll S, int q, vector<int> A, vector<int> B,
    vector<ll> L, vector<ll> C, vector<int> U, vector<int> V, vector<ll> T){

    vector<ll> ans(q, inf);
    
    for(int i = 0; i < m; i++){
        v[A[i]].pb({B[i], L[i], C[i]});
        v[B[i]].pb({A[i], L[i], C[i]});
    }

    for(int s = 0; s < n; s++){
        fill(dis2[s], dis2[s]+n, inf);
        for(int i = 0; i < n; i++) bl[i] = 0;

        dis2[s][s] = 0;
        pq.push({0, s});

        while(!pq.empty()){
            int nd = pq.top().sc; pq.pop();
            if(bl[nd]) continue;
            bl[nd] = 1;

            for(auto [x, l, c]: v[nd])
                if((dis2[s][nd]%S+l<=c ? dis2[s][nd]+l:((dis2[s][nd]+S-1)/S)*S+l) < dis2[s][x]){
                dis2[s][x] = (dis2[s][nd]%S+l<=c ? dis2[s][nd]+l:((dis2[s][nd]+S-1)/S)*S+l);
                pq.push({-dis2[s][x], x});
            }
        }
    }

    for(int s = 0; s < n; s++){
        ll t = 0;
        while(1){
            ll nxt = inf;
            sth[s].pb(t);
            for(int i = 0; i < n; i++) dis1[s][i].pb(inf);
            for(int i = 0; i < n; i++) bl[i] = 0;
            fill(change, change+n, inf);

            dis1[s][s].back() = t;
            pq.push({0, s});

            while(!pq.empty()){
                int nd = pq.top().sc; pq.pop();
                if(bl[nd]) continue;
                bl[nd] = 1;

                for(auto [x, l, c]: v[nd]) if(dis1[s][nd].back()+l <= min(dis1[s][x].back()-1, c)){
                    dis1[s][x].back() = dis1[s][nd].back()+l;
                    change[x] = c-dis1[s][x].back()+1;
                    pq.push({-dis1[s][x].back(), x});
                }
            }

            for(int i = 0; i < n; i++) if(dis1[s][i].back() != inf) dis1[s][i].back()-=t;
            nxt = min(nxt, *min_element(change, change+n));

            if(nxt == inf) break;
            t+=nxt;
        }
    }

    for(int i = 0; i < q; i++){
        int p = (upper_bound(sth[U[i]].begin(), sth[U[i]].end(), T[i])-sth[U[i]].begin())-1;
        if(dis1[U[i]][V[i]][p] != inf){
            ans[i] = min(ans[i], dis1[U[i]][V[i]][p]);
            continue;
        }
        for(int j = 0; j < n; j++) if(dis1[U[i]][j][p] != inf)
            ans[i] = min(ans[i], S-T[i]+dis2[j][V[i]]);
    }

    return ans;
}
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...