Problem:
イクタ君の住む町にはN個の塔がある。
それぞれの塔には0からN−1までの異なる番号が与えられており、番号iの塔を塔iと呼ぶ。
好奇心旺盛なイクタ君はN個の塔の高さに興味を持ち、それらの大小関係を表す表Tを作ることにした。
TはN × N個の要素を持ち、各要素Ti, j (0 ≤ i, j ≤ N − 1)は次のように定義される。

Ti, j = −1 ⇔ 塔iの高さが塔jの高さより小さい

Ti, j = 0 ⇔ 塔iの高さと塔jの高さが等しい

Ti, j = 1 ⇔ 塔iの高さが塔jの高さより大きい

イクタ君は表Tを作成するための調査として、2つの塔を選びそれらの高さを比較するということをN−1回繰り返した。

イクタ君の調査に関して以下のことがわかっている。

i回目の比較(1 ≤ i ≤ \ N − 1)で塔ai, 塔biが選ばれたとすると塔aiの高さが塔biの高さより大きかった。つまりTai, bi = 1, Tbi, ai = −1であった。

それぞれの塔は自分自身よりも大きな塔とは高々一回しか比較されていない。

残念ながらイクタ君の調査による情報だけで表Tの内容を一意に決めることができるとは限らない。
表Tがイクタ君の調査と矛盾せず、Tが定義されるような塔の高さの組み合わせが存在するときTを正しい表と呼ぶことにする。
正しい表として何通りのものが考えられるかを計算してイクタ君に教えて欲しい。

ただし、比較された２つの塔の高さは互いに異なるがすべての塔の高さが互いに異なるとは限らない。

Input

入力は以下の形式で与えられる。

N
a1 b1
...
aN−1 bN−1

Nは塔の数を表す。
ai, bi (1 ≤ i ≤ N − 1)は塔aiが塔biより高いことを表す。

Constraints

入力中の各変数は以下の制約を満たす。

1 ≤ N ≤ 200

0 ≤ ai, bi < N

ai ≠ bi

異なる塔が同じ高さであることもある。

イクタ君の調査結果自体が矛盾していることはなく、少なくとも1つイクタ君の調査結果と矛盾しない表Tが存在する。

Output

考えられる正しい表Tの個数の数を1,000,000,007で割ったあまりを出力せよ。

Sample Input 1

3
0 1
1 2

Output for the Sample Input 1

1

正しい表Tは上図の1通りである。

Sample Input 2

3
0 1
0 2

Output for the Sample Input 2

3

正しい表Tは上図の3通りである。

Sample Input 3

1

Output for the Sample Input 3

1

正しい表Tは上図の1通りである。

Sample Input 4

7
0 1
1 2
2 3
3 4
0 5
0 6

Output for the Sample Input 4

91

Solution in C++11 (judge verdict: Accepted)

#include "bits/stdc++.h"
#include<unordered_map>
#pragma warning(disable:4996)
using namespace std;
const int mod = 1000000007;
struct Mod {
public:
	int num;
	Mod() : num(0) { ; }
	Mod(long long int n) : num((n % mod + mod) % mod) { ; }
	Mod(int n) : num((n % mod + mod) % mod) { ; }
	operator int() { return num; }
};

