結果
問題 | No.1728 [Cherry 3rd Tune] Bullet |
ユーザー |
![]() |
提出日時 | 2021-10-29 22:31:28 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 135 ms / 2,000 ms |
コード長 | 837 bytes |
コンパイル時間 | 4,954 ms |
コンパイル使用メモリ | 255,204 KB |
最終ジャッジ日時 | 2025-01-25 09:11:49 |
ジャッジサーバーID (参考情報) |
judge3 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 27 |
ソースコード
#include <stdio.h> #include <bits/stdc++.h> #include <atcoder/all> using namespace atcoder; using mint = modint1000000007; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf (long long)1000000000000000000 void solve(){ long long N,C; cin>>N>>C; vector<long long> y; for(long long i=1;i*i<=N;i++){ if(N%i==0){ y.push_back(i); if(i*i!=N)y.push_back(N/i); } } sort(y.rbegin(),y.rend()); vector<long long> cnt(y.size()); rep(i,y.size()){ cnt[i] = N / y[i]; rep(j,i){ if(y[j] % y[i] == 0)cnt[i] -= cnt[j]; } } mint ans = 0LL; ans += mint(C).pow(N); ans *= N; rep(i,cnt.size()){ mint temp = mint(C).pow(y[i]*2); temp *= cnt[i]; ans += temp; } ans /= N*2; cout<<ans.val()<<endl; } int main(){ int _t; cin>>_t; rep(_,_t){ solve(); } return 0; }