結果
問題 | No.3000 Optimal Run Length Encoding |
ユーザー | 👑 rin204 |
提出日時 | 2024-12-25 20:30:03 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 21,928 bytes |
コンパイル時間 | 4,003 ms |
コンパイル使用メモリ | 282,588 KB |
実行使用メモリ | 110,300 KB |
最終ジャッジ日時 | 2024-12-25 20:31:55 |
合計ジャッジ時間 | 89,495 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 202 ms
5,248 KB |
testcase_02 | AC | 178 ms
5,248 KB |
testcase_03 | AC | 170 ms
5,248 KB |
testcase_04 | AC | 181 ms
5,248 KB |
testcase_05 | AC | 159 ms
5,248 KB |
testcase_06 | AC | 161 ms
5,248 KB |
testcase_07 | AC | 195 ms
5,248 KB |
testcase_08 | AC | 177 ms
5,248 KB |
testcase_09 | AC | 162 ms
5,248 KB |
testcase_10 | AC | 152 ms
5,248 KB |
testcase_11 | AC | 158 ms
5,248 KB |
testcase_12 | AC | 193 ms
5,248 KB |
testcase_13 | AC | 162 ms
5,248 KB |
testcase_14 | AC | 152 ms
5,248 KB |
testcase_15 | AC | 150 ms
5,248 KB |
testcase_16 | AC | 163 ms
5,248 KB |
testcase_17 | AC | 148 ms
5,248 KB |
testcase_18 | AC | 124 ms
5,248 KB |
testcase_19 | AC | 265 ms
5,248 KB |
testcase_20 | AC | 267 ms
5,248 KB |
testcase_21 | AC | 293 ms
5,248 KB |
testcase_22 | AC | 289 ms
5,248 KB |
testcase_23 | AC | 267 ms
5,248 KB |
testcase_24 | AC | 299 ms
5,248 KB |
testcase_25 | AC | 262 ms
5,248 KB |
testcase_26 | AC | 282 ms
5,248 KB |
testcase_27 | AC | 236 ms
5,248 KB |
testcase_28 | AC | 275 ms
5,248 KB |
testcase_29 | AC | 302 ms
5,248 KB |
testcase_30 | AC | 365 ms
5,248 KB |
testcase_31 | AC | 369 ms
5,248 KB |
testcase_32 | AC | 422 ms
5,248 KB |
testcase_33 | AC | 515 ms
5,248 KB |
testcase_34 | AC | 583 ms
9,636 KB |
testcase_35 | AC | 322 ms
22,156 KB |
testcase_36 | AC | 296 ms
22,116 KB |
testcase_37 | AC | 300 ms
22,308 KB |
testcase_38 | AC | 294 ms
22,164 KB |
testcase_39 | AC | 309 ms
22,124 KB |
testcase_40 | AC | 294 ms
22,240 KB |
testcase_41 | AC | 292 ms
22,248 KB |
testcase_42 | AC | 298 ms
22,196 KB |
testcase_43 | AC | 319 ms
22,252 KB |
testcase_44 | AC | 298 ms
22,116 KB |
testcase_45 | AC | 691 ms
47,108 KB |
testcase_46 | AC | 395 ms
24,224 KB |
testcase_47 | AC | 695 ms
43,200 KB |
testcase_48 | AC | 468 ms
30,116 KB |
testcase_49 | AC | 658 ms
38,348 KB |
testcase_50 | AC | 424 ms
28,716 KB |
testcase_51 | AC | 654 ms
43,112 KB |
testcase_52 | AC | 480 ms
30,528 KB |
testcase_53 | AC | 661 ms
43,188 KB |
testcase_54 | AC | 434 ms
27,888 KB |
testcase_55 | AC | 765 ms
47,108 KB |
testcase_56 | AC | 479 ms
33,500 KB |
testcase_57 | AC | 628 ms
37,996 KB |
testcase_58 | AC | 383 ms
24,192 KB |
testcase_59 | AC | 677 ms
43,096 KB |
testcase_60 | AC | 691 ms
47,236 KB |
testcase_61 | AC | 695 ms
44,128 KB |
testcase_62 | AC | 500 ms
24,584 KB |
testcase_63 | AC | 660 ms
39,436 KB |
testcase_64 | AC | 387 ms
24,040 KB |
testcase_65 | AC | 685 ms
47,468 KB |
testcase_66 | AC | 389 ms
23,988 KB |
testcase_67 | AC | 657 ms
39,664 KB |
testcase_68 | AC | 393 ms
22,812 KB |
testcase_69 | AC | 663 ms
43,812 KB |
testcase_70 | AC | 1,372 ms
110,300 KB |
testcase_71 | AC | 290 ms
22,240 KB |
testcase_72 | AC | 614 ms
35,160 KB |
testcase_73 | AC | 318 ms
23,812 KB |
testcase_74 | AC | 722 ms
35,972 KB |
testcase_75 | AC | 290 ms
22,244 KB |
testcase_76 | AC | 686 ms
35,604 KB |
testcase_77 | AC | 284 ms
22,712 KB |
testcase_78 | AC | 636 ms
36,532 KB |
testcase_79 | AC | 309 ms
23,424 KB |
testcase_80 | AC | 518 ms
35,776 KB |
testcase_81 | AC | 58 ms
5,248 KB |
testcase_82 | AC | 959 ms
67,836 KB |
testcase_83 | AC | 919 ms
67,460 KB |
testcase_84 | AC | 947 ms
67,364 KB |
testcase_85 | AC | 884 ms
66,592 KB |
testcase_86 | AC | 912 ms
63,888 KB |
testcase_87 | AC | 941 ms
71,484 KB |
testcase_88 | AC | 913 ms
64,916 KB |
testcase_89 | AC | 937 ms
71,356 KB |
testcase_90 | AC | 850 ms
61,932 KB |
testcase_91 | AC | 939 ms
63,176 KB |
testcase_92 | WA | - |
testcase_93 | WA | - |
testcase_94 | AC | 934 ms
66,936 KB |
testcase_95 | AC | 1,069 ms
67,048 KB |
testcase_96 | AC | 942 ms
66,960 KB |
testcase_97 | AC | 376 ms
25,928 KB |
testcase_98 | AC | 274 ms
21,224 KB |
testcase_99 | AC | 260 ms
19,356 KB |
testcase_100 | AC | 312 ms
23,568 KB |
testcase_101 | AC | 313 ms
24,864 KB |
testcase_102 | AC | 400 ms
26,804 KB |
testcase_103 | AC | 413 ms
28,224 KB |
testcase_104 | AC | 440 ms
30,084 KB |
testcase_105 | AC | 354 ms
27,232 KB |
testcase_106 | AC | 483 ms
30,136 KB |
testcase_107 | AC | 356 ms
26,784 KB |
testcase_108 | AC | 418 ms
29,392 KB |
testcase_109 | AC | 326 ms
26,920 KB |
testcase_110 | AC | 413 ms
27,920 KB |
testcase_111 | AC | 482 ms
31,840 KB |
testcase_112 | AC | 761 ms
41,116 KB |
testcase_113 | AC | 772 ms
40,988 KB |
testcase_114 | AC | 924 ms
40,980 KB |
testcase_115 | AC | 783 ms
40,852 KB |
testcase_116 | AC | 736 ms
39,892 KB |
testcase_117 | AC | 722 ms
40,596 KB |
testcase_118 | AC | 743 ms
43,488 KB |
testcase_119 | AC | 717 ms
40,020 KB |
testcase_120 | AC | 768 ms
40,276 KB |
testcase_121 | AC | 718 ms
40,392 KB |
testcase_122 | AC | 731 ms
40,312 KB |
testcase_123 | AC | 760 ms
40,660 KB |
testcase_124 | AC | 753 ms
40,540 KB |
testcase_125 | AC | 787 ms
47,536 KB |
testcase_126 | AC | 759 ms
55,664 KB |
testcase_127 | AC | 725 ms
40,792 KB |
testcase_128 | AC | 780 ms
47,072 KB |
testcase_129 | AC | 747 ms
47,728 KB |
testcase_130 | AC | 772 ms
43,740 KB |
testcase_131 | AC | 789 ms
47,440 KB |
testcase_132 | AC | 1,117 ms
86,560 KB |
testcase_133 | AC | 1,044 ms
84,696 KB |
testcase_134 | AC | 1,064 ms
83,020 KB |
testcase_135 | AC | 1,052 ms
80,736 KB |
testcase_136 | AC | 988 ms
77,008 KB |
testcase_137 | AC | 1,032 ms
76,248 KB |
testcase_138 | AC | 966 ms
75,616 KB |
testcase_139 | AC | 1,014 ms
75,016 KB |
testcase_140 | AC | 984 ms
74,224 KB |
testcase_141 | AC | 992 ms
73,120 KB |
testcase_142 | AC | 980 ms
73,232 KB |
ソースコード
// #pragma GCC target("avx2") // #pragma GCC optimize("O3") // #pragma GCC optimize("unroll-loops") // #define INTERACTIVE #include <bits/stdc++.h> using namespace std; namespace templates { // type using ll = long long; using ull = unsigned long long; using Pii = pair<int, int>; using Pil = pair<int, ll>; using Pli = pair<ll, int>; using Pll = pair<ll, ll>; template <class T> using pq = priority_queue<T>; template <class T> using qp = priority_queue<T, vector<T>, greater<T>>; // clang-format off #define vec(T, A, ...) vector<T> A(__VA_ARGS__); #define vvec(T, A, h, ...) vector<vector<T>> A(h, vector<T>(__VA_ARGS__)); #define vvvec(T, A, h1, h2, ...) vector<vector<vector<T>>> A(h1, vector<vector<T>>(h2, vector<T>(__VA_ARGS__))); // clang-format on // for loop #define fori1(a) for (ll _ = 0; _ < (a); _++) #define fori2(i, a) for (ll i = 0; i < (a); i++) #define fori3(i, a, b) for (ll i = (a); i < (b); i++) #define fori4(i, a, b, c) for (ll i = (a); ((c) > 0 || i > (b)) && ((c) < 0 || i < (b)); i += (c)) #define overload4(a, b, c, d, e, ...) e #define fori(...) overload4(__VA_ARGS__, fori4, fori3, fori2, fori1)(__VA_ARGS__) // declare and input // clang-format off #define INT(...) int __VA_ARGS__; inp(__VA_ARGS__); #define LL(...) ll __VA_ARGS__; inp(__VA_ARGS__); #define STRING(...) string __VA_ARGS__; inp(__VA_ARGS__); #define CHAR(...) char __VA_ARGS__; inp(__VA_ARGS__); #define DOUBLE(...) double __VA_ARGS__; STRING(str___); __VA_ARGS__ = stod(str___); #define VEC(T, A, n) vector<T> A(n); inp(A); #define VVEC(T, A, n, m) vector<vector<T>> A(n, vector<T>(m)); inp(A); // clang-format on // const value const ll MOD1 = 1000000007; const ll MOD9 = 998244353; const double PI = acos(-1); // other macro #if !defined(RIN__LOCAL) && !defined(INTERACTIVE) #define endl "\n" #endif #define spa ' ' #define len(A) ll(A.size()) #define all(A) begin(A), end(A) // function vector<char> stoc(string &S) { int n = S.size(); vector<char> ret(n); for (int i = 0; i < n; i++) ret[i] = S[i]; return ret; } string ctos(vector<char> &S) { int n = S.size(); string ret = ""; for (int i = 0; i < n; i++) ret += S[i]; return ret; } template <class T> auto min(const T &a) { return *min_element(all(a)); } template <class T> auto max(const T &a) { return *max_element(all(a)); } template <class T, class S> auto clamp(T &a, const S &l, const S &r) { return (a > r ? r : a < l ? l : a); } template <class T, class S> inline bool chmax(T &a, const S &b) { return (a < b ? a = b, 1 : 0); } template <class T, class S> inline bool chmin(T &a, const S &b) { return (a > b ? a = b, 1 : 0); } template <class T, class S> inline bool chclamp(T &a, const S &l, const S &r) { auto b = clamp(a, l, r); return (a != b ? a = b, 1 : 0); } template <typename T> T sum(vector<T> &A) { T tot = 0; for (auto a : A) tot += a; return tot; } template <typename T> vector<T> compression(vector<T> X) { sort(all(X)); X.erase(unique(all(X)), X.end()); return X; } // input and output namespace io { // __int128_t std::istream &operator>>(std::istream &is, __int128_t &value) { std::string str; is >> str; value = 0; int sign = 1; for (size_t i = 0; i < str.size(); i++) { if (i == 0 && str[i] == '-') { sign = -1; continue; } value = value * 10 + str[i] - '0'; } value *= sign; return is; } std::ostream &operator<<(std::ostream &dest, __int128_t value) { std::ostream::sentry s(dest); if (s) { __uint128_t tmp = value < 0 ? -value : value; char buffer[128]; char *d = std::end(buffer); do { --d; *d = "0123456789"[tmp % 10]; tmp /= 10; } while (tmp != 0); if (value < 0) { --d; *d = '-'; } int len = std::end(buffer) - d; if (dest.rdbuf()->sputn(d, len) != len) { dest.setstate(std::ios_base::badbit); } } return dest; } // vector<T> template <typename T> istream &operator>>(istream &is, vector<T> &A) { for (auto &a : A) is >> a; return is; } template <typename T> ostream &operator<<(ostream &os, vector<T> &A) { for (size_t i = 0; i < A.size(); i++) { os << A[i]; if (i != A.size() - 1) os << ' '; } return os; } // vector<vector<T>> template <typename T> istream &operator>>(istream &is, vector<vector<T>> &A) { for (auto &a : A) is >> a; return is; } template <typename T> ostream &operator<<(ostream &os, vector<vector<T>> &A) { for (size_t i = 0; i < A.size(); i++) { os << A[i]; if (i != A.size() - 1) os << endl; } return os; } // pair<S, T> template <typename S, typename T> istream &operator>>(istream &is, pair<S, T> &A) { is >> A.first >> A.second; return is; } template <typename S, typename T> ostream &operator<<(ostream &os, pair<S, T> &A) { os << A.first << ' ' << A.second; return os; } // vector<pair<S, T>> template <typename S, typename T> istream &operator>>(istream &is, vector<pair<S, T>> &A) { for (size_t i = 0; i < A.size(); i++) { is >> A[i]; } return is; } template <typename S, typename T> ostream &operator<<(ostream &os, vector<pair<S, T>> &A) { for (size_t i = 0; i < A.size(); i++) { os << A[i]; if (i != A.size() - 1) os << endl; } return os; } // tuple template <typename T, size_t N> struct TuplePrint { static ostream &print(ostream &os, const T &t) { TuplePrint<T, N - 1>::print(os, t); os << ' ' << get<N - 1>(t); return os; } }; template <typename T> struct TuplePrint<T, 1> { static ostream &print(ostream &os, const T &t) { os << get<0>(t); return os; } }; template <typename... Args> ostream &operator<<(ostream &os, const tuple<Args...> &t) { TuplePrint<decltype(t), sizeof...(Args)>::print(os, t); return os; } // io functions void FLUSH() { cout << flush; } void print() { cout << endl; } template <class Head, class... Tail> void print(Head &&head, Tail &&...tail) { cout << head; if (sizeof...(Tail)) cout << spa; print(std::forward<Tail>(tail)...); } template <typename T, typename S> void prisep(vector<T> &A, S sep) { int n = A.size(); for (int i = 0; i < n; i++) { cout << A[i]; if (i != n - 1) cout << sep; } cout << endl; } template <typename T, typename S> void priend(T A, S end) { cout << A << end; } template <typename T> void prispa(T A) { priend(A, spa); } template <typename T, typename S> bool printif(bool f, T A, S B) { if (f) print(A); else print(B); return f; } template <class... T> void inp(T &...a) { (cin >> ... >> a); } } // namespace io using namespace io; // read graph vector<vector<int>> read_edges(int n, int m, bool direct = false, int indexed = 1) { vector<vector<int>> edges(n, vector<int>()); for (int i = 0; i < m; i++) { INT(u, v); u -= indexed; v -= indexed; edges[u].push_back(v); if (!direct) edges[v].push_back(u); } return edges; } vector<vector<int>> read_tree(int n, int indexed = 1) { return read_edges(n, n - 1, false, indexed); } template <typename T = long long> vector<vector<pair<int, T>>> read_wedges(int n, int m, bool direct = false, int indexed = 1) { vector<vector<pair<int, T>>> edges(n, vector<pair<int, T>>()); for (int i = 0; i < m; i++) { INT(u, v); T w; inp(w); u -= indexed; v -= indexed; edges[u].push_back({v, w}); if (!direct) edges[v].push_back({u, w}); } return edges; } template <typename T = long long> vector<vector<pair<int, T>>> read_wtree(int n, int indexed = 1) { return read_wedges<T>(n, n - 1, false, indexed); } // yes / no namespace yesno { // yes inline bool yes(bool f = true) { cout << (f ? "yes" : "no") << endl; return f; } inline bool Yes(bool f = true) { cout << (f ? "Yes" : "No") << endl; return f; } inline bool YES(bool f = true) { cout << (f ? "YES" : "NO") << endl; return f; } // no inline bool no(bool f = true) { cout << (!f ? "yes" : "no") << endl; return f; } inline bool No(bool f = true) { cout << (!f ? "Yes" : "No") << endl; return f; } inline bool NO(bool f = true) { cout << (!f ? "YES" : "NO") << endl; return f; } // possible inline bool possible(bool f = true) { cout << (f ? "possible" : "impossible") << endl; return f; } inline bool Possible(bool f = true) { cout << (f ? "Possible" : "Impossible") << endl; return f; } inline bool POSSIBLE(bool f = true) { cout << (f ? "POSSIBLE" : "IMPOSSIBLE") << endl; return f; } // impossible inline bool impossible(bool f = true) { cout << (!f ? "possible" : "impossible") << endl; return f; } inline bool Impossible(bool f = true) { cout << (!f ? "Possible" : "Impossible") << endl; return f; } inline bool IMPOSSIBLE(bool f = true) { cout << (!f ? "POSSIBLE" : "IMPOSSIBLE") << endl; return f; } // Alice Bob inline bool Alice(bool f = true) { cout << (f ? "Alice" : "Bob") << endl; return f; } inline bool Bob(bool f = true) { cout << (f ? "Bob" : "Alice") << endl; return f; } // Takahashi Aoki inline bool Takahashi(bool f = true) { cout << (f ? "Takahashi" : "Aoki") << endl; return f; } inline bool Aoki(bool f = true) { cout << (f ? "Aoki" : "Takahashi") << endl; return f; } } // namespace yesno using namespace yesno; } // namespace templates using namespace templates; ////// https://nyaannyaan.github.io/library/string/run-enumerate.hpp // #line 2 "string/run-enumerate.hpp" #include <algorithm> #include <set> #include <utility> #include <vector> using namespace std; // #line 2 "string/z-algorithm.hpp" template <typename Container> vector<int> z_algorithm(const Container &s) { int n = s.size(); if (n == 0) return {}; vector<int> a(n); a[0] = n; int i = 1, j = 0; while (i < n) { while (i + j < n && s[j] == s[i + j]) j++; a[i] = j; if (j == 0) { i++; continue; } int k = 1; while (i + k < n && k + a[k] < j) a[i + k] = a[k], k++; i += k, j -= k; } return a; } /** * @brief Z algorithm */ // #line 10 "string/run-enumerate.hpp" // (p, l, r) : S[l, r) は周期 p かつ極大 // sum_{(p,l,r)} 1 <= n // sum_{(p,l,r)} (r-l)/p <= 3n // sum_{(p,l,r)} (r-l+1-2*p) = O(n log n) template <typename C> vector<tuple<int, int, int>> run_enumerate(const C &S) { int N = S.size(); using T = tuple<int, int, int>; vector<vector<pair<int, int>>> by_p(N + 1); auto solve_sub = [&](const C &l, const C &r) { vector<T> res; int n = l.size(), m = r.size(); C s = l, t = r; t.insert(end(t), begin(l), end(l)); t.insert(end(t), begin(r), end(r)); reverse(begin(s), end(s)); auto ZS = z_algorithm(s), ZT = z_algorithm(t); for (int p = 1; p <= n; p++) { int a = p == n ? p : min(ZS[p] + p, n); int b = min(ZT[n + m - p], m); if (a + b < 2 * p) continue; res.emplace_back(p, a, b); } return res; }; auto dfs = [&](auto rc, int L, int R) -> void { if (R - L <= 1) return; int M = (L + R) / 2; rc(rc, L, M), rc(rc, M, R); C SL{begin(S) + L, begin(S) + M}; C SR{begin(S) + M, begin(S) + R}; auto sub_res1 = solve_sub(SL, SR); for (auto &[p, a, b] : sub_res1) by_p[p].emplace_back(M - a, M + b); reverse(begin(SL), end(SL)); reverse(begin(SR), end(SR)); auto sub_res2 = solve_sub(SR, SL); for (auto &[p, a, b] : sub_res2) by_p[p].emplace_back(M - b, M + a); }; dfs(dfs, 0, N); vector<T> res; set<pair<int, int>> done; for (int p = 0; p <= N; p++) { auto &LR = by_p[p]; sort(begin(LR), end(LR), [](auto &x, auto &y) { if (x.first == y.first) return x.second > y.second; return x.first < y.first; }); int r = -1; for (auto &lr : LR) { if (r >= lr.second) continue; r = lr.second; if (!done.count(lr)) { done.insert(lr); res.emplace_back(p, lr.first, lr.second); } } } return res; } void solve() { STRING(S); int n = len(S); auto runs = run_enumerate(S); const int inf = 1 << 30; vec(int, dp, n + 1, inf); vec(Pii, bef, n + 1, {-1, -1}); dp[0] = 0; auto dc = [&](auto &&self, int l, int r, vector<tuple<int, int, int>> &runs) -> void { if (r - l == 1) { if (chmin(dp[r], dp[l] + 2)) { bef[r] = {l, 1}; } return; } int mid = (l + r) / 2; vector<tuple<int, int, int>> L; vector<tuple<int, int, int>> R; for (auto [p, l, r] : runs) { if (l <= mid) { int r_ = min(r, mid); if (r_ - l >= 2 * p) { L.emplace_back(p, l, r_); } } if (mid < r) { int l_ = max(l, mid); if (r - l_ >= 2 * p) { R.emplace_back(p, l_, r); } } } self(self, l, mid, L); for (auto [p, l, r] : runs) { if (r < mid or mid <= l) continue; fori(s, l, min(mid, l + p)) { int add = 1; int nex = 10; int fr = s; for (int j = 1; fr + j * p <= r; j++) { if (j == nex) { add++; nex *= 10; } if (dp[fr + j * p] > dp[fr] + add + p) { dp[fr + j * p] = dp[fr] + add + p; bef[fr + j * p] = {fr, p}; } else if (dp[fr + j * p] < dp[fr] + add) { add = 1; nex = 10; fr = fr + j * p; j = 0; } } } } { int mi = dp[mid]; int fr = mid; for (int i = mid - 1; i >= l; i--) { if (chmin(mi, dp[i] + mid - i)) { fr = i; } } int add = 1; int nex = 10; for (int j = 1; fr + j <= r; j++) { if (chmin(dp[fr + j], dp[fr] + j + 1)) { bef[fr + j] = {fr, j}; } } } self(self, mid, r, R); }; dc(dc, 0, n, runs); string ans; int i = n; while (i > 0) { auto [fr, p] = bef[i]; int c = (i - fr) / p; i = fr; string S_ = S.substr(i, p) + to_string(c); reverse(all(S_)); ans += S_; } reverse(all(ans)); print(ans); } int main() { #ifndef INTERACTIVE std::cin.tie(0)->sync_with_stdio(0); #endif // std::cout << std::fixed << std::setprecision(12); int t; t = 1; std::cin >> t; while (t--) solve(); return 0; } // // #pragma GCC target("avx2") // // #pragma GCC optimize("O3") // // #pragma GCC optimize("unroll-loops") // // #define INTERACTIVE // // #include "kyopro-cpp/template.hpp" // // ////// https://nyaannyaan.github.io/library/string/run-enumerate.hpp // // // #line 2 "string/run-enumerate.hpp" // // #include <algorithm> // #include <set> // #include <utility> // #include <vector> // using namespace std; // // // #line 2 "string/z-algorithm.hpp" // // template <typename Container> // vector<int> z_algorithm(const Container &s) { // int n = s.size(); // if (n == 0) return {}; // vector<int> a(n); // a[0] = n; // int i = 1, j = 0; // while (i < n) { // while (i + j < n && s[j] == s[i + j]) j++; // a[i] = j; // if (j == 0) { // i++; // continue; // } // int k = 1; // while (i + k < n && k + a[k] < j) a[i + k] = a[k], k++; // i += k, j -= k; // } // return a; // } // // /** // * @brief Z algorithm // */ // // #line 10 "string/run-enumerate.hpp" // // // (p, l, r) : S[l, r) は周期 p かつ極大 // // sum_{(p,l,r)} 1 <= n // // sum_{(p,l,r)} (r-l)/p <= 3n // // sum_{(p,l,r)} (r-l+1-2*p) = O(n log n) // template <typename C> // vector<tuple<int, int, int>> run_enumerate(const C &S) { // int N = S.size(); // using T = tuple<int, int, int>; // vector<vector<pair<int, int>>> by_p(N + 1); // // auto solve_sub = [&](const C &l, const C &r) { // vector<T> res; // int n = l.size(), m = r.size(); // C s = l, t = r; // t.insert(end(t), begin(l), end(l)); // t.insert(end(t), begin(r), end(r)); // reverse(begin(s), end(s)); // auto ZS = z_algorithm(s), ZT = z_algorithm(t); // for (int p = 1; p <= n; p++) { // int a = p == n ? p : min(ZS[p] + p, n); // int b = min(ZT[n + m - p], m); // if (a + b < 2 * p) continue; // res.emplace_back(p, a, b); // } // return res; // }; // // auto dfs = [&](auto rc, int L, int R) -> void { // if (R - L <= 1) return; // int M = (L + R) / 2; // rc(rc, L, M), rc(rc, M, R); // C SL{begin(S) + L, begin(S) + M}; // C SR{begin(S) + M, begin(S) + R}; // auto sub_res1 = solve_sub(SL, SR); // for (auto &[p, a, b] : sub_res1) by_p[p].emplace_back(M - a, M + b); // reverse(begin(SL), end(SL)); // reverse(begin(SR), end(SR)); // auto sub_res2 = solve_sub(SR, SL); // for (auto &[p, a, b] : sub_res2) by_p[p].emplace_back(M - b, M + a); // }; // dfs(dfs, 0, N); // // vector<T> res; // set<pair<int, int>> done; // for (int p = 0; p <= N; p++) { // auto &LR = by_p[p]; // sort(begin(LR), end(LR), [](auto &x, auto &y) { // if (x.first == y.first) return x.second > y.second; // return x.first < y.first; // }); // int r = -1; // for (auto &lr : LR) { // if (r >= lr.second) continue; // r = lr.second; // if (!done.count(lr)) { // done.insert(lr); // res.emplace_back(p, lr.first, lr.second); // } // } // } // return res; // } // // void solve() { // STRING(S); // int n = len(S); // auto runs = run_enumerate(S); // const int inf = 1 << 30; // vec(int, dp, n + 1, inf); // vec(Pii, bef, n + 1, {-1, -1}); // dp[0] = 0; // // auto dc = [&](auto &&self, int l, int r, vector<tuple<int, int, int>> &runs) -> void { // if (r - l == 1) { // if (chmin(dp[r], dp[l] + 2)) { // bef[r] = {l, 1}; // } // return; // } // int mid = (l + r) / 2; // vector<tuple<int, int, int>> L; // vector<tuple<int, int, int>> R; // for (auto [p, l, r] : runs) { // if (l <= mid) { // int r_ = min(r, mid); // if (r_ - l >= 2 * p) { // L.emplace_back(p, l, r_); // } // } // if (mid < r) { // int l_ = max(l, mid); // if (r - l_ >= 2 * p) { // R.emplace_back(p, l_, r); // } // } // } // // self(self, l, mid, L); // // for (auto [p, l, r] : runs) { // if (r < mid or mid <= l) continue; // // fori(s, l, min(mid, l + p)) { // int add = 1; // int nex = 10; // int fr = s; // for (int j = 1; fr + j * p <= r; j++) { // if (j == nex) { // add++; // nex *= 10; // } // if (dp[fr + j * p] > dp[fr] + add + p) { // dp[fr + j * p] = dp[fr] + add + p; // bef[fr + j * p] = {fr, p}; // } else if (dp[fr + j * p] < dp[fr] + add) { // add = 1; // nex = 10; // fr = fr + j * p; // j = 0; // } // } // } // } // // { // int mi = dp[mid]; // int fr = mid; // for (int i = mid - 1; i >= l; i--) { // if (chmin(mi, dp[i] + mid - i)) { // fr = i; // } // } // int add = 1; // int nex = 10; // for (int j = 1; fr + j <= r; j++) { // if (chmin(dp[fr + j], dp[fr] + j + 1)) { // bef[fr + j] = {fr, j}; // } // } // } // // self(self, mid, r, R); // }; // // dc(dc, 0, n, runs); // // string ans; // int i = n; // while (i > 0) { // auto [fr, p] = bef[i]; // int c = (i - fr) / p; // i = fr; // string S_ = S.substr(i, p) + to_string(c); // reverse(all(S_)); // ans += S_; // } // reverse(all(ans)); // // print(ans); // } // // int main() { // #ifndef INTERACTIVE // std::cin.tie(0)->sync_with_stdio(0); // #endif // // std::cout << std::fixed << std::setprecision(12); // int t; // t = 1; // std::cin >> t; // while (t--) solve(); // return 0; // }