結果

問題 No.1164 GCD Products hard
ユーザー titia
提出日時 2025-08-20 02:10:17
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,066 ms / 2,500 ms
コード長 630 bytes
コンパイル時間 841 ms
コンパイル使用メモリ 82,032 KB
実行使用メモリ 248,728 KB
最終ジャッジ日時 2025-08-20 02:10:46
合計ジャッジ時間 26,288 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

A,B,N=map(int,input().split())

x=10**7+1
L=round(x**(1/2))

Primelist=[i for i in range(x+1)]
Primelist[1]=0 # 1は素数でないので0にする.
 
for i in Primelist:
    if i>L:
        break
    if i==0:
        continue
    for j in range(2*i,x+1,i):
        Primelist[j]=0

Primes=[Primelist[j] for j in range(x+1) if Primelist[j]!=0]

mod=10**9+7
mod2=mod-1

X=[0]*(B+1)
ANS=1

for p in Primes:

    now=p    
    for i in range(B+1):
        if now>B:
            break

        ko=B//now-(A-1)//now

        X[now]=pow(ko,N,mod2)

        ANS=ANS*pow(p,X[now],mod)%mod

        now=now*p

print(ANS)
0