#pragma GCC optimize(2,3,"Ofast","inline")
#include <bits/stdc++.h>
#define int unsigned long long
#define double long double
#define lowbit(i) (i&(-i))
#define mid ((l+r)>>1)
#define add(i,j) (((i+j)>=mod)?i+j-mod:i+j)
using namespace std;
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int mod=1e9+7;
int siz[300005];
vector<int> vc[300005];
int cnt[300005];
int dp[300005][35][3];
int f[2400005][2],g[2400005][2];
int tot=0;
void dfs(int now,int fa){
	for(auto v:vc[now]){
		if(v==fa) continue;
		dfs(v,now);
	}
	for(auto v:vc[now]) if(v!=fa) cnt[siz[v]]++;
	for(int i=20;i>=0;i--){
		cnt[i]+=cnt[i+1];
		if(i==0) break;
	}
	for(int i=20;i>=0;i--){
		bool ok=1;
		for(int j=0;j<i;j++){
			if(cnt[j]<i-j) ok=0;
		}
		if(ok){
			siz[now]=i;
			break;
		}
		if(i==0) break;
	}
	tot+=(1<<(siz[now]+1));
//	if(siz[now]>20) exit(2);
//	if(tot>=16000000) exit(2);
	for(int i=0;i<=20;i++) cnt[i]=0;
	for(int i=0;i<(1<<(siz[now]+1));i++) f[i][0]=f[i][1]=0;
	f[0][0]=1;
	for(auto v:vc[now]){
		if(v==fa) continue;
//		for(int i=0;i<(1<<(siz[now]+2));i++) g[i][0]=g[i][1]=0;
		int tot=0;
		for(int j=0;j<=siz[v];j++) (tot+=dp[v][j][0]+dp[v][j][2])%=mod;
		for(int j=0;j<=min(siz[v],siz[now]);j++){
			for(int i=0;i<(1<<(siz[now]+1));i++){
				if(!((i>>j)&1)){
					g[i|(1<<j)][0]=add(g[i|(1<<j)][0],f[i][0]*dp[v][j][0]);
					if((1<<j)<i) g[i|(1<<j)][1]=add(g[i|(1<<j)][1],f[i][1]*dp[v][j][0]);
					if((1<<j)>i) g[i|(1<<j)][1]=add(g[i|(1<<j)][1],f[i][0]*dp[v][j][1]);
				}
			}
		}
		for(int i=0;i<(1<<(siz[now]+1));i++) f[i][0]=(f[i][0]*tot+g[i][0])%mod,f[i][1]=(f[i][1]*tot+g[i][1])%mod,g[i][0]=0,g[i][1]=0;
	}
	for(int i=0;i<=siz[now];i++){
		dp[now][i][0]=f[(1<<i)-1][0];
		dp[now][i][2]=f[(1<<i)-1+(1<<i)][1];
	}
	for(int i=0;i<=siz[now]+1;i++){
		for(int j=0;j<i;j++){
			dp[now][j][1]=add(dp[now][j][1],add(f[(1<<i)-1-(1<<j)][0],f[(1<<i)-1-(1<<j)+(1<<i)][1]));
//			if(now==4&&i==1&&j==0) cout<<"OK"<<f[(1<<i)-1-(1<<j)][0]<<" "<<f[(1<<i)-1-(1<<j)+(1<<i)][1]<<" "<<dp[now][j][1]<<" ";
		}
	}
//	cout<<now<<" "<<siz[now]<<"\n";
//	for(int i=0;i<=siz[now]+1;i++) cout<<dp[now][i][0]<<" "<<dp[now][i][1]<<" "<<dp[now][i][2]<<"\n";
}
signed main(){
	int n=read(); //cin>>n;
	for(int i=1;i<n;i++){
		int u=read(),v=read(); //cin>>u>>v;
		vc[u].push_back(v);
		vc[v].push_back(u);
	}
	dfs(1,0);
	int ans=0;
	for(int i=0;i<=siz[1];i++) (ans+=dp[1][i][0]+dp[1][i][2])%=mod;
	cout<<ans;
	return 0;
}