結果

問題 No.196 典型DP (1)
コンテスト
ユーザー 259_Momone
提出日時 2019-02-26 11:38:21
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 28 ms / 2,000 ms
コード長 1,496 bytes
コンパイル時間 1,409 ms
コンパイル使用メモリ 166,436 KB
実行使用メモリ 27,136 KB
最終ジャッジ日時 2024-06-11 22:03:03
合計ジャッジ時間 2,884 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

long dynamic_programming(const unsigned long K, const long now, const long prev, vector<vector<long>>& dp, const vector<vector<long>>& edge){
    if(edge[now].size() == 1 && edge[now][0] == prev){
        dp[now].push_back(1);
        return 1;
    }
    auto tapi = [&K](vector<long>& to, const vector<long>& src) -> void{
        size_t N = min(to.size() + src.size(), K + 1);
        vector<long> hoge(N);
        for(int i = 0; i < to.size(); ++i){
            if(!to[i])continue;
            for(int j = 0; j < src.size(); ++j){
                if(i + j > K)break;
                hoge[i + j] += to[i] * src[j] % 1000000007;
            }
        }
        for(auto& i : hoge)i %= 1000000007;
        to = vector<long>(hoge);
    };
    long ret = 0;
    dp[now][0] = 1;
    for(const auto i : edge[now])
        if(i != prev) ret += dynamic_programming(K, i, now, dp, edge);
    for(const auto i : edge[now])
        if(i != prev){
            tapi(dp[now], dp[i]);
        }

    if(ret + 1 <= K)dp[now].resize(ret + 2);
    if(ret + 1 < dp[now].size())dp[now][ret + 1] = 1;
    return ret + 1;
}

int main(){
    unsigned long N, K;
    cin >> N >> K;

    vector<vector<long>> edge(N);
    for(int i = 1, a, b; i < N; ++i){
        cin >> a >> b;
        edge[a].push_back(b);
        edge[b].push_back(a);
    }

    vector<vector<long>> dp(N, vector<long>(1, 1));
    dynamic_programming(K, 0, 0, dp, edge);
    cout << dp[0][K] << endl;
}
0