#755746 (C++17) No.1909 Detect from Substrings

提出ソース

結果

問題	No.1909 Detect from Substrings
ユーザー	keijak
提出日時	2022-04-23 10:19:04
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	MLE
実行時間	-
コード長	6,398 bytes
コンパイル時間	3,059 ms
コンパイル使用メモリ	228,056 KB
実行使用メモリ	813,916 KB
最終ジャッジ日時	2024-06-24 21:14:01
合計ジャッジ時間	6,070 ms
ジャッジサーバーID （参考情報）	judge2 / judge5

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 6,812 KB
testcase_01	AC	2 ms 6,940 KB
testcase_02	AC	2 ms 6,944 KB
testcase_03	AC	11 ms 6,944 KB
testcase_04	AC	12 ms 6,944 KB
testcase_05	AC	19 ms 9,216 KB
testcase_06	AC	28 ms 14,080 KB
testcase_07	AC	19 ms 12,800 KB
testcase_08	MLE	-
testcase_09	--	-
testcase_10	--	-
testcase_11	--	-
testcase_12	--	-
testcase_13	--	-
testcase_14	--	-
testcase_15	--	-
testcase_16	--	-
testcase_17	--	-
testcase_18	--	-
testcase_19	--	-
testcase_20	--	-
testcase_21	--	-
testcase_22	--	-
testcase_23	--	-
testcase_24	--	-
testcase_25	--	-
testcase_26	--	-
testcase_27	--	-
testcase_28	--	-
testcase_29	--	-
testcase_30	--	-
testcase_31	--	-
testcase_32	--	-
testcase_33	--	-
testcase_34	--	-
testcase_35	--	-
testcase_36	--	-
testcase_37	--	-
testcase_38	--	-

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

ソースコード

raw source code

#include <bits/stdc++.h>
#define REP_(i, a_, b_, a, b, ...) for (int i = (a), END_##i = (b); i < END_##i; ++i)
#define REP(i, ...) REP_(i, __VA_ARGS__, __VA_ARGS__, 0, __VA_ARGS__)
#define ALL(x) std::begin(x), std::end(x)
using Int = long long;
using Uint = unsigned long long;
using Real = long double;

template<typename T, typename U>
inline bool chmax(T &a, U b) { return a < b and ((a = b), true); }
template<typename T, typename U>
inline bool chmin(T &a, U b) { return a > b and ((a = b), true); }
template<typename T>
inline int ssize(const T &a) { return (int) a.size(); }
template<typename T>
constexpr T kBigVal = std::numeric_limits<T>::max() / 2;

struct CastInput {
  template<typename T>
  operator T() const {
    T x;
    std::cin >> x;
    return x;
  }
  struct Sized {
    int n;
    template<typename T>
    operator T() const {
      T xs(n);
      for (auto &x: xs) std::cin >> x;
      return xs;
    }
  };
  Sized operator()(int n) const { return {n}; }
} in;

template<typename Container>
std::ostream &print_seq(const Container &seq,
                        const char *sep = " ",
                        const char *ends = "\n",
                        std::ostream &os = std::cout) {
  const auto itl = std::begin(seq), itr = std::end(seq);
  for (auto it = itl; it != itr; ++it) {
    if (it != itl) os << sep;
    os << *it;
  }
  return os << ends;
}

template<typename T>
inline std::ostream &print_one(const T &x, char endc) {
  if constexpr (std::is_same<T, bool>::value) {
    return std::cout << (x ? "Yes" : "No") << endc;
  } else {
    return std::cout << x << endc;
  }
}
template<typename T>
inline std::ostream &print(const T &x) { return print_one(x, '\n'); }
template<typename T, typename... Ts>
std::ostream &print(const T &head, Ts... tail) {
  return print_one(head, ' '), print(tail...);
}
inline std::ostream &print() { return std::cout << '\n'; }
void exit_() { std::cout.flush(), std::cerr.flush(), std::_Exit(0); }

#ifdef MY_DEBUG
#include "debug_dump.hpp"
#include "backward.hpp"
backward::SignalHandling kSignalHandling;
#else
#define DUMP(...)
#define cerr if(false)cerr
#endif

using namespace std;

struct RollingHash {
  using u64 = std::uint64_t;
  using u128 = __uint128_t;
  static const u64 kMod = (1ULL << 61) - 1;

  static u64 base() {
    static const u64 kBase = []() {
      std::mt19937_64 engine(std::random_device{}());
      std::uniform_int_distribution<u64> rand(1000, kMod - 1);
      return rand(engine);
    }();
    return kBase;
  }

  static inline u64 add(u64 a, u64 b) {
    a += b;
    return (a >= kMod) ? (a - kMod) : a;
  }

  static inline u64 sub(u64 a, u64 b) { return add(a, kMod - b); }

