結果

問題 No.2179 Planet Traveler
ユーザー H20H20
提出日時 2023-01-06 22:54:32
言語 PyPy3
(7.3.15)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 2,361 bytes
コンパイル時間 290 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 166,256 KB
最終ジャッジ日時 2024-11-30 19:28:53
合計ジャッジ時間 44,687 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 25 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff #

# UnionFind 参考は以下のサイト
# https://note.nkmk.me/python-union-find/
from collections import defaultdict
from decimal import Decimal

class UnionFind():
    def __init__(self, n):
        self.n = n
        self.parents = [-1] * n

    def find(self, x):
        if self.parents[x] < 0:
            return x
        else:
            self.parents[x] = self.find(self.parents[x])
            return self.parents[x]

    def union(self, x, y):
        x = self.find(x)
        y = self.find(y)

        if x == y:
            return

        if self.parents[x] > self.parents[y]:
            x, y = y, x

        self.parents[x] += self.parents[y]
        self.parents[y] = x

    def size(self, x):
        return -self.parents[self.find(x)]

    def same(self, x, y):
        return self.find(x) == self.find(y)

    def members(self, x):
        root = self.find(x)
        return [i for i in range(self.n) if self.find(i) == root]

    def roots(self):
        return [i for i, x in enumerate(self.parents) if x < 0]

    def group_count(self):
        return len(self.roots())

    def all_group_members(self):
        group_members = defaultdict(list)
        for member in range(self.n):
            group_members[self.find(member)].append(member)
        return group_members

    def __str__(self):
        return '\n'.join(f'{r}: {m}' for r, m in self.all_group_members().items())


def is_ok(s):
    uf = UnionFind(N)
    hs = Decimal(s)**Decimal('0.5')
    for i in range(N):
        for j in range(i+1,N):
            if not uf.same(i,j):
                xi,yi,ti = XYT[i]
                xj,yj,tj = XYT[j]
                if ti==tj:
                    if (xi-xj)**2+(yi-yj)**2<=s:
                        uf.union(i,j)
                else:
                    if XY[i][j]<=hs:
                        uf.union(i,j)
    return uf.same(0,N-1)

def meguru_bisect(ng, ok):
    while (abs(ok - ng) > 1):
        mid = (ok + ng) // 2
        if is_ok(mid):
            ok = mid
        else:
            ng = mid
    return ok

N = int(input())
XYT = [list(map(int, input().split())) for _ in range(N)]
D = []
for i in range(N):
    D.append(Decimal(XYT[i][0]**2+XYT[i][1]**2)**Decimal('0.5'))

XY = [[0]*N for _ in range(N)]
for i in range(N):
    for j in range(i+1,N):
        XY[i][j]=abs(D[i]-D[j])

print(meguru_bisect(-1,3200000001))
0