#1033271 (C++23) No.3000 Optimal Run Length Encoding

提出ソース

結果

問題	No.3000 Optimal Run Length Encoding
ユーザー	👑 rin204
提出日時	2024-12-25 21:01:08
言語	C++23 (gcc 13.3.0 + boost 1.87.0)
結果	AC
実行時間	3,603 ms / 10,000 ms
コード長	25,459 bytes
コンパイル時間	4,534 ms
コンパイル使用メモリ	291,996 KB
実行使用メモリ	110,560 KB
最終ジャッジ日時	2024-12-25 21:03:49
合計ジャッジ時間	148,089 ms
ジャッジサーバーID （参考情報）	judge3 / judge4

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	AC	269 ms 5,248 KB
testcase_02	AC	219 ms 5,248 KB
testcase_03	AC	217 ms 5,248 KB
testcase_04	AC	211 ms 5,248 KB
testcase_05	AC	179 ms 5,248 KB
testcase_06	AC	192 ms 5,248 KB
testcase_07	AC	202 ms 5,248 KB
testcase_08	AC	229 ms 5,248 KB
testcase_09	AC	183 ms 5,248 KB
testcase_10	AC	175 ms 5,248 KB
testcase_11	AC	186 ms 5,248 KB
testcase_12	AC	196 ms 5,248 KB
testcase_13	AC	191 ms 5,248 KB
testcase_14	AC	183 ms 5,248 KB
testcase_15	AC	163 ms 5,248 KB
testcase_16	AC	192 ms 5,248 KB
testcase_17	AC	158 ms 5,248 KB
testcase_18	AC	110 ms 5,248 KB
testcase_19	AC	334 ms 5,248 KB
testcase_20	AC	368 ms 5,248 KB
testcase_21	AC	362 ms 5,248 KB
testcase_22	AC	397 ms 5,248 KB
testcase_23	AC	354 ms 5,248 KB
testcase_24	AC	342 ms 5,248 KB
testcase_25	AC	367 ms 5,248 KB
testcase_26	AC	363 ms 5,248 KB
testcase_27	AC	311 ms 5,248 KB
testcase_28	AC	348 ms 5,248 KB
testcase_29	AC	379 ms 5,248 KB
testcase_30	AC	408 ms 5,248 KB
testcase_31	AC	444 ms 5,248 KB
testcase_32	AC	484 ms 5,248 KB
testcase_33	AC	563 ms 5,248 KB
testcase_34	AC	636 ms 9,844 KB
testcase_35	AC	314 ms 22,188 KB
testcase_36	AC	349 ms 22,112 KB
testcase_37	AC	317 ms 22,244 KB
testcase_38	AC	322 ms 22,336 KB
testcase_39	AC	325 ms 22,128 KB
testcase_40	AC	322 ms 22,244 KB
testcase_41	AC	320 ms 22,248 KB
testcase_42	AC	350 ms 22,120 KB
testcase_43	AC	312 ms 22,280 KB
testcase_44	AC	338 ms 22,348 KB
testcase_45	AC	1,773 ms 47,108 KB
testcase_46	AC	1,291 ms 24,100 KB
testcase_47	AC	1,632 ms 43,200 KB
testcase_48	AC	1,358 ms 29,856 KB
testcase_49	AC	1,594 ms 38,348 KB
testcase_50	AC	1,481 ms 28,724 KB
testcase_51	AC	1,649 ms 42,988 KB
testcase_52	AC	1,431 ms 30,656 KB
testcase_53	AC	1,678 ms 43,316 KB
testcase_54	AC	1,268 ms 27,888 KB
testcase_55	AC	1,728 ms 47,236 KB
testcase_56	AC	1,552 ms 33,252 KB
testcase_57	AC	1,593 ms 37,992 KB
testcase_58	AC	1,333 ms 24,112 KB
testcase_59	AC	1,742 ms 43,092 KB
testcase_60	AC	1,777 ms 47,364 KB
testcase_61	AC	1,791 ms 44,260 KB
testcase_62	AC	1,355 ms 24,584 KB
testcase_63	AC	1,739 ms 39,312 KB
testcase_64	AC	1,320 ms 24,040 KB
testcase_65	AC	1,083 ms 47,456 KB
testcase_66	AC	737 ms 23,884 KB
testcase_67	AC	1,067 ms 39,660 KB
testcase_68	AC	679 ms 22,808 KB
testcase_69	AC	1,168 ms 43,936 KB
testcase_70	AC	1,890 ms 110,560 KB
testcase_71	AC	364 ms 22,292 KB
testcase_72	AC	981 ms 35,200 KB
testcase_73	AC	595 ms 23,788 KB
testcase_74	AC	915 ms 36,120 KB
testcase_75	AC	410 ms 22,244 KB
testcase_76	AC	888 ms 35,480 KB
testcase_77	AC	540 ms 22,680 KB
testcase_78	AC	1,002 ms 36,356 KB
testcase_79	AC	544 ms 23,424 KB
testcase_80	AC	966 ms 35,776 KB
testcase_81	AC	124 ms 5,248 KB
testcase_82	AC	1,372 ms 67,960 KB
testcase_83	AC	1,352 ms 67,588 KB
testcase_84	AC	1,340 ms 67,504 KB
testcase_85	AC	1,317 ms 66,508 KB
testcase_86	AC	1,303 ms 63,892 KB
testcase_87	AC	1,432 ms 71,488 KB
testcase_88	AC	1,259 ms 65,040 KB
testcase_89	AC	1,354 ms 71,360 KB
testcase_90	AC	1,283 ms 61,832 KB
testcase_91	AC	1,303 ms 63,164 KB
testcase_92	AC	1,316 ms 67,272 KB
testcase_93	AC	1,337 ms 67,432 KB
testcase_94	AC	1,322 ms 66,940 KB
testcase_95	AC	1,288 ms 67,180 KB
testcase_96	AC	1,295 ms 67,088 KB
testcase_97	AC	485 ms 25,936 KB
testcase_98	AC	440 ms 21,224 KB
testcase_99	AC	358 ms 19,484 KB
testcase_100	AC	425 ms 23,564 KB
testcase_101	AC	494 ms 24,688 KB
testcase_102	AC	2,532 ms 26,800 KB
testcase_103	AC	2,012 ms 28,220 KB
testcase_104	AC	2,393 ms 30,120 KB
testcase_105	AC	972 ms 27,344 KB
testcase_106	AC	3,603 ms 30,136 KB
testcase_107	AC	716 ms 26,920 KB
testcase_108	AC	1,932 ms 29,392 KB
testcase_109	AC	962 ms 27,040 KB
testcase_110	AC	950 ms 28,048 KB
testcase_111	AC	3,286 ms 31,716 KB
testcase_112	AC	906 ms 41,240 KB
testcase_113	AC	963 ms 41,116 KB
testcase_114	AC	966 ms 40,980 KB
testcase_115	AC	936 ms 40,856 KB
testcase_116	AC	933 ms 39,888 KB
testcase_117	AC	922 ms 40,728 KB
testcase_118	AC	962 ms 43,664 KB
testcase_119	AC	886 ms 40,016 KB
testcase_120	AC	944 ms 40,284 KB
testcase_121	AC	909 ms 40,392 KB
testcase_122	AC	918 ms 40,440 KB
testcase_123	AC	887 ms 40,528 KB
testcase_124	AC	920 ms 40,668 KB
testcase_125	AC	932 ms 47,664 KB
testcase_126	AC	949 ms 55,792 KB
testcase_127	AC	872 ms 41,048 KB
testcase_128	AC	961 ms 47,328 KB
testcase_129	AC	935 ms 47,860 KB
testcase_130	AC	978 ms 43,996 KB
testcase_131	AC	982 ms 47,448 KB
testcase_132	AC	1,566 ms 86,568 KB
testcase_133	AC	1,478 ms 84,696 KB
testcase_134	AC	1,448 ms 83,152 KB
testcase_135	AC	1,402 ms 80,740 KB
testcase_136	AC	1,407 ms 77,012 KB
testcase_137	AC	1,390 ms 76,376 KB
testcase_138	AC	1,366 ms 75,632 KB
testcase_139	AC	1,412 ms 75,044 KB
testcase_140	AC	1,442 ms 74,136 KB
testcase_141	AC	1,385 ms 72,988 KB
testcase_142	AC	1,363 ms 73,292 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

