#387080 (C++17) No.899 γatheree

提出ソース

結果

問題	No.899 γatheree
ユーザー	👑 emthrm
提出日時	2019-10-05 02:50:42
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	532 ms / 2,000 ms
コード長	5,618 bytes
コンパイル時間	1,836 ms
コンパイル使用メモリ	138,928 KB
実行使用メモリ	16,768 KB
最終ジャッジ日時	2024-10-04 08:50:19
合計ジャッジ時間	12,666 ms
ジャッジサーバーID （参考情報）	judge1 / judge5

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	AC	2 ms 5,248 KB
testcase_02	AC	2 ms 5,248 KB
testcase_03	AC	2 ms 5,248 KB
testcase_04	AC	2 ms 5,248 KB
testcase_05	AC	3 ms 5,248 KB
testcase_06	AC	507 ms 16,384 KB
testcase_07	AC	479 ms 16,512 KB
testcase_08	AC	524 ms 16,384 KB
testcase_09	AC	501 ms 16,396 KB
testcase_10	AC	513 ms 16,384 KB
testcase_11	AC	512 ms 16,384 KB
testcase_12	AC	516 ms 16,512 KB
testcase_13	AC	518 ms 16,512 KB
testcase_14	AC	532 ms 16,512 KB
testcase_15	AC	519 ms 16,512 KB
testcase_16	AC	521 ms 16,512 KB
testcase_17	AC	518 ms 16,512 KB
testcase_18	AC	512 ms 16,512 KB
testcase_19	AC	512 ms 16,512 KB
testcase_20	AC	495 ms 16,512 KB
testcase_21	AC	486 ms 16,768 KB
testcase_22	AC	495 ms 16,768 KB
testcase_23	AC	495 ms 16,768 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <algorithm>
#include <bitset>
#include <cassert>
#include <cctype>
#include <chrono>
#define _USE_MATH_DEFINES
#include <cmath>
#include <cstdio>
#include <cstring>
#include <ctime>
#include <deque>
#include <functional>
#include <iostream>
#include <iterator>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <sstream>
#include <stack>
#include <string>
#include <tuple>
#include <utility>
#include <vector>
using namespace std;

#define FOR(i,m,n) for(int i=(m);i<(n);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(v) (v).begin(),(v).end()

const int INF = 0x3f3f3f3f;
const long long LINF = 0x3f3f3f3f3f3f3f3fLL;
const double EPS = 1e-8;
const int MOD = 1000000007; // 998244353;
const int dy[] = {1, 0, -1, 0}, dx[] = {0, -1, 0, 1};
/*-------------------------------------------------*/
template <typename Monoid>
struct RSQandRUQ {
  RSQandRUQ(int sz, const Monoid &UNITY = 0) : UNITY(UNITY) {
    init(sz);
    dat.assign((n << 1) - 1, UNITY);
  }

  RSQandRUQ(const vector<Monoid> &a, const Monoid &UNITY = 0) : UNITY(UNITY) {
    int a_sz = a.size();
    init(a_sz);
    dat.resize((n << 1) - 1);
    REP(i, a_sz) dat[n - 1 + i] = a[i];
    for (int i = n - 2; i >= 0; --i) dat[i] = dat[(i << 1) + 1] + dat[(i << 1) + 2];
  }

  void update(int a, int b, const Monoid &value) { update(a, b, value, 0, 0, n); }

  Monoid sum(int a, int b) { return sum(a, b, 0, 0, n); }

  Monoid operator[](const int idx) { return sum(idx, idx + 1); }

  int find(int a, int b, const Monoid &value) { return find(a, b, value, 0, 0, n); }

private:
  int n = 1;
  const Monoid UNITY;
  vector<Monoid> dat, lazy;
  vector<bool> need_to_be_eval;

  void init(int sz) {
    while (n < sz) n <<= 1;
    lazy.assign((n << 1) - 1, UNITY);
    need_to_be_eval.assign((n << 1) - 1, false);
  }

  inline void evaluate(int node, int left, int right) {
    if (need_to_be_eval[node]) {
      dat[node] = (right - left) * lazy[node];
      if (node < n - 1) {
        lazy[(node << 1) + 1] = lazy[(node << 1) + 2] = lazy[node];
        need_to_be_eval[(node << 1) + 1] = need_to_be_eval[(node << 1) + 2] = true;
      }
      lazy[node] = UNITY;
      need_to_be_eval[node] = false;
    }
  }

