結果

問題 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
コンパイル時間 220 ms
コンパイル使用メモリ 12,032 KB
実行使用メモリ 11,008 KB
最終ジャッジ日時 2023-12-20 12:15:18
合計ジャッジ時間 8,524 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 32 ms
10,880 KB
testcase_01 AC 32 ms
10,880 KB
testcase_02 AC 31 ms
10,880 KB
testcase_03 AC 32 ms
10,880 KB
testcase_04 AC 32 ms
10,880 KB
testcase_05 WA -
testcase_06 AC 32 ms
10,880 KB
testcase_07 AC 31 ms
10,880 KB
testcase_08 RE -
testcase_09 WA -
testcase_10 WA -
testcase_11 AC 32 ms
10,880 KB
testcase_12 AC 32 ms
10,880 KB
testcase_13 AC 32 ms
10,880 KB
testcase_14 AC 32 ms
10,880 KB
testcase_15 AC 32 ms
10,880 KB
testcase_16 AC 31 ms
10,880 KB
testcase_17 AC 31 ms
10,880 KB
testcase_18 AC 32 ms
10,880 KB
testcase_19 AC 32 ms
10,880 KB
testcase_20 AC 32 ms
10,880 KB
testcase_21 AC 31 ms
10,880 KB
testcase_22 WA -
testcase_23 AC 33 ms
10,880 KB
testcase_24 AC 31 ms
10,880 KB
testcase_25 AC 31 ms
10,880 KB
testcase_26 WA -
testcase_27 AC 31 ms
10,880 KB
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 AC 32 ms
10,880 KB
testcase_33 AC 33 ms
10,880 KB
testcase_34 AC 59 ms
11,008 KB
testcase_35 WA -
testcase_36 AC 62 ms
11,008 KB
testcase_37 WA -
testcase_38 AC 60 ms
11,008 KB
testcase_39 WA -
testcase_40 AC 60 ms
11,008 KB
testcase_41 WA -
testcase_42 AC 60 ms
11,008 KB
testcase_43 WA -
testcase_44 AC 60 ms
11,008 KB
testcase_45 AC 54 ms
11,008 KB
testcase_46 AC 60 ms
11,008 KB
testcase_47 WA -
testcase_48 AC 60 ms
11,008 KB
testcase_49 WA -
testcase_50 AC 60 ms
11,008 KB
testcase_51 AC 54 ms
11,008 KB
testcase_52 AC 62 ms
11,008 KB
testcase_53 WA -
testcase_54 AC 34 ms
10,880 KB
testcase_55 AC 32 ms
10,880 KB
testcase_56 WA -
testcase_57 AC 31 ms
10,880 KB
testcase_58 AC 30 ms
10,880 KB
testcase_59 AC 31 ms
10,880 KB
testcase_60 AC 31 ms
10,880 KB
testcase_61 AC 31 ms
10,880 KB
testcase_62 AC 31 ms
10,880 KB
testcase_63 AC 31 ms
10,880 KB
testcase_64 AC 31 ms
10,880 KB
testcase_65 AC 32 ms
10,880 KB
testcase_66 AC 39 ms
10,880 KB
testcase_67 AC 32 ms
10,880 KB
testcase_68 AC 39 ms
10,880 KB
testcase_69 AC 32 ms
10,880 KB
testcase_70 AC 36 ms
10,880 KB
testcase_71 AC 32 ms
10,880 KB
testcase_72 AC 31 ms
10,880 KB
testcase_73 AC 31 ms
10,880 KB
testcase_74 AC 31 ms
10,880 KB
testcase_75 AC 31 ms
10,880 KB
testcase_76 AC 31 ms
10,880 KB
testcase_77 AC 31 ms
10,880 KB
testcase_78 AC 31 ms
10,880 KB
testcase_79 AC 31 ms
10,880 KB
testcase_80 AC 42 ms
10,880 KB
testcase_81 AC 32 ms
10,880 KB
testcase_82 AC 35 ms
10,880 KB
testcase_83 AC 32 ms
10,880 KB
testcase_84 AC 33 ms
10,880 KB
testcase_85 AC 31 ms
10,880 KB
testcase_86 AC 32 ms
10,880 KB
testcase_87 AC 32 ms
10,880 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