#780647 (C++17) No.2020 Sum of Common Prefix Length

提出ソース

結果

問題	No.2020 Sum of Common Prefix Length
ユーザー	👑 emthrm
提出日時	2022-07-22 22:32:58
言語	C++17 (gcc 13.2.0 + boost 1.83.0)
結果	AC
実行時間	204 ms / 2,000 ms
コード長	7,696 bytes
コンパイル時間	2,523 ms
コンパイル使用メモリ	218,396 KB
実行使用メモリ	82,132 KB
最終ジャッジ日時	2023-09-17 11:02:08
合計ジャッジ時間	8,291 ms
ジャッジサーバーID （参考情報）	judge14 / judge12

このコードへのチャレンジ（β）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 4,376 KB
testcase_01	AC	2 ms 4,376 KB
testcase_02	AC	1 ms 4,380 KB
testcase_03	AC	3 ms 4,380 KB
testcase_04	AC	3 ms 4,376 KB
testcase_05	AC	2 ms 4,376 KB
testcase_06	AC	3 ms 4,376 KB
testcase_07	AC	48 ms 5,252 KB
testcase_08	AC	48 ms 5,208 KB
testcase_09	AC	48 ms 5,252 KB
testcase_10	AC	48 ms 5,224 KB
testcase_11	AC	50 ms 5,224 KB
testcase_12	AC	49 ms 5,304 KB
testcase_13	AC	48 ms 5,204 KB
testcase_14	AC	50 ms 5,280 KB
testcase_15	AC	60 ms 11,160 KB
testcase_16	AC	60 ms 11,120 KB
testcase_17	AC	117 ms 36,792 KB
testcase_18	AC	97 ms 25,376 KB
testcase_19	AC	101 ms 27,596 KB
testcase_20	AC	163 ms 75,928 KB
testcase_21	AC	190 ms 67,484 KB
testcase_22	AC	199 ms 68,656 KB
testcase_23	AC	192 ms 66,072 KB
testcase_24	AC	204 ms 73,528 KB
testcase_25	AC	201 ms 72,800 KB
testcase_26	AC	99 ms 50,340 KB
testcase_27	AC	1 ms 4,380 KB
testcase_28	AC	52 ms 7,228 KB
testcase_29	AC	54 ms 7,432 KB
testcase_30	AC	52 ms 6,996 KB
testcase_31	AC	144 ms 74,712 KB
testcase_32	AC	149 ms 82,132 KB
testcase_33	AC	106 ms 50,344 KB
testcase_34	AC	76 ms 22,472 KB
testcase_35	AC	73 ms 22,100 KB
testcase_36	AC	111 ms 49,160 KB
testcase_37	AC	90 ms 32,836 KB

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

ソースコード

raw source code

#define _USE_MATH_DEFINES
#include <bits/stdc++.h>
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()
using ll = long long;
constexpr int INF = 0x3f3f3f3f;
constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL;
constexpr double EPS = 1e-8;
constexpr int MOD = 1000000007;
// constexpr int MOD = 998244353;
constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1};
constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1};
constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1};
template <typename T, typename U>
inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; }
template <typename T, typename U>
inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; }
struct IOSetup {
  IOSetup() {
    std::cin.tie(nullptr);
    std::ios_base::sync_with_stdio(false);
    std::cout << fixed << setprecision(20);
  }
} iosetup;

struct HeavyLightDecomposition {
  std::vector<int> parent, subtree, id, inv, head;

  explicit HeavyLightDecomposition(const std::vector<std::vector<int>>& graph,
                                   const int root = 0)
      : graph(graph) {
    const int n = graph.size();
    parent.assign(n, -1);
    subtree.assign(n, 1);
    dfs1(root);
    id.resize(n);
    inv.resize(n);
    head.assign(n, root);
    int cur_id = 0;
    dfs2(root, &cur_id);
  }

  template <typename Fn>
  void update_v(int u, int v, const Fn f) const {
    while (true) {
      if (id[u] > id[v]) std::swap(u, v);
      f(std::max(id[head[v]], id[u]), id[v] + 1);
      if (head[u] == head[v]) break;
      v = parent[head[v]];
    }
  }

  template <typename F, typename G, typename T>
  T query_v(int u, int v, const F f, const G g, const T id_t) const {
    T left = id_t, right = id_t;
    while (true) {
      if (id[u] > id[v]) {
        std::swap(u, v);
        std::swap(left, right);
      }
      left = g(left, f(std::max(id[head[v]], id[u]), id[v] + 1));
      if (head[u] == head[v]) break;
      v = parent[head[v]];
    }
    return g(left, right);
  }

  template <typename Fn>
  void update_subtree_v(const int ver, const Fn f) const {
    f(id[ver], id[ver] + subtree[ver]);
  }

  template <typename T, typename Fn>
  T query_subtree_v(const int ver, const Fn f) const {
    return f(id[ver], id[ver] + subtree[ver]);
  }

  template <typename Fn>
  void update_e(int u, int v, const Fn f) const {
    while (true) {
      if (id[u] > id[v]) std::swap(u, v);
      if (head[u] == head[v]) {
        f(id[u], id[v]);
        break;
      } else {
        f(id[head[v]] - 1, id[v]);
        v = parent[head[v]];
      }
    }
  }

  template <typename F, typename G, typename T>
  T query_e(int u, int v, const F f, const G g, const T id_t) const {
    T left = id_t, right = id_t;
    while (true) {
      if (id[u] > id[v]) {
        std::swap(u, v);
        std::swap(left, right);
      }
      if (head[u] == head[v]) {
        left = g(left, f(id[u], id[v]));
        break;
      } else {
        left = g(left, f(id[head[v]] - 1, id[v]));
        v = parent[head[v]];
      }
    }
    return g(left, right);
  }

