結果

問題 No.2179 Planet Traveler
ユーザー taiga0629kyopro
提出日時 2022-09-26 19:16:22
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 416 ms / 3,000 ms
コード長 1,598 bytes
コンパイル時間 182 ms
コンパイル使用メモリ 82,472 KB
実行使用メモリ 118,292 KB
最終ジャッジ日時 2024-11-30 16:29:43
合計ジャッジ時間 6,373 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 




########################################
from heapq import heappush, heappop
def dijkstra( G, s, INF=10 ** 18):
    """
    https://tjkendev.github.io/procon-library/python/graph/dijkstra.html
    O((|E|+|V|)log|V|)
    V: 頂点数
    G[v] = [(nod, cost)]:
        頂点vから遷移可能な頂点(nod)とそのコスト(cost)
    s: 始点の頂点"""

    N=len(G)
    N+=1
    dist = [INF] * N
    hp = [(0, s)]  # (c, v)
    dist[s] = 0
    while hp:
        c, v = heappop(hp)  #vまで行くコストがc
        if dist[v] < c:
            continue
        for u, cost in G[v]:
            if max(dist[v] , cost) < dist[u]:
                dist[u] = max(dist[v] , cost)
                heappush(hp, (dist[u], u))
    return dist
##################################################



n=int(input())
x=[0]
y=[0]
t=[0]
for i in range(n):
    a,b,c=map(int,input().split())
    x.append(a)
    y.append(b)
    t.append(c)
def cost(i,j):
    if t[i]==t[j]:
        return (x[i]-x[j])**2+(y[i]-y[j])**2
    else:
        ri2 = x[i] ** 2 + y[i] ** 2
        rj2 = x[j] ** 2 + y[j] ** 2
        def f(s):
            if s<0:return 0
            if ri2+rj2-s<0:return 1
            return (ri2+rj2-s)**2<=4*ri2*rj2
        ng=int(abs(ri2**0.5-rj2**0.5)**2)-2
        ok=ng+4
        while ok-ng>1:
            mid=(ok+ng)//2
            if f(mid):ok=mid
            else:ng=mid
        return ok


root=[[] for i in range(n+5)]
for i in range(1,n+1):
    for j in range(i+1,n+1):
        c=cost(i,j)
        root[i].append((j,c))
        root[j].append((i,c))

print(dijkstra(root,1)[n])
0