結果

問題 No.1333 Squared Sum
ユーザー とりゐ
提出日時 2022-11-07 16:37:38
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,977 ms / 2,000 ms
コード長 1,873 bytes
コンパイル時間 252 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 247,632 KB
最終ジャッジ日時 2024-07-21 00:13:37
合計ジャッジ時間 52,175 ms
ジャッジサーバーID
(参考情報)
judge1 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 44
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

def rerooting():
dp=[[E]*len(edge[v]) for v in range(n)]
# dfs1
memo=[E]*n
for v in order[::-1]:
res=E
for i in range(len(edge[v])):
if edge[v][i]==par[v]:
continue
dp[v][i]=memo[edge[v][i]]
res=merge(res,f(dp[v][i],edge[v][i],v))
memo[v]=g(res,v)
# dfs2
memo2=[E]*n
for v in order:
for i in range(len(edge[v])):
if edge[v][i]==par[v]:
dp[v][i]=memo2[v]
s=len(edge[v])
cumR=[E]*(s+1)
cumR[s]=E
for i in range(s,0,-1):
cumR[i-1]=merge(cumR[i],f(dp[v][i-1],v,edge[v][i-1]))
cumL=E
for i in range(s):
if edge[v][i]!=par[v]:
val=val=merge(cumL,cumR[i+1])
memo2[edge[v][i]]=g(val,v)
cumL=merge(cumL,f(dp[v][i],edge[v][i],v))
ans=[E for i in range(n)]
for v in range(n):
res=E
for i in range(len(edge[v])):
res=merge(res,f(dp[v][i],edge[v][i],v))
ans[v]=g(res,v)
#ans[v]=calc_ans(res,v)
return ans
E=(0,0,0)
mod=10**9+7
def f(res,v,u):
a,b,c=res
w=weight[v+(u<<20)]
a+=2*w*b+c*w%mod*w
a%=mod
b+=c*w
b%=mod
return (a,b,c)
def g(res,v):
a,b,c=res
return (a,b,c+1)
def merge(a,b):
xa,ya,za=a
xb,yb,zb=b
return ((xa+xb)%mod,(ya+yb)%mod,za+zb)
'''
def calc_ans(res,v):
return
'''
from sys import stdin
input=lambda :stdin.readline()[:-1]
weight={}
n=int(input())
edge=[[] for i in range(n)]
for i in range(n-1):
x,y,w=map(int,input().split())
x,y=x-1,y-1
edge[x].append(y)
edge[y].append(x)
weight[x+(y<<20)]=w
weight[y+(x<<20)]=w
# make order table
# root = 0
from collections import deque
order=[]
par=[-1]*n
todo=deque([0])
while todo:
v=todo.popleft()
order.append(v)
for u in edge[v]:
if u!=par[v]:
par[u]=v
todo.append(u)
ans=rerooting()
res=0
for i in ans:
res+=i[0]
res%=mod
res*=pow(2,mod-2,mod)
print(res%mod)
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0