結果

問題 No.599 回文かい
ユーザー stoqstoq
提出日時 2021-10-26 06:50:37
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 2,737 bytes
コンパイル時間 2,519 ms
コンパイル使用メモリ 214,244 KB
実行使用メモリ 11,040 KB
最終ジャッジ日時 2024-04-15 04:03:55
合計ジャッジ時間 8,273 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 3 ms
6,940 KB
testcase_05 AC 7 ms
6,940 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 2 ms
6,940 KB
testcase_08 AC 2 ms
6,940 KB
testcase_09 AC 4 ms
6,940 KB
testcase_10 AC 13 ms
6,940 KB
testcase_11 AC 12 ms
6,940 KB
testcase_12 AC 37 ms
6,944 KB
testcase_13 AC 68 ms
6,944 KB
testcase_14 TLE -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
evil_0.txt -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#include <atcoder/fenwicktree>
#include <atcoder/modint>
using namespace std;
using namespace atcoder;

#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

using M1 = modint1000000007;
using M2 = modint998244353;
using hash_t = pair<M1, M2>;
namespace atcoder {
bool operator<(const M1 &t1, const M1 &t2) { return t1.val() < t2.val(); }
bool operator<(const M2 &t1, const M2 &t2) { return t1.val() < t2.val(); }
}  // namespace atcoder

struct Hash {
  int n;
  string s;
  long long B;
  bool calc_rev;
  fenwick_tree<M1> bit1, bit1_rev;
  fenwick_tree<M2> bit2, bit2_rev;
  vector<M1> Bp1;
  vector<M2> Bp2;

  Hash() {}
  Hash(string &s, bool calc_rev = true)
      : s(s), n(s.length()), calc_rev(calc_rev), Bp1(n + 1), Bp2(n + 1) {
    random_device seed;
    mt19937 engine(seed());
    B = (long long)(engine()) + 1000;
    bit1 = fenwick_tree<M1>(n);
    bit2 = fenwick_tree<M2>(n);
    M1 p1 = 1;
    M2 p2 = 1;
    Bp1[0] = 1, Bp2[0] = 1;
    for (int i = 0; i < n; i++) {
      bit1.add(i, p1 * s[i]);
      bit2.add(i, p2 * s[i]);
      p1 *= B, p2 *= B;
      Bp1[i + 1] = p1;
      Bp2[i + 1] = p2;
    }
    if (not calc_rev) return;
    bit1_rev = fenwick_tree<M1>(n);
    bit2_rev = fenwick_tree<M2>(n);
    for (int i = 0; i < n; i++) {
      bit1_rev.add(i, Bp1[n - i - 1] * s[i]);
      bit2_rev.add(i, Bp2[n - i - 1] * s[i]);
    }
  }

  void update(int i, char c) {
    bit1.add(i, Bp1[i] * c - bit1.sum(i, i + 1));
    bit2.add(i, Bp2[i] * c - bit2.sum(i, i + 1));
    if (not calc_rev) return;
    bit1_rev.add(i, Bp1[n - i - 1] * c - bit1_rev.sum(i, i + 1));
    bit2_rev.add(i, Bp2[n - i - 1] * c - bit2_rev.sum(i, i + 1));
  }

  hash_t query(int l, int r, bool rev = false) {
    assert(l <= r);
    if (rev) {
      assert(calc_rev);
      hash_t res = {bit1_rev.sum(l, r), bit2_rev.sum(l, r)};
      res.first /= Bp1[n - r];
      res.second /= Bp2[n - r];
      return res;
    }
    hash_t res = {bit1.sum(l, r), bit2.sum(l, r)};
    res.first /= Bp1[l];
    res.second /= Bp2[l];
    return res;
  }

  bool is_palindrome(int l, int r) {
    assert(calc_rev);
    return query(l, (l + r) / 2, false) == query((l + r + 1) / 2, r, true);
  }
};

int main() {
  string s;
  cin >> s;
  int n = s.length();
  Hash hash(s, false);
  using mint = M1;
  vector<mint> dp(n);
  vector<int> used(n);
  auto f = [&](auto self, int l, int r) -> mint {
    if (used[l]) return dp[l];
    mint res = 1;
    for (int l2 = l + 1, r2 = r - 1; l2 <= r2; l2++, r2--) {
      if (hash.query(l, l2) == hash.query(r2, r)) {
        res += self(self, l2, r2);
      }
    }
    used[l] = true;
    return dp[l] = res;
  };
  cout << f(f, 0, n).val() << "\n";
}
0