結果
問題 | No.1333 Squared Sum |
ユーザー |
|
提出日時 | 2024-06-27 18:13:21 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,831 ms / 2,000 ms |
コード長 | 5,521 bytes |
コンパイル時間 | 1,315 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 290,372 KB |
最終ジャッジ日時 | 2024-06-27 18:14:03 |
合計ジャッジ時間 | 39,140 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 44 |
ソースコード
import os,sys,random,threading#sys.exit() 退出程序from random import randint,choice,shuffle#randint(a,b)从[a,b]范围随机选择一个数#choice(seq)seq可以是一个列表,元组或字符串,从seq中随机选取一个元素#shuffle(x)将一个可变的序列x中的元素打乱from copy import deepcopyfrom io import BytesIO,IOBasefrom types import GeneratorTypefrom functools import lru_cache,reduce#reduce(op,迭代对象)from bisect import bisect_left,bisect_right#bisect_left(x) 大于等于x的第一个下标#bisect_right(x) 大于x的第一个下标from collections import Counter,defaultdict,dequefrom itertools import accumulate,combinations,permutations#accumulate(a)用a序列生成一个累积迭代器,一般list化前面放个[0]做前缀和用#combinations(a,k)a序列选k个 组合迭代器#permutations(a,k)a序列选k个 排列迭代器from heapq import heapify,heappop,heappush#heapify将列表转为堆from typing import Generic,Iterable,Iterator,TypeVar,Union,Listfrom string import ascii_lowercase,ascii_uppercase,digits#小写字母,大写字母,十进制数字from math import ceil,floor,sqrt,pi,factorial,gcd,log,log10,log2,inf#ceil向上取整,floor向下取整 ,sqrt开方 ,factorial阶乘from decimal import Decimal,getcontext#Decimal(s) 实例化Decimal对象,一般使用字符串#getcontext().prec=100 修改精度from sys import stdin, stdout, setrecursionlimitinput = lambda: sys.stdin.readline().rstrip("\r\n")MI = lambda :map(int,input().split())li = lambda :list(MI())ii = lambda :int(input())mod = int(1e9 + 7) #998244353inf = int(1e20)py = lambda :print("YES")pn = lambda :print("NO")DIRS = [(0, 1), (1, 0), (0, -1), (-1, 0)] # 右下左上DIRS8 = [(0, 1), (1, 1), (1, 0), (1, -1), (0, -1), (-1, -1), (-1, 0),(-1, 1)] # →↘↓↙←↖↑↗import os,sys,random,threading#sys.exit() 退出程序#sys.setrecursionlimit(10**6) #调整栈空间from random import randint,choice,shuffle#randint(a,b)从[a,b]范围随机选择一个数#choice(seq)seq可以是一个列表,元组或字符串,从seq中随机选取一个元素#shuffle(x)将一个可变的序列x中的元素打乱from copy import deepcopyfrom io import BytesIO,IOBasefrom types import GeneratorTypefrom functools import lru_cache,reduce#reduce(op,迭代对象)from bisect import bisect_left,bisect_right#bisect_left(x) 大于等于x的第一个下标#bisect_right(x) 大于x的第一个下标from collections import Counter,defaultdict,dequefrom itertools import accumulate,combinations,permutations#accumulate(a)用a序列生成一个累积迭代器,一般list化前面放个[0]做前缀和用#combinations(a,k)a序列选k个 组合迭代器#permutations(a,k)a序列选k个 排列迭代器from heapq import heapify,heappop,heappush#heapify将列表转为堆from typing import Generic,Iterable,Iterator,TypeVar,Union,Listfrom string import ascii_lowercase,ascii_uppercase,digits#小写字母,大写字母,十进制数字from math import ceil,floor,sqrt,pi,factorial,gcd,log,log10,log2,inf#ceil向上取整,floor向下取整 ,sqrt开方 ,factorial阶乘from decimal import Decimal,getcontext#Decimal(s) 实例化Decimal对象,一般使用字符串#getcontext().prec=100 修改精度from sys import stdin, stdout, setrecursionlimitinput = lambda: sys.stdin.readline().rstrip("\r\n")MI = lambda :map(int,input().split())li = lambda :list(MI())ii = lambda :int(input())mod = int(1e9 + 7) #998244353inf = 1<<60py = lambda :print("YES")pn = lambda :print("NO")DIRS = [(0, 1), (1, 0), (0, -1), (-1, 0)] # 右下左上DIRS8 = [(0, 1), (1, 1), (1, 0), (1, -1), (0, -1), (-1, -1), (-1, 0),(-1, 1)] # →↘↓↙←↖↑↗def bootstrap(f, stack=[]): #yielddef wrappedfunc(*args, **kwargs):if stack:return f(*args, **kwargs)else:to = f(*args, **kwargs)while True:if type(to) is GeneratorType:stack.append(to)to = next(to)else:stack.pop()if not stack:breakto = stack[-1].send(to)return toreturn wrappedfuncn=ii()g=[[] for _ in range(n+1)]for _ in range(n-1):u,v,w=li()w%=modg[u].append((v,w))g[v].append((u,w))f1=[0]*(n+1) #i子树内所有点到i点的 距离 平方和f2=[0]*(n+1) #i子树内所有点到i点的 距离 和f3=[1]*(n+1) #i子树内所有点的 个数 之和#u是父节点,v是子节点,v的的信息合并到u中#f1[u]+=f1[v]+2*f2[v]+f3[v]#f2[u]+=f2[v]+f3[v]#f3[u]+=f3[v]@bootstrapdef dfs(u,fa):for v,w in g[u]:if v==fa:continueyield dfs(v,u)f1[u]+=f1[v]+(2*f2[v]*w)%mod+(f3[v]*w*w)%modf1[u]%=modf2[u]+=f2[v]+(f3[v]*w)%modf2[u]%=modf3[u]+=f3[v]f3[u]%=modyield Nonedfs(1,-1)res=[0]*(n+1)@bootstrapdef dfs2(u,fa):t1=t2=t3=0t3+=1for v,w in g[u]:t1=(t1+f1[v]+(2*f2[v]*w)%mod+(f3[v]*w*w)%mod)%modt2=(t2+f2[v]+(f3[v]*w)%mod)%modt3=(t3+f3[v])%modres[u]=t1for v,w in g[u]:if v==fa:continuef1[u]=(t1-(f1[v]+(2*f2[v]*w)%mod+(f3[v]*w*w)%mod))%modf2[u]=(t2-(f2[v]+(f3[v]*w)%mod))%modf3[u]=(t3-f3[v])%modyield dfs2(v,u)yield Nonedfs2(1,-1)tot=sum(res)%modprint(tot*pow(2,mod-2,mod)%mod)