結果
| 問題 |
No.2337 Equidistant
|
| コンテスト | |
| ユーザー |
yupooh
|
| 提出日時 | 2023-06-02 22:59:53 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 1,519 ms / 4,000 ms |
| コード長 | 2,539 bytes |
| コンパイル時間 | 348 ms |
| コンパイル使用メモリ | 82,512 KB |
| 実行使用メモリ | 169,076 KB |
| 最終ジャッジ日時 | 2024-12-28 21:29:22 |
| 合計ジャッジ時間 | 22,895 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 28 |
ソースコード
import sys
input = sys.stdin.readline
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=[1]*n
pa=[-1]*n
from collections import deque
que=deque()
que.append((0,1))
seen[0]=1
while que:
v,task=que.pop()
if task==0:
for w in adj[v]:
if w!=pa[v]:
size[v]+=size[w]
else:
que.append((v,0))
seen[v]=1
for w in adj[v]:
if seen[w]==0:
que.append((w,1))
pa[w]=v
depth[w]=depth[v]+1
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 depth[s]<depth[t]:
s,t=t,s
if d%2==1:
print(0)
continue
if depth[s]==depth[t]:
u1=up(d//2-1,s)
u2=up(d//2-1,t)
print(n-size[u1]-size[u2])
continue
u1=up(d//2,s)
u2=up(d//2-1,s)
print(size[u1]-size[u2])
yupooh