結果

問題 No.515 典型LCP
ユーザー pekempeypekempey
提出日時 2017-09-29 15:00:10
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
(最新)
AC  
(最初)
実行時間 -
コード長 3,676 bytes
コンパイル時間 2,062 ms
コンパイル使用メモリ 183,440 KB
実行使用メモリ 66,752 KB
最終ジャッジ日時 2024-11-15 19:53:55
合計ジャッジ時間 14,251 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

inline unsigned lg(unsigned n) {
  asm("lzcnt %1,%0" : "=r"(n) : "r"(n));
  // why...?
  // return 31 - n;
  return n;
}

class SuffixDataStructure {
public:
  // a[i]>=0
  SuffixDataStructure(const std::vector<int> &a) : a(a) {
    sa = suffix_array(a);
    sa.insert(sa.begin(), a.size());
    build_lcp();
    build_rmq();
  }

  int get_lcp(int i, int j) {
    if (i == j) {
      return a.size() - std::max(i, j);
    }
    i = isa[i];
    j = isa[j];
    if (i > j) {
      std::swap(i, j);
    }
    int k = lg(j - i);
    return std::min(rmq[k][i], rmq[k][j - (1 << k)]);
  }

private:
  const std::vector<int> &a;
  std::vector<int> sa;
  std::vector<int> isa;
  std::vector<int> lcp;
  std::vector<std::vector<int>> rmq;

  vector<int> suffix_array(vector<int> a) {
    const int n = a.size();
    if (n == 0) return {};
    const int H = *max_element(begin(a), end(a)) + 1;
    vector<int> s(n), ss, next(n), rank(n);
    for (int i = n - 2; i >= 0; i--) {
      s[i] = a[i] == a[i + 1] ? s[i + 1] : a[i] < a[i + 1];
      if (!s[i] && s[i + 1]) {
        next[i + 1] = ss.empty() ? n : ss.back();
        ss.push_back(i + 1);
      }
    }
    auto mp = [&](vector<int> f, vector<int> a) { for (int &x : a) x = f[x]; return a; };
    auto induced_sort = [&]() {
      vector<int> sa(n, -1), L(H + 1);
      for (int x : a) L[x + 1]++;
      for (int i = 0; i < H; i++) L[i + 1] += L[i];
      auto S = L, SS = L;
      sa[L[a[n - 1]]++] = n - 1;
      for_each(ss.rbegin(), ss.rend(), [&](int i) { sa[--SS[a[i] + 1]] = i; });
      for_each(sa.begin(), sa.end(), [&](int i) { i--; if (i >= 0 && !s[i]) sa[L[a[i]]++] = i; });
      for_each(sa.rbegin(), sa.rend(), [&](int i) { i--; if (i >= 0 && s[i]) sa[--S[a[i] + 1]] = i; });
      return sa;
    };
    reverse(ss.begin(), ss.end());
    int j = -1;
    for (int i : induced_sort()) if (0 < i && s[i] && !s[i - 1]) {
      if (j != -1) rank[i] = rank[j] + lexicographical_compare(a.begin() + j, a.begin() + next[j], a.begin() + i, a.begin() + next[i]);
      j = i;
    }
    ss = mp(ss, suffix_array(mp(rank, ss)));
    return induced_sort();
  }

  void build_lcp() {
    const int n = a.size();
    isa.resize(n + 1);
    lcp.assign(n + 1, 0);
    for (int i = 0; i <= n; i++) {
      isa[sa[i]] = i;
    }
    int k = 0;
    for (int i = 0; i < n; i++) {
      int j = sa[isa[i] - 1];
      k = std::max(0, k - 1);
      while (i + k < n && j + k < n && a[i + k] == a[j + k]) {
        k++;
      }
      lcp[isa[i] - 1] = k;
    }
  }

  void build_rmq() {
    const int n = lcp.size();
    const int m = lg(n);
    rmq.assign(m + 1, std::vector<int>(n));
    for (int i = 0; i < n; i++) {
      rmq[0][i] = lcp[i];
    }
    for (int i = 0; i < m; i++) {
      for (int j = 0; j + (1 << i) < n; j++) {
        rmq[i + 1][j] = std::min(rmq[i][j], rmq[i][j + (1 << i)]);
      }
    }
  }
};

int main() {
  int n;
  cin >> n;

  std::vector<int> start(n);
  std::vector<int> t;
  std::vector<int> len(n);
  for (int i = 0; i < n; i++) {
    start[i] = t.size();
    static char buf[800001];
    scanf("%s", buf);
    int j;
    for (j = 0; buf[j] != 0; j++) {
      t.push_back(buf[j]);
    }
    len[i] = j;
  }

  SuffixDataStructure suf(t);

  int m;
  long long x, d;
  cin >> m >> x >> d;

  long long ans = 0;
  const long long N = 1LL * n * (n - 1);
  while (m--) {
    int i = x / (n - 1);
    int j = x % (n - 1);
    if (i > j) {
      std::swap(i, j);
    } else {
      j = j + 1;
    }
    x += d;
    if (x >= N) {
      x -= N;
    }
    ans += min(suf.get_lcp(start[i], start[j]), min(len[i], len[j]));
  }
  std::cout << ans << std::endl;
}
0