結果
| 問題 |
No.196 典型DP (1)
|
| コンテスト | |
| ユーザー |
startcpp
|
| 提出日時 | 2018-03-07 21:35:17 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,670 bytes |
| コンパイル時間 | 510 ms |
| コンパイル使用メモリ | 60,948 KB |
| 実行使用メモリ | 41,540 KB |
| 最終ジャッジ日時 | 2024-10-04 12:10:18 |
| 合計ジャッジ時間 | 6,140 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | -- * 3 |
| other | AC * 23 TLE * 1 -- * 17 |
ソースコード
//与えられた木から部分木を0個以上カットしてできる大きさN-Kの木は何通り?という問題.
#include <iostream>
#include <vector>
#define int long long
#define rep(i, n) for(i = 0; i < n; i++)
using namespace std;
int mod = 1000000007;
int N, K, u, v;
vector<int> et[2000];
int tsz[2000]; //tsz[i] = 頂点i以下の頂点の個数
int tdp[2000][2000]; //tdp[i][j] = 頂点i以下でj頂点残す方法の数
int cdp[2000][2000]; //cdp[i][j] = 子i個で合計j頂点残す方法の数
bool solvedTdp[2000];
int dfs(int p, int v) {
int i, ret = 1;
rep(i, et[v].size()) {
if (et[v][i] == p) continue;
ret += dfs(v, et[v][i]);
}
return tsz[v] = ret;
}
void dfs2(int p, int v) {
int i, j, k;
vector<int> childs;
rep(i, et[v].size()) {
if (et[v][i] == p) continue;
if (!solvedTdp[et[v][i]]) dfs2(v, et[v][i]);
childs.push_back(et[v][i]);
}
if (childs.size() == 0) { //葉の場合
tdp[v][0] = tdp[v][1] = 1;
solvedTdp[v] = true;
return;
}
rep(i, childs.size() + 1) rep(j, tsz[v]) cdp[i][j] = 0;
cdp[0][0] = 1;
rep(i, childs.size()) { //今まで子をi個見た
rep(j, tsz[v]) { //今までj頂点残した
rep(k, tsz[childs[i]] + 1) { //i番目の子でk頂点残す
cdp[i + 1][j + k] += cdp[i][j] * tdp[childs[i]][k];
cdp[i + 1][j + k] %= mod;
}
}
}
rep(i, tsz[v] + 1) {
if (i == 0) tdp[v][0] = 1;
else tdp[v][i] = cdp[childs.size()][i - 1];
}
solvedTdp[v] = true;
return;
}
signed main() {
int i;
cin >> N >> K;
rep(i, N - 1) {
cin >> u >> v;
et[u].push_back(v);
et[v].push_back(u);
}
dfs(0, 0);
dfs2(0, 0);
cout << tdp[0][N - K] << endl;
return 0;
}
startcpp