  void update(int a, int b, const Monoid &value, int node, int left, int right) {
    evaluate(node, left, right);
    if (right <= a || b <= left) return;
    if (a <= left && right <= b) {
      lazy[node] = value;
      need_to_be_eval[node] = true;
      evaluate(node, left, right);
    } else {
      update(a, b, value, (node << 1) + 1, left, (left + right) >> 1);
      update(a, b, value, (node << 1) + 2, (left + right) >> 1, right);
      dat[node] = dat[(node << 1) + 1] + dat[(node << 1) + 2];
    }
  }

  Monoid sum(int a, int b, int node, int left, int right) {
    evaluate(node, left, right);
    if (right <= a || b <= left) return UNITY;
    if (a <= left && right <= b) return dat[node];
    return sum(a, b, (node << 1) + 1, left, (left + right) >> 1) + sum(a, b, (node << 1) + 2, (left + right) >> 1, right);
  }

  int find(int a, int b, const Monoid &value, int node, int left, int right) {
    evaluate(node, left, right);
    if (dat[node] < value || right <= a || b <= left) return -1;
    if (right - left == 1) return node - (n - 1);
    int res_l = find(a, b, value, (node << 1) + 1, left, (left + right) >> 1);
    if (res_l != -1) return res_l;
    return find(a, b, value, (node << 1) + 2, (left + right) >> 1, right);
  }
};

int main() {
  cin.tie(nullptr); ios::sync_with_stdio(false);
  // freopen("input.txt", "r", stdin);

  int n; cin >> n;
  vector<vector<int> > graph(n);
  REP(_, n - 1) {
    int u, v; cin >> u >> v;
    graph[u].emplace_back(v);
    graph[v].emplace_back(u);
  }
  vector<int> idx(n, -1), par(n, -1), left(n, INF), right(n, -1);
  int now = 0;
  queue<int> que;
  idx[0] = now++;
  que.emplace(0);
  while (!que.empty()) {
    int ver = que.front(); que.pop();
    for (int e : graph[ver]) {
      if (idx[e] == -1) {
        idx[e] = now++;
        que.emplace(e);
        par[idx[e]] = idx[ver];
        left[idx[ver]] = min(left[idx[ver]], idx[e]);
        right[idx[ver]] = max(right[idx[ver]], idx[e]);
      }
    }
  }
  vector<int> child_left(n, INF), child_right(n, -1);
  REP(i, n) {
    for (int e : graph[i]) if (idx[e] != par[idx[i]]) {
      child_left[idx[i]] = min(child_left[idx[i]], left[idx[e]]);
      child_right[idx[i]] = max(child_right[idx[i]], right[idx[e]]);
    }
  }
  vector<int> a(n); REP(i, n) cin >> a[i];
  vector<long long> fairy(n);
  REP(i, n) fairy[idx[i]] = a[i];
  RSQandRUQ<long long> rsq(fairy);
  int q; cin >> q;
  while (q--) {
    int x; cin >> x;
    x = idx[x];
    long long ans = rsq[x];
    rsq.update(x, x + 1, 0);
    if (par[x] != -1) {
      ans += rsq[par[x]];
      rsq.update(par[x], par[x] + 1, 0);
      ans += rsq.sum(left[par[x]], right[par[x]] + 1);
      rsq.update(left[par[x]], right[par[x]] + 1, 0);
      if (par[par[x]] != -1) {
        ans += rsq[par[par[x]]];
        rsq.update(par[par[x]], par[par[x]] + 1, 0);
      }
    }
    if (left[x] <= right[x]) {
      ans += rsq.sum(left[x], right[x] + 1);
      rsq.update(left[x], right[x] + 1, 0);
      if (child_left[x] <= child_right[x]) {
        ans += rsq.sum(child_left[x], child_right[x] + 1);
        rsq.update(child_left[x], child_right[x] + 1, 0);
      }
    }
    cout << ans << '\n';
    rsq.update(x, x + 1, ans);
  }
  return 0;
}

yukicoder

結果

テストケース

ソースコード