結果

問題 No.2594 Mix shake!!
ユーザー hitonanodehitonanode
提出日時 2023-12-20 12:15:08
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
WA  
実行時間 -
コード長 2,021 bytes
コンパイル時間 323 ms
コンパイル使用メモリ 13,056 KB
実行使用メモリ 11,904 KB
最終ジャッジ日時 2024-09-27 09:46:17
合計ジャッジ時間 6,588 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 36 ms
11,520 KB
testcase_01 AC 37 ms
11,520 KB
testcase_02 AC 37 ms
11,648 KB
testcase_03 AC 37 ms
11,520 KB
testcase_04 AC 38 ms
11,520 KB
testcase_05 WA -
testcase_06 AC 37 ms
11,520 KB
testcase_07 AC 40 ms
11,648 KB
testcase_08 RE -
testcase_09 WA -
testcase_10 WA -
testcase_11 AC 37 ms
11,520 KB
testcase_12 AC 36 ms
11,648 KB
testcase_13 AC 37 ms
11,520 KB
testcase_14 AC 37 ms
11,520 KB
testcase_15 AC 37 ms
11,520 KB
testcase_16 AC 37 ms
11,520 KB
testcase_17 AC 37 ms
11,520 KB
testcase_18 AC 37 ms
11,648 KB
testcase_19 AC 38 ms
11,520 KB
testcase_20 AC 37 ms
11,648 KB
testcase_21 AC 36 ms
11,520 KB
testcase_22 WA -
testcase_23 AC 39 ms
11,520 KB
testcase_24 AC 37 ms
11,520 KB
testcase_25 AC 36 ms
11,648 KB
testcase_26 WA -
testcase_27 AC 36 ms
11,648 KB
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 AC 38 ms
11,648 KB
testcase_33 AC 39 ms
11,648 KB
testcase_34 AC 65 ms
11,904 KB
testcase_35 WA -
testcase_36 AC 66 ms
11,776 KB
testcase_37 WA -
testcase_38 AC 67 ms
11,776 KB
testcase_39 WA -
testcase_40 AC 65 ms
11,776 KB
testcase_41 WA -
testcase_42 AC 67 ms
11,776 KB
testcase_43 WA -
testcase_44 AC 66 ms
11,776 KB
testcase_45 AC 59 ms
11,776 KB
testcase_46 AC 65 ms
11,904 KB
testcase_47 WA -
testcase_48 AC 66 ms
11,904 KB
testcase_49 WA -
testcase_50 AC 65 ms
11,776 KB
testcase_51 AC 60 ms
11,776 KB
testcase_52 AC 65 ms
11,904 KB
testcase_53 WA -
testcase_54 AC 38 ms
11,648 KB
testcase_55 AC 37 ms
11,520 KB
testcase_56 WA -
testcase_57 AC 36 ms
11,648 KB
testcase_58 AC 36 ms
11,520 KB
testcase_59 AC 36 ms
11,520 KB
testcase_60 AC 37 ms
11,520 KB
testcase_61 AC 36 ms
11,648 KB
testcase_62 AC 36 ms
11,648 KB
testcase_63 AC 38 ms
11,520 KB
testcase_64 AC 36 ms
11,520 KB
testcase_65 AC 37 ms
11,520 KB
testcase_66 AC 45 ms
11,648 KB
testcase_67 AC 38 ms
11,520 KB
testcase_68 AC 44 ms
11,648 KB
testcase_69 AC 37 ms
11,648 KB
testcase_70 AC 41 ms
11,648 KB
testcase_71 AC 37 ms
11,648 KB
testcase_72 AC 37 ms
11,520 KB
testcase_73 AC 36 ms
11,520 KB
testcase_74 AC 37 ms
11,648 KB
testcase_75 AC 37 ms
11,520 KB
testcase_76 AC 37 ms
11,520 KB
testcase_77 AC 36 ms
11,520 KB
testcase_78 AC 36 ms
11,520 KB
testcase_79 AC 36 ms
11,520 KB
testcase_80 AC 47 ms
11,648 KB
testcase_81 AC 37 ms
11,776 KB
testcase_82 AC 40 ms
11,648 KB
testcase_83 AC 38 ms
11,648 KB
testcase_84 AC 38 ms
11,648 KB
testcase_85 AC 39 ms
11,648 KB
testcase_86 AC 38 ms
11,776 KB
testcase_87 AC 37 ms
11,648 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

# バグりまくり
from fractions import Fraction


def to_slope(vec: list[tuple[Fraction, int, int, int]]) -> list[tuple[int, int]]:
    x, y = 0, 0
    ret: list[tuple[int, int]] = []
    for _, a, b, _ in vec:
        x += a
        y += b
        ret.append((x, y))
    return ret


def solve() -> bool:

    n = int(input())

    A = list(map(int, input().split()))
    B = list(map(int, input().split()))
    C = list(map(int, input().split()))
    D = list(map(int, input().split()))

    AB = [(Fraction(b, a), a, b, i) for i, (a, b) in enumerate(zip(A, B))]
    CD = [(Fraction(d, c), c, d, i) for i, (c, d) in enumerate(zip(C, D))]

    AB.sort()
    CD.sort()

    lb_xys = to_slope(AB)  # 初期状態の y = f(x) の折れ線
    ub_xys = to_slope(CD)  # 所望の最終状態の y = f(x) の折れ線

    # 最終状態の折れ線が初期状態の折れ線の下側を通っていたら明らかに No
    for (x0, y0), (x1, y1) in zip(lb_xys[:-1], lb_xys[1:], strict=True):
        for x, y in ub_xys:
            if (x1 - x0) * (y - y0) - (y1 - y0) * (x - x0) < 0:
                return False

    # 初期状態と最終状態で濃度の大小関係が一致していたら Yes
    if [i for _, _, _, i in AB] == [i for _, _, _, i in CD]:
        return True

    # n - 1 本以下の線分で (0, 0) から右上まで到達できたら Yes
    vx, vy = Fraction(0, 1), Fraction(0, 1)  # 現在位置
    for i in range(n - 1):
        # (vx, vy) から出て (px, py) を通る直線を考える
        px, py = ub_xys[i]

        # if px <= vx or (1 - vx) * (py - vy) - (1 - vy) * (px - vx) >= 0:
        if px <= vx:
            return True

        # y = a1 x + b1
        a1 = (py - vy) / (px - vx)
        b1 = vy - a1 * vx

        # y = a2 x + b2
        ax, ay = lb_xys[i]
        bx, by = lb_xys[i + 1]

        a2 = Fraction(by - ay, bx - ax)
        b2 = ay - a2 * ax

        vx = (b2 - b1) / (a1 - a2)
        vy = a1 * vx + b1

    return False

print('Yes' if solve() else 'No')
0