結果

問題 No.196 典型DP (1)
ユーザー Mao-betaMao-beta
提出日時 2024-04-04 02:17:35
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 1,195 ms / 2,000 ms
コード長 1,963 bytes
コンパイル時間 163 ms
コンパイル使用メモリ 13,056 KB
実行使用メモリ 64,768 KB
最終ジャッジ日時 2024-10-01 00:13:17
合計ジャッジ時間 24,394 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
import math
import bisect
from heapq import heapify, heappop, heappush
from collections import deque, defaultdict, Counter
from functools import lru_cache
from itertools import accumulate, combinations, permutations, product

sys.setrecursionlimit(1000000)
MOD = 10 ** 9 + 7
MOD99 = 998244353

input = lambda: sys.stdin.readline().strip()
NI = lambda: int(input())
NMI = lambda: map(int, input().split())
NLI = lambda: list(NMI())
SI = lambda: input()
SMI = lambda: input().split()
SLI = lambda: list(SMI())
EI = lambda m: [NLI() for _ in range(m)]


def adjlist(n, edges, directed=False, in_origin=1):
    if len(edges) == 0:
        return [[] for _ in range(n)]

    weighted = True if len(edges[0]) > 2 else False
    if in_origin == 1:
        if weighted:
            edges = [[x-1, y-1, w] for x, y, w in edges]
        else:
            edges = [[x-1, y-1] for x, y in edges]

    res = [[] for _ in range(n)]

    if weighted:
        for u, v, c in edges:
            res[u].append([v, c])
            if not directed:
                res[v].append([u, c])

    else:
        for u, v in edges:
            res[u].append(v)
            if not directed:
                res[v].append(u)

    return res


def main():
    N, K = NMI()
    AB = EI(N-1)
    # i以下の部分木でj個塗る
    dp = [[0]*(N+1) for _ in range(N)]
    P = [N] * N
    G = adjlist(N, AB, in_origin=0)
    # i以下の部分木のサイズ
    SZ = [0] * N

    def dfs(now):
        SZ[now] = 1
        dp[now][0] = 1
        for v in G[now]:
            if v == P[now]:
                continue
            P[v] = now
            dfs(v)
            for i in range(SZ[now], -1, -1):
                for j in range(SZ[v], 0, -1):
                    dp[now][i+j] += dp[now][i] * dp[v][j] % MOD
                    dp[now][i+j] %= MOD
            SZ[now] += SZ[v]
        dp[now][SZ[now]] = 1

    dfs(0)
    print(dp[0][K])


if __name__ == "__main__":
    main()
0