結果
問題 | No.3005 トレミーの問題 |
ユーザー | lif4635 |
提出日時 | 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 |
ソースコード
# 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")