結果
問題 |
No.1164 GCD Products hard
|
ユーザー |
![]() |
提出日時 | 2023-09-10 19:35:24 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,423 bytes |
コンパイル時間 | 4,308 ms |
コンパイル使用メモリ | 252,384 KB |
最終ジャッジ日時 | 2025-02-16 21:35:41 |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 26 TLE * 1 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using namespace atcoder; using ll = long long; using P = pair<int, int>; using mint = modint1000000007; using mint2 = modint; #define rep(i, n) for (int i = 0; i < (int)(n); i++) vector<int> Sieve(int n){ vector<int> lst(n + 1, 1); int i = 2; while ((ll)i * i <= n){ if (lst[i]){ for (int j = 2 * i; j <= n; j += i){ lst[j] = 0; } } i++; } vector<int> p; for (int i = 2; i <= n; i++){ if (lst[i]){ p.push_back(i); } } return p; } ll modpow(ll v, ll N, int mod){ ll val = 1; ll now = v; while (N){ if (N & 1){ val *= now; val %= mod; } now = now * now; now %= mod; N >>= 1; } return val; } int main() { int A, B, N; cin >> A >> B >> N; int mod = 1000000007; mint2::set_mod(mod - 1); vector<mint2> dp(B + 1); for (int i = 1; i <= B; i++){ mint2 v = B/i - (A - 1)/i; dp[i] = v.pow(N); } vector<int> lst = Sieve(B + 1); for (int p : lst){ for (int j = 2 * p; j <= B; j += p){ int n = j/p; dp[n] -= dp[j]; } } mint ans = 1; for (int i = 1; i <= B; i++){ ans *= mint(i).pow((ll)dp[i].val()); } cout << ans.val() << endl; }