#pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #ifndef HIDDEN_IN_VS // 折りたたみ用 // 警告の抑制 #define _CRT_SECURE_NO_WARNINGS // ライブラリの読み込み #include using namespace std; // 型名の短縮 using ll = long long; // -2^63 ~ 2^63 = 9 * 10^18(int は -2^31 ~ 2^31 = 2 * 10^9) using pii = pair; using pll = pair; using pil = pair; using pli = pair; using vi = vector; using vvi = vector; using vvvi = vector; using vl = vector; using vvl = vector; using vvvl = vector; using vb = vector; using vvb = vector; using vvvb = vector; using vc = vector; using vvc = vector; using vvvc = vector; using vd = vector; using vvd = vector; using vvvd = vector; template using priority_queue_rev = priority_queue, greater>; using Graph = vvi; // 定数の定義 const double PI = acos(-1); const vi DX = { 1, 0, -1, 0 }; // 4 近傍(下,右,上,左) const vi DY = { 0, 1, 0, -1 }; int INF = 1001001001; ll INFL = 4004004004004004004LL; double EPS = 1e-12; // 入出力高速化 struct fast_io { fast_io() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(18); } } fastIOtmp; // 汎用マクロの定義 #define all(a) (a).begin(), (a).end() #define sz(x) ((int)(x).size()) #define lbpos(a, x) (int)distance((a).begin(), std::lower_bound(all(a), x)) #define ubpos(a, x) (int)distance((a).begin(), std::upper_bound(all(a), x)) #define Yes(b) {cout << ((b) ? "Yes\n" : "No\n");} #define rep(i, n) for(int i = 0, i##_len = int(n); i < i##_len; ++i) // 0 から n-1 まで昇順 #define repi(i, s, t) for(int i = int(s), i##_end = int(t); i <= i##_end; ++i) // s から t まで昇順 #define repir(i, s, t) for(int i = int(s), i##_end = int(t); i >= i##_end; --i) // s から t まで降順 #define repe(v, a) for(const auto& v : (a)) // a の全要素(変更不可能) #define repea(v, a) for(auto& v : (a)) // a の全要素(変更可能) #define repb(set, d) for(int set = 0; set < (1 << int(d)); ++set) // d ビット全探索(昇順) #define repp(a) sort(all(a)); for(bool a##_perm = true; a##_perm; a##_perm = next_permutation(all(a))) // a の順列全て(昇順) #define smod(n, m) ((((n) % (m)) + (m)) % (m)) // 非負mod #define uniq(a) {sort(all(a)); (a).erase(unique(all(a)), (a).end());} // 重複除去 #define EXIT(a) {cout << (a) << endl; exit(0);} // 強制終了 // 汎用関数の定義 template inline ll pow(T n, int k) { ll v = 1; rep(i, k) v *= n; return v; } template inline bool chmax(T& M, const T& x) { if (M < x) { M = x; return true; } return false; } // 最大値を更新(更新されたら true を返す) template inline bool chmin(T& m, const T& x) { if (m > x) { m = x; return true; } return false; } // 最小値を更新(更新されたら true を返す) // 演算子オーバーロード template inline istream& operator>>(istream& is, pair& p) { is >> p.first >> p.second; return is; } template inline istream& operator>>(istream& is, vector& v) { repea(x, v) is >> x; return is; } template inline vector& operator--(vector& v) { repea(x, v) --x; return v; } template inline vector& operator++(vector& v) { repea(x, v) ++x; return v; } // 手元環境(Visual Studio) #ifdef _MSC_VER #include "local.hpp" // 提出用(gcc) #else inline int popcount(int n) { return __builtin_popcount(n); } inline int popcount(ll n) { return __builtin_popcountll(n); } inline int lsb(int n) { return n != 0 ? __builtin_ctz(n) : -1; } inline int lsb(ll n) { return n != 0 ? __builtin_ctzll(n) : -1; } inline int msb(int n) { return n != 0 ? (31 - __builtin_clz(n)) : -1; } inline int msb(ll n) { return n != 0 ? (63 - __builtin_clzll(n)) : -1; } #define gcd __gcd #define dump(...) #define dumpel(v) #define dump_list(v) #define dump_list2D(v) #define input_from_file(f) #define output_to_file(f) #define Assert(b) { if (!(b)) while (1) cout << "OLE"; } #endif #endif // 折りたたみ用 //--------------AtCoder 専用-------------- #include using namespace atcoder; //using mint = modint1000000007; using mint = modint998244353; //using mint = modint; // mint::set_mod(m); istream& operator>>(istream& is, mint& x) { ll x_; is >> x_; x = x_; return is; } ostream& operator<<(ostream& os, const mint& x) { os << x.val(); return os; } using vm = vector; using vvm = vector; using vvvm = vector; //---------------------------------------- //【幅優先探索(動的)】O(|V| + |E|) /* * st から到達可能な t までの最短距離を dist[t] に格納したリストを返す. * nxt(s) は s の次に訪れることのできる頂点のリストを返す. */ template map dynamic_BFS(T st, const function(T)>& nxt) { // verify : https://atcoder.jp/contests/abc241/tasks/abc241_f map dist; // st からの最短距離を保持するテーブル dist[st] = 0; queue que; // 次に探索する頂点を入れておくキュー que.push(st); while (!que.empty()) { // 未探索の頂点 s を得る. auto s = que.front(); que.pop(); repe(t, nxt(s)) { // t が発見済みの頂点なら何もしない. if (dist.count(t)) continue; // スタートからの最短距離を確定する. dist[t] = dist[s] + 1; // 未探索の頂点として t を追加する. que.push(t); } } return dist; /* nxt の定義の雛形 using T = ll; function(T)> nxt = [&](T s) { vector res; return res; }; */ } void zikken() { int n = 6; using T = pair; function(T)> nxt = [&](T s) { vector res; string t, u; tie(t, u) = s; if (sz(t) > n) return res; res.push_back({ t + 'a', u }); res.push_back({ t + 'b', u }); res.push_back({ t, t }); res.push_back({ t + u, t }); return res; }; auto dist = dynamic_BFS({ "", "" }, nxt); map res; repe(tmp, dist) { if (sz(tmp.first.first) > n) continue; if (res.count(tmp.first.first)) { chmin(res[tmp.first.first], tmp.second); } else { res[tmp.first.first] = tmp.second; } } dumpel(res); exit(0); } //【ローリングハッシュ(列)】 /* * 列 s[0..n) の連続部分列 s[l, r) のハッシュ値を計算する. * * Rolling_hash(STR s) : O(n) * 列 s[0..n) で初期化する. * STR は string,vector など. * * ll get(int l, int r) : O(1) * 連続部分列 s[l, r) のハッシュ値を返す(空なら 0) * * ll join(ll hs, ll ht, int ls) : O(1) * ハッシュ値 hs をもつ s[0..ls) とハッシュ値 ht をもつ t を連結した s+t のハッシュ値を返す. */ template struct Rolling_hash_sub { using mint = static_modint; using vm = vector; mint B = BASE; // 適当な基数 mint S = SHIFT; // 適当なシフト // 列の長さ int n; // powB[i] : B^i, powB_inv[i] : B^(-i) vm powB, powB_inv; // v[i] : s[0, i) のハッシュ値 vm v; // コンストラクタ(列 s で初期化) Rolling_hash_sub() : n(0) {} Rolling_hash_sub(const STR& s) : n(sz(s)), powB(n + 1), powB_inv(n + 1), v(n + 1) { // ハッシュ値計算用の B の累乗 powB[0] = 1; rep(i, n) powB[i + 1] = powB[i] * B; // ハッシュ値計算用の B の逆元の累乗 mint invB = B.inv(); powB_inv[0] = 1; rep(i, n) powB_inv[i + 1] = powB_inv[i] * invB; // s[0, i) のハッシュ値 v[i] の計算 rep(i, n) v[i + 1] = v[i] + (s[i] + S) * powB[i]; } // 代入 Rolling_hash_sub(const Rolling_hash_sub& rh) = default; Rolling_hash_sub& operator=(const Rolling_hash_sub& rh) = default; // s[l, r) のハッシュ値の取得 int get(int l, int r) { // ハッシュ値は Σi=[0..r-l) (s[l+i] + S) * B^i (mod MOD) return ((v[r] - v[l]) * powB_inv[l]).val(); } // ハッシュ値 hs をもつ s[0..ls) とハッシュ値 ht をもつ t を連結した s+t のハッシュ値を返す. int join(int hs, int ht, int ls) { return (hs + ht * powB[ls]).val(); } }; template class Rolling_hash { int n; // 列の長さ // 衝突の可能性を減らすため,二つのハッシュ値を統合する. Rolling_hash_sub rh1; Rolling_hash_sub rh2; public: // コンストラクタ(文字列 s で初期化) Rolling_hash() : n(0) {} Rolling_hash(const STR& s) : n(sz(s)), rh1(s), rh2(s) {} // 代入 Rolling_hash(const Rolling_hash& rh) = default; Rolling_hash& operator=(const Rolling_hash& rh) = default; // 列の長さの取得 int size() { return n; } // s[l, r) のハッシュ値の取得 ll get(int l, int r) { // verify : https://atcoder.jp/contests/abc175/tasks/abc175_f if (r <= 0 || l >= n || l >= r) return 0; chmax(l, 0); chmin(r, n); return (ll(rh1.get(l, r)) << 32) + ll(rh2.get(l, r)); } // ハッシュ値 hs をもつ s[0..ls) とハッシュ値 ht をもつ t を連結した s+t のハッシュ値を返す. ll join(ll hs, ll ht, int ls) { // verify : https://atcoder.jp/contests/arc050/tasks/arc050_d int hs1 = (int)(hs >> 32), hs2 = (int)(hs % (1LL << 32)); int ht1 = (int)(ht >> 32), ht2 = (int)(ht % (1LL << 32)); return (ll(rh1.join(hs1, ht1, ls)) << 32) + ll(rh2.join(hs2, ht2, ls)); } }; int naive(string s) { int n = sz(s); Rolling_hash S(s); INF = 99; vvi dp(n + 1, vi(n, INF)); dp[0][0] = 0; rep(i, n) { rep(j, n) { if (S.get(i, i + j) == S.get(0, j)) { chmin(dp[i + j][j], dp[i][j] + 1); } chmin(dp[i + 1][j], dp[i][j] + 1); chmin(dp[i][i], dp[i][j] + 1); } } dumpel(dp); int res = INF; rep(j, n) chmin(res, dp[n][j]); return res; } /* babbababa 0: 0 99 99 99 99 99 99 99 99 1: 1 2 99 99 99 99 99 99 99 2: 2 3 3 99 99 99 99 99 99 3: 3 4 4 4 99 99 99 99 99 4: 4 5 5 5 5 99 99 99 99 5: 5 6 5 6 6 6 99 99 99 6: 6 7 6 5 7 7 6 99 99 7: 7 8 6 6 8 8 7 7 99 8: 8 9 7 7 9 9 8 8 8 9: 9 10 7 8 10 10 9 9 9 7 */ int TLE(string s) { int n = sz(s); Rolling_hash S(s); vi diag(n, INF); diag[0] = 0; int res = INF; rep(j, n) { vi dp(n + 1 - j, INF); dp[j - j] = diag[j]; repi(i, j, n - 1) { if (S.get(i, i + j) == S.get(0, j)) { chmin(dp[i + j - j], dp[i - j] + 1); } chmin(dp[i + 1 - j], dp[i - j] + 1); chmin(diag[i], dp[i - j] + 1); } chmin(res, dp[n - j]); dump(dp); } return res; } int main() { input_from_file("input.txt"); // output_to_file("output.txt"); // zikken(); string s; cin >> s; dump(naive(s)); cout << TLE(s) << endl; }