結果

問題 No.2602 Real Collider
ユーザー cleantted
提出日時 2023-11-28 20:21:33
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 337 ms / 2,000 ms
コード長 2,445 bytes
コンパイル時間 285 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 121,728 KB
最終ジャッジ日時 2024-09-26 13:04:44
合計ジャッジ時間 24,964 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 78
権限があれば一括ダウンロードができます

ソースコード

diff # 

import copy
import heapq
import itertools
import math
import operator
import sys
from bisect import bisect, bisect_left, bisect_right, insort
from collections import Counter, deque
from fractions import Fraction
from functools import cmp_to_key, lru_cache, partial
from inspect import currentframe
from math import ceil, gcd, log10, pi, sqrt

# import pypyjit
# pypyjit.set_param('max_unroll_recursion=-1')
input = sys.stdin.readline
sys.setrecursionlimit(10000000)
# mod = 10 ** 9 + 7
mod = 998244353
# mod = 1 << 128
# mod = 10 ** 30 + 1
INF = 1 << 61
DIFF = 10 ** -9
DX = [1, 0, -1, 0, 1, 1, -1, -1]
DY = [0, 1, 0, -1, 1, -1, 1, -1]

def read_values(): return tuple(map(int, input().split()))
def read_index(): return tuple(map(lambda x: int(x) - 1, input().split()))
def read_list(): return list(read_values())
def read_lists(N): return [read_list() for _ in range(N)]
def dprint(*values): print(*values, file=sys.stderr)
def dprint2(*values):
    names = {id(v): k for k, v in currentframe().f_back.f_locals.items()}
    dprint(", ".join(f"{names.get(id(value), '???')}={repr(value)}" for value in values))


def main():
    Q = int(input())
    XA, YA, XB, YB, XC, YC = read_list()
    L = read_lists(Q)

    def dot(x1, y1, x2, y2):
        return x1 * x2 + y1 * y2

    def rot(x1, y1, x2, y2):
        return x1 * y2 - x2 * y1

    def dist2(x1, y1, x2, y2):
        return (x1 - x2) ** 2 + (y1 - y2) ** 2
    
    def g(x1, y1, x2, y2, x3, y3):
        return dist2(x1, y1, x2, y2), dot(x3 - x1, y3 - y1, x3 - x2, y3 - y2), (x1, y1), (x2, y2), (x3, y3)

    _, c, A, B, C = max(
        g(XA, YA, XB, YB, XC, YC),
        g(XB, YB, XC, YC, XA, YA),
        g(XC, YC, XA, YA, XB, YB),
    )

    k1 = rot(*A, *B) + rot(*B, *C) + rot(*C, *A)

    if c <= 0:
        t1 = 0
        t2 = 1
    else:
        t1 = c ** 2
        t2 = dist2(*A, *C) * dist2(*B, *C)

    res = []
    for X, Y in L:
        s1 = dot(X - A[0], Y - A[1], X - B[0], Y - B[1])
        s2 = dist2(X, Y, *A) * dist2(X, Y, *B)
        k2 = rot(*A, *B) + rot(*B, X, Y) + rot(X, Y, *A)
        if k1 * k2 >= 0:
            if s1 < 0 or t1 * s2 >= (s1 ** 2) * t2:
                res.append("Yes")
            else:
                res.append("No")
        else:
            if s1 <= 0 and t1 * s2 <= (s1 ** 2) * t2:
                res.append("Yes")
            else:
                res.append("No")

    print(*res, sep="\n")


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