  static inline u64 mul(u64 a, u64 b) {
    u128 t = u128(a) * b;
    u64 na = u64(t >> 61);
    u64 nb = u64(t & kMod);
    na += nb;
    return (na >= kMod) ? (na - kMod) : na;
  }

  static u64 pow_base(int i) {
    static std::vector<u64> kPowBase(1, 1);
    while (int(kPowBase.size()) <= i) {
      u64 val = mul(kPowBase.back(), base());
      kPowBase.push_back(val);
    }
    return kPowBase[i];
  }

  // Calculates hash(a || b) from hash(a), hash(b) and length(b).
  static u64 concat(u64 a_hash, u64 b_hash, int b_length) {
    return add(mul(a_hash, pow_base(b_length)), b_hash);
  }
};

// Computes hash for any substring in O(1).
struct SpanHash : public RollingHash {
  std::vector<u64> cum_hash;

  // Constructionn: O(n).
  // All elements must be non-zero. Otherwise we won't be able to distinguish
  // between [1] and [0, 1].
  template<class Sequence>
  explicit SpanHash(const Sequence &s) : cum_hash(s.size() + 1) {
    int i = 0;
    for (const auto &x: s) {
      u64 val = static_cast<u64>(x);
      assert(val != 0);  // Ensure all elements are non-zero!
      cum_hash[i + 1] = add(mul(cum_hash[i], base()), val);
      ++i;
    }
  }

  // Returns the hash value of the substring s[l:r]. O(1).
  u64 get(int l, int r) const {
    // Compute `hash(s[0:r]) - hash(s[0:l]) * B^(r-l) (mod M)`
    return sub(cum_hash.at(r), mul(cum_hash.at(l), pow_base(r - l)));
  }
};

// Binary search over integers
template<class T, class F>
auto bisect(T truthy, T falsy, F pred) -> T {
  static_assert(std::is_integral_v<T>);
  static_assert(std::is_invocable_r_v<bool, F, T>);
  while (std::max(truthy, falsy) - std::min(truthy, falsy) > T(1)) {
    auto mid = (truthy & falsy) + (truthy ^ falsy) / T(2);
    auto ok = pred(mid);
    (ok ? truthy : falsy) = std::move(mid);
  }
  return truthy;
}

Int solve() {
  int n = in, m = in;
  vector<string> S = in(n);
  vector<SpanHash> SH;
  if (m > 20) {
    for (const auto &s: S) {
      assert(ssize(s) == m);
      SH.emplace_back(s);
    }
  }

  auto sub_eq = [&](int i, int j, int l, int r) {
    if (m <= 20) {
      return S[i].substr(l, r - l) == S[j].substr(l, r - l);
    } else {
      return SH[i].get(l, r) == SH[j].get(l, r);
    }
  };

  std::set < string > candidates;
  for (int j = 1; j < 2; ++j) {
    for (int i = 0; i <= m; ++i) {
      if (sub_eq(0, j, 0, i)) {
        candidates.insert(S[j].substr(0, i) + S[0][i] + S[j].substr(i));
        candidates.insert(S[j].substr(0, i) + S[j][i] + S[0].substr(i));
      }
      if (sub_eq(0, j, m - i, m)) {
        candidates.insert(S[j].substr(0, m - i) + S[0][m - 1 - i] + S[j].substr(m - i));
        candidates.insert(S[j].substr(0, m - i) + S[j][m - 1 - i] + S[0].substr(m - i));
      }
    }
  }

  auto is_subseq = [&](int j, string_view T, SpanHash &th) -> bool {
    int plen = 0, slen = 0;
    if (m <= 20) {
      while (plen < m and S[j][plen] == T[plen]) ++plen;
      while (slen < m and S[j][m - 1 - slen] == T[m - slen]) ++slen;
    } else {
      plen = bisect(0, m + 1, [&](int i) {
        return SH[j].get(0, i) == th.get(0, i);
      });
      slen = bisect(0, m + 1, [&](int i) {
        return SH[j].get(m - i, m) == th.get(m + 1 - i, m + 1);
      });
    }
    return (plen + slen >= m);
  };

  DUMP(candidates);
  int count = 0;
  for (const string &T: candidates) {
    SpanHash th(T);
    bool ok = true;
    for (int j = 0; j < n; ++j) {
      if (not is_subseq(j, T, th)) {
        ok = false;
        break;
      }
    }
    if (ok) ++count;
  }
  return count;
}

int main() {
  std::ios::sync_with_stdio(false), cin.tie(nullptr);
  cout << std::fixed << std::setprecision(18);
  const int T = 1;//in;
  REP(t, T) { print(solve()); }
  exit_();
}

yukicoder

結果

テストケース

ソースコード