結果
問題 |
No.900 aδδitivee
|
ユーザー |
![]() |
提出日時 | 2020-01-03 00:58:18 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 984 ms / 2,000 ms |
コード長 | 4,001 bytes |
コンパイル時間 | 1,565 ms |
コンパイル使用メモリ | 165,148 KB |
実行使用メモリ | 43,648 KB |
最終ジャッジ日時 | 2024-11-22 18:50:03 |
合計ジャッジ時間 | 14,564 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 27 |
ソースコード
#include "bits/stdc++.h" using namespace std; #define int long long #define FOR(i, a, b) for(int i=(a);i<(b);i++) #define RFOR(i, a, b) for(int i=(b-1);i>=(a);i--) #define REP(i, n) for(int i=0; i<(n); i++) #define RREP(i, n) for(int i=(n-1); i>=0; i--) #define ALL(a) (a).begin(),(a).end() #define UNIQUE_SORT(l) sort(ALL(l)); l.erase(unique(ALL(l)), l.end()); #define CONTAIN(a, b) find(ALL(a), (b)) != (a).end() #define array2(type, x, y) array<array<type, y>, x> #define vector2(type) vector<vector<type> > #define out(...) printf(__VA_ARGS__) int dxy[] = {0, 1, 0, -1, 0}; /*================================*/ int N,Q; /********************************** * * Graph * **********************************/ struct Edge { int src, to, w; Edge(int u, int v, int w): src(u), to(v), w(w) {}; }; vector< vector<Edge> > E(111111); int W[111111]; int D[111111]; int P[111111][20]; // 平方分割: 各番号へ、追加の重み情報が入る int S[2222]; int T[222222]; // 深さに掛ける分 int SB[2222]; int TB[222222]; // オイラーツアー: int et = 0; int LS[222222]; // iに対するTの添字(スタート時) int RS[222222]; // iに対するTの添字(終了時) void doubling() { FOR(j, 1, 20) { REP(i, N) { P[i][j] = (P[i][j-1] == -1) ? -1 : P[P[i][j-1]][j-1]; } } } void add_edge(int u, int v, int w) { E[u].push_back(*new Edge(u,v,w)); E[v].push_back(*new Edge(v,u,w)); } void dfs(int parent=-1, int i=0, int weight=0, int depth=0) { P[i][0] = parent; W[i] = weight; D[i] = depth; LS[i] = et++; for (auto e : E[i]) { if (e.to == parent) continue; dfs(i, e.to, weight + e.w, depth+1); } RS[i] = et++; } int lca(int u, int v) { // depth調整 if (D[u] > D[v]) swap(u, v); RFOR(i, 0, 20) { if ((D[v] - D[u])&(1<<i)) v = P[v][i]; } if (u == v) return u; // lca RFOR(i, 0, 20) { if (P[v][i] != P[u][i]) { v = P[v][i]; u = P[u][i]; } } return P[u][0]; } int route_weight(int u, int v) { return W[u] + W[v] - W[lca(u, v)] * 2; } void setup_tree() { REP(i,N)REP(j,20) P[i][j]=-1; dfs(); doubling(); } void add_weight(int i, int w) { int l = LS[i] + 1; // 範囲内 int r = RS[i]; // 範囲外 int k = sqrt(N); while(l < r) { int lk = l / k; int rk = r / k; if (lk == rk) { FOR(j, l, r) { T[j] += w; } } else if (l % k == 0) { S[lk] += w; } else { FOR(j, l, (lk+1)*k) { T[j] += w; } } l = (lk+1)*k; } } void add_weightb(int i, int w) { int l = LS[i] + 1; // 範囲内 int r = RS[i]; // 範囲外 int k = sqrt(N); while(l < r) { int lk = l / k; int rk = r / k; if (lk == rk) { FOR(j, l, r) { TB[j] += w; } } else if (l % k == 0) { SB[lk] += w; } else { FOR(j, l, (lk+1)*k) { TB[j] += w; } } l = (lk+1)*k; } } int query(int i) { int k = sqrt(N); int l = LS[i]; return W[i] + (T[l] + S[l / k]) + (TB[l] + SB[l/k]) * D[i]; } /**********************************************************/ signed main() { #if DEBUG std::ifstream in("input.txt"); std::cin.rdbuf(in.rdbuf()); #endif cin>>N; int u,v,w; REP(i,N-1) { cin >> u >> v >> w; add_edge(u, v, w); } REP(i, N*2) T[i] = TB[i] = 0; REP(i, sqrt(N)*2+2) S[i] = SB[i] = 0; setup_tree(); cin >> Q; int c,a,b,x; REP(i,Q) { cin >> c; if (c == 1) { cin >> a >> x; add_weight(a, -x * D[a]); add_weightb(a, x); } else { cin >> b; out("%lld\n", query(b)); } } return 0; }