結果

問題 No.2589 Prepare Integers
コンテスト
ユーザー fdironia
提出日時 2023-12-26 03:19:09
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,310 ms / 2,000 ms
コード長 2,641 bytes
コンパイル時間 1,481 ms
コンパイル使用メモリ 121,312 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-27 14:53:52
合計ジャッジ時間 7,375 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 24
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <cmath>
#include <vector>
#include <tuple>
#include <algorithm>

using namespace std;

using ll = long long;
using vi = vector<int>;
using v2i = vector<vi>;
using t4i = tuple<int, int, int, int>;

constexpr int X = 1e9;

t4i exgcd(int a, int b) {
    if (b == 0) {
        return {1, 0, 0, 1};
    }

    auto [k1, k2, k3, k4] = exgcd(b, a % b);
    return {k2, k1 - a / b * k2, k4, k3 - a / b * k4};
}

ll solve(const v2i &base, int k, ll x) {
    x--;
    int n = base.size();

    int xd[n];
    for (int i = 0; i < n; i++) {
        if (!base[i].empty()) {
            xd[i] = x % (k / base[i][i]);
            x /= k / base[i][i];
        }
    }

    if (x > 0) {
        return -1;
    }

    int ans[n];
    fill(ans, ans + n, 0);
    for (int i = n - 1; i >= 0; i--) {
        if (!base[i].empty()) {
            int t = ((xd[i] - ans[i] / base[i][i]) % k + k) % k;
            for (int j = 0; j < n; j++) {
                ans[j] = (ans[j] + 1ll * t * base[i][j]) % k;
            }
        }
    }

    ll t = 0;
    for (int i = n - 1; i >= 0; i--) {
        t = t * k + ans[i];
    }

    return t;
}

int main() {
    int k, q;
    cin >> k >> q;
    int n = log(X) / log(k) + 1;

    v2i base(n);

    while (q--) {
        int op;
        ll x;
        cin >> op >> x;

        if (op == 1) {
            int xd[n];
            for (int i = 0; i < n; i++) {
                xd[i] = x % k;
                x /= k;
            }

            for (int i = n - 1; i >= 0; i--) {
                if (xd[i] == 0) {
                    continue;
                }

                auto [k1, k2, k3, k4] = exgcd(xd[i], base[i].empty() ? k : base[i][i]);

                if (base[i].empty()) {
                    base[i] = vi(n);
                }

                int t[n];
                copy(base[i].begin(), base[i].end(), t);

                for (int j = 0; j < n; j++) {
                    base[i][j] = ((1ll * k1 * xd[j] + 1ll * k2 * t[j]) % k + k) % k;
                }
                for (int j = 0; j < n; j++) {
                    xd[j] = ((1ll * k3 * xd[j] + 1ll * k4 * t[j]) % k + k) % k;
                }
            }
        } else if (op == 2) {
            cout << solve(base, k, x) << endl;
        } else {
            ll l = 0, r = x + 2;
            while (r - l > 1) {
                ll mid = (l + r) / 2;

                ll t = solve(base, k, mid);
                if (t <= x && t != -1) {
                    l = mid;
                } else {
                    r = mid;
                }
            }

            cout << l << endl;
        }
    }

    return 0;
}
0