結果

問題 No.1333 Squared Sum
ユーザー Kiri8128Kiri8128
提出日時 2021-01-08 22:54:05
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,149 ms / 2,000 ms
コード長 1,580 bytes
コンパイル時間 335 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 197,152 KB
最終ジャッジ日時 2024-11-16 14:37:13
合計ジャッジ時間 31,060 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 44 ms
53,376 KB
testcase_01 AC 43 ms
53,888 KB
testcase_02 AC 48 ms
53,760 KB
testcase_03 AC 1,149 ms
178,332 KB
testcase_04 AC 1,119 ms
181,268 KB
testcase_05 AC 1,148 ms
181,396 KB
testcase_06 AC 1,144 ms
181,068 KB
testcase_07 AC 1,120 ms
181,428 KB
testcase_08 AC 1,120 ms
178,296 KB
testcase_09 AC 1,128 ms
178,132 KB
testcase_10 AC 1,137 ms
180,956 KB
testcase_11 AC 1,108 ms
178,900 KB
testcase_12 AC 1,119 ms
177,332 KB
testcase_13 AC 719 ms
182,496 KB
testcase_14 AC 1,098 ms
184,860 KB
testcase_15 AC 1,108 ms
182,024 KB
testcase_16 AC 42 ms
53,504 KB
testcase_17 AC 43 ms
53,760 KB
testcase_18 AC 42 ms
53,760 KB
testcase_19 AC 43 ms
53,632 KB
testcase_20 AC 42 ms
53,760 KB
testcase_21 AC 43 ms
53,760 KB
testcase_22 AC 43 ms
54,016 KB
testcase_23 AC 52 ms
54,016 KB
testcase_24 AC 42 ms
53,888 KB
testcase_25 AC 44 ms
53,632 KB
testcase_26 AC 1,095 ms
182,564 KB
testcase_27 AC 1,124 ms
182,844 KB
testcase_28 AC 1,118 ms
182,580 KB
testcase_29 AC 733 ms
181,920 KB
testcase_30 AC 494 ms
119,044 KB
testcase_31 AC 301 ms
98,176 KB
testcase_32 AC 661 ms
138,772 KB
testcase_33 AC 538 ms
126,064 KB
testcase_34 AC 872 ms
158,824 KB
testcase_35 AC 679 ms
137,608 KB
testcase_36 AC 437 ms
113,520 KB
testcase_37 AC 429 ms
115,056 KB
testcase_38 AC 534 ms
119,088 KB
testcase_39 AC 767 ms
145,656 KB
testcase_40 AC 965 ms
197,152 KB
testcase_41 AC 990 ms
195,520 KB
testcase_42 AC 978 ms
195,792 KB
testcase_43 AC 968 ms
194,780 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

# 全方位木DP
import sys
input = sys.stdin.readline
from collections import deque

mod = 10 ** 9 + 7
N = int(input())
X = [[] for i in range(N)]
for i in range(N-1):
    x, y, w = map(int, input().split())
    X[x-1].append((y-1, w))
    X[y-1].append((x-1, w))

P = [-1] * N
Q = deque([0])
R = []
D = [0] * N
while Q:
    i = deque.popleft(Q)
    R.append(i)
    for a, w in X[i]:
        if a != P[i]:
            P[a] = i
            D[a] = w
            X[a].remove((i, w))
            deque.append(Q, a)

# 
# a0 += b0
# a1 += b1 + b0 * d
# a2 += b0 * d ** 2 + 2 * b1 * d + b2

##### Settings
unit = (0, 0, 0)
merge = lambda a, b: ((a[0] + b[0]) % mod, (a[1] + b[1]) % mod, (a[2] + b[2]) % mod)

## ↓ この2つは必ず対称的にする（じゃないとバグる）
def adj_bu(a, i, p):
    a0, a1, a2 = a
    a0 += 1
    d = D[i]
    return (a0, (a1 + a0 * d) % mod, (a0 * d ** 2 + 2 * a1 * d + a2) % mod)
def adj_td(a, i, p):
    a0, a1, a2 = a
    a0 += 1
    d = D[i]
    return (a0, (a1 + a0 * d) % mod, (a0 * d ** 2 + 2 * a1 * d + a2) % mod)

adj_fin = lambda a, i: a[2]
#####

ME = [unit] * N
XX = [0] * N
TD = [unit] * N
for i in R[1:][::-1]:
    p = P[i]
    XX[i] = adj_bu(ME[i], i, p)
    ME[p] = merge(ME[p], XX[i])
XX[R[0]] = adj_fin(ME[R[0]], R[0])

for i in R:
    ac = TD[i]
    for j, _ in X[i]:
        TD[j] = ac
        ac = merge(ac, XX[j])
    ac = unit
    for j, _ in X[i][::-1]:
        TD[j] = adj_td(merge(TD[j], ac), j, i)
        ac = merge(ac, XX[j])
        XX[j] = adj_fin(merge(ME[j], TD[j]), j)

print(sum(XX) * (mod + 1) // 2 % mod)
0