#include #include #include using namespace std; class query { public: int tp, t, x, y; query() : tp(-1), t(-1), x(-1), y(-1) {} query(int tp_, int t_, int x_, int y_) : tp(tp_), t(t_), x(x_), y(y_) {} bool operator<(const query& q) const { return t < q.t; } }; vector solve(int N, int M, int Q, const vector& TA, const vector& TP, const vector& X, const vector& Y, const vector& S, const vector& I) { // step #1. preparation vector qs(M + Q); for (int i = 0; i < M; i++) { qs[i] = query(TP[i], i * 2 + 1, X[i], Y[i]); } for (int i = 0; i < Q; i++) { qs[M + i] = query(5, S[i] * 2, I[i], i); } sort(qs.begin(), qs.end()); // step #2. renumber vertices vector > group(N); vector root(N), id(N), idl(N), idr(N); for (int i = 0; i < N; i++) { group[i] = { i }; root[i] = i; id[i] = i; } for (query q : qs) { if (q.tp == 1) { int rx = root[q.x]; int ry = root[q.y]; if (rx != ry) { if (group[rx].size() < group[ry].size()) { swap(rx, ry); } group[rx].insert(group[rx].end(), group[ry].begin(), group[ry].end()); int delta = idr[rx] - idl[rx]; for (int j : group[ry]) { id[j] += delta; } idl[ry] += delta; idr[ry] += delta; for (int j : group[ry]) { root[j] = rx; } } } } for (query &q : qs) { if (q.tp == 1) { q.x = id[q.x]; q.y = id[q.y]; } else { q.x = id[q.x]; } } vector A(N); for (int i = 0; i < N; i++) { A[id[i]] = TA[i]; } // step #3. find left/right of update queries for (int i = 0; i < N; i++) { group[i] = { i }; root[i] = i; idl[i] = i; idr[i] = i + 1; } vector l(M + Q, -1), r(M + Q, -1); for (int i = 0; i < M + Q; i++) { query q = qs[i]; if (q.tp == 1) { int rx = root[q.x]; int ry = root[q.y]; if (rx != ry) { if (group[rx].size() < group[ry].size()) { swap(rx, ry); } group[rx].insert(group[rx].end(), group[ry].begin(), group[ry].end()); idl[rx] = min(idl[rx], idl[ry]); idr[rx] = max(idr[rx], idr[ry]); for (int j : group[ry]) { root[j] = rx; } } } if (q.tp == 4) { l[i] = idl[root[q.x]]; r[i] = idr[root[q.x]]; } } // step #4. find all high-degree vertices vector deg(N, 0); for (query q : qs) { if (q.tp == 1) { deg[q.x] += 1; deg[q.y] += 1; } } int Z = 0; vector vert_id(N, -1); for (int i = 0; i < N; i++) { if (1LL * deg[i] * deg[i] >= N) { vert_id[i] = Z; Z += 1; } } vector > flag(N, vector(Z, -1)); // step #5. binary indexed tree vector bit(N + 1); auto add = [&](int pos, int val) { for (int i = pos + 1; i <= N; i += i & (-i)) { bit[i] += val; } }; auto getsum = [&](int pos) { long long answer = 0; for (int i = pos; i >= 1; i -= i & (-i)) { answer += bit[i]; } return answer; }; // step #6. process queries vector > g(N); vector val(N, 0), answer(Q, -1); vector > bigval(M + Q + 1); bigval[0] = vector(Z, 0); for (int i = 0; i < M + Q; i++) { bigval[i + 1] = bigval[i]; query q = qs[i]; if (q.tp == 1) { if (vert_id[q.x] == -1) { g[q.x].push_back(q.y); } else if (flag[q.y][vert_id[q.x]] == -1) { flag[q.y][vert_id[q.x]] = i; } if (vert_id[q.y] == -1) { g[q.y].push_back(q.x); } else if (flag[q.x][vert_id[q.y]] == -1) { flag[q.x][vert_id[q.y]] = i; } } if (q.tp == 2) { val[q.x] += q.y; } if (q.tp == 3) { if (vert_id[q.x] == -1) { for (int j : g[q.x]) { val[j] += q.y; } } else { bigval[i + 1][vert_id[q.x]] += q.y; } val[q.x] += q.y; } if (q.tp == 4) { add(l[i], q.y); add(r[i], -q.y); } if (q.tp == 5) { long long subval = val[q.x]; for (int j = 0; j < Z; j++) { if (flag[q.x][j] != -1) { subval += bigval[i + 1][j] - bigval[flag[q.x][j]][j]; } } subval += getsum(q.x + 1); answer[q.y] = max(A[q.x] - subval, 0LL); } } return answer; } vector solve_easy(int N, int M, int Q, const vector& TA, const vector& TP, const vector& X, const vector& Y, const vector& S, const vector& I) { // step #1. preparation vector qs(M + Q); for (int i = 0; i < M; i++) { qs[i] = query(TP[i], i * 2 + 1, X[i], Y[i]); } for (int i = 0; i < Q; i++) { qs[M + i] = query(5, S[i] * 2, I[i], i); } sort(qs.begin(), qs.end()); // step #2. renumber vertices vector A = TA; vector answer(Q, -1); vector > g(N); for (query q : qs) { if (q.tp == 1) { g[q.x].push_back(q.y); g[q.y].push_back(q.x); } if (q.tp == 2) { A[q.x] = max(A[q.x] - q.y, 0LL); } if (q.tp == 3) { A[q.x] = max(A[q.x] - q.y, 0LL); for (int i : g[q.x]) { A[i] = max(A[i] - q.y, 0LL); } } if (q.tp == 4) { vector vis(N, false); auto dfs = [&](auto& self, int pos) -> void { vis[pos] = true; A[pos] = max(A[pos] - q.y, 0LL); for (int i : g[pos]) { if (!vis[i]) { self(self, i); } } }; dfs(dfs, q.x); } if (q.tp == 5) { answer[q.y] = A[q.x]; } } return answer; } int main() { cin.tie(0); ios::sync_with_stdio(false); int N; cin >> N; vector A(N); for (int i = 0; i < N; i++) { cin >> A[i]; } int M; cin >> M; vector TP(M), X(M), Y(M); for (int i = 0; i < M; i++) { cin >> TP[i] >> X[i] >> Y[i]; X[i] -= 1; Y[i] -= (TP[i] == 1 ? 1 : 0); } int Q; cin >> Q; vector S(Q), I(Q); for (int i = 0; i < Q; i++) { cin >> S[i] >> I[i]; I[i] -= 1; } vector answer = solve_easy(N, M, Q, A, TP, X, Y, S, I); for (int i = 0; i < Q; i++) { cout << answer[i] << '\n'; } return 0; }