結果
問題 | No.2296 Union Path Query (Hard) |
ユーザー |
|
提出日時 | 2023-02-13 08:02:35 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,292 ms / 7,000 ms |
コード長 | 4,423 bytes |
コンパイル時間 | 2,139 ms |
コンパイル使用メモリ | 192,856 KB |
実行使用メモリ | 68,172 KB |
最終ジャッジ日時 | 2024-11-23 04:06:48 |
合計ジャッジ時間 | 35,075 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 45 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; int main(){ ios::sync_with_stdio(false); cin.tie(0); ll N, X, Q; cin >> N >> X >> Q; vector<vector<pair<int, ll>>> G(N); //重み付きグラフ vector<vector<int>> parent(18, vector<int>(N, -1)); //親の頂点のダブリング配列 vector<vector<ll>> dp(18, vector<ll>(N)); //辺の重みのダブリング配列 vector<int> depth(N); //根からの深さを入れる vector<int> parent_or_size(N, -1); //Union Find vector<tuple<ll, int, int>> diameter_and_vertex(N); //直径と頂点を入れる for(int i = 0; i < N; i++) diameter_and_vertex[i] = make_tuple(0, i, i); //leaderはUnion Find用の経路圧縮用関数 function<int(int)> leader = [&](int v){ if(parent_or_size[v] < 0) return v; return parent_or_size[v] = leader(parent_or_size[v]); }; //dfsでダブリングを更新する function<void(int)> dfs = [&](int v) { //ダブリングの更新ループ for(int i = 0; i + 1 < 18; i++){ if(parent[i][v] == -1) { parent[i + 1][v] = -1; dp[i + 1][v] = dp[i][v]; } else { parent[i + 1][v] = parent[i][parent[i][v]]; dp[i + 1][v] = dp[i][parent[i][v]] + dp[i][v]; } } int u; ll w; for(auto &&edge : G[v]){ tie(u, w) = edge; if(parent[0][v] == u) continue; depth[u] = depth[v] + 1; parent[0][u] = v; dp[0][u] = w; dfs(u); } }; //頂点u, 頂点v間の距離を求める関数 auto dist = [&](int u, int v){ if(leader(u) != leader(v)) return -1ll; if(depth[u] > depth[v]) swap(u, v); ll result = 0; //頂点 v を頂点 u と高さが同じになるようにする for(int i = 0; i < 18; i++){ if((depth[v] - depth[u]) >> i & 1){ result += dp[i][v]; v = parent[i][v]; } } if(u == v) return result; for(int i = 17; i >= 0; i--){ if(parent[i][u] != parent[i][v]){ result += dp[i][u]; result += dp[i][v]; u = parent[i][u]; v = parent[i][v]; } } result += dp[0][u]; result += dp[0][v]; return result; }; //頂点u, 頂点vを重みwの辺でマージ auto merge = [&](int u, int v, ll w){ int x = leader(u), y = leader(v); if(x == y) return true; assert(x != y); //非連結なことが保証されているはず if(-parent_or_size[x] < -parent_or_size[y]) { swap(x, y); swap(u, v); } //x - u 側の連結成分が親となるようにする //Union Findの更新 parent_or_size[x] += parent_or_size[y]; parent_or_size[y] = x; //辺を追加してDFSでダブリングを更新 G[u].emplace_back(v, w); G[v].emplace_back(u, w); depth[v] = depth[u] + 1; parent[0][v] = u; dp[0][v] = w; dfs(v); //直径を更新する tuple<ll, int, int> mx = diameter_and_vertex[x]; ll temp; array<int, 4> candidate_vertex{}; tie(temp, candidate_vertex[0], candidate_vertex[1]) = diameter_and_vertex[x]; tie(temp, candidate_vertex[2], candidate_vertex[3]) = diameter_and_vertex[y]; for(int i = 0; i < 4; i++) { int v1 = candidate_vertex[i]; for(int j = 0; j < i; j++) { int v2 = candidate_vertex[j]; ll length = dist(v1, v2); mx = max(mx, make_tuple(length, v1, v2)); } } diameter_and_vertex[x] = mx; return false; }; ll type, u, v, w; int cnt = 1; while(Q--){ cin >> type; if(type == 1){ cin >> v >> w; merge(v, X, w); }else if(type == 2){ cin >> u >> v; ll ans = dist(u, v); cout << ans << '\n'; if(ans != -1) ( X += ans ) %= N; }else if(type == 3){ cin >> v; cout << get<0>( diameter_and_vertex[ leader(v) ]) << '\n'; }else{ cin >> v; (X += v) %= N; } } }