AtCoder Regular Contest 055 C: ABCAC
解法
ABCACの5つの区間に分けるとき、末尾のA+Cの部分の長さ s を全探索する。
1 つめの A の始点と、2つめの C の終点は決まっているが、A+Cの長さを定めることで、2つめの A の始点と 1 つめの C の終点が同時に求まる。Aの2つの始点からの Longest Common Prefix を求め、Cの2つの終点からの Longest Common Suffix を求める。これで、A+C の範囲内で A と C の分け方を数えることができる。
コード
#include <bits/stdc++.h> using namespace std; class RollingHash { private: vector<long long> mod; void init() { // mod.push_back(999999937LL); mod.push_back(1000000009LL); } const long long base = 1000000007LL; int n; vector<vector<long long>> hs, pw; RollingHash() {} public: int cnt = 0; RollingHash(const string &s) : n(s.size()) { init(); hs.resize(mod.size()), pw.resize(mod.size()); for (int i = 0; i < mod.size(); i++) { hs[i].assign(n + 1, 0); pw[i].assign(n + 1, 0); hs[i][0] = 0; pw[i][0] = 1; for (int j = 0; j < n; j++) { // pw[i][j + 1] = pw[i][j] * base % mod[i]; // hs[i][j + 1] = (hs[i][j] * base + s[j]) % mod[i]; pw[i][j + 1] = pw[i][j] * base; hs[i][j + 1] = (hs[i][j] * base + s[j]); } } } inline long long hash(int l, int r, int i) { // return ((hs[i][r] - hs[i][l] * pw[i][r - l]) % mod[i] + mod[i]) % mod[i]; return hs[i][r] - hs[i][l] * pw[i][r - l]; } inline bool match(int l1, int r1, int l2, int r2) { bool ret = 1; for (int i = 0; i < mod.size(); i++) ret &= hash(l1, r1, i) == hash(l2, r2, i); return ret; } inline bool match(int l1, int l2, int k) { return match(l1, l1 + k, l2, l2 + k); } int lcp(int i, int j) { int l = 0, r = min(n - i, n - j) + 1; while (l + 1 < r) { int m = (l + r) / 2; if (match(i, i + m, j, j + m)) l = m; else r = m; } return l; } }; int main() { cin.tie(0); ios::sync_with_stdio(false); string S; cin >> S; RollingHash rolling(S); reverse(S.begin(), S.end()); RollingHash invroll(S); int N = S.size(); int64_t ans = 0; for (int suffix = 2; suffix * 2 < N; ++suffix) { int lcp = rolling.lcp(0, N - suffix); int inv = invroll.lcp(0, suffix); lcp = min(lcp, suffix - 1); inv = min(inv, suffix - 1); if (lcp == 0 || inv == 0) continue; if (lcp + inv < suffix) continue; ans += lcp - (suffix - inv) + 1; } cout << ans << endl; }