結果
問題 | No.1339 循環小数 |
ユーザー | 👑 SPD_9X2 |
提出日時 | 2021-01-23 17:36:22 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 470 ms / 2,000 ms |
コード長 | 1,570 bytes |
コンパイル時間 | 735 ms |
コンパイル使用メモリ | 87,260 KB |
実行使用メモリ | 131,664 KB |
最終ジャッジ日時 | 2023-08-30 02:09:28 |
合計ジャッジ時間 | 9,693 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 72 ms
71,812 KB |
testcase_01 | AC | 75 ms
76,212 KB |
testcase_02 | AC | 76 ms
76,092 KB |
testcase_03 | AC | 76 ms
76,092 KB |
testcase_04 | AC | 77 ms
76,144 KB |
testcase_05 | AC | 75 ms
76,288 KB |
testcase_06 | AC | 74 ms
76,220 KB |
testcase_07 | AC | 76 ms
76,212 KB |
testcase_08 | AC | 76 ms
76,232 KB |
testcase_09 | AC | 75 ms
76,212 KB |
testcase_10 | AC | 77 ms
75,848 KB |
testcase_11 | AC | 80 ms
76,152 KB |
testcase_12 | AC | 82 ms
76,516 KB |
testcase_13 | AC | 79 ms
76,728 KB |
testcase_14 | AC | 82 ms
76,552 KB |
testcase_15 | AC | 80 ms
76,736 KB |
testcase_16 | AC | 81 ms
76,376 KB |
testcase_17 | AC | 81 ms
76,576 KB |
testcase_18 | AC | 82 ms
76,712 KB |
testcase_19 | AC | 81 ms
76,604 KB |
testcase_20 | AC | 79 ms
76,488 KB |
testcase_21 | AC | 291 ms
108,524 KB |
testcase_22 | AC | 321 ms
120,328 KB |
testcase_23 | AC | 312 ms
112,584 KB |
testcase_24 | AC | 300 ms
115,000 KB |
testcase_25 | AC | 312 ms
114,708 KB |
testcase_26 | AC | 319 ms
117,340 KB |
testcase_27 | AC | 318 ms
115,032 KB |
testcase_28 | AC | 326 ms
118,876 KB |
testcase_29 | AC | 279 ms
114,456 KB |
testcase_30 | AC | 304 ms
112,528 KB |
testcase_31 | AC | 426 ms
131,664 KB |
testcase_32 | AC | 439 ms
130,340 KB |
testcase_33 | AC | 315 ms
111,180 KB |
testcase_34 | AC | 222 ms
102,432 KB |
testcase_35 | AC | 470 ms
130,560 KB |
testcase_36 | AC | 308 ms
114,844 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)