結果

問題 No.199 星を描こう
ユーザー はむ吉🐹はむ吉🐹
提出日時 2016-09-04 12:11:38
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
WA  
実行時間 -
コード長 1,796 bytes
コンパイル時間 137 ms
コンパイル使用メモリ 13,056 KB
実行使用メモリ 11,008 KB
最終ジャッジ日時 2024-04-27 16:27:45
合計ジャッジ時間 1,932 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 27 ms
10,880 KB
testcase_01 AC 26 ms
10,880 KB
testcase_02 AC 29 ms
11,008 KB
testcase_03 AC 26 ms
10,880 KB
testcase_04 AC 28 ms
10,880 KB
testcase_05 AC 28 ms
10,880 KB
testcase_06 AC 26 ms
10,880 KB
testcase_07 AC 25 ms
10,880 KB
testcase_08 AC 26 ms
10,880 KB
testcase_09 AC 26 ms
10,880 KB
testcase_10 AC 26 ms
10,880 KB
testcase_11 AC 26 ms
10,880 KB
testcase_12 AC 25 ms
10,880 KB
testcase_13 AC 26 ms
10,880 KB
testcase_14 AC 26 ms
10,880 KB
testcase_15 WA -
testcase_16 AC 27 ms
10,880 KB
testcase_17 WA -
testcase_18 AC 26 ms
10,880 KB
testcase_19 WA -
testcase_20 AC 27 ms
10,880 KB
testcase_21 AC 27 ms
10,880 KB
testcase_22 AC 26 ms
10,880 KB
testcase_23 AC 26 ms
10,880 KB
testcase_24 AC 26 ms
10,880 KB
testcase_25 AC 27 ms
10,880 KB
testcase_26 AC 28 ms
10,880 KB
testcase_27 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#!/usr/bin/env python3


import enum


NUM_POINTS = 5
EPS = 1e-10


class PointsRelation(enum.Enum):
    counter_clockwise = 1
    clockwise = 2
    online_back = 3
    on_segment = 4
    online_front = 5


def inner_product(v1, v2):
    return v1.real * v2.real + v1.imag * v2.imag


def outer_product(v1, v2):
    return v1.real * v2.imag - v1.imag * v2.real


def project(a, b):
    return a * inner_product(a, b) / (abs(a) ** 2)


def points_relation(p0, p1, p2):
    v1 = p1 - p0
    v2 = p2 - p0
    op = outer_product(v1, v2)
    if op > EPS:
        return PointsRelation.counter_clockwise
    elif op < -EPS:
        return PointsRelation.clockwise
    elif inner_product(v1, v2) < -EPS:
        return PointsRelation.online_back
    elif abs(v1) < abs(v2):
        return PointsRelation.online_front
    else:
        return PointsRelation.on_segment

        
def in_place_convex_hull_andrew(ps):
    def judge(p0, p1, p2):
        b = points_relation(p0, p1, p2) != PointsRelation.counter_clockwise
        b &= points_relation(p0, p1, p2) != PointsRelation.online_front
        return b
    if len(ps) < 3:
        return ps
    ps.sort(key=lambda p: p.imag)
    n = len(ps)
    k = 0
    ch = [None for _ in range(2 * n)]
    for i in range(n):
        while k >= 2 and judge(ch[k - 2], ch[k - 1], ps[i]):
            k -= 1
        ch[k] = ps[i]
        k += 1
    t = k + 1
    for i in range(n - 1)[::-1]:
        while k >= t and judge(ch[k - 2], ch[k - 1], ps[i]):
            k -= 1
        ch[k] = ps[i]
        k += 1
    ch = ch[:k - 1]
    return ch


def main():
    ps = [complex(*map(float, input().split())) for _ in range(NUM_POINTS)]
    ch = in_place_convex_hull_andrew(ps)
    print("YES" if len(ch) == NUM_POINTS else "NO")


if __name__ == '__main__':
    main()
0