  template <typename Fn>
  void update_subtree_e(const int ver, const Fn f) const {
    f(id[ver], id[ver] + subtree[ver] - 1);
  }

  template <typename T, typename Fn>
  T query_subtree_e(const int ver, const Fn f) const {
    return f(id[ver], id[ver] + subtree[ver] - 1);
  }

  int lowest_common_ancestor(int u, int v) const {
    while (true) {
      if (id[u] > id[v]) std::swap(u, v);
      if (head[u] == head[v]) break;
      v = parent[head[v]];
    }
    return u;
  }

 private:
  std::vector<std::vector<int>> graph;

  void dfs1(const int ver) {
    for (int i = 0; i < graph[ver].size(); ++i) {
      if (graph[ver][i] != parent[ver]) {
        parent[graph[ver][i]] = ver;
        dfs1(graph[ver][i]);
        subtree[ver] += subtree[graph[ver][i]];
        if (subtree[graph[ver][i]] > subtree[graph[ver].front()]) {
          std::swap(graph[ver][i], graph[ver].front());
        }
      }
    }
  }

  void dfs2(const int ver, int* cur_id) {
    id[ver] = (*cur_id)++;
    inv[id[ver]] = ver;
    for (const int e : graph[ver]) {
      if (e != parent[ver]) {
        head[e] = (e == graph[ver].front() ? head[ver] : e);
        dfs2(e, cur_id);
      }
    }
  }
};

template <typename Abelian>
struct FenwickTree {
  explicit FenwickTree(const int n, const Abelian ID = 0)
      : n(n), ID(ID), data(n, ID) {}

  void add(int idx, const Abelian val) {
    for (; idx < n; idx |= idx + 1) {
      data[idx] += val;
    }
  }

  Abelian sum(int idx) const {
    Abelian res = ID;
    for (--idx; idx >= 0; idx = (idx & (idx + 1)) - 1) {
      res += data[idx];
    }
    return res;
  }

  Abelian sum(const int left, const int right) const {
    return left < right ? sum(right) - sum(left) : ID;
  }

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

  int lower_bound(Abelian val) const {
    if (val <= ID) return 0;
    int res = 0, exponent = 1;
    while (exponent <= n) exponent <<= 1;
    for (int mask = exponent >> 1; mask > 0; mask >>= 1) {
      const int idx = res + mask - 1;
      if (idx < n && data[idx] < val) {
        val -= data[idx];
        res += mask;
      }
    }
    return res;
  }

 private:
  const int n;
  const Abelian ID;
  std::vector<Abelian> data;
};

struct Xor128 {
  int rand() {
    unsigned int t = x ^ (x << 11);
    x = y; y = z; z = w;
    return w = (w ^ (w >> 19)) ^ (t ^ (t >> 8));
  }
  int rand(const int ub) {
    const int res = rand() % ub;
    return res < 0 ? res + ub : res;
  }
  int rand(const int lb, const int ub) { return lb + rand(ub - lb); }
  long long randll() {
    return (static_cast<unsigned long long>(rand()) << 32) | rand();
  }
  long long randll(const long long ub) {
    const long long res = randll() % ub;
    return res < 0 ? res + ub : res;
  }
  long long randll(const long long lb, const long long ub) {
    return lb + randll(ub - lb);
  }
 private:
  unsigned int x = 123456789, y = 362436069, z = 521288629;
  unsigned int w = std::time(nullptr);
} xor128;

int main() {
  constexpr int C = 26;
  int n; cin >> n;
  vector<string> s(n); REP(i, n) cin >> s[i];
  int q; cin >> q;
  vector<int> x(q);
  vector<char> c(q, '$');
  REP(i, q) {
    int type; cin >> type >> x[i]; --x[i];
    if (type == 1) cin >> c[i];
  }
  int num = 0;
  REP(i, n) num += s[i].length();
  REP(i, q) num += c[i] != '$';
  vector<array<int, C>> trie(1);
  fill(ALL(trie.front()), -1);
  trie.reserve(num + 1);
  auto add = [&](int index, char c) -> int {
    if (trie[index][c - 'a'] == -1) {
      trie[index][c - 'a'] = trie.size();
      trie.emplace_back();
      fill(ALL(trie.back()), -1);
    }
    return trie[index][c - 'a'];
  };
  vector<int> pos(n, 0);
  REP(i, n) for (char c : s[i]) pos[i] = add(pos[i], c);
  REP(i, q) {
    if (c[i] != '$') pos[x[i]] = add(pos[x[i]], c[i]);
  }
  fill(ALL(pos), 0);
  const int m = trie.size();
  vector<vector<int>> graph(m);
  REP(i, m) REP(c, C) {
    if (trie[i][c] != -1) graph[i].emplace_back(trie[i][c]);
  }
  HeavyLightDecomposition hld(graph, 0);
  FenwickTree<ll> weight(m);
  REP(i, n) for (char c : s[i]) {
    pos[i] = trie[pos[i]][c - 'a'];
    weight.add(hld.id[pos[i]], 1);
  }
  REP(i, q) {
    if (c[i] == '$') {
      cout << hld.query_v(0, pos[x[i]],
                          [&weight](const int l, const int r) -> ll { return weight.sum(l, r); },
                          [](const ll a, const ll b) -> ll { return a + b; },
                          0LL)
           << '\n';
    } else {
      pos[x[i]] = trie[pos[x[i]]][c[i] - 'a'];
      weight.add(hld.id[pos[x[i]]], 1);
    }
  }
  return 0;
}

yukicoder

結果

テストケース

ソースコード