結果
問題 | No.1333 Squared Sum |
ユーザー |
![]() |
提出日時 | 2022-11-07 16:42:11 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,648 ms / 2,000 ms |
コード長 | 1,845 bytes |
コンパイル時間 | 351 ms |
コンパイル使用メモリ | 82,316 KB |
実行使用メモリ | 214,232 KB |
最終ジャッジ日時 | 2024-07-21 00:17:42 |
合計ジャッジ時間 | 43,596 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 44 |
ソースコード
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],weight[v][i])) 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,weight[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],weight[v][i])) ans=0 for v in range(n): res=E for i in range(len(edge[v])): res=merge(res,f(dp[v][i],edge[v][i],weight[v][i])) ans+=g(res,v)[0] #ans[v]=calc_ans(res,v) return ans E=(0,0,0) mod=10**9+7 def f(res,v,w): a,b,c=res 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] n=int(input()) edge=[[] for i in range(n)] weight=[[] 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].append(w) weight[y].append(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() ans*=pow(2,mod-2,mod) print(ans%mod)