結果

問題 No.3651 K-th Sum of Divisors
コンテスト
ユーザー Rumain831
提出日時 2026-08-29 01:56:25
言語 C++23(gcc16)
(gcc 16.1.0 + boost 1.90.0)
コンパイル:
g++-16 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 29 ms / 2,000 ms
+ 703µs
コード長 698 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,792 ms
コンパイル使用メモリ 179,628 KB
実行使用メモリ 28,160 KB
最終ジャッジ日時 2026-08-29 01:56:42
合計ジャッジ時間 6,440 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 55
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
using ll = long long;
ll mod = 100003;

int main(void){
  ll n, k; cin >> n >> k;
  if(k==1){
    cout << n << endl; return 0;
  }
  vector<ll> sum(mod+1);
  for(int i=1; i<=mod; i++){
    int now=i;
    while(now<=mod){
      sum[now]+=i, sum[now]%=mod;
      now+=i;
    }
  }
  ll x=0;
  for(int i=1; i<=n; i++)if(n%i==0) x+=i, x%=mod;
  k-=2;
  vector db(60, vector<int>(mod+1));
  for(int i=1; i<=mod; i++) db[0][i]=sum[i];
  for(int i=1; i<60; i++){
    for(int j=1; j<=mod; j++) db[i][j]=db[i-1][db[i-1][j]];
  }
  int c=0;
  while(k){
    if(k&1) x=db[c][x];
    c++, k/=2;
  }
  cout << x << endl;
  return 0; 
}
0