結果

問題 No.1396 Giri
ユーザー merom686
提出日時 2021-02-14 21:57:33
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 121 ms / 2,000 ms
コード長 2,122 bytes
コンパイル時間 1,009 ms
コンパイル使用メモリ 99,072 KB
最終ジャッジ日時 2025-01-18 20:58:57
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <string>
#include <algorithm>
#include <array>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
using ll = long long;

constexpr int P = 998244353;

ll powmod(ll n, ll k) {
    ll r = 1, t = n % P;
    for (; k != 0; k /= 2) {
        if (k & 1) r = r * t % P;
        t = t * t % P;
    }
    return r;
}
ll modinv(ll n) {
    return powmod(n, P - 2);
}

struct Sieve {
    Sieve(int n) : l(n + 1, 0) {
        for (int i = 2; i <= n; i++) {
            if (!l[i]) p.push_back(l[i] = i);
            for (int q : p) {
                if (q > l[i]) break;
                if (int j = i * q; j <= n) l[j] = q; else break;
            }
        }
    }
    bool is_prime(int i) {
        return l[i] == i;
    }
    vector<int> l, p;
};

int main() {
    int n;
    cin >> n;

    Sieve si(n);

    vector<array<int, 2>> x(n + 1);
    vector<array<array<int, 2>, 8>> y(n + 1);

    for (int i = 1; i <= n; i++) {
        int k = 0, l = 0, m = 0;
        for (int t = i; t != 1;) {
            int p = si.l[t];
            t /= p;
            if (p != l) {
                if (k) y[i][m++] = { l, k };
                k = 1;
                l = p;
            } else {
                k++;
            }
        }
        if (k) y[i][m++] = { l, k };

        for (int j = 0; j < m; j++) {
            auto [l, k] = y[i][j];
            if (k > x[l][1]) {
                if (k > x[l][0]) {
                    x[l][1] = x[l][0];
                    x[l][0] = k;
                } else {
                    x[l][1] = k;
                }
            }
        }
    }

    ll a0 = 1, a1 = 1;
    for (int i = 1; i <= n; i++) {
        if (x[i][0]) a0 = a0 * powmod(i, x[i][0]) % P;
    }
    for (int i = 1; i <= n; i++) {
        ll t = 1;
        for (int j = 0; j < 8; j++) {
            auto [l, k] = y[i][j];
            if (!k) break;
            if (k == x[l][0]) {
                t *= powmod(l, x[l][0] - x[l][1]);
            }
        }
        a1 = max(a1, t);
    }

    ll r = a0 * modinv(a1) % P;
    cout << r << endl;

    return 0;
}
0