結果

問題 No.1075 木の上の山
ユーザー mkawa2
提出日時 2021-10-12 23:55:35
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 245 ms / 2,000 ms
コード長 1,974 bytes
コンパイル時間 591 ms
コンパイル使用メモリ 82,508 KB
実行使用メモリ 110,532 KB
最終ジャッジ日時 2024-09-17 11:32:52
合計ジャッジ時間 5,839 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge6
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys

# sys.setrecursionlimit(200005)
int1 = lambda x: int(x)-1
p2D = lambda x: print(*x, sep="\n")
def II(): return int(sys.stdin.readline())
def LI(): return list(map(int, sys.stdin.readline().split()))
def LI1(): return list(map(int1, sys.stdin.readline().split()))
def LLI(rows_number): return [LI() for _ in range(rows_number)]
def LLI1(rows_number): return [LI1() for _ in range(rows_number)]
def SI(): return sys.stdin.readline().rstrip()
# dij = [(0, 1), (-1, 0), (0, -1), (1, 0)]
dij = [(0, 1), (-1, 0), (0, -1), (1, 0), (1, 1), (1, -1), (-1, 1), (-1, -1)]
inf = 10**16
# md = 998244353
md = 10**9+7

n, k = LI()
to = [[] for _ in range(n)]
for _ in range(n-1):
    u, v = LI1()
    to[u].append(v)
    to[v].append(u)

# ***************ここから本体*****************

# 設定
root = 0
e = lambda: [1]*k  # 単位元

def op(a, b):
    return [a*b%md for a, b in zip(a, b)]

def trans_up(u):
    for i in range(len(dp[u])-1):
        dp[u][i+1] += dp[u][i]
        dp[u][i+1] %= md
    return

def trans_down(u):
    for i in range(len(rdp[u])-1):
        rdp[u][i+1] += rdp[u][i]
        rdp[u][i+1] %= md
    return

# dfs順を作る
uu = []
stack = [root]
par = [-1]*n
while stack:
    u = stack.pop()
    uu.append(u)
    vv = []
    # ↓変数注意
    for v in to[u]:
        if v == par[u]: continue
        vv.append(v)
        par[v] = u
        stack.append(v)
    to[u] = vv

dp = [e() for _ in range(n)]
rdp = [e() for _ in range(n)]
# 下から更新
for u in uu[::-1]:
    p = par[u]
    if p != -1:
        trans_up(u)
        dp[p] = op(dp[p], dp[u])
# print(dp)

# 上から更新
ans = 1
for u in uu:
    cs = op(e(), rdp[u])
    # ↓変数注意
    for v in to[u]:
        rdp[v] = cs
        cs = op(cs, dp[v])
    cs = e()
    for v in to[u][::-1]:
        rdp[v] = op(rdp[v], cs)
        trans_down(v)
        cs = op(cs, dp[v])
    cs = op(cs, [0]+rdp[u][:-1])
    # print(u,cs)
    for a in cs: ans = (ans+a)%md

print(ans)
0