結果
| 問題 | 
                            No.1164 GCD Products hard
                             | 
                    
| ユーザー | 
                             rlangevin
                         | 
                    
| 提出日時 | 2023-09-10 14:54:57 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                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