結果

問題 No.187 中華風 (Hard)
ユーザー akasia_midoriakasia_midori
提出日時 2023-11-23 03:02:31
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 157 ms / 3,000 ms
コード長 1,856 bytes
コンパイル時間 347 ms
コンパイル使用メモリ 81,700 KB
実行使用メモリ 76,648 KB
最終ジャッジ日時 2023-11-23 03:02:35
合計ジャッジ時間 3,809 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 64 ms
68,040 KB
testcase_01 AC 62 ms
68,040 KB
testcase_02 AC 99 ms
76,164 KB
testcase_03 AC 109 ms
76,648 KB
testcase_04 AC 157 ms
76,648 KB
testcase_05 AC 121 ms
76,164 KB
testcase_06 AC 112 ms
76,172 KB
testcase_07 AC 120 ms
76,164 KB
testcase_08 AC 133 ms
76,408 KB
testcase_09 AC 133 ms
76,536 KB
testcase_10 AC 133 ms
76,408 KB
testcase_11 AC 120 ms
76,164 KB
testcase_12 AC 111 ms
76,172 KB
testcase_13 AC 66 ms
68,040 KB
testcase_14 AC 65 ms
68,040 KB
testcase_15 AC 76 ms
72,268 KB
testcase_16 AC 77 ms
72,396 KB
testcase_17 AC 44 ms
55,728 KB
testcase_18 AC 58 ms
65,844 KB
testcase_19 AC 42 ms
55,728 KB
testcase_20 AC 107 ms
76,164 KB
testcase_21 AC 43 ms
55,728 KB
testcase_22 AC 118 ms
76,164 KB
testcase_23 AC 43 ms
55,728 KB
testcase_24 AC 42 ms
55,728 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import functools
import math

def lcm2(x, y):
    gcd = math.gcd(x, y)
    return gcd * (x // gcd) * (y // gcd)

gcd = lambda l: functools.reduce(math.gcd, l)
lcm = lambda l: functools.reduce(lcm2, l)
import math

# Extended Euclidean Algorithm
def extgcd(a, b):
    if b:
        d, y, x = extgcd(b, a % b)
        y -= (a // b)*x
        return d, x, y
    return a, 1, 0

def Chineese_mod_theory(b, m):
    # 求めたいxがあって m1とm2が互いに疎であり、x%m1≡b1, x%m2≡b2であるとき  
    # x = b1m2 + b2m1 である
    # 例えば b1m2%m1≡b1, b2m1%m1≡0である
    # これは拡張ユークリッドの互除法で考えると 
    # m1p1 + m2p2 = d(=gcd(m1, m2))
    # s = (b1-b2)/d　と置くと
    # m1p1s + m2p2s = b1-b2
    # x = b2 + m2p2s,  x = b1 - m1p1s
    # のどっちかを解けばいい
    MOD = 10**9 + 7
    M = 1
    r = 0
    for i in range(len(b)):
        d, p1, p2 = extgcd(M, m[i])
        if (b[i]-r)%d != 0:
            return (-1,0)
        s = (b[i]-r)//d
        tmp = (s * p1) % (m[i]//d)

        r = (r + M*tmp)
        M = (M * (m[i]//d))
    return (r%M)%MOD, M%MOD



def oi(): return int(input())
def os(): return input()
def mi(): return list(map(int, input().split()))

# import sys
# input = sys.stdin.readline
# import sys
# sys.setrecursionlimit(10**8)
# import pypyjit
# pypyjit.set_param('max_unroll_recursion=-1')
# input_count = 0

input_count = 0
N = oi()
amari = []
wari = []
all_zero = True
for _ in range(N):
    x,y = mi()
    wari.append(y)
    amari.append(x)
    if x != 0:
        all_zero = False
if all_zero:
    # amari, wari = pre_garner(amari, wari)
    out = lcm(wari)
    print(out%(10**9+7))
else:
    # amari, wari = pre_garner(amari, wari)
    # if amari is None:
        # print(-1)
    # else:
    out = Chineese_mod_theory(amari, wari)
    print(out[0])
0