結果
| 問題 | 
                            No.1164 GCD Products hard
                             | 
                    
| ユーザー | 
                             rlangevin
                         | 
                    
| 提出日時 | 2024-04-14 11:34:57 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 1,060 ms / 2,500 ms | 
| コード長 | 626 bytes | 
| コンパイル時間 | 289 ms | 
| コンパイル使用メモリ | 82,428 KB | 
| 実行使用メモリ | 244,792 KB | 
| 最終ジャッジ日時 | 2024-10-03 08:34:46 | 
| 合計ジャッジ時間 | 18,125 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge3 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 2 | 
| other | AC * 27 | 
ソースコード
from math import ceil, sqrt
def Sieve(n):
    lst = [True] * (n + 1)
    lst[0] = lst[1] = False
    S = set()
    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.add(i)
    return S
A, B, N = map(int, input().split())
mod = 10**9 + 7
ans = 1
for p in Sieve(B):
    now = p
    cnt = 0
    while now <= B:
        v = B // now - (A - 1) // now
        cnt += pow(v, N, mod - 1)
        cnt %= mod - 1
        now *= p
    ans *= pow(p, cnt, mod)
    ans %= mod
    
print(ans)
            
            
            
        
            
rlangevin