結果

問題 No.196 典型DP (1)
ユーザー startcppstartcpp
提出日時 2018-03-07 21:40:40
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 65 ms / 2,000 ms
コード長 1,716 bytes
コンパイル時間 604 ms
コンパイル使用メモリ 61,668 KB
実行使用メモリ 53,732 KB
最終ジャッジ日時 2024-10-04 12:39:19
合計ジャッジ時間 3,046 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff #

//与えられた木から部分木を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;
	
	int chnum = 0;
	rep(i, childs.size()) {	//今まで子をi個見た
		rep(j, chnum + 1) {	//今まで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;
			}
		}
		chnum += tsz[childs[i]];
	}
	
	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;
}
0