結果

問題 No.3651 K-th Sum of Divisors
コンテスト
ユーザー shingo0909
提出日時 2026-08-28 22:38:33
言語 C++23
(gcc 15.2.0 + boost 1.90.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 56 ms / 2,000 ms
+ 957µs
コード長 2,927 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,257 ms
コンパイル使用メモリ 343,980 KB
実行使用メモリ 40,636 KB
最終ジャッジ日時 2026-08-28 22:38:53
合計ジャッジ時間 7,830 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 55
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define rep(i, n) for (int i = 0; i < (int)(n); i++)

// https://github.com/shingo0909/kyopro/blob/39029fab39c1735c38e5e9e4494d9934f81d6918/Library/doubling.cpp

struct doubling {
    vector<vector<int>> d;
    int n;
    doubling(vector<int> p) : d(60, vector<int>(p.size())) {
        n = p.size();
        d[0] = p;
        for (int i = 1; i < 60; i++) {
            for (int j = 0; j < n; j++) {
                d[i][j] = d[i - 1][d[i - 1][j]];
            }
        }
    }
    int get_ans(int x, long long t) {
        for (int i = 0; i < 60; i++) {
            if ((t >> i) % 2) {
                x = d[i][x];
            }
        }
        return x;
    }
};

// https://github.com/shingo0909/kyopro/blob/39029fab39c1735c38e5e9e4494d9934f81d6918/Library/PrimeTable.cpp

struct PrimeTable {
  private:
    int n;
    vector<int> spf;
    vector<int> primes;
    void build() {
        spf[0] = -1;
        spf[1] = -1;

        for (int i = 2; i <= n; i++) {
            if (spf[i] == 0) {
                spf[i] = i;
                primes.push_back(i);
            }
            for (int p : primes) {
                if (p > spf[i] || 1LL * p * i > n)
                    break;
                spf[p * i] = p;
            }
        }
        return;
    }

  public:
    PrimeTable(int n) : n(n), spf(n + 1, 0) {
        build();
    }

    bool is_prime(int x) const {
        return x >= 2 && spf[x] == x;
    }

    // 素因数分解
    vector<pair<int, int>> factorize(int x) const {
        vector<pair<int, int>> res;
        while (x > 1) {
            int p = spf[x];
            int cnt = 0;
            while (x % p == 0) {
                x /= p;
                cnt++;
            }
            res.emplace_back(p, cnt);
        }
        return res;
    }

    // 約数列挙
    // ソートされていないことに注意
    vector<int> divisors(int n) const {
        vector<int> res{1};
        auto pf = factorize(n);
        for (auto p : pf) {
            int s = (int)res.size();
            for (int i = 0; i < s; i++) {
                int v = 1;
                for (int j = 0; j < p.second; j++) {
                    v *= p.first;
                    res.push_back(res[i] * v);
                }
            }
        }
        return res;
    }
};

int main() {
    cin.tie(nullptr);
    ios_base::sync_with_stdio(false);
    ll n, k;
    cin >> n >> k;
    if (k == 1) {
        cout << n << endl;
        return 0;
    }
    k -= 2;
    int m = 100003;
    PrimeTable p(3000010);
    vector<int> g(m, 0);
    rep(i, m) {
        if (i == 0)
            continue;
        ll s = 0;
        for (auto j : p.divisors(i))
            s += j;
        g[i] = s % m;
    }
    doubling d(g);
    ll f = 0;
    for (auto i : p.divisors(n))
        f += i;
    f %= m;
    cout << d.get_ans(f, k) << endl;
    return 0;
}
0