// #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;
}

template <class S, S (*op)(S, S), S (*e)()>
struct SWAG {
    std::vector<S> L, R, Lcum;
    S Rall;

    SWAG() {
        Lcum = {e()};
        Rall = e();
    }

    void push(S x) {
        R.push_back(x);
        Rall = op(Rall, x);
    }

    void pop() {
        if (L.empty()) {
            assert(!R.empty());
            while (!R.empty()) {
                S x = R.back();
                R.pop_back();
                L.push_back(x);
                Lcum.push_back(op(x, Lcum[Lcum.size() - 1]));
            }
            Rall = e();
        }
        L.pop_back();
        Lcum.pop_back();
    }

    S prod() {
        return op(Lcum[Lcum.size() - 1], Rall);
    }

    int size() {
        return L.size() + R.size();
    }

    void clear() {
        L.clear();
        R.clear();
        Lcum = {e()};
        Rall = e();
    }
};
struct SS {
    int x;
    int idx;
};
SS op(SS l, SS r) {
    return l.x < r.x ? l : r;
}
SS e() {
    return {1 << 29, -1};
}

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;
    using SW = SWAG<SS, op, e>;

    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;
                vec(SW, sw, 1);
                sw[0].push({dp[s], fr});
                for (int j = 1; fr + j * p <= r; j++) {
                    for (int ii = 0; ii < sw.size(); ii++) {
                        if (chmin(dp[fr + j * p], sw[ii].prod().x + ii + 1 + p)) {
                            bef[fr + j * p] = {sw[ii].prod().idx, p};
                        }
                    }
                    // 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] + 1 <= dp[fr] + add) {
                    //     add = 1;
                    //     nex = 10;
                    //     fr  = fr + j * p;
                    //     j   = 0;
                    // }
                    {
                        int ii = 0;
                        SS add = {dp[fr + j * p], fr + j * p};
                        int ma = 9;
                        while (1) {
                            sw[ii].push(add);
                            if (sw[ii].size() > ma) {
                                if (!sw[ii].L.empty()) {
                                    add = sw[ii].L.back();
                                } else {
                                    add = sw[ii].R[0];
                                }
                                sw[ii].pop();
                                ma *= 10;
                                ii++;
                                if (ii == sw.size()) sw.emplace_back();
                            } else {
                                break;
                            }
                        }
                    }
                }
            }
        }

        {
            int mi = dp[mid];
            int fr = mid;
            for (int i = mid - 1; i >= l; i--) {
                if (chmin(mi, dp[i] + mid - i)) {
                    fr = i;
                }
            }
            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);
    // print(len(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;
