結果
問題 | No.1339 循環小数 |
ユーザー | 👑 SPD_9X2 |
提出日時 | 2021-01-23 17:36:22 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 423 ms / 2,000 ms |
コード長 | 1,570 bytes |
コンパイル時間 | 374 ms |
コンパイル使用メモリ | 82,124 KB |
実行使用メモリ | 129,276 KB |
最終ジャッジ日時 | 2024-06-10 01:16:40 |
合計ジャッジ時間 | 6,744 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 37 ms
53,328 KB |
testcase_01 | AC | 53 ms
58,936 KB |
testcase_02 | AC | 38 ms
60,128 KB |
testcase_03 | AC | 37 ms
59,244 KB |
testcase_04 | AC | 37 ms
60,888 KB |
testcase_05 | AC | 38 ms
59,552 KB |
testcase_06 | AC | 38 ms
59,956 KB |
testcase_07 | AC | 38 ms
59,752 KB |
testcase_08 | AC | 38 ms
59,836 KB |
testcase_09 | AC | 39 ms
59,000 KB |
testcase_10 | AC | 38 ms
60,032 KB |
testcase_11 | AC | 42 ms
62,152 KB |
testcase_12 | AC | 43 ms
63,500 KB |
testcase_13 | AC | 42 ms
63,116 KB |
testcase_14 | AC | 43 ms
63,100 KB |
testcase_15 | AC | 43 ms
63,636 KB |
testcase_16 | AC | 43 ms
63,600 KB |
testcase_17 | AC | 43 ms
63,844 KB |
testcase_18 | AC | 44 ms
63,152 KB |
testcase_19 | AC | 43 ms
63,580 KB |
testcase_20 | AC | 42 ms
62,904 KB |
testcase_21 | AC | 251 ms
109,964 KB |
testcase_22 | AC | 279 ms
116,300 KB |
testcase_23 | AC | 277 ms
115,280 KB |
testcase_24 | AC | 251 ms
113,132 KB |
testcase_25 | AC | 271 ms
112,496 KB |
testcase_26 | AC | 272 ms
113,020 KB |
testcase_27 | AC | 262 ms
113,112 KB |
testcase_28 | AC | 275 ms
115,356 KB |
testcase_29 | AC | 241 ms
112,204 KB |
testcase_30 | AC | 259 ms
110,000 KB |
testcase_31 | AC | 388 ms
129,276 KB |
testcase_32 | AC | 384 ms
129,144 KB |
testcase_33 | AC | 269 ms
112,272 KB |
testcase_34 | AC | 174 ms
99,424 KB |
testcase_35 | AC | 423 ms
128,908 KB |
testcase_36 | AC | 261 ms
112,616 KB |
ソースコード
""" https://yukicoder.me/problems/no/1339 循環節は、 10^K = 1 (mod N) となるKを求めれば行ける K == 0 が満たしてしまうので 10^(K-1) == 10^(-1) (mod N) となるK-1を求めてあげればよい """ import sys import math from sys import stdin #aのmodを法にした逆元を返す #法が素数でなくても利用可能 #aとmodが互いに素でないならば、-1を返す import math def extGCD(a,b): g = math.gcd(a,b) x, y, u, v = 1, 0, 0, 1 while b: k = a // b x -= k * u y -= k * v x, u = u, x y, v = v, y a, b = b, a % b return g ,x, y def exinv(a,mod): g,x,y = extGCD(a,mod) if g != 1: return -1 if x > 0: return x else: return x + mod # X^K = Y (mod M) from math import ceil, sqrt def BsGs(X,Y,M): dic = {} dic[1] = 0 sq = ceil(M**0.5) #baby-step Z = 1 for i in range(sq): Z = Z * X % M if Z not in dic: dic[Z] = i+1 if Y in dic: return dic[Y] #giant-step R = exinv(Z,M) for i in range(1,sq+1): Y = Y * R % M if Y in dic: return dic[Y] + i * sq return -1 #print (BsGs(10,247,2469)) TT = int(stdin.readline()) for loop in range(TT): N = int(stdin.readline()) while N % 2 == 0: N //= 2 while N % 5 == 0: N //= 5 if N == 1: print (1) else: r = 10 #print (exinv(r,N),N,file=sys.stderr) k = BsGs(r,exinv(r,N),N) print (k+1)