結果

問題 No.1442 I-wate Shortest Path Problem
ユーザー chocorusk
提出日時 2021-03-26 23:11:14
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,439 ms / 3,000 ms
コード長 3,516 bytes
コンパイル時間 5,245 ms
コンパイル使用メモリ 195,316 KB
最終ジャッジ日時 2025-01-19 23:30:38
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#include <utility>
#include <functional>
#include <time.h>
#include <stack>
#include <array>
#include <list>
#include <atcoder/all>
#define popcount __builtin_popcount
using namespace std;
using namespace atcoder;
typedef long long ll;
typedef pair<ll, int> P;
struct LCA{
    vector<vector<int>> g;
    vector<int> d;
    vector<vector<int>> p;
    int log;
    int n;
    LCA(const vector<vector<int>> &g):n(g.size()), g(g), d(g.size()){
        log=0;
        while(1<<log<=n) log++;
        p.resize(log, vector<int>(n));
    }
    void dfs(int x, int prev){
        for(auto y:g[x]){
            if(y==prev) continue;
            d[y]=d[x]+1;
            p[0][y]=x;
            dfs(y, x);
        }
    }
    void build(){
        dfs(0, -1);
        for(int i=1; i<log; i++){
            for(int j=0; j<n; j++){
                p[i][j]=p[i-1][p[i-1][j]];
            }
        }
    }
    int lca(int a, int b){
        if(d[a]>d[b]) swap(a, b);
        int dd=d[b]-d[a], i=0;
        int a1=a, b1=b;
        while(dd){
            if(dd&1) b1=p[i][b1];
            dd>>=1;
            i++;
        }
        if(a1==b1) return a1;
        for(int j=log-1; j>=0; j--){
            if(p[j][a1]!=p[j][b1]){
                a1=p[j][a1], b1=p[j][b1];
            }
        }
        return p[0][a1];
    }
    int dist(int a, int b){
        return d[a]+d[b]-2*d[lca(a, b)];
    }
};
int a[100010], b[100010];
ll c[100010];
vector<P> g[100010];
vector<P> g1[100020];
ll d1[10][100020];
ll d[100010];
void dfs(int x, int p){
    for(auto q:g[x]){
        int y=q.first;
        if(y==p) continue;
        d[y]=d[x]+q.second;
        dfs(y, x);
    }
}
int main()
{
    int n, k;
    cin>>n>>k;
    vector<vector<int>> gt(n);
    for(int i=0; i<n-1; i++){
        cin>>a[i]>>b[i]>>c[i];
        a[i]--; b[i]--;
        gt[a[i]].push_back(b[i]);
        gt[b[i]].push_back(a[i]);
        g[a[i]].push_back({b[i], c[i]*2});
        g[b[i]].push_back({a[i], c[i]*2});
        g1[a[i]].push_back({b[i], c[i]*2});
        g1[b[i]].push_back({a[i], c[i]*2});
    }
    for(int i=0; i<k; i++){
        int m; ll p; cin>>m>>p;
        for(int j=0; j<m; j++){
            int x; cin>>x; x--;
            g1[x].push_back({i+n, p});
            g1[i+n].push_back({x, p});
        }
    }
    const ll INF=1e18;
    for(int i=0; i<k; i++){
        priority_queue<P, vector<P>, greater<P>> que;
        fill(d1[i], d1[i]+n+k, INF);
        d1[i][i+n]=0;
        que.push({0, i+n});
        while(!que.empty()){
            auto p=que.top(); que.pop();
            int x=p.second;
            if(d1[i][x]<p.first) continue;
            for(auto q:g1[x]){
                int y=q.first;
                if(d1[i][y]>d1[i][x]+q.second){
                    d1[i][y]=d1[i][x]+q.second;
                    que.push({d1[i][y], y});
                }
            }
        }
    }
    LCA lca(gt);
    lca.build();
    dfs(0, -1);
    int q; cin>>q;
    while(q--){
        int u, v; cin>>u>>v; u--; v--;
        ll ans=d[u]+d[v]-2*d[lca.lca(u, v)];
        for(int i=0; i<k; i++){
            ans=min(ans, d1[i][u]+d1[i][v]);
        }
        ans/=2;
        cout<<ans<<endl;
    }
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0