// }
//
// #include "data_structure/SWAG.hpp"
// struct SS {
//     int x;
//     int idx;
// };
// SS op(SS l, SS r) {
//     return l.x < r.x ? l : r;
// }
// SS e() {
//     return {1 << 29, -1};
// }
//
// 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;
//     using SW = SWAG<SS, op, e>;
//
//     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;
//                 vec(SW, sw, 1);
//                 sw[0].push({dp[s], fr});
//                 for (int j = 1; fr + j * p <= r; j++) {
//                     for (int ii = 0; ii < sw.size(); ii++) {
//                         if (chmin(dp[fr + j * p], sw[ii].prod().x + ii + 1 + p)) {
//                             bef[fr + j * p] = {sw[ii].prod().idx, p};
//                         }
//                     }
//                     // 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] + 1 <= dp[fr] + add) {
//                     //     add = 1;
//                     //     nex = 10;
//                     //     fr  = fr + j * p;
//                     //     j   = 0;
//                     // }
//                     {
//                         int ii = 0;
//                         SS add = {dp[fr + j * p], fr + j * p};
//                         int ma = 9;
//                         while (1) {
//                             sw[ii].push(add);
//                             if (sw[ii].size() > ma) {
//                                 if (!sw[ii].L.empty()) {
//                                     add = sw[ii].L.back();
//                                 } else {
//                                     add = sw[ii].R[0];
//                                 }
//                                 sw[ii].pop();
//                                 ma *= 10;
//                                 ii++;
//                                 if (ii == sw.size()) sw.emplace_back();
//                             } else {
//                                 break;
//                             }
//                         }
//                     }
//                 }
//             }
//         }
//
//         {
//             int mi = dp[mid];
//             int fr = mid;
//             for (int i = mid - 1; i >= l; i--) {
//                 if (chmin(mi, dp[i] + mid - i)) {
//                     fr = i;
//                 }
//             }
//             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);
//     // print(len(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;
// }

yukicoder

結果

テストケース

ソースコード