結果

問題 No.430 文字列検索
ユーザー yuppe19 😺yuppe19 😺
提出日時 2019-11-26 18:03:45
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 208 ms / 2,000 ms
コード長 2,110 bytes
コンパイル時間 1,000 ms
コンパイル使用メモリ 94,272 KB
実行使用メモリ 23,544 KB
最終ジャッジ日時 2024-04-27 03:00:19
合計ジャッジ時間 2,346 ms
ジャッジサーバーID
(参考情報)
judge1 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 208 ms
23,544 KB
testcase_02 AC 15 ms
6,940 KB
testcase_03 AC 16 ms
6,944 KB
testcase_04 AC 2 ms
6,944 KB
testcase_05 AC 2 ms
6,944 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 2 ms
6,944 KB
testcase_08 AC 201 ms
23,524 KB
testcase_09 AC 2 ms
6,940 KB
testcase_10 AC 9 ms
6,944 KB
testcase_11 AC 37 ms
6,948 KB
testcase_12 AC 39 ms
6,944 KB
testcase_13 AC 39 ms
6,944 KB
testcase_14 AC 30 ms
6,944 KB
testcase_15 AC 26 ms
6,944 KB
testcase_16 AC 17 ms
6,944 KB
testcase_17 AC 16 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <cassert>
#include <iostream>
#include <unordered_map>
#include <vector>
using namespace std;
using u64 = uint64_t;

u64 uz = time(NULL);
inline u64 xorshift64() {
  uz ^= uz << 18;
  uz ^= uz >> 19;
  uz ^= uz << 19;
  return uz;
}

class rolling_hash {
  static constexpr u64 M = (1ULL << 63) - (1ULL << 41) + 1,
                       M2 = M - 2,
                       R2 = static_cast<u64>(-__uint128_t(M) % M);
  static constexpr u64 MR(__uint128_t t) {
    u64 r = u64((__uint128_t(u64(t) * M2) * M + t) >> 64);
    return r < M ? r : r - M;
  }
  static constexpr u64 init(__uint128_t x) {
    return MR(x * R2);
  }
  static constexpr u64 mod_mul(__uint128_t x, u64 y) {
    return MR(x * y);
  }

  string s;
  size_t n;
  u64 b61;
  vector<u64> B, H;
 public:
  rolling_hash() {}
  explicit rolling_hash(const string &_s, u64 seed=xorshift64()) : s(_s), b61(seed % (M-128) + 127), n(s.size()) {
    B.resize(n+1);
    H.resize(n+1);
    B[0] = init(1);
    for(size_t i=0; i<n; ++i) {
      B[i+1] = mod_mul(B[i], b61);
      H[i+1] = mod_mul(H[i], b61) + s[i];
      if(H[i+1] >= M) { H[i+1] -= M; }
    }
  }
  u64 calc_hash(size_t l, size_t r) const {    // [l, r]. つまり l も r も含む。1-indexed.
    assert(l <= r);
    assert(l <= n && r <= n); // if(l > n || r > n) { return -1; }
    return calc_hash2(l, r-l+1);
  }
  u64 calc_hash2(size_t l, size_t len) const { // l から len 文字。1-indexed.
    assert(1 <= len);
    assert(l <= n && l+len <= n+1);
    u64 x = H[l+len-1], yz = mod_mul(B[len], H[l-1]);
    return x < yz ? M - yz + x : x - yz;
  }
};

int main(void) {
  cin.tie(nullptr); ios::sync_with_stdio(false);
  string s; int m; cin >> s >> m;
  u64 seed = xorshift64();
  rolling_hash rh(s, seed);
  unordered_map<u64, int> mp;
  for(size_t i=1, n=s.size(); i<=n; ++i) {
    for(size_t len=1; len<=10 && i+len<=n+1; ++len) {
      ++mp[rh.calc_hash2(i, len)];
    }
  }
  int res = 0;
  for(int i=0; i<m; ++i) {
    string ci; cin >> ci;
    rolling_hash rh2(ci, seed);
    res += mp[rh2.calc_hash(1, ci.size())];
  }
  cout << res << '\n';
  return 0;
}
0