結果

問題 No.3005 トレミーの問題
ユーザー lif4635lif4635
提出日時 2025-01-17 21:33:19
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 65 ms / 2,000 ms
コード長 3,371 bytes
コンパイル時間 281 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 53,248 KB
最終ジャッジ日時 2025-01-17 21:33:31
合計ジャッジ時間 2,726 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #

# input
import sys
input = sys.stdin.readline
II = lambda : int(input())
MI = lambda : map(int, input().split())
LI = lambda : list(map(int, input().split()))
SI = lambda : input().rstrip()
LLI = lambda n : [list(map(int, input().split())) for _ in range(n)]
LSI = lambda n : [input().rstrip() for _ in range(n)]
MI_1 = lambda : map(lambda x:int(x)-1, input().split())
LI_1 = lambda : list(map(lambda x:int(x)-1, input().split()))

def graph(n:int, m:int, dir:bool=False, index:int=-1) -> list[set[int]]:
    edge = [set() for i in range(n+1+index)]
    for _ in range(m):
        a,b = map(int, input().split())
        a += index
        b += index
        edge[a].add(b)
        if not dir:
            edge[b].add(a)
    return edge

def graph_w(n:int, m:int, dir:bool=False, index:int=-1) -> list[set[tuple]]:
    edge = [set() for i in range(n+1+index)]
    for _ in range(m):
        a,b,c = map(int, input().split())
        a += index
        b += index
        edge[a].add((b,c))
        if not dir:
            edge[b].add((a,c))
    return edge

mod = 998244353
inf = 1001001001001001001
ordalp = lambda s : ord(s)-65 if s.isupper() else ord(s)-97
ordallalp = lambda s : ord(s)-39 if s.isupper() else ord(s)-97
yes = lambda : print("Yes")
no = lambda : print("No")
yn = lambda flag : print("Yes" if flag else "No")
def acc(a:list[int]):
    sa = [0]*(len(a)+1)
    for i in range(len(a)):
        sa[i+1] = a[i] + sa[i]
    return sa

prinf = lambda ans : print(ans if ans < 1000001001001001001 else -1)
alplow = "abcdefghijklmnopqrstuvwxyz"
alpup = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
alpall = "abcdefghijklmnopqrstuvwxyzABCDEFGHIJKLMNOPQRSTUVWXYZ"
URDL = {'U':(-1,0), 'R':(0,1), 'D':(1,0), 'L':(0,-1)}
DIR_4 = [[-1,0],[0,1],[1,0],[0,-1]]
DIR_8 = [[-1,0],[-1,1],[0,1],[1,1],[1,0],[1,-1],[0,-1],[-1,-1]]
DIR_BISHOP = [[-1,1],[1,1],[1,-1],[-1,-1]]
prime60 = [2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59]
import pypyjit
pypyjit.set_param('max_unroll_recursion=-1')
sys.setrecursionlimit(10**7)
sys.set_int_max_str_digits(0)


# import numpy as np
p = []
mat = []
for i in range(4):
    x,y = MI()
    p.append((x,y))
    mat.append([x,y,x**2+y**2,1])


from itertools import combinations

def line(p1,p2,p3):
    x1, y1 = p1
    x2, y2 = p2
    x3, y3 = p3
    return (x2 - x1) * (y3 - y1) == (x3 - x1) * (y2 - y1)

for i,j,k in combinations(p,3):
    if line(i,j,k):
        print("NO")
        exit()

def determinant(a, mod = mod): #行列式
    n = len(a)
    res = 1
    for i in range(n):
        #0でない行をもってくる
        for j in range(i+1,n):
            if a[j][i] == 0: continue
            a[i],a[j] = a[j],a[i]
            res *= -1
            break
        # else:
        #     # 行列式が0になる
        #     return 0
        
        res *= a[i][i]
        res %= mod
        # 行列式自体は0ではないがmodで0になる
        # 逆元がとれないのでreturnしておく
        if a[i][i]%mod == 0:
            return 0
        
        inv = pow(a[i][i],-1,mod)
        for j in range(n):
            a[i][j] *= inv
            a[i][j] %= mod
        for j in range(i+1,n):
            tmp = a[j][i]
            for k in range(n):
                a[j][k] -= a[i][k]*tmp
                a[j][k] %= mod
        # print(res)
    return res%mod



if determinant(mat) == 0:
    print("YES")
else:
    print("NO")
0