結果

問題 No.196 典型DP (1)
ユーザー maspy
提出日時 2020-01-01 14:21:44
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 660 ms / 2,000 ms
コード長 1,100 bytes
コンパイル時間 383 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 57,068 KB
最終ジャッジ日時 2024-11-22 03:23:11
合計ジャッジ時間 29,395 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines

import numpy as np

N,K = map(int,readline().split())
m = map(int,read().split())
AB = zip(m,m)

MOD = 10 ** 9 + 7

graph = [[] for _ in range(N+1)]
for a,b in AB:
    a+=1; b+=1
    graph[a].append(b)
    graph[b].append(a)

root = 1
parent = [0] * (N+1)
order = []
stack = [root]
while stack:
    x = stack.pop()
    order.append(x)
    for y in graph[x]:
        if y == parent[x]:
            continue
        parent[y] = x
        stack.append(y)

size = [1] * (N+1) # 部分木の大きさ
for v in order[::-1]:
    p = parent[v]
    size[p] += size[v]

def convolve(A,B):
    LA = len(A); LB = len(B)
    if LA > LB:
        A,B = B,A
        LA,LB = LB,LA
    C = np.zeros(LA+LB-1,np.int64)
    for i,x in enumerate(A):
        C[i:i+LB] += A[i] * B % MOD
    C %= MOD
    return C

dp = [np.array([1],np.int64) for _ in range(N+1)]

for v in order[::-1]:
    p = parent[v]
    dp[v] = np.append(dp[v],1)
    dp[p] = convolve(dp[p], dp[v])

answer = dp[root][K]
print(answer)
0