結果
| 問題 | No.827 総神童数 |
| ユーザー |
👑  Kazun
|
| 提出日時 | 2021-05-19 13:31:16 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 538 ms / 2,000 ms |
| コード長 | 12,922 bytes |
| コンパイル時間 | 321 ms |
| コンパイル使用メモリ | 82,264 KB |
| 実行使用メモリ | 161,468 KB |
| 最終ジャッジ日時 | 2024-10-09 17:38:49 |
| 合計ジャッジ時間 | 12,366 ms |
|
ジャッジサーバーID （参考情報） |
judge4 / judge3 |
（要ログイン）
テストケース
テストケース表示| 入力 | 結果 | 実行時間 実行使用メモリ |
|---|---|---|
| testcase_00 | AC | 45 ms
56,576 KB |
| testcase_01 | AC | 46 ms
56,576 KB |
| testcase_02 | AC | 46 ms
56,704 KB |
| testcase_03 | AC | 45 ms
56,576 KB |
| testcase_04 | AC | 45 ms
56,448 KB |
| testcase_05 | AC | 45 ms
56,576 KB |
| testcase_06 | AC | 45 ms
56,576 KB |
| testcase_07 | AC | 48 ms
56,320 KB |
| testcase_08 | AC | 46 ms
56,320 KB |
| testcase_09 | AC | 531 ms
161,468 KB |
| testcase_10 | AC | 191 ms
98,580 KB |
| testcase_11 | AC | 84 ms
77,964 KB |
| testcase_12 | AC | 118 ms
84,352 KB |
| testcase_13 | AC | 408 ms
151,868 KB |
| testcase_14 | AC | 92 ms
79,112 KB |
| testcase_15 | AC | 199 ms
102,260 KB |
| testcase_16 | AC | 434 ms
148,904 KB |
| testcase_17 | AC | 520 ms
161,144 KB |
| testcase_18 | AC | 259 ms
113,464 KB |
| testcase_19 | AC | 508 ms
150,228 KB |
| testcase_20 | AC | 384 ms
130,284 KB |
| testcase_21 | AC | 291 ms
115,784 KB |
| testcase_22 | AC | 449 ms
136,360 KB |
| testcase_23 | AC | 79 ms
75,164 KB |
| testcase_24 | AC | 501 ms
141,388 KB |
| testcase_25 | AC | 365 ms
125,432 KB |
| testcase_26 | AC | 506 ms
140,944 KB |
| testcase_27 | AC | 327 ms
120,572 KB |
| testcase_28 | AC | 310 ms
117,368 KB |
| testcase_29 | AC | 286 ms
107,608 KB |
| testcase_30 | AC | 134 ms
81,756 KB |
| testcase_31 | AC | 222 ms
103,024 KB |
| testcase_32 | AC | 234 ms
102,804 KB |
| testcase_33 | AC | 538 ms
148,244 KB |
| testcase_34 | AC | 504 ms
141,864 KB |
| testcase_35 | AC | 247 ms
103,868 KB |
| testcase_36 | AC | 295 ms
115,084 KB |
| testcase_37 | AC | 378 ms
125,640 KB |
| testcase_38 | AC | 488 ms
145,232 KB |
ソースコード
class Tree:
def __init__(self,N,index=0):
"""N頂点(index, index+1, ..., N-1+index)の根付き木を生成する.
"""
self.N=N
self.index=index
self.parent=[-1]*(N+index)
self.__mutable=True
def vertex_exist(self,x):
return self.index<=x<self.index+self.N
def __after_seal_check(self,*vertexes):
if self.__mutable:
return False
for v in vertexes:
if not self.vertex_exist(v):
return False
return True
def is_mutable(self):
return self.__mutable
#設定パート
def root_set(self,root):
"""頂点xを根に設定する.
"""
assert self.vertex_exist(root)
assert self.__mutable
self.root=root
def parent_set(self,x,y):
"""頂点xの親をyに設定する.
"""
assert self.vertex_exist(x)
assert self.vertex_exist(y)
assert self.__mutable
self.parent[x]=y
def child_set(self,x,y):
"""頂点xの子の一つにyを設定する.
"""
assert self.vertex_exist(x)
assert self.vertex_exist(y)
assert self.__mutable
self.parent[y]=x
def seal(self):
"""木の情報を確定させる.
"""
assert self.__mutable
assert hasattr(self,"root")
a=self.index
b=self.index+self.N
C=[[] for _ in range(b)]
p=self.parent
ve=self.vertex_exist
for i in range(a,b):
if i!=self.root:
assert ve(p[i])
C[p[i]].append(i)
self.__mutable=False
self.children=C
#データを求める.
def depth_search(self,Mode=True):
"""木の深さを求める.
"""
assert self.__after_seal_check()
if hasattr(self,"depth"):
return self.depth
from collections import deque
C=self.children
D=[-1]*(self.index+self.N)
E=[[] for _ in range(self.N)]
Q=deque([self.root])
D[self.root]=0
E[0]=[self.root]
while Q:
x=Q.popleft()
d=D[x]
for y in C[x]:
D[y]=d+1
E[d+1].append(y)
Q.append(y)
self.depth=D
self.tower=E
if Mode:
return D
def vertex_depth(self,x):
"""頂点xの深さを求める.
"""
assert self.__after_seal_check(x)
if not hasattr(self,"depth"):
self.depth_search(Mode=False)
return self.depth[x]
def __upper_list(self):
assert self.__after_seal_check()
if hasattr(self,"upper_list"):
return
if not hasattr(self,"depth"):
self.depth_search(False)
b=max(self.depth).bit_length()
X=[[-1]*(self.index+self.N) for _ in range(b)]
Y=X[0]
p=self.parent
rg=range(self.index,self.index+self.N)
for x in rg:
if x!=self.root:
Y[x]=p[x]
else:
Y[x]=self.root
for k in range(1,b):
Y=X[k-1]
Z=X[k]
for x in rg:
Z[x]=Y[Y[x]]
self.upper_list=X
def upper(self,x,k,over=True):
"""頂点xから見てk個親の頂点を求める.
over:(頂点xの深さ)<kのときにTrueならば根を返し, Falseならばエラーを吐く.
"""
assert self.__after_seal_check(x)
assert 0<=k
if not hasattr(self,"upper_list"):
self.__upper_list()
if self.vertex_depth(x)<k:
if over:
return self.root
else:
raise ValueError
i=0
while k:
if k&1:
x=self.upper_list[i][x]
k>>=1
i+=1
return x
def lowest_common_ancestor(self,x,y):
"""頂点x,yの最小共通先祖(x,yに共通する先祖で最も深いもの)を求める.
"""
assert self.__after_seal_check(x,y)
dd=self.vertex_depth(y)-self.vertex_depth(x)
if dd<0:
x,y=y,x
dd=-dd
y=self.upper(y,dd)
if x==self.root:
return x
if x==y:
return x
d=self.vertex_depth(x)
b=d.bit_length()
X=self.upper_list
for k in range(b-1,-1,-1):
px=X[k][x];py=X[k][y]
if px!=py:
x=px;y=py
return self.upper(x,1)
def __degree_count(self):
assert self.__after_seal_check()
if hasattr(self,"deg"):
return
self.deg=[0]*(self.index+self.N)
for v in range(self.index,self.index+self.N):
d=len(self.children[v])+1
if d!=self.root:
d-=1
self.deg[v]=d
return
def degree(self,v):
"""頂点vの次数を求める.
"""
assert self.__after_seal_check(v)
if not hasattr(self,"deg"):
self.__degree_count()
return self.deg[v]
def diameter(self):
"""木の直径を求める.
"""
assert self.__after_seal_check()
from collections import deque
def bfs(start):
X=[-1]*(self.index+self.N)
Q=deque([start])
X[start]=0
pa=self.parent
ch=self.children
while Q:
x=Q.popleft()
if X[pa[x]]==-1:
Q.append(pa[x])
X[pa[x]]=X[x]+1
for y in ch[x]:
if X[y]==-1:
Q.append(y)
X[y]=X[x]+1
y=max(range(self.index,self.index+self.N),key=lambda x:X[x])
return y,X[y]
y,_=bfs(self.root)
z,d=bfs(y)
return y,z,d
def path(self,u,v):
"""頂点u,v間のパスを求める.
"""
assert self.__after_seal_check(u,v)
w=self.lowest_common_ancestor(u,v)
pa=self.parent
X=[u]
while u!=w:
u=pa[u]
X.append(u)
Y=[v]
while v!=w:
v=pa[v]
Y.append(v)
return X+Y[-2::-1]
def is_brother(self,u,v):
"""2つの頂点u,vは兄弟 (親が同じ) か?
"""
assert self.__after_seal_check(u,v)
if u==self.root or v==self.root:
return False
return self.parent[u]==self.parent[v]
def is_ancestor(self,u,v):
"""頂点uは頂点vの先祖か?
"""
assert self.__after_seal_check(u,v)
dd=self.vertex_depth(v)-self.vertex_depth(u)
if dd<0:
return False
v=self.upper(v,dd)
return u==v
def is_descendant(self,u,v):
"""頂点uは頂点vの子孫か?
"""
assert self.__after_seal_check(u,v)
return self.is_ancestor(v,u)
def is_leaf(self,v):
"""頂点vは葉?
"""
return not bool(self.children[v])
def distance(self,u,v):
"""2頂点u,v間の距離を求める.
"""
assert self.__after_seal_check(u,v)
dep=self.vertex_depth
return dep(u)+dep(v)-2*dep(self.lowest_common_ancestor(u,v))
def __descendant_count(self):
assert self.__after_seal_check()
if hasattr(self,"des_count"):
return
if not hasattr(self,"tower"):
self.depth_search(False)
self.des_count=[1]*(self.index+self.N)
pa=self.parent
for T in self.tower[:0:-1]:
for x in T:
self.des_count[pa[x]]+=self.des_count[x]
return
def descendant_count(self,v):
"""頂点vの子孫の数を求める.
"""
assert self.__after_seal_check(v)
self.__descendant_count()
return self.des_count[v]
def subtree_size(self,v):
"""頂点vを根とした部分根付き木のサイズを求める.
"""
return self.descendant_count(v)
def preorder(self,v):
"""頂点vの行きがけ順を求める.
"""
assert self.__after_seal_check(v)
if hasattr(self,"preorder_number"):
self.preorder_number[v]
from collections import deque
Q=deque([self.root])
T=[-1]*(self.N+self.index)
p=1
while Q:
x=Q.popleft()
T[x]=p
p+=1
C=self.children[x]
for y in C:
Q.append(y)
self.preorder_number=T
return T[v]
def dfs_yielder(self):
"""DFSにおける頂点の出入りをyieldする.
(v,1): 頂点vに入る
(v,0): 頂点vを出る
"""
assert self.__after_seal_check()
#最初
yield (self.root,1)
v=self.root
ch=self.children
pa=self.parent
R=[-1]*self.index+[len(ch[x]) for x in range(self.index,self.index+self.N)]
S=[0]*(self.index+self.N)
while True:
if R[v]==S[v]: #もし,進めないならば
yield (v,0) #頂点vを出る
if v==self.root:
break
else:
v=pa[v]
else: #進める
w=v
v=ch[v][S[v]]
S[w]+=1
yield (v,1)
def top_down(self):
assert self.__after_seal_check()
if not hasattr(self,"tower"):
self.depth_search(False)
for E in self.tower:
for v in E:
yield v
def bottom_up(self):
assert self.__after_seal_check()
if not hasattr(self,"tower"):
self.depth_search(False)
for E in self.tower[::-1]:
for v in E:
yield v
def tree_dp(self,calc,unit,f,g,Mode=False):
"""葉から木DPを行う.
[input]
calc:モノイドを成す2項演算 M x M -> M
unit:Mの単位元
f,g: M -> M
Mode: False->根の値のみ, True->全ての値
[補足]
頂点 v の子が x,y,z,...のとき, 更新式は
dp[v]=g(f(x)*f(y)*f(z)*...)
になる.
"""
data=[unit]*(self.index+self.N)
pa=self.parent
for x in self.bottom_up():
if x==self.root:
break
data[x]=g(data[x])
data[pa[x]]=calc(data[pa[x]],f(data[x]))
if Mode:
return data
else:
return data[self.root]
def euler_tour(self):
""" オイラーツアーに関する計算を行う.
"""
assert self.__after_seal_check()
#最初
X=[]; X_append=X.append #X: Euler Tour のリスト
v=self.root
ch=self.children
pa=self.parent
R=[-1]*self.index+[len(ch[x]) for x in range(self.index,self.index+self.N)]
S=[0]*(self.index+self.N)
while True:
X_append(v)
if R[v]==S[v]: #もし,進めないならば
if v==self.root:
break
else:
v=pa[v]
else: #進める
w=v
v=ch[v][S[v]]
S[w]+=1
self.euler=X
self.in_time=[-1]*(self.index+self.N)
self.out_time=[-1]*(self.index+self.N)
for i in range(len(X)):
v=X[i]
if self.in_time[v]==-1:
self.in_time[v]=self.out_time[v]=i
else:
self.out_time[v]=i
#=================================================
def Making_Tree(N,E,root,index=0):
"""木を作る.
N:頂点数
E:辺のリスト
root:根
"""
from collections import deque
F=[[] for _ in range(index+N)]
for u,v in E:
assert index<=u<index+N
assert index<=v<index+N
assert u!=v
F[u].append(v)
F[v].append(u)
X=[-1]*(index+N)
X[root]=root
C=[[] for _ in range(index+N)]
Q=deque([root])
while Q:
x=Q.popleft()
for y in F[x]:
if X[y]==-1:
X[y]=x
Q.append(y)
C[x].append(y)
T=Tree(N,index)
T.root_set(root)
T.parent=X
T.children=C
T.seal()
return T
#==================================================
import sys
input=sys.stdin.readline
N=int(input())
E=[]
for _ in range(N-1):
u,v=map(int,input().split())
E.append((u,v))
T=Making_Tree(N,E,1,1)
Mod=10**9+7
N_fact=1
for i in range(1,N+1):
N_fact*=i
N_fact%=Mod
Inv=[0]*(N+1); Inv[1]=1
for i in range(2,N+1):
q,r=divmod(Mod,i)
Inv[i]=(-q*Inv[r])%Mod
T.depth_search(Mode=False)
dep=T.depth
X=0
for i in range(1,N+1):
X+=N_fact*Inv[dep[i]+1]
X%=Mod
print(X)