結果
問題 | No.3000 Optimal Run Length Encoding |
ユーザー | NyaanNyaan |
提出日時 | 2024-12-13 16:06:52 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 6,003 ms / 10,000 ms |
コード長 | 18,447 bytes |
コンパイル時間 | 4,933 ms |
コンパイル使用メモリ | 297,104 KB |
実行使用メモリ | 368,948 KB |
最終ジャッジ日時 | 2024-12-25 19:48:31 |
合計ジャッジ時間 | 220,371 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 3 ms
5,248 KB |
testcase_01 | AC | 290 ms
5,248 KB |
testcase_02 | AC | 303 ms
5,248 KB |
testcase_03 | AC | 270 ms
5,248 KB |
testcase_04 | AC | 281 ms
5,248 KB |
testcase_05 | AC | 239 ms
5,248 KB |
testcase_06 | AC | 244 ms
5,248 KB |
testcase_07 | AC | 276 ms
5,248 KB |
testcase_08 | AC | 265 ms
5,248 KB |
testcase_09 | AC | 242 ms
5,248 KB |
testcase_10 | AC | 229 ms
5,248 KB |
testcase_11 | AC | 257 ms
5,248 KB |
testcase_12 | AC | 267 ms
5,248 KB |
testcase_13 | AC | 259 ms
5,248 KB |
testcase_14 | AC | 230 ms
5,248 KB |
testcase_15 | AC | 232 ms
5,248 KB |
testcase_16 | AC | 250 ms
5,248 KB |
testcase_17 | AC | 215 ms
5,248 KB |
testcase_18 | AC | 158 ms
5,248 KB |
testcase_19 | AC | 453 ms
5,248 KB |
testcase_20 | AC | 457 ms
5,248 KB |
testcase_21 | AC | 455 ms
5,248 KB |
testcase_22 | AC | 458 ms
5,248 KB |
testcase_23 | AC | 456 ms
5,248 KB |
testcase_24 | AC | 462 ms
5,248 KB |
testcase_25 | AC | 462 ms
5,248 KB |
testcase_26 | AC | 462 ms
5,248 KB |
testcase_27 | AC | 405 ms
5,248 KB |
testcase_28 | AC | 443 ms
5,248 KB |
testcase_29 | AC | 473 ms
5,248 KB |
testcase_30 | AC | 529 ms
5,248 KB |
testcase_31 | AC | 556 ms
5,248 KB |
testcase_32 | AC | 598 ms
5,248 KB |
testcase_33 | AC | 646 ms
5,248 KB |
testcase_34 | AC | 771 ms
15,232 KB |
testcase_35 | AC | 383 ms
34,752 KB |
testcase_36 | AC | 395 ms
34,240 KB |
testcase_37 | AC | 397 ms
34,624 KB |
testcase_38 | AC | 392 ms
34,504 KB |
testcase_39 | AC | 394 ms
34,636 KB |
testcase_40 | AC | 388 ms
34,496 KB |
testcase_41 | AC | 382 ms
34,764 KB |
testcase_42 | AC | 386 ms
34,616 KB |
testcase_43 | AC | 403 ms
34,376 KB |
testcase_44 | AC | 388 ms
34,368 KB |
testcase_45 | AC | 6,003 ms
368,820 KB |
testcase_46 | AC | 4,899 ms
331,692 KB |
testcase_47 | AC | 5,288 ms
354,244 KB |
testcase_48 | AC | 4,824 ms
337,992 KB |
testcase_49 | AC | 5,095 ms
346,908 KB |
testcase_50 | AC | 4,984 ms
343,764 KB |
testcase_51 | AC | 5,304 ms
351,036 KB |
testcase_52 | AC | 5,108 ms
353,496 KB |
testcase_53 | AC | 5,163 ms
351,040 KB |
testcase_54 | AC | 4,942 ms
328,824 KB |
testcase_55 | AC | 5,698 ms
368,948 KB |
testcase_56 | AC | 5,205 ms
351,488 KB |
testcase_57 | AC | 5,030 ms
342,964 KB |
testcase_58 | AC | 4,659 ms
335,124 KB |
testcase_59 | AC | 5,215 ms
354,856 KB |
testcase_60 | AC | 5,680 ms
368,944 KB |
testcase_61 | AC | 5,594 ms
363,500 KB |
testcase_62 | AC | 4,813 ms
342,936 KB |
testcase_63 | AC | 5,054 ms
355,192 KB |
testcase_64 | AC | 4,672 ms
331,340 KB |
testcase_65 | AC | 923 ms
69,800 KB |
testcase_66 | AC | 505 ms
43,024 KB |
testcase_67 | AC | 936 ms
62,352 KB |
testcase_68 | AC | 470 ms
41,752 KB |
testcase_69 | AC | 1,070 ms
77,592 KB |
testcase_70 | AC | 1,650 ms
110,924 KB |
testcase_71 | AC | 448 ms
38,724 KB |
testcase_72 | AC | 1,064 ms
91,868 KB |
testcase_73 | AC | 559 ms
76,152 KB |
testcase_74 | AC | 1,223 ms
95,296 KB |
testcase_75 | AC | 599 ms
77,508 KB |
testcase_76 | AC | 1,135 ms
96,480 KB |
testcase_77 | AC | 470 ms
66,356 KB |
testcase_78 | AC | 987 ms
89,180 KB |
testcase_79 | AC | 497 ms
66,912 KB |
testcase_80 | AC | 880 ms
84,600 KB |
testcase_81 | AC | 83 ms
5,248 KB |
testcase_82 | AC | 1,333 ms
89,204 KB |
testcase_83 | AC | 1,330 ms
90,024 KB |
testcase_84 | AC | 1,325 ms
89,836 KB |
testcase_85 | AC | 1,277 ms
86,656 KB |
testcase_86 | AC | 1,259 ms
83,992 KB |
testcase_87 | AC | 1,355 ms
93,256 KB |
testcase_88 | AC | 1,278 ms
86,744 KB |
testcase_89 | AC | 1,347 ms
93,132 KB |
testcase_90 | AC | 1,249 ms
85,600 KB |
testcase_91 | AC | 1,269 ms
83,032 KB |
testcase_92 | AC | 1,290 ms
87,796 KB |
testcase_93 | AC | 1,328 ms
87,704 KB |
testcase_94 | AC | 1,341 ms
89,616 KB |
testcase_95 | AC | 1,304 ms
88,984 KB |
testcase_96 | AC | 1,317 ms
88,592 KB |
testcase_97 | AC | 452 ms
36,760 KB |
testcase_98 | AC | 388 ms
36,500 KB |
testcase_99 | AC | 343 ms
32,208 KB |
testcase_100 | AC | 430 ms
33,728 KB |
testcase_101 | AC | 459 ms
39,448 KB |
testcase_102 | AC | 622 ms
59,436 KB |
testcase_103 | AC | 650 ms
59,396 KB |
testcase_104 | AC | 722 ms
66,656 KB |
testcase_105 | AC | 625 ms
64,012 KB |
testcase_106 | AC | 743 ms
64,444 KB |
testcase_107 | AC | 628 ms
62,444 KB |
testcase_108 | AC | 663 ms
65,320 KB |
testcase_109 | AC | 594 ms
58,792 KB |
testcase_110 | AC | 648 ms
65,260 KB |
testcase_111 | AC | 798 ms
71,572 KB |
testcase_112 | AC | 1,045 ms
61,408 KB |
testcase_113 | AC | 1,047 ms
61,512 KB |
testcase_114 | AC | 1,023 ms
61,504 KB |
testcase_115 | AC | 1,022 ms
61,616 KB |
testcase_116 | AC | 1,019 ms
61,836 KB |
testcase_117 | AC | 1,020 ms
61,760 KB |
testcase_118 | AC | 1,034 ms
61,540 KB |
testcase_119 | AC | 1,011 ms
61,708 KB |
testcase_120 | AC | 1,024 ms
61,668 KB |
testcase_121 | AC | 1,085 ms
61,732 KB |
testcase_122 | AC | 1,039 ms
61,724 KB |
testcase_123 | AC | 1,057 ms
61,720 KB |
testcase_124 | AC | 999 ms
61,700 KB |
testcase_125 | AC | 1,042 ms
60,960 KB |
testcase_126 | AC | 1,053 ms
60,700 KB |
testcase_127 | AC | 1,020 ms
61,800 KB |
testcase_128 | AC | 1,014 ms
61,076 KB |
testcase_129 | AC | 1,043 ms
60,824 KB |
testcase_130 | AC | 1,028 ms
61,500 KB |
testcase_131 | AC | 1,199 ms
60,796 KB |
testcase_132 | AC | 1,490 ms
90,720 KB |
testcase_133 | AC | 1,518 ms
110,028 KB |
testcase_134 | AC | 1,508 ms
104,008 KB |
testcase_135 | AC | 1,489 ms
95,332 KB |
testcase_136 | AC | 1,462 ms
90,840 KB |
testcase_137 | AC | 1,471 ms
88,100 KB |
testcase_138 | AC | 1,419 ms
86,328 KB |
testcase_139 | AC | 1,417 ms
84,432 KB |
testcase_140 | AC | 1,404 ms
88,428 KB |
testcase_141 | AC | 1,368 ms
85,532 KB |
testcase_142 | AC | 1,322 ms
76,648 KB |
ソースコード
/** * date : 2024-12-12 23:46:05 * author : Nyaan */ using namespace std; // intrinstic #include <immintrin.h> #include <algorithm> #include <array> #include <bitset> #include <cassert> #include <cctype> #include <cfenv> #include <cfloat> #include <chrono> #include <cinttypes> #include <climits> #include <cmath> #include <complex> #include <cstdarg> #include <cstddef> #include <cstdint> #include <cstdio> #include <cstdlib> #include <cstring> #include <deque> #include <fstream> #include <functional> #include <initializer_list> #include <iomanip> #include <ios> #include <iostream> #include <istream> #include <iterator> #include <limits> #include <list> #include <map> #include <memory> #include <new> #include <numeric> #include <ostream> #include <queue> #include <random> #include <set> #include <sstream> #include <stack> #include <streambuf> #include <string> #include <tr2/dynamic_bitset> #include <tuple> #include <type_traits> #include <typeinfo> #include <unordered_map> #include <unordered_set> #include <utility> #include <vector> // utility namespace Nyaan { using ll = long long; using i64 = long long; using u64 = unsigned long long; using i128 = __int128_t; using u128 = __uint128_t; template <typename T> using V = vector<T>; template <typename T> using VV = vector<vector<T>>; using vi = vector<int>; using vl = vector<long long>; using vd = V<double>; using vs = V<string>; using vvi = vector<vector<int>>; using vvl = vector<vector<long long>>; template <typename T> using minpq = priority_queue<T, vector<T>, greater<T>>; template <typename T, typename U> struct P : pair<T, U> { template <typename... Args> constexpr P(Args... args) : pair<T, U>(args...) {} using pair<T, U>::first; using pair<T, U>::second; P &operator+=(const P &r) { first += r.first; second += r.second; return *this; } P &operator-=(const P &r) { first -= r.first; second -= r.second; return *this; } P &operator*=(const P &r) { first *= r.first; second *= r.second; return *this; } template <typename S> P &operator*=(const S &r) { first *= r, second *= r; return *this; } P operator+(const P &r) const { return P(*this) += r; } P operator-(const P &r) const { return P(*this) -= r; } P operator*(const P &r) const { return P(*this) *= r; } template <typename S> P operator*(const S &r) const { return P(*this) *= r; } P operator-() const { return P{-first, -second}; } }; using pl = P<ll, ll>; using pi = P<int, int>; using vp = V<pl>; constexpr int inf = 1001001001; constexpr long long infLL = 4004004004004004004LL; template <typename T> int sz(const T &t) { return t.size(); } template <typename T, typename U> inline bool amin(T &x, U y) { return (y < x) ? (x = y, true) : false; } template <typename T, typename U> inline bool amax(T &x, U y) { return (x < y) ? (x = y, true) : false; } template <typename T> inline T Max(const vector<T> &v) { return *max_element(begin(v), end(v)); } template <typename T> inline T Min(const vector<T> &v) { return *min_element(begin(v), end(v)); } template <typename T> inline long long Sum(const vector<T> &v) { return accumulate(begin(v), end(v), 0LL); } template <typename T> int lb(const vector<T> &v, const T &a) { return lower_bound(begin(v), end(v), a) - begin(v); } template <typename T> int ub(const vector<T> &v, const T &a) { return upper_bound(begin(v), end(v), a) - begin(v); } constexpr long long TEN(int n) { long long ret = 1, x = 10; for (; n; x *= x, n >>= 1) ret *= (n & 1 ? x : 1); return ret; } template <typename T, typename U> pair<T, U> mkp(const T &t, const U &u) { return make_pair(t, u); } template <typename T> vector<T> mkrui(const vector<T> &v, bool rev = false) { vector<T> ret(v.size() + 1); if (rev) { for (int i = int(v.size()) - 1; i >= 0; i--) ret[i] = v[i] + ret[i + 1]; } else { for (int i = 0; i < int(v.size()); i++) ret[i + 1] = ret[i] + v[i]; } return ret; }; template <typename T> vector<T> mkuni(const vector<T> &v) { vector<T> ret(v); sort(ret.begin(), ret.end()); ret.erase(unique(ret.begin(), ret.end()), ret.end()); return ret; } template <typename F> vector<int> mkord(int N, F f) { vector<int> ord(N); iota(begin(ord), end(ord), 0); sort(begin(ord), end(ord), f); return ord; } template <typename T> vector<int> mkinv(vector<T> &v) { int max_val = *max_element(begin(v), end(v)); vector<int> inv(max_val + 1, -1); for (int i = 0; i < (int)v.size(); i++) inv[v[i]] = i; return inv; } vector<int> mkiota(int n) { vector<int> ret(n); iota(begin(ret), end(ret), 0); return ret; } template <typename T> T mkrev(const T &v) { T w{v}; reverse(begin(w), end(w)); return w; } template <typename T> bool nxp(T &v) { return next_permutation(begin(v), end(v)); } // 返り値の型は入力の T に依存 // i 要素目 : [0, a[i]) template <typename T> vector<vector<T>> product(const vector<T> &a) { vector<vector<T>> ret; vector<T> v; auto dfs = [&](auto rc, int i) -> void { if (i == (int)a.size()) { ret.push_back(v); return; } for (int j = 0; j < a[i]; j++) v.push_back(j), rc(rc, i + 1), v.pop_back(); }; dfs(dfs, 0); return ret; } // F : function(void(T&)), mod を取る操作 // T : 整数型のときはオーバーフローに注意する template <typename T> T Power(T a, long long n, const T &I, const function<void(T &)> &f) { T res = I; for (; n; f(a = a * a), n >>= 1) { if (n & 1) f(res = res * a); } return res; } // T : 整数型のときはオーバーフローに注意する template <typename T> T Power(T a, long long n, const T &I = T{1}) { return Power(a, n, I, function<void(T &)>{[](T &) -> void {}}); } template <typename T> T Rev(const T &v) { T res = v; reverse(begin(res), end(res)); return res; } template <typename T> vector<T> Transpose(const vector<T> &v) { using U = typename T::value_type; if (v.empty()) return {}; int H = v.size(), W = v[0].size(); vector res(W, T(H, U{})); for (int i = 0; i < H; i++) { for (int j = 0; j < W; j++) { res[j][i] = v[i][j]; } } return res; } template <typename T> vector<T> Rotate(const vector<T> &v, int clockwise = true) { using U = typename T::value_type; int H = v.size(), W = v[0].size(); vector res(W, T(H, U{})); for (int i = 0; i < H; i++) { for (int j = 0; j < W; j++) { if (clockwise) { res[W - 1 - j][i] = v[i][j]; } else { res[j][H - 1 - i] = v[i][j]; } } } return res; } } // namespace Nyaan // bit operation namespace Nyaan { __attribute__((target("popcnt"))) inline int popcnt(const u64 &a) { return __builtin_popcountll(a); } inline int lsb(const u64 &a) { return a ? __builtin_ctzll(a) : 64; } inline int ctz(const u64 &a) { return a ? __builtin_ctzll(a) : 64; } inline int msb(const u64 &a) { return a ? 63 - __builtin_clzll(a) : -1; } template <typename T> inline int gbit(const T &a, int i) { return (a >> i) & 1; } template <typename T> inline void sbit(T &a, int i, bool b) { if (gbit(a, i) != b) a ^= T(1) << i; } constexpr long long PW(int n) { return 1LL << n; } constexpr long long MSK(int n) { return (1LL << n) - 1; } } // namespace Nyaan // inout namespace Nyaan { template <typename T, typename U> ostream &operator<<(ostream &os, const pair<T, U> &p) { os << p.first << " " << p.second; return os; } template <typename T, typename U> istream &operator>>(istream &is, pair<T, U> &p) { is >> p.first >> p.second; return is; } template <typename T> ostream &operator<<(ostream &os, const vector<T> &v) { int s = (int)v.size(); for (int i = 0; i < s; i++) os << (i ? " " : "") << v[i]; return os; } template <typename T> istream &operator>>(istream &is, vector<T> &v) { for (auto &x : v) is >> x; return is; } istream &operator>>(istream &is, __int128_t &x) { string S; is >> S; x = 0; int flag = 0; for (auto &c : S) { if (c == '-') { flag = true; continue; } x *= 10; x += c - '0'; } if (flag) x = -x; return is; } istream &operator>>(istream &is, __uint128_t &x) { string S; is >> S; x = 0; for (auto &c : S) { x *= 10; x += c - '0'; } return is; } ostream &operator<<(ostream &os, __int128_t x) { if (x == 0) return os << 0; if (x < 0) os << '-', x = -x; string S; while (x) S.push_back('0' + x % 10), x /= 10; reverse(begin(S), end(S)); return os << S; } ostream &operator<<(ostream &os, __uint128_t x) { if (x == 0) return os << 0; string S; while (x) S.push_back('0' + x % 10), x /= 10; reverse(begin(S), end(S)); return os << S; } void in() {} template <typename T, class... U> void in(T &t, U &...u) { cin >> t; in(u...); } void out() { cout << "\n"; } template <typename T, class... U, char sep = ' '> void out(const T &t, const U &...u) { cout << t; if (sizeof...(u)) cout << sep; out(u...); } struct IoSetupNya { IoSetupNya() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(15); cerr << fixed << setprecision(7); } } iosetupnya; } // namespace Nyaan // debug #ifdef NyaanDebug #define trc(...) (void(0)) #endif #ifndef NyaanDebug #define trc(...) (void(0)) #endif #ifndef NyaanLocal #define trc2(...) (void(0)) #endif // macro #define each(x, v) for (auto &&x : v) #define each2(x, y, v) for (auto &&[x, y] : v) #define all(v) (v).begin(), (v).end() #define rep(i, N) for (long long i = 0; i < (long long)(N); i++) #define repr(i, N) for (long long i = (long long)(N) - 1; i >= 0; i--) #define rep1(i, N) for (long long i = 1; i <= (long long)(N); i++) #define repr1(i, N) for (long long i = (N); (long long)(i) > 0; i--) #define reg(i, a, b) for (long long i = (a); i < (b); i++) #define regr(i, a, b) for (long long i = (b) - 1; i >= (a); i--) #define fi first #define se second #define ini(...) \ int __VA_ARGS__; \ in(__VA_ARGS__) #define inl(...) \ long long __VA_ARGS__; \ in(__VA_ARGS__) #define ins(...) \ string __VA_ARGS__; \ in(__VA_ARGS__) #define in2(s, t) \ for (int i = 0; i < (int)s.size(); i++) { \ in(s[i], t[i]); \ } #define in3(s, t, u) \ for (int i = 0; i < (int)s.size(); i++) { \ in(s[i], t[i], u[i]); \ } #define in4(s, t, u, v) \ for (int i = 0; i < (int)s.size(); i++) { \ in(s[i], t[i], u[i], v[i]); \ } #define die(...) \ do { \ Nyaan::out(__VA_ARGS__); \ return; \ } while (0) namespace Nyaan { void solve(); } int main() { Nyaan::solve(); } // using namespace std; 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 */ // (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; } using namespace Nyaan; constexpr int INF = 1001001001; template <typename T> struct SWAG { vector<T> a0, a1, r0, r1; T f0, f1; T f(const T &a, const T &b) { return a.first < b.first ? a : b; } static constexpr T I{INF, -1}; SWAG() : f0(I), f1(I) {} private: void push_s0(const T &x) { a0.push_back(x); r0.push_back(f0 = f(x, f0)); } void push_s1(const T &x) { a1.push_back(x); r1.push_back(f1 = f(f1, x)); } void transfer() { while (!a1.empty()) { push_s0(a1.back()); a1.pop_back(); } while (!r1.empty()) r1.pop_back(); f1 = I; } public: void push(const T &x) { if (a0.empty()) { push_s0(x); transfer(); } else { push_s1(x); } } T pop() { if (a0.empty()) transfer(); T res = a0.back(); a0.pop_back(); r0.pop_back(); f0 = r0.empty() ? I : r0.back(); return res; } T query() { return f(f0, f1); } int size() { return a0.size() + a1.size(); } }; struct Data { vector<SWAG<pair<int, int>>> dat; static constexpr int size_table[9] = { 9, 90, 900, 9000, 90000, 900000, 9000000, 90000000, 900000000, }; Data() { dat.push_back({}); } void push(int val, int pos) { dat[0].push({val, pos}); for (int i = 0; i < (int)dat.size(); i++) { if ((int)dat[i].size() <= size_table[i]) break; if (i + 1 == (int)dat.size()) dat.push_back({}); dat[i + 1].push(dat[i].pop()); } } pair<int, int> query() { int res = INF, pos = -1; for (int i = 0; i < (int)dat.size(); i++) { auto [val, idx] = dat[i].query(); val += i + 1; if (val < res) res = val, pos = idx; } return {res, pos}; } }; int min_period(string S) { int N = sz(S); vi z = z_algorithm(S); rep1(i, N - 1) { if (N % i != 0) continue; if (z[i] == N - i) return i; } return N; } int num_len(int x) { return to_string(x).size(); } string naive(string S) { int N = sz(S); vector<int> dp(N + 1, inf); vector<int> prev(N + 1, -1); dp[0] = 0; rep(i, N) reg(j, i + 1, N + 1) { int p = min_period(S.substr(i, j - i)); int val = dp[i] + p + num_len((j - i) / p); if (amin(dp[j], val)) prev[j] = i; } string ans; for (int j = N; j != 0;) { int i = prev[j]; int p = min_period(S.substr(i, j - i)); string cur = S.substr(i, p) + to_string((j - i) / p); reverse(all(cur)); ans += cur; j = i; } reverse(all(ans)); assert(sz(ans) == dp[N]); return ans; } string calc(string S) { int N = sz(S); auto run = run_enumerate(S); VV<pi> rundata(N + 1); for (auto &[p, l, r] : run) { rep(i, p) if (l + i + 2 * p <= r) rundata[l + i + 2 * p].emplace_back(r, p); } vector<int> dp(N + 1, inf); vector<int> prev(N + 1, -1); dp[0] = 0; // dp[i] - i を保持する int prev_best = inf, prev_pos = -1; struct Data2 { int p, r; Data dat; }; VV<Data2> datas(N + 1); rep(i, N + 1) { // dp[j] (j<i) からの遷移 if (amin(dp[i], prev_best + i + 1)) prev[i] = prev_pos; // 新規 run の追加 for (auto &[r, p] : rundata[i]) { Data2 dat; dat.r = r, dat.p = p; dat.dat.push(dp[i - 2 * p], i - 2 * p); dat.dat.push(dp[i - 1 * p], i - 1 * p); datas[i].push_back(dat); } // run からの遷移 each(dat, datas[i]) { auto [val, pos] = dat.dat.query(); // trc(i, dat.p, val, pos); if (amin(dp[i], val + dat.p)) prev[i] = pos; } // dp[i] の値, 決定 // segtree の更新 if (amin(prev_best, dp[i] - i)) prev_pos = i; // run の更新 each(dat, datas[i]) { if (i + dat.p <= dat.r) { dat.dat.push(dp[i], i); int nxt = i + dat.p; datas[nxt].push_back(std::move(dat)); } else { Data2 d = std::move(dat); } } } string ans; for (int j = N; j != 0;) { int i = prev[j]; int p = min_period(S.substr(i, j - i)); string cur = S.substr(i, p) + to_string((j - i) / p); reverse(all(cur)); ans += cur; j = i; } reverse(all(ans)); assert(sz(ans) == dp[N]); return ans; } string restore(string T, int limit) { string res; for (int i = 0; i < (int)T.size();) { int j = i; while (j < (int)T.size() and 'a' <= T[j] and T[j] <= 'z') j++; int k = j; while (k < (int)T.size() and '0' <= T[k] and T[k] <= '9') k++; if (i == j or j == k) return {}; if (k - j > 15) return {}; string sub = T.substr(i, j - i); int num = stoll(T.substr(j, k - j)); while (num--) { res += sub; if ((int)res.size() > limit) return res; } i = k; } return res; } void check(string S) { auto an = naive(S); auto ac = calc(S); assert(restore(an, sz(S)) == S); assert(restore(ac, sz(S)) == S); if (sz(an) != sz(ac)) { trc2(S, an, ac); exit(1); } } void q() { /* check("baaaaaaaaa"); check("aabababc"); check("aaaaaaaaaa"); rep1(len, 11) { rep(b, Power(3, len)) { string S; for (int x = b; sz(S) < len; x /= 3) S.push_back(x % 3 + 'a'); check(S); } trc2(len, "OK"); } rep1(len, 18) { rep(b, PW(len)) { string S; rep(i, len) S.push_back('a' + gbit(b, i)); check(S); } trc2(len, "OK"); } */ ins(S); auto ans = calc(S); assert(restore(ans, sz(S)) == S); out(ans); } void Nyaan::solve() { int t = 1; in(t); while (t--) q(); }