Mod operator+(const Mod a, const Mod b) { return Mod((a.num + b.num) % mod); }
Mod operator+(const long long int a, const Mod b) { return Mod(a) + b; }
Mod operator++(Mod &a) { return a + Mod(1); }
Mod operator-(const Mod a, const Mod b) { return Mod((mod + a.num - b.num) % mod); }
Mod operator-(const long long int a, const Mod b) { return Mod(a) - b; }
Mod operator--(Mod &a) { return a - Mod(1); }
Mod operator*(const Mod a, const Mod b) { return Mod(((long long)a.num * b.num) % mod); }
Mod operator*(const long long int a, const Mod b) { return Mod(a)*b; }
Mod operator+=(Mod &a, const Mod b) { return a = a + b; }
Mod operator+=(long long int &a, const Mod b) { return a = a + b; }
Mod operator-=(Mod &a, const Mod b) { return a = a - b; }
Mod operator-=(long long int &a, const Mod b) { return a = a - b; }
Mod operator*=(Mod &a, const Mod b) { return a = a * b; }
Mod operator*=(long long int &a, const Mod b) { return a = a * b; }
Mod operator^(const Mod a, const int n) {
	if (n == 0) return Mod(1);
	Mod res = (a * a) ^ (n / 2);
	if (n % 2) res = res * a;
	return res;
}
Mod mod_pow(const Mod a, const int n) {
	if (n == 0) return Mod(1);
	Mod res = mod_pow((a * a), (n / 2));
	if (n % 2) res = res * a;
	return res;
}
Mod inv(const Mod a) { return a ^ (mod - 2); }
Mod operator/(const Mod a, const Mod b) {
	assert(b.num != 0);
	return a * inv(b);
}
Mod operator/(const long long int a, const Mod b) {
	assert(b.num != 0);
	return Mod(a) * inv(b);
}
Mod operator/=(Mod &a, const Mod b) {
	assert(b.num != 0);
	return a = a * inv(b);
}

#define MAX_MOD_N 1024000

Mod fact[MAX_MOD_N], factinv[MAX_MOD_N];
void init() {
	fact[0] = Mod(1); factinv[0] = 1;
	for (int i = 0; i < MAX_MOD_N - 1; ++i) {
		fact[i + 1] = fact[i] * Mod(i + 1);
		factinv[i + 1] = factinv[i] / Mod(i + 1);
	}
}
Mod comb(const int a, const int b) {
	return fact[a] * factinv[b] * factinv[a - b];
}
struct Node {
	int id;
	int paid;
	vector<int>chid;
	Node(int a) :chid() {
		id = a;
		paid = -1;
	}
};
vector<Node>nodes;
map<int, Mod>getans(const int nid) {
	if (nodes[nid].chid.size() == 0) {
		const map<int, Mod>mp = { {1,1 } };
		return mp;
	}
	else if (nodes[nid].chid.size() == 1) {
		const map<int, Mod>mp(getans(nodes[nid].chid[0]));
		map<int, Mod>nmp;
		for (auto i : mp) {
			nmp[i.first + 1] = i.second;
		}
		return nmp;
	}
	else {
		map<int, Mod> now = getans(nodes[nid].chid[0]);
		for (int i = 1; i < nodes[nid].chid.size(); ++i) {
			const map<int, Mod>next = getans(nodes[nid].chid[i]);
			map<int, Mod>newmap;
			for (const auto j : now) {
				for (const auto k : next) {
					pair< int, Mod>l(j);
					pair< int, Mod>r(k);
					if (l.first < r.first) {
						swap(l, r);
					}
					assert(l.first >= r.first);
					for (int sum = l.first; sum <= l.first + r.first; ++sum) {
						const int kaburi = l.first + r.first - sum;
						const int kaburazu = sum - l.first;

						Mod ans = comb(l.first, kaburi);
						ans *= comb(l.first + kaburazu, kaburazu);
						ans *= l.second*r.second;
						newmap[sum] += ans;
					}

				}
			}
			now = newmap;
		}
		map<int, Mod>nmp;
		for (auto i : now) {
			nmp[i.first + 1] = i.second;
		}
		return nmp;
	}
}
int main() {
	init();
	int N; cin >> N;
	for (int i = 0; i < N; ++i) {
		nodes.push_back(Node(i));
	}
	for (int i = 0; i < N - 1; ++i) {
		int a, b; cin >> a >> b;
		nodes[a].chid.push_back(b);
		nodes[b].paid = a;
	}
	int root;
	for (int i = 0; i < N; ++i) {
		if (nodes[i].paid == -1) {
			root = i;
		}
	}
	map<int, Mod>ans(getans(root));
	Mod sum = 0;
	for (auto i : ans) {
		sum += i.second;
	}
	cout << sum << endl;
	return 0;
}