結果
問題 | No.2337 Equidistant |
ユーザー | yupooh |
提出日時 | 2023-06-02 22:20:41 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 2,574 bytes |
コンパイル時間 | 271 ms |
コンパイル使用メモリ | 13,184 KB |
実行使用メモリ | 133,248 KB |
最終ジャッジ日時 | 2024-06-08 23:42:56 |
合計ジャッジ時間 | 57,408 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 32 ms
11,136 KB |
testcase_01 | RE | - |
testcase_02 | RE | - |
testcase_03 | RE | - |
testcase_04 | RE | - |
testcase_05 | RE | - |
testcase_06 | RE | - |
testcase_07 | RE | - |
testcase_08 | RE | - |
testcase_09 | RE | - |
testcase_10 | RE | - |
testcase_11 | RE | - |
testcase_12 | RE | - |
testcase_13 | RE | - |
testcase_14 | RE | - |
testcase_15 | RE | - |
testcase_16 | RE | - |
testcase_17 | RE | - |
testcase_18 | RE | - |
testcase_19 | RE | - |
testcase_20 | RE | - |
testcase_21 | RE | - |
testcase_22 | WA | - |
testcase_23 | RE | - |
testcase_24 | TLE | - |
testcase_25 | RE | - |
testcase_26 | TLE | - |
testcase_27 | RE | - |
testcase_28 | RE | - |
ソースコード
import sys input = sys.stdin.readline sys.setrecursionlimit(500000) n,q=map(int,input().split()) adj=[[] for _ in range(n)] S=[] for _ in range(n-1): a,b=map(int,input().split()) adj[a-1].append(b-1) adj[b-1].append(a-1) S.append((a,b)) depth = [0]*n seen = [0]*n size=[0]*n def dfs(v, e): depth[v] = e seen[v]=1 tot=1 for w in adj[v]: if seen[w]==0: tot+=dfs(w, e+1) size[v]=tot return tot dfs(0, 0) N=n INF = 10 ** 16 depth = [INF] * N dfs_seq = [INF] * N # 0を根にして深さdepthとDFS順序dfs_seqを求める todo = [(0, 0)] count = 0 while todo: d, pos = todo.pop() if dfs_seq[pos] < INF: continue dfs_seq[pos] = count count += 1 depth[pos] = d for next_ in adj[pos]: todo.append((d + 1, next_)) # 2 ** k 個上の頂点を前計算する(k < log2 N) import math K = int(math.log2(N)) + 1 up2 = [[INF] * N for _ in range(K)] for a, b in S: if depth[a-1] < depth[b-1]: up2[0][b-1] = a-1 else: up2[0][a-1] = b-1 for k in range(1, K): for p in range(N): if up2[k-1][p] < INF: up2[k][p] = up2[k-1][up2[k-1][p]] else: up2[k][p] = INF # p個上の頂点をダブリングで求める(p <= depth[n]) def up(n, p): k = 0 while n > 0: if n % 2: p = up2[k][p] if p == INF: break n //= 2 k += 1 return p # 最小共通祖先LCAを求める def lca(p, q): # 深い方のdepthを、浅い方に揃える if depth[p] > depth[q]: p, q = q, p d = depth[p] q = up(depth[q] - d, q) # depthをそろえた頂点が同じであれば、それが答え if p == q: return p # さもなければ、前計算した2**k個上の頂点を順に見比べて、答えを探索する k = K - 1 while k >= 0: p1, q1 = up2[k][p], up2[k][q] if p1 != q1: p, q = p1, q1 else: ans = p1 k -= 1 return ans for _ in range(q): s,t=map(int,input().split()) s-=1 t-=1 x=lca(s,t) d=depth[s]+depth[t]-2*depth[x] if d%2==1: print(0) continue if x==s: u1=up(t,d//2) u2=up(t,d//2-1) print(size[u1]-size[u2]) elif x==t: u1=up(s,d//2) u2=up(s,d//2-1) print(size[u1]-size[u2]) elif depth[s]-depth[x]==d//2: print(n-size[x]+1) else: if depth[s]>depth[t]: u1=up(s,d//2) u2=up(s,d//2-1) print(size[u1]-size[u2]) else: u1=up(t,d//2) u2=up(t,d//2-1) print(size[u1]-size[u2])