結果
| 問題 | No.1164 GCD Products hard |
| ユーザー |
rlangevin
|
| 提出日時 | 2023-09-10 14:54:57 |
| 言語 | PyPy3 (7.3.17) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 669 bytes |
| 記録 | |
| コンパイル時間 | 288 ms |
| コンパイル使用メモリ | 82,304 KB |
| 実行使用メモリ | 212,992 KB |
| 最終ジャッジ日時 | 2024-06-28 06:07:26 |
| 合計ジャッジ時間 | 7,736 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | -- * 2 |
| other | TLE * 1 -- * 26 |
ソースコード
from math import ceil, sqrt
def Sieve(n):
lst = [True] * (n + 1)
S = []
for i in range(2, ceil(sqrt(n)) + 1):
if lst[i]:
for j in range(2 * i, n + 1, i):
lst[j] = False
for i in range(2, n + 1):
if lst[i]:
S.append(i)
return S
A, B, N = map(int, input().split())
mod = 10**9 + 7
dp = [0] * (B + 1)
for i in range(1, B + 1):
v = B//i - (A - 1)//i
dp[i] = pow(v, N, mod - 1)
for p in Sieve(B + 1):
for j in range(2 * p, B + 1, p):
dp[j//p] -= dp[j]
ans = 1
for i in range(1, B + 1):
ans *= pow(i, dp[i]%(mod - 1), mod)
ans %= mod
print(ans)
rlangevin