結果
問題 |
No.2180 Comprehensive Line Segments
|
ユーザー |
![]() |
提出日時 | 2022-12-17 13:39:51 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 2,247 bytes |
コンパイル時間 | 460 ms |
コンパイル使用メモリ | 82,104 KB |
実行使用メモリ | 214,828 KB |
最終ジャッジ日時 | 2024-11-17 01:37:14 |
合計ジャッジ時間 | 77,776 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 TLE * 1 |
other | AC * 15 WA * 2 TLE * 8 |
ソースコード
from __future__ import annotations import sys from itertools import permutations from fractions import Fraction as frac input = sys.stdin.readline class Vector2: def __init__(self, x, y): self.x = x self.y = y def __str__(self): return str((self.x, self.y)) __repr__ = __str__ def __eq__(self, other): return (self.x == other.x and self.y == other.y) def __hash__(self): return hash((self.x, self.y)) def __add__(self, other): return Vector2(self.x + other.x, self.y + other.y) def __sub__(self, other): return Vector2(self.x - other.x, self.y - other.y) def __mul__(self, other): return self.x * other.y - self.y * other.x def normalize(self) -> Vector2: assert(self.x != 0 or self.y != 0) norm = self.x ** 2 + self.y ** 2 self.x *= abs(self.x) / norm self.y *= abs(self.y) / norm return self def sgn(x: frac) -> int: if(x > 0): return 1 if(x < 0): return -1 return 0 """ Main Code """ N = int(input()) P = [Vector2(*map(frac, input().split())) for _ in [0] * N] if(N == 1): print(1) exit(0) vec_lis = [[None]*N for _ in [0]*N] for i in range(N - 1): for j in range(i + 1, N): vec_lis[i][j] = (P[j] - P[i]).normalize() vec_lis[j][i] = (P[i] - P[j]).normalize() ans = N dp = [[[N]*N for _ in [0]*N] for _ in [0]*(1 << N)] num = [*range(N)] for i, j in permutations(num, 2): dp[(1 << i) | (1 << j)][i][j] = 1 goal = (1 << N) - 1 for b_now in range(3, goal): for v_prev, v_now in permutations(num, 2): if(dp[b_now][v_prev][v_now] == N or not((b_now >> v_prev) & 1) or not((b_now >> v_now) & 1)): continue c_next = c_now = dp[b_now][v_prev][v_now] v1 = vec_lis[v_prev][v_now] for v_next1, v_next2 in permutations(num, 2): if((((b_now >> v_next1) & 1) and v_now != v_next1) or ((b_now >> v_next2) & 1)): continue b_next = b_now | (1 << v_next1) | (1 << v_next2) v2, v3 = vec_lis[v_now][v_next1], vec_lis[v_next1][v_next2] if(v_now == v_next1): c_next = c_now + (v1 != v3) else: d = (v1 == v2) + (v2 == v3) if(d == 0): d = (sgn(v1 * v2) == sgn(v2 * v3) == sgn(v1 * v3)) c_next = c_now + 2 - d if(dp[b_next][v_next1][v_next2] <= c_next): continue dp[b_next][v_next1][v_next2] = c_next ans = min(min(dp[goal])) print(ans)