結果
問題 | No.186 中華風 (Easy) |
ユーザー | Ueeek |
提出日時 | 2021-03-02 13:57:29 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,558 bytes |
コンパイル時間 | 115 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 10,752 KB |
最終ジャッジ日時 | 2024-10-03 01:43:15 |
合計ジャッジ時間 | 1,678 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 30 ms
10,624 KB |
testcase_01 | AC | 30 ms
10,624 KB |
testcase_02 | AC | 30 ms
10,752 KB |
testcase_03 | AC | 30 ms
10,624 KB |
testcase_04 | AC | 30 ms
10,624 KB |
testcase_05 | AC | 30 ms
10,624 KB |
testcase_06 | AC | 30 ms
10,624 KB |
testcase_07 | AC | 30 ms
10,624 KB |
testcase_08 | AC | 30 ms
10,624 KB |
testcase_09 | AC | 30 ms
10,624 KB |
testcase_10 | AC | 30 ms
10,624 KB |
testcase_11 | AC | 30 ms
10,752 KB |
testcase_12 | AC | 30 ms
10,624 KB |
testcase_13 | AC | 31 ms
10,624 KB |
testcase_14 | AC | 30 ms
10,624 KB |
testcase_15 | AC | 29 ms
10,624 KB |
testcase_16 | WA | - |
testcase_17 | WA | - |
testcase_18 | AC | 29 ms
10,624 KB |
testcase_19 | AC | 29 ms
10,752 KB |
testcase_20 | AC | 30 ms
10,496 KB |
testcase_21 | AC | 30 ms
10,624 KB |
testcase_22 | AC | 30 ms
10,752 KB |
ソースコード
#/ref https://atcoder.jp/contests/abc193/submissions/20516429 def _inv_gcd(a: int, b: int) -> tuple: a %= b if a == 0: return(b, 0) s = b t = a m0 = 0 m1 = 1 while t: u = s//t s -= t*u m0 -= m1*u s, t = t, s m0, m1 = m1, m0 if m0 < 0: m0 += b//s return(s, m0) def CRT(R: list, M: list): """ Chinese Remainder Theorem Arguments: R: list of remainders M: list of modulo all_same(x%ri(mod mi) for i in range(len(r))) Return: x: s.t x==ri(mod mi) for all i in range(len(r))(0<x<lcm) lcm: lcm of list m(if 0, anser is not exisis) """ assert len(R) == len(M) N = len(R) #2-var CRTを繰り返し適応していく。一つ目を(r,m)=(0,1)としておく。(単位元的な) r0 = 0 m0 = 1 #Contracts: 0<=r0<m0 for i in range(N): assert 1 <= M[i] r1 = R[i] % M[i] m1 = M[i] if m0 < m1: r0, r1 = r1, r0 m0, m1 = m1, m0 if m0 % m1 == 0: if r0 % m1 != r1: return (0, 0) continue g, im = _inv_gcd(m0, m1) u1 = m1//g if (r1-r0) % g: return (0, 0) x = (r1-r0)//g % u1*im % u1 r0 += x*m0 m0 *= u1 if r0 < 0: r0 += m0 return (r0, m0) X=[] Y=[] for _ in range(3): x,y = map(int,input().split()) X.append(x) Y.append(y) x,lcm = CRT(X,Y) if lcm==0: print(-1) else: print(x)