結果

問題 No.3607 Sum of Powers of GCDs
コンテスト
ユーザー tau1235
提出日時 2026-07-31 23:18:11
言語 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  
実行時間 623 ms / 2,500 ms
+ 952µs
コード長 1,167 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,417 ms
コンパイル使用メモリ 345,412 KB
実行使用メモリ 93,248 KB
最終ジャッジ日時 2026-07-31 23:18:21
合計ジャッジ時間 9,281 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 11
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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

int mm=1e6;
int kk=10;
using mint=atcoder::modint998244353;
vector<vector<mint>> a;
vector<vector<mint>> sa;

template< typename T >
vector< pair< pair< T, T >, T > > quotient_range(T N) {
  T M;
  vector< pair< pair< T, T >, T > > ret;
  for(M = 1; M * M <= N; M++) {
    ret.emplace_back(make_pair(M, M), N / M);
  }
  for(T i = M; i >= 1; i--) {
    T L = N / (i + 1) + 1;
    T R = N / i;
    if(L <= R && ret.back().first.second < L) ret.emplace_back(make_pair(L, R), N / L);
  }
  return ret;
}

void f(int k){
  vector<mint> dp(mm+1);
  for (int i=1;i<=mm;i++){
    dp[i]+=mint(i).pow(k);
    for (int j=i*2;j<=mm;j+=i){
      dp[j]-=dp[i];
    }
  }
  a[k]=dp;
  for (int i=1;i<=mm;i++) sa[k][i]=sa[k][i-1]+dp[i];
}

void solve(){
  int n,m,k;
  cin>>n>>m>>k;
  auto vp=quotient_range(m);
  mint ans=0;
  for (auto [p,q]:vp){
    auto [l,r]=p;
    ans+=(sa[k][r]-sa[k][l-1])*mint(q).pow(n);
  }
  cout<<ans.val()<<endl;
}

int main(){
  a=vector(kk+1,vector<mint>(mm+1));
  sa=vector(kk+1,vector<mint>(mm+1));
  for (int k=1;k<=kk;k++) f(k);
  int t=1;
  cin>>t;
  while (t--) solve();
}
0