結果

問題 No.1094 木登り / Climbing tree
ユーザー lloyzlloyz
提出日時 2022-10-14 00:13:44
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,770 ms / 2,000 ms
コード長 3,539 bytes
コンパイル時間 178 ms
コンパイル使用メモリ 82,156 KB
実行使用メモリ 161,280 KB
最終ジャッジ日時 2024-06-26 12:15:01
合計ジャッジ時間 40,092 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 43 ms
54,656 KB
testcase_01 AC 1,672 ms
142,900 KB
testcase_02 AC 414 ms
161,280 KB
testcase_03 AC 452 ms
82,464 KB
testcase_04 AC 596 ms
107,656 KB
testcase_05 AC 977 ms
134,132 KB
testcase_06 AC 844 ms
98,280 KB
testcase_07 AC 1,670 ms
143,092 KB
testcase_08 AC 1,674 ms
143,768 KB
testcase_09 AC 1,770 ms
142,824 KB
testcase_10 AC 1,699 ms
143,012 KB
testcase_11 AC 1,674 ms
143,356 KB
testcase_12 AC 1,683 ms
143,108 KB
testcase_13 AC 1,709 ms
143,340 KB
testcase_14 AC 1,657 ms
143,076 KB
testcase_15 AC 635 ms
93,952 KB
testcase_16 AC 839 ms
141,568 KB
testcase_17 AC 728 ms
114,100 KB
testcase_18 AC 688 ms
104,276 KB
testcase_19 AC 797 ms
130,176 KB
testcase_20 AC 1,674 ms
143,448 KB
testcase_21 AC 751 ms
117,760 KB
testcase_22 AC 1,643 ms
143,600 KB
testcase_23 AC 1,657 ms
143,012 KB
testcase_24 AC 1,681 ms
142,600 KB
testcase_25 AC 1,685 ms
142,848 KB
testcase_26 AC 1,682 ms
142,588 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

class HeavyLightDecomposition:
    # HL 分解 (HLD)
    def __init__(self, E, root=1):
        # E は双方向に辺を張った木で、破壊的に有向グラフに変更される
        # O(N)
        self.E = E
        self.root = root
        self.N = N = len(E)  # 頂点数
        self.Parent = [-1] * N  # 頂点番号 v -> 親ノード
        self.Size = [-1] * N  # 頂点番号 v -> 部分木のサイズ
        self.dfs1()

        self.Mapping = [-1] * N  # 頂点番号 v -> 内部インデックス
        self.Head = list(range(N))  # 頂点番号 v -> v を含む heavy path の左端の頂点番号
        self.Depth = [0] * N  # 頂点番号 v -> 深さ（root から v までの距離）
        self.dfs2()

    def dfs1(self):
        E = self.E
        Parent, Size = self.Parent, self.Size
        Path = [self.root]
        Idx_edge = [0]
        while Path:
            v = Path[-1]
            idx_edge = Idx_edge[-1]
            Ev = E[v]
            if idx_edge != len(Ev):
                # 行きがけ・通りがけ 辺の数だけ実行される
                u = Ev[idx_edge]
                Idx_edge[-1] += 1
                E[u].remove(v)  # 有向グラフならここをコメントアウトする
                Parent[u] = v
                Path.append(u)
                Idx_edge.append(0)
            else:
                # 帰りがけ 頂点の数だけ実行される
                if len(Ev) >= 2:
                    ma = -1
                    argmax = None
                    for i, u in enumerate(Ev):
                        if Size[u] > ma:
                            ma = Size[u]
                            argmax = i
                    Ev[0], Ev[argmax] = Ev[argmax], Ev[0]
                Size[v] = sum(Size[u] for u in Ev) + 1
                Path.pop()
                Idx_edge.pop()
    
    def dfs2(self):
        E = self.E
        Mapping = self.Mapping
        Head = self.Head
        Depth = self.Depth
        k = 0
        St = [self.root]
        while St:
            v = St.pop()
            Mapping[v] = k
            k += 1
            Ev = E[v]
            if Ev:
                Head[Ev[0]] = Head[v]
                St += Ev[::-1]
                for u in Ev:
                    Depth[u] = Depth[v] + 1  # distance を使わないのならここをコメントアウトする
    
    def lca(self, v, u):
        # O(logN)
        Parent = self.Parent
        Mapping = self.Mapping
        Head = self.Head
        while True:
            if Mapping[v] > Mapping[u]:
                v, u = u, v  # v の方を根に近くする
            if Head[v] == Head[u]:
                return v
            u = Parent[Head[u]]
    
    def distance(self, v, u):
        # O(logN)
        Depth = self.Depth
        return Depth[v] + Depth[u] - 2 * Depth[self.lca(v, u)]

from collections import deque

n = int(input())
E = [[] for _ in range(n + 1)]
EC = [[] for _ in range(n + 1)]
for _ in range(n - 1):
    a, b, c = map(int, input().split())
    E[a].append(b)
    EC[a].append((b, c))
    E[b].append(a)
    EC[b].append((a, c))

INF = 10**18
C = [INF for _ in range(n + 1)]
C[1] = 0
Que = deque([1])
while Que:
    cp = Que.popleft()
    cc = C[cp]
    for np, dc in EC[cp]:
        if cc + dc >= C[np]:
            continue
        C[np] = cc + dc
        Que.append(np)
HL = HeavyLightDecomposition(E)

q = int(input())
for _ in range(q):
    s, t = map(int, input().split())
    st_lca = HL.lca(s, t)
    print(C[s] + C[t] - 2 * C[st_lca])
0