結果
問題 | No.901 K-ary εxtrεεmε |
ユーザー | terasa |
提出日時 | 2023-01-15 17:13:10 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 433 ms / 3,000 ms |
コード長 | 4,628 bytes |
コンパイル時間 | 307 ms |
コンパイル使用メモリ | 82,372 KB |
実行使用メモリ | 119,632 KB |
最終ジャッジ日時 | 2024-06-09 04:59:26 |
合計ジャッジ時間 | 14,027 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 277 ms
119,632 KB |
testcase_01 | AC | 57 ms
69,596 KB |
testcase_02 | AC | 91 ms
78,612 KB |
testcase_03 | AC | 91 ms
79,032 KB |
testcase_04 | AC | 89 ms
79,060 KB |
testcase_05 | AC | 86 ms
78,924 KB |
testcase_06 | AC | 93 ms
79,108 KB |
testcase_07 | AC | 433 ms
115,880 KB |
testcase_08 | AC | 428 ms
115,908 KB |
testcase_09 | AC | 411 ms
116,256 KB |
testcase_10 | AC | 419 ms
116,184 KB |
testcase_11 | AC | 411 ms
116,040 KB |
testcase_12 | AC | 421 ms
115,840 KB |
testcase_13 | AC | 416 ms
116,056 KB |
testcase_14 | AC | 423 ms
115,856 KB |
testcase_15 | AC | 433 ms
116,228 KB |
testcase_16 | AC | 405 ms
115,880 KB |
testcase_17 | AC | 430 ms
115,752 KB |
testcase_18 | AC | 426 ms
115,880 KB |
testcase_19 | AC | 421 ms
116,352 KB |
testcase_20 | AC | 421 ms
115,864 KB |
testcase_21 | AC | 420 ms
116,052 KB |
testcase_22 | AC | 424 ms
115,596 KB |
testcase_23 | AC | 417 ms
116,384 KB |
testcase_24 | AC | 422 ms
116,224 KB |
testcase_25 | AC | 406 ms
116,152 KB |
testcase_26 | AC | 413 ms
116,232 KB |
testcase_27 | AC | 399 ms
116,032 KB |
testcase_28 | AC | 403 ms
115,760 KB |
testcase_29 | AC | 418 ms
116,172 KB |
ソースコード
from typing import List, Tuple, Callable, TypeVar, Optional import sys import itertools import heapq import bisect import math from collections import deque, defaultdict from functools import lru_cache, cmp_to_key input = sys.stdin.readline if __file__ != 'prog.py': sys.setrecursionlimit(10 ** 6) def readints(): return map(int, input().split()) def readlist(): return list(readints()) def readstr(): return input()[:-1] def readlist1(): return list(map(lambda x: int(x) - 1, input().split())) class CumulativeSum: def __init__(self, A): self.S = [0] acc = 0 for a in A: acc += a self.S.append(acc) def get(self, l, r): """return sum(A[l:r]), i.e. sum of A[x] (l <= x < r) """ return self.S[r] - self.S[l] class HLD: # reference: https://codeforces.com/blog/entry/53170 def __init__(self, N, E, root: int = 0): self.N = N self.E = E self.root = root self.D = [0] * self.N self.par = [-1] * self.N self.sz = [0] * self.N self.top = [0] * self.N self.ord = [None] * self.N self._dfs_sz() self._dfs_hld() def path_query_range(self, u: int, v: int, is_edge_query: bool = False) -> List[Tuple[int, int]]: """return list of [l, r) ranges that cover u-v path""" ret = [] while True: if self.ord[u] > self.ord[v]: u, v = v, u if self.top[u] == self.top[v]: ret.append((self.ord[u] + is_edge_query, self.ord[v] + 1)) return ret ret.append((self.ord[self.top[v]], self.ord[v] + 1)) v = self.par[self.top[v]] def subtree_query_range(self, v: int, is_edge_query: bool = False) -> Tuple[int, int]: """return [l, r) range that cover vertices of subtree v""" return (self.ord[v] + is_edge_query, self.ord[v] + self.sz[v]) def get_index(self, v: int) -> int: """return euler tour order of given vertex""" return self.ord[v] def lca(self, u, v): while True: if self.ord[u] > self.ord[v]: u, v = v, u if self.top[u] == self.top[v]: return u v = self.par[self.top[v]] def dist(self, u, v): return self.D[u] + self.D[v] - 2 * self.D[self.lca(u, v)] def _dfs_sz(self): stack = [(self.root, -1)] while stack: v, p = stack.pop() if v < 0: v = ~v self.sz[v] = 1 for i, dst in enumerate(self.E[v]): if dst == p: continue self.sz[v] += self.sz[dst] # v -> E[v][0] will be heavy path if self.sz[self.E[v][0]] < self.sz[dst]: self.E[v][0], self.E[v][i] = self.E[v][i], self.E[v][0] else: if ~p: self.D[v] = self.D[p] + 1 self.par[v] = p # avoid first element of E[v] is parent of vertex v if v has some children if len(self.E[v]) >= 2 and self.E[v][0] == p: self.E[v][0], self.E[v][1] = self.E[v][1], self.E[v][0] stack.append((~v, p)) for dst in self.E[v]: if dst == p: continue stack.append((dst, v)) def _dfs_hld(self): stack = [(self.root, -1)] cnt = 0 while stack: v, p = stack.pop() self.ord[v] = cnt cnt += 1 heavy_path_idx = len(self.E[v]) - 1 for i, dst in enumerate(self.E[v][::-1]): if dst == p: continue # top[dst] is top[v] if v -> dst is heavy path otherwise dst itself self.top[dst] = self.top[v] if i == heavy_path_idx else dst stack.append((dst, v)) N = int(input()) E = [[] for _ in range(N)] edges = [] for _ in range(N - 1): u, v, w = readints() E[u].append(v) E[v].append(u) edges.append((u, v, w)) solver = HLD(N, E) W = [0] * N for u, v, w in edges: if solver.ord[u] > solver.ord[v]: u, v = v, u W[solver.get_index(v)] = w S = CumulativeSum(W) Q = int(input()) for _ in range(Q): k, *X = readlist() X = sorted(X, key=lambda x: solver.ord[x]) ans = 0 for i in range(k): s = X[i] t = X[(i + 1) % k] for l, r in solver.path_query_range(s, t, is_edge_query=True): ans += S.get(l, r) print(ans // 2)