結果

問題 No.187 中華風 (Hard)
ユーザー いい
提出日時 2023-09-21 14:20:18
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
TLE  
実行時間 -
コード長 2,828 bytes
コンパイル時間 208 ms
コンパイル使用メモリ 11,136 KB
実行使用メモリ 10,480 KB
最終ジャッジ日時 2023-09-21 14:20:49
合計ジャッジ時間 30,538 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 33 ms
10,016 KB
testcase_01 AC 36 ms
10,152 KB
testcase_02 AC 540 ms
10,272 KB
testcase_03 AC 511 ms
10,400 KB
testcase_04 AC 1,740 ms
10,288 KB
testcase_05 AC 1,663 ms
10,412 KB
testcase_06 AC 1,727 ms
10,332 KB
testcase_07 AC 1,718 ms
10,408 KB
testcase_08 TLE -
testcase_09 TLE -
testcase_10 TLE -
testcase_11 AC 1,693 ms
10,444 KB
testcase_12 AC 1,737 ms
10,416 KB
testcase_13 AC 2,339 ms
10,200 KB
testcase_14 AC 1,906 ms
10,072 KB
testcase_15 AC 352 ms
10,276 KB
testcase_16 AC 379 ms
10,284 KB
testcase_17 AC 31 ms
10,020 KB
testcase_18 AC 32 ms
10,068 KB
testcase_19 AC 31 ms
10,016 KB
testcase_20 AC 1,452 ms
10,392 KB
testcase_21 AC 33 ms
10,016 KB
testcase_22 AC 1,726 ms
10,480 KB
testcase_23 AC 31 ms
10,088 KB
testcase_24 AC 31 ms
10,016 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys

V = sys.version_info
_39 = False
_310 = False
_311 = False
if V.major == 3:
    _39 = V.minor >= 9
    _310 = V.minor >= 10
    _311 = V.minor >= 11

if _39:
    li = list
    tup = tuple
    dic = dict
    st = set
    ty = type
else:
    from typing import (
        List,
        Tuple,
        Type,
        Dict,
        Set,
    )

    li = List
    tup = Tuple
    dic = Dict
    st = Set
    ty = Type

from sys import stdin

read = stdin.buffer.read
rl = stdin.buffer.readline
rb = lambda: rl().split()
rls = stdin.buffer.readlines


from typing import (
    TypeVar,
    Sequence as Seq,
    Iterable as Iter,
    Protocol as Proto,
    Generic as Gen,
)
from typing import Callable as Fn


from typing import Iterable


def prints(
    a: Iterable[object],
    sep: str = "\n",
) -> None:
    print(sep.join(map(str, a)))


# import numpy as np


def egcd(
    a: int,
    b: int,
) -> tup[int, int, int]:
    if not b:
        if a < 0:
            return -a, -1, 0
        return a, 1, 0
    q, r = divmod(a, b)
    g, s, t = egcd(b, r)
    return g, t, s - q * t


def crt_mod(
    m: Seq[int],
    r: Seq[int],
    mod: int,
) -> int:
    l = len(m)
    assert len(r) == l
    m = list(m) + [mod]
    p = [1] * (l + 1)
    x = [0] * (l + 1)
    for i, (n, b) in enumerate(zip(m, r)):
        assert n > 0
        g, ip, _ = egcd(p[i], n)
        assert g == 1
        t = (b - x[i]) % n * ip % n
        for j in range(i + 1, l + 1):
            x[j] += t * p[j]
            x[j] %= m[j]
            p[j] = p[j] * n % m[j]
    return x[-1]


def factorize(n: int) -> li[tup[int, int]]:
    f = []
    for i in range(2, n):
        if i * i > n:
            break
        if n % i:
            continue
        c = 0
        while n % i == 0:
            n //= i
            c += 1
        f.append((i, c))
    if n > 1:
        f.append((n, 1))
    return f


def to_coprime(mr: li[tup[int, int]]):
    f = dict()
    for m, r in mr:
        for p, e in factorize(m):
            n = pow(p, e)
            b = r % n
            l, a = f.get(p, (1, 0))
            if l < n:
                l, n = n, l
                a, b = b, a
            if a % n != b:
                mr.clear()
                return
            f[p] = (l, a)
    mr.clear()
    for x in f.values():
        mr.append(x)


def solve() -> None:
    n = int(rl())
    mr = []
    for _ in range(n):
        r, m = map(int, rb())
        mr.append((m, r))
    to_coprime(mr)
    n = len(mr)
    if n == 0:
        print(-1)
        return
    m, r = tuple(zip(*mr))
    mod = 10**9 + 7
    if max(r) > 0:
        x = crt_mod(m, r, mod)
        print(x)
        return
    l = 1
    for x in m:
        l = l * x % mod
    print(l)


def main() -> None:
    t = 1
    # t = int(rl())
    for _ in range(t):
        solve()


main()
0