結果
問題 | No.901 K-ary εxtrεεmε |
ユーザー | kwm_t |
提出日時 | 2023-05-06 20:23:00 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 125 ms / 3,000 ms |
コード長 | 4,283 bytes |
コンパイル時間 | 2,964 ms |
コンパイル使用メモリ | 217,872 KB |
実行使用メモリ | 24,320 KB |
最終ジャッジ日時 | 2024-11-24 01:48:41 |
合計ジャッジ時間 | 8,328 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 73 ms
24,320 KB |
testcase_01 | AC | 2 ms
5,248 KB |
testcase_02 | AC | 3 ms
5,248 KB |
testcase_03 | AC | 3 ms
5,248 KB |
testcase_04 | AC | 3 ms
5,248 KB |
testcase_05 | AC | 3 ms
5,248 KB |
testcase_06 | AC | 3 ms
5,248 KB |
testcase_07 | AC | 125 ms
21,888 KB |
testcase_08 | AC | 118 ms
21,888 KB |
testcase_09 | AC | 119 ms
22,016 KB |
testcase_10 | AC | 119 ms
21,888 KB |
testcase_11 | AC | 118 ms
22,144 KB |
testcase_12 | AC | 122 ms
22,016 KB |
testcase_13 | AC | 121 ms
22,016 KB |
testcase_14 | AC | 123 ms
22,144 KB |
testcase_15 | AC | 117 ms
22,144 KB |
testcase_16 | AC | 119 ms
22,016 KB |
testcase_17 | AC | 124 ms
22,016 KB |
testcase_18 | AC | 117 ms
22,016 KB |
testcase_19 | AC | 121 ms
22,016 KB |
testcase_20 | AC | 119 ms
22,144 KB |
testcase_21 | AC | 120 ms
21,888 KB |
testcase_22 | AC | 112 ms
22,144 KB |
testcase_23 | AC | 112 ms
22,144 KB |
testcase_24 | AC | 112 ms
22,144 KB |
testcase_25 | AC | 113 ms
22,016 KB |
testcase_26 | AC | 115 ms
22,144 KB |
testcase_27 | AC | 94 ms
22,016 KB |
testcase_28 | AC | 96 ms
22,144 KB |
testcase_29 | AC | 97 ms
22,144 KB |
ソースコード
#include <bits/stdc++.h> //#include <atcoder/all> using namespace std; //using namespace atcoder; //using mint = modint1000000007; //const int mod = 1000000007; //using mint = modint998244353; //const int mod = 998244353; //const int INF = 1e9; //const long long LINF = 1e18; #define rep(i, n) for (int i = 0; i < (n); ++i) #define rep2(i,l,r)for(int i=(l);i<(r);++i) #define rrep(i, n) for (int i = (n-1); i >= 0; --i) #define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i) #define all(x) (x).begin(),(x).end() #define allR(x) (x).rbegin(),(x).rend() #define endl "\n" #define P pair<int,int> template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; } template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; } class Tree { public: Tree(int n, int root) : n(n), root(root) { edge.resize(n); for (int i = 0; i < MAXLOGV; i++) parent[i].resize(n); depth.resize(n); dist_.resize(n); idx.resize(n); cnt = 0; } // uとvをつなぐ // lcaを求めることが主目的なので無向グラフとしている void unite(int u, int v, long long w) { edge[u].push_back({ v, w }); edge[v].push_back({ u, w }); //edge[u].emplace_back(v, w); //edge[v].emplace_back(u, w); } // initする // コンストラクタだけじゃなくてこれも呼ばないとlcaが求められないぞ void init() { dfs(root, -1, 0); for (int k = 0; k + 1 < MAXLOGV; k++) { for (int v = 0; v < n; v++) { if (parent[k][v] < 0) parent[k + 1][v] = -1; else parent[k + 1][v] = parent[k][parent[k][v]]; } } } // uとvのlcaを求める int lca(int u, int v) const { if (depth[u] > depth[v]) swap(u, v); for (int k = 0; k < MAXLOGV; k++) { if ((depth[v] - depth[u]) >> k & 1) { v = parent[k][v]; } } if (u == v) return u; for (int k = MAXLOGV - 1; k >= 0; k--) { if (parent[k][u] != parent[k][v]) { u = parent[k][u]; v = parent[k][v]; } } return parent[0][u]; } // uのn個親を求める int pare(int v, int n) { if (depth[v] < n)return -1; n = min(n, depth[v]); int idx = MAXLOGV; while (n) { for (int i = idx - 1; i >= 0; --i) { if (n < (1 << i))continue; if (-1 == parent[i][v])continue; n -= (1 << i); v = parent[i][v]; idx = i; break; } } return v; } // uからvに向かってd進んだ頂点を返す int JumpOnTree(int u, int v, int d) { if (0 == d)return u; int distuv = dist(u, v); if (distuv < d)return -1; int l = lca(u, v); if (l == u)return pare(v, distuv - d); if (l == v)return pare(u, d); int distlu = dist(l, u); if (distlu >= d)return pare(u, d); return pare(v, distuv - d); } // uとvの距離を求める // edgeを定義しないといけない時はこれじゃダメ long long dist(int u, int v) const { int p = lca(u, v); return dist_[u] + dist_[v] - dist_[p] * 2; //return (depth[u] - depth[p]) + (depth[v] - depth[p]); } //頂点wが頂点u,vのパス上に存在するか bool on_path(int u, int v, int w) { return (dist(u, w) + dist(v, w) == dist(u, v)); } int dfs(int v, int p, int d) { idx[v] = cnt; cnt++; int ret = 1; parent[0][v] = p; depth[v] = d; for (auto e : edge[v]) { if (e.first == p) continue; dist_[e.first] = dist_[v] + e.second; auto get = dfs(e.first, v, d + 1); ret += get; } return ret; } static const int MAXLOGV = 25; // グラフの隣接リスト表現 vector<vector<pair<int, long long>>>edge; // 頂点の数 int n; // 根ノードの番号 int root; // 親ノード vector<int> parent[MAXLOGV]; // 根からの深さ vector<int> depth; vector<long long>dist_; vector<int>idx; int cnt; }; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int n; cin >> n; Tree tree(n, 0); for (int i = 0; i < n - 1; i++) { int x, y, w; cin >> x >> y >> w; tree.unite(x, y, w); } tree.init(); int q; cin >> q; while (q--) { int k; cin >> k; vector<P>v(k); rep(i, k) { int x; cin >> x; v[i] = { tree.idx[x],x }; } sort(all(v)); long long ans = 0; rep(i, k) { int v0 = v[i].second; int v1 = v[(i + 1) % k].second; ans += tree.dist(v0, v1); } cout << ans / 2 << endl; } return 0; }