結果
問題 | No.899 γatheree |
ユーザー | Haar |
提出日時 | 2020-04-06 09:12:42 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 699 ms / 2,000 ms |
コード長 | 7,465 bytes |
コンパイル時間 | 2,219 ms |
コンパイル使用メモリ | 225,316 KB |
実行使用メモリ | 19,456 KB |
最終ジャッジ日時 | 2024-07-05 19:01:37 |
合計ジャッジ時間 | 16,131 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 2 ms
5,376 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 699 ms
19,328 KB |
testcase_07 | AC | 637 ms
19,328 KB |
testcase_08 | AC | 633 ms
19,328 KB |
testcase_09 | AC | 635 ms
19,328 KB |
testcase_10 | AC | 643 ms
19,456 KB |
testcase_11 | AC | 645 ms
19,328 KB |
testcase_12 | AC | 643 ms
19,200 KB |
testcase_13 | AC | 643 ms
19,328 KB |
testcase_14 | AC | 630 ms
19,328 KB |
testcase_15 | AC | 645 ms
19,200 KB |
testcase_16 | AC | 646 ms
19,456 KB |
testcase_17 | AC | 643 ms
19,328 KB |
testcase_18 | AC | 643 ms
19,200 KB |
testcase_19 | AC | 650 ms
19,456 KB |
testcase_20 | AC | 679 ms
19,328 KB |
testcase_21 | AC | 567 ms
19,296 KB |
testcase_22 | AC | 550 ms
19,416 KB |
testcase_23 | AC | 547 ms
19,320 KB |
ソースコード
#include <bits/stdc++.h> #ifdef DEBUG #include <Mylib/Debug/debug.cpp> #else #define dump(...) #endif template <typename Cost = int> class Edge{ public: int from,to; Cost cost; Edge() {} Edge(int to, Cost cost): to(to), cost(cost){} Edge(int from, int to, Cost cost): from(from), to(to), cost(cost){} Edge rev() const {return Edge(to,from,cost);} friend std::ostream& operator<<(std::ostream &os, const Edge &e){ os << "(FROM: " << e.from << "," << "TO: " << e.to << "," << "COST: " << e.cost << ")"; return os; } }; template <typename T> using Graph = std::vector<std::vector<Edge<T>>>; template <typename T> using Tree = std::vector<std::vector<Edge<T>>>; template <typename C, typename T> void add_edge(C &g, int from, int to, T w){ g[from].emplace_back(from, to, w); } template <typename C, typename T> void add_undirected(C &g, int a, int b, T w){ add_edge(g, a, b, w); add_edge(g, b, a, w); } template <typename Monoid> class LazySegmentTree{ using value_type_get = typename Monoid::value_type_get; using value_type_update = typename Monoid::value_type_update; using monoid_get = typename Monoid::monoid_get; using monoid_update = typename Monoid::monoid_update; const int depth, size, hsize; std::vector<value_type_get> data; std::vector<value_type_update> lazy; inline void propagate(int i){ if(lazy[i] == monoid_update::id()) return; if(i < hsize){ lazy[i << 1 | 0] = monoid_update::op(lazy[i], lazy[i << 1 | 0]); lazy[i << 1 | 1] = monoid_update::op(lazy[i], lazy[i << 1 | 1]); } int len = hsize >> (31 - __builtin_clz(i)); data[i] = Monoid::op(data[i], lazy[i], len); lazy[i] = monoid_update::id(); } inline value_type_get update_aux(int i, int l, int r, int s, int t, const value_type_update &x){ propagate(i); if(r <= s || t <= l) return data[i]; else if(s <= l && r <= t){ lazy[i] = monoid_update::op(x, lazy[i]); propagate(i); return data[i]; } else return data[i] = monoid_get::op(update_aux(i << 1 | 0, l, (l+r) / 2, s, t, x), update_aux(i << 1 | 1, (l+r) / 2, r, s, t, x)); } inline value_type_get get_aux(int i, int l, int r, int x, int y){ propagate(i); if(r <= x || y <= l) return monoid_get::id(); else if(x <= l && r <= y) return data[i]; else return monoid_get::op(get_aux(i << 1 | 0, l, (l+r) / 2, x, y), get_aux(i << 1 | 1, (l+r) / 2, r, x, y)); } public: LazySegmentTree(){} LazySegmentTree(int n): depth(n > 1 ? 32-__builtin_clz(n-1) + 1 : 1), size(1 << depth), hsize(size / 2), data(size, monoid_get::id()), lazy(size, monoid_update::id()) {} inline void update(int l, int r, const value_type_update &x){update_aux(1, 0, hsize, l, r, x);} inline void update_at(int i, const value_type_update &x){update(i, i+1, x);} inline value_type_get get(int l, int r){return get_aux(1, 0, hsize, l, r);} inline value_type_get at(int i){return get(i, i+1);} template <typename T> inline void init(const T &val){ init_with_vector(std::vector<T>(hsize, val)); } template <typename T> inline void init_with_vector(const std::vector<T> &val){ data.assign(size, monoid_get::id()); lazy.assign(size, monoid_update::id()); for(int i = 0; i < (int)val.size(); ++i) data[hsize + i] = val[i]; for(int i = hsize-1; i > 0; --i) data[i] = monoid_get::op(data[i << 1 | 0], data[i << 1 | 1]); } }; template <typename T> struct SumMonoid{ using value_type = T; constexpr inline static value_type id(){return 0;} constexpr inline static value_type op(const value_type &a, const value_type &b){return a + b;} }; template <typename T> struct UpdateMonoid{ using value_type = std::optional<T>; constexpr inline static value_type id(){return std::nullopt;} constexpr inline static value_type op(const value_type &a, const value_type &b){return (a ? a : b);} }; template <typename T, typename U> struct UpdateSum{ using monoid_get = SumMonoid<T>; using monoid_update = UpdateMonoid<U>; using value_type_get = typename monoid_get::value_type; using value_type_update = typename monoid_update::value_type; constexpr inline static value_type_get op(const value_type_get &a, const value_type_update &b, int len){ return b ? *b * len : a; } }; template <typename T> struct EulerTourBFS{ int N; std::vector<int> parent; std::vector<int> depth; std::vector<std::vector<int>> bfs_order; std::vector<std::vector<int>> dfs_order; std::vector<int> left, right; EulerTourBFS(const Tree<T> &tree, int root): N(tree.size()), parent(N), depth(N), left(N), right(N) { { int ord = 0; dfs(tree, root, -1, 0, ord); } { std::queue<std::pair<int,int>> q; q.emplace(root, 0); int ord = 0; while(not q.empty()){ auto [i, d] = q.front(); q.pop(); if((int)bfs_order.size() <= d) bfs_order.push_back(std::vector<int>()); bfs_order[d].push_back(ord); ++ord; for(auto &e : tree[i]){ if(e.to == parent[i]) continue; q.emplace(e.to, d + 1); } } } } void dfs(const Tree<T> &tree, int cur, int par, int d, int &ord){ parent[cur] = par; depth[cur] = d; if((int)dfs_order.size() <= d) dfs_order.push_back(std::vector<int>()); dfs_order[d].push_back(ord); left[cur] = ord; ++ord; for(auto &e : tree[cur]){ if(e.to == par) continue; dfs(tree, e.to, cur, d + 1, ord); } right[cur] = ord; } public: template <typename Func> inline void query_children(int i, int d, const Func &f) const { if(i != -1){ d += depth[i]; if((int)bfs_order.size() > d){ int l = std::lower_bound(dfs_order[d].begin(), dfs_order[d].end(), left[i]) - dfs_order[d].begin(); int r = std::lower_bound(dfs_order[d].begin(), dfs_order[d].end(), right[i]) - dfs_order[d].begin(); if(l >= (int)bfs_order[d].size()) return; if(r == l) return; f(bfs_order[d][l], bfs_order[d][r-1]+1); } } } template <typename Func> inline void query_at(int i, const Func &f) const { query_children(i, 0, f); } inline int get_parent(int i) const { if(i == -1) return -1; return parent[i]; } }; int main(){ std::cin.tie(0); std::ios::sync_with_stdio(false); int N; std::cin >> N; Tree<int> tree(N); for(int i = 0; i < N-1; ++i){ int u, v; std::cin >> u >> v; add_undirected(tree, u, v, 1); } auto res = EulerTourBFS<int>(tree, 0); std::vector<int64_t> A(N); for(int i = 0; i < N; ++i) std::cin >> A[i]; LazySegmentTree<UpdateSum<int64_t, int64_t>> seg(N); for(int i = 0; i < N; ++i){ res.query_at(i, [&](int l, int r){seg.update(l, r, A[i]);}); } int Q; std::cin >> Q; while(Q--){ int x; std::cin >> x; int64_t ans = 0; auto f = [&](int l, int r){ ans += seg.get(l, r); seg.update(l, r, 0); }; // 親の親 res.query_at(res.get_parent(res.get_parent(x)), f); // 親 res.query_at(res.get_parent(x), f); // 親の子 res.query_children(res.get_parent(x), 1, f); // 自分 res.query_at(x, f); // 子 res.query_children(x, 1, f); // 子の子 res.query_children(x, 2, f); res.query_at(x, [&](int l, int r){seg.update(l, r, ans);}); std::cout << ans << std::endl; } return 0; }