結果
問題 | No.1287 えぬけー |
ユーザー | persimmon-persimmon |
提出日時 | 2021-02-16 21:16:00 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 737 bytes |
コンパイル時間 | 236 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 365,688 KB |
最終ジャッジ日時 | 2024-09-13 10:41:12 |
合計ジャッジ時間 | 11,880 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 30 ms
10,752 KB |
testcase_01 | AC | 30 ms
10,624 KB |
testcase_02 | AC | 30 ms
10,624 KB |
testcase_03 | AC | 30 ms
10,752 KB |
testcase_04 | AC | 30 ms
10,624 KB |
testcase_05 | TLE | - |
testcase_06 | TLE | - |
testcase_07 | TLE | - |
testcase_08 | TLE | - |
ソースコード
import sys sys.setrecursionlimit(10**7) t=int(input()) cases=[list(map(int,input().split())) for _ in range(t)] mod=10**9+7 memo={} def ext_gcd(a, b): # a*x+b*y==gcd(a,b)たるgcd(a,b),x,y if (a,b) in memo:return memo[(a,b)] if b > 0: d,x,y = ext_gcd(b,a % b) memo[(a,b)]=d,y,x-(a//b)*y return d,y,x-(a//b)*y return a,1,0 for x,k in cases: # x=n^k (mod)となるnを求める。k:奇数 <=5*10**8 # x^i=n^(i*k) (mod) # ここでi*k%(mod-1)=1なら、フェルマーの小定理より # x^i=n (mod) # となる。 # i*k%(mod-1)=1となるiを探す。 # kの制約よりkとmod-1は互いに素 # ある(i,j)があり、i*k=j*(mod-1)+1 # i*k+j*(mod-1)=1 g,i,j=ext_gcd(k,mod-1) print(pow(x,i,mod))