結果

問題 No.1396 Giri
ユーザー sten_sansten_san
提出日時 2021-02-14 23:16:46
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 439 ms / 2,000 ms
コード長 2,861 bytes
コンパイル時間 2,444 ms
コンパイル使用メモリ 215,460 KB
実行使用メモリ 198,104 KB
最終ジャッジ日時 2024-07-22 11:13:46
合計ジャッジ時間 6,431 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 439 ms
198,104 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,944 KB
testcase_05 AC 423 ms
197,976 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 2 ms
6,944 KB
testcase_08 AC 2 ms
6,940 KB
testcase_09 AC 2 ms
6,940 KB
testcase_10 AC 2 ms
6,940 KB
testcase_11 AC 2 ms
6,940 KB
testcase_12 AC 2 ms
6,944 KB
testcase_13 AC 2 ms
6,940 KB
testcase_14 AC 1 ms
6,944 KB
testcase_15 AC 2 ms
6,944 KB
testcase_16 AC 2 ms
6,944 KB
testcase_17 AC 6 ms
6,944 KB
testcase_18 AC 30 ms
18,368 KB
testcase_19 AC 213 ms
105,156 KB
testcase_20 AC 300 ms
143,460 KB
testcase_21 AC 387 ms
178,460 KB
testcase_22 AC 423 ms
195,568 KB
testcase_23 AC 427 ms
197,944 KB
testcase_24 AC 420 ms
198,104 KB
testcase_25 AC 418 ms
197,980 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

struct iofast_t {
    iofast_t() {
        ios::sync_with_stdio(false);
        cin.tie(nullptr); cout.tie(nullptr);
    }
} iofast;

struct uns_t {} uns;
template <typename Element, typename Head, typename ...Args>
auto vec(Element init, Head arg, Args ...args) {
    if constexpr (sizeof...(Args) == 0) return std::vector(arg, init);
    else return std::vector(arg, vec(init, args...));
}
template <typename Element, typename Head, typename ...Args>
auto vec(uns_t, Head arg, Args ...args) {
    return vec(Element(), arg, args...);
}

int main() {
    constexpr int64_t mod = 998244353;

    auto mpow = [&](int64_t a, int64_t b, int64_t m) {
        int64_t ans = 1;
        while (b) {
            if (b & 1) {
                ans = (ans * a) % m;
            }
            a = (a * a) % m;
            b >>= 1;
        }
        return ans;
    };

    auto minv = [&](int64_t a, int64_t m) {
        return mpow(a, m - 2, m);
    };

    int n; cin >> n;

    auto sieve = vec<int>(uns, n + 1);
    iota(begin(sieve), end(sieve), 0);

    for (int64_t i = 2; i * i <= n; ++i) {
        if (sieve[i] != i) continue;
        for (int64_t j = i * i; j <= n; j += i) {
            if (sieve[j] == j) {
                sieve[j] = i;
            }
        }
    }

    auto prime_factor = [&](int64_t a) {
        auto pf = vec<tuple<int, int>>(uns, 0);
        while (a != 1) {
            int d = sieve[a], cnt = 0;
            while (a % d == 0) {
                a /= d;
                ++cnt;
            }
            pf.emplace_back(d, cnt);
        }
        return pf;
    };

    auto primes = vec<tuple<int, int>>(uns, n + 1, 0);
    for (int i = 2; i <= n; ++i) {
        primes[i] = prime_factor(i);
    }

    auto lcm = [&]() {
        auto p = vec<int>(0, n + 1);
        for (int i = 2; i <= n; ++i) {
            for (auto [a, b] : primes[i]) {
                p[a] = max(p[a], b);
            }
        }

        int64_t ans = 1;
        for (int i = 0; i < size(p); ++i) {
            ans = (ans * mpow(i, p[i], mod)) % mod;
        }

        return ans;
    }();

    auto p = vec<array<int, 32>>(uns, n + 1);
    for (int i = 1; i <= n; ++i) {
        for (auto [a, b] : primes[i]) {
            ++p[a][b];
        }
    }

    int64_t ans = lcm, removed = 1;
    for (int i = 2; i <= n; ++i) {
        int64_t cnt = lcm, prod = 1;
        for (auto [a, b] : primes[i]) {
            if (1 < p[a][b] || accumulate(begin(p[a]) + b + 1, end(p[a]), 0)) continue;
            auto c = b - 1;
            while (0 < c && p[a][c] == 0) --c;
            prod *= mpow(a, b - c, mod);
            cnt *= minv(mpow(a, b - c, mod), mod);
            cnt %= mod;
        }

        if (removed < prod) {
            removed = prod;
            ans = cnt;
        }
    }

    cout << ans << endl;
}
0