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

const int maxn = 310000;
const int maxb = 22;
const int mod  = 1e9+7;
inline void add(int &a,const int &b){ a+=b;if(a>=mod)a-=mod; }
inline void sub(int &a,const int &b){ a-=b;if(a<0)a+=mod; }

int n;
vector<int>E[maxn];
int f[maxn][maxb],g[maxn][maxb],mxb[maxn];
int fx[2][maxn],gx[2][maxn];

inline bool cmp(const int x,const int y)
{
	return mxb[x]<mxb[y];
}
void dp(const int x,const int fa)
{
	vector<int>pi;
	for(auto y:E[x]) if(y!=fa)
	{
		dp(y,x);
		pi.push_back(y);
	}
	sort(pi.begin(),pi.end(),cmp);
	
	if(pi.size()>0)
	{
		int now=0,nowsiz= 1<<mxb[pi[0]];
		for(int i=0;i<nowsiz;i++) fx[now][i]=gx[now][i]=0;
		fx[now][0]=1;
		
		for(int sz=pi.size(),i=0;i<sz;i++)
		{
			int y=pi[i];
			
			now=!now;
			for(int j=0;j<(1<<mxb[y]);j++) fx[now][j]=gx[now][j]=0;
			
			for(int j=0;j<nowsiz;j++) if(fx[!now][j]) for(int by=0;by<=mxb[y];by++)
			{
				add(fx[now][j],(ll)fx[!now][j]*f[y][by]%mod);
				if( by>0 && !(j>>(by-1)&1))
					add(fx[now][j|1<<(by-1)],(ll)fx[!now][j]*f[y][by]%mod);
				
				if( by>0 && (1+j) <= (1<<(by-1)) )
					add(gx[now][j|1<<(by-1)],(ll)fx[!now][j]*g[y][by]%mod);
			}
			for(int j=0;j<nowsiz;j++) if(gx[!now][j])
			{
				int k= j==0?-1: 31-__builtin_clz(j);
				for(int by=0;by<=mxb[y];by++)
				{
					add(gx[now][j],(ll)gx[!now][j]*f[y][by]%mod);
					if( by>0 && by-1<k && !(j>>(by-1)&1) )
						add(gx[now][j|1<<(by-1)],(ll)gx[!now][j]*f[y][by]%mod);
				}
			}
			
			nowsiz=1<<mxb[y];
		}
		
		for(int j=0;j<nowsiz;j++)
		{
			/*cerr<<j<<endl;
			cerr<<__builtin_clz((int)j)<<endl;
			cerr<<__builtin_clz((unsigned int)j)<<endl;
			cerr<<__builtin_clz(0)<<endl;*/
			
			int k= j==0?-1: 31-__builtin_clz(j);
			int num= __builtin_popcount(j);
			
			//f
			if(k+1==num)
			{
				add(f[x][k+2],fx[now][j]);
				add(f[x][0],gx[now][j]);
			}
			//g
			if(k+1==num+1)
			{
				for(int l=0;l<k;k++) if( !(j>>l&1) )
				{
					add(g[x][l+1],gx[now][j]);
					break;
				}
			}
		}
		for(int b=1;b<maxb;b++) if(f[x][b])
			add(g[x][b],f[x][b]);
	}
	else f[x][1]=g[x][1]=1;
	for(int b=0;b<maxb;b++) if(f[x][b]||g[x][b])
		mxb[x]=b;
}

signed main()
{
	/////////////////////////////////////////////////////////////////////////
	//freopen("tmp.in","r",stdin);
	ios_base::sync_with_stdio(false);
	cin.tie(0);
	
	//cerr<< ( 31-__builtin_clz(0) ) <<endl;
	
	cin>>n;
	for(int i=1;i<n;i++)
	{
		int x,y; cin>>x>>y;
		E[x].push_back(y);
		E[y].push_back(x);
	}
	dp(1,0);
	
	int ans=0;
	for(int i=0;i<=mxb[1];i++) 
		add(ans,f[1][i]);
	cout<<ans<<endl;
	
	return 0;
}