結果
| 問題 |
No.898 tri-βutree
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2019-10-26 22:15:53 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
AC
|
| 実行時間 | 3,237 ms / 4,000 ms |
| コード長 | 2,242 bytes |
| コンパイル時間 | 210 ms |
| コンパイル使用メモリ | 12,672 KB |
| 実行使用メモリ | 70,272 KB |
| 最終ジャッジ日時 | 2024-09-14 19:11:01 |
| 合計ジャッジ時間 | 55,689 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 21 |
ソースコード
from collections import deque
class LCA:
def __init__(self, size):
self.size = size
self.level = (size - 1).bit_length()
self.edges = [[] for _ in range(size)]
self.depth = [0] * size
self.length = [float('inf')] * size
self.prev = [None] * size
self.kprev = None
def addEdge(self, fr, to, cost=1):
self.edges[fr].append((to, cost))
self.edges[to].append((fr, cost))
def dfs(self):
st = deque([(0, 0)])
self.depth[0] = 0
self.length[0] = 0
while st:
now, prev = st.pop()
for to, cost in self.edges[now]:
if to == prev:
continue
st.append((to, now))
self.depth[to] = self.depth[now] + 1
self.length[to] = self.length[now] + cost
self.prev[to] = now
def construct(self):
self.dfs()
kprev = [self.prev]
S = self.prev
for _ in range(self.level):
T = [0] * self.size
for i in range(self.size):
if S[i] is None:
continue
T[i] = S[S[i]]
kprev.append(T)
S = T
self.kprev = kprev
def lca(self, u, v):
dist = self.depth[v] - self.depth[u]
if dist < 0:
u, v = v, u
dist = -dist
for k in range(self.level + 1):
if dist & 1:
v = self.kprev[k][v]
dist //= 2
if u == v:
return u
for k in range(self.level)[:: -1]:
prevU = self.kprev[k][u]
prevV = self.kprev[k][v]
if prevU != prevV:
u = prevU
v = prevV
return self.kprev[0][u]
def dist(self, u, v):
return self.length[u] + self.length[v] - self.length[self.lca(u, v)] * 2
N = int(input())
tree = LCA(N)
for _ in range(N - 1):
fr, to, cost = map(int, input().split())
tree.addEdge(fr, to, cost)
tree.construct()
Q = int(input())
ans = []
for _ in range(Q):
x, y, z = map(int, input().split())
ans.append((tree.dist(x, y) + tree.dist(y, z) + tree.dist(z, x)) // 2)
print(*ans, sep='\n')