結果
| 問題 |
No.1164 GCD Products hard
|
| ユーザー |
rlangevin
|
| 提出日時 | 2023-09-10 14:56:25 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 672 bytes |
| コンパイル時間 | 362 ms |
| コンパイル使用メモリ | 82,304 KB |
| 実行使用メモリ | 171,392 KB |
| 最終ジャッジ日時 | 2024-06-28 06:08:45 |
| 合計ジャッジ時間 | 7,825 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / 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//2 + 5):
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