#include<bits/stdc++.h>
using namespace std;
namespace IO_ER{
    #define LL long long
    #define db double
    #define U unsigned
    #define ULL U LL
    #define In inline
    #define Re register
    #define f(a,b,i) for(int i=a,MAXN##i=b;i<=MAXN##i;i++)
    #define ff(a,b,i) for(int i=a;i<b;i++)
    #define f_(a,b,i) for(int i=a;i>=b;i--)
    #define ff_(a,b,i) for(int i=a;i>b;i--)
    typedef pair<LL,int> Pi;
    int inf=0x3f3f3f3f;
    int INF=0x7fffffff;
    LL infll=0x3f3f3f3f3f3f3f3fll;
    LL INFll=0x7fffffffffffffffll;
    template<typename T>void read(T &x){
        x=0;
        int fl=0;
        char ch=getchar();
        while(ch<'0'||'9'<ch){
            if(ch=='-')fl=1;
            ch=getchar();
        }
        while('0'<=ch&&ch<='9'){
            x=(x<<1)+(x<<3)+(ch^48);
            ch=getchar();
        }
        x=fl?-x:x;
    }
    template<typename T,typename ...Args>void read(T &x,Args &...args){
        read(x);
        read(args...);
    }
}
using namespace IO_ER;

#define N 300050
#define mod 1000000007

template<typename T>In T mul(const T &a,const T &b){return 1ll*a*b%mod;}
template<typename T,typename ...Args>In T mul(const T &a,const T &b,const Args &...args){return mul(a,mul(b,args...));}

template<typename T>In T add(const T &a,const T &b){return a+b>=mod?a+b-mod:a+b;}
template<typename T,typename ...Args>In T add(const T &a,const T &b,const Args &...args){return add(a,add(b,args...));}

template<typename T>In void cadd(T &a,const T &b){a+=b;a>=mod?a-=mod:1;}
template<typename T,typename ...Args>In void cadd(T &a,const T &b,const Args &...args){cadd(a,b);cadd(a,args...);}

typedef vector<int> Vi;
typedef vector<array<int,2> > Vai;

int n;

Vi e[N];

int sz[N];

In bool cmp(int x,int y){return sz[x]<sz[y];}

int dp[N][20][2];

int sum[N];

In void dfs(int u,int fa){
    sz[u]=1;
    if(fa)e[u].erase(find(e[u].begin(),e[u].end(),fa));
    for(auto v:e[u]){
        dfs(v,u);
        sz[u]+=sz[v];
    }
    sort(e[u].begin(),e[u].end(),cmp);
    Vi g(4),h(4);
    g[0]=1;
    int lim=0;
    for(auto v:e[u]){
        if(v==e[u].back())lim++;
        else lim=__lg(sz[v]);
        int mx=(1<<(lim+2))-1;
        g.resize(mx+2);
        h.resize(mx+2);
        f(0,mx,s)h[s]=mul(g[s],sum[v]);
        f(0,lim+2,j){
            if(!dp[v][j][0])continue;
            f(0,mx,s){
                if(s>>j&1)continue;
                cadd(h[s|(1<<j)],mul(g[s],dp[v][j][0]));
            }
        }  
        swap(g,h);
    }
    f(0,lim+2,i){
        dp[u][i][0]=g[(1<<i)-1];
        cadd(sum[u],dp[u][i][0]);
        ff(0,i,j)cadd(dp[u][j][1],g[(1<<i)-(1<<j)-1]);
    }

    Vai gg(4);
    Vai hh(4);
    gg[0][0]=1;
    lim=0;
    for(auto v:e[u]){
        if(v==e[u].back())lim++;
        else lim=__lg(sz[v]);
        int mx=(1<<(lim+2))-1;
        gg.resize(mx+2);
        hh.resize(mx+2);
        f(0,mx,s){
            hh[s][0]=mul(gg[s][0],sum[v]);
            hh[s][1]=mul(gg[s][1],sum[v]);
        }
        f(0,lim+2,j){
            if(dp[v][j][0]){
                f(0,mx,s){
                    if(s>>j&1)continue;
                    cadd(hh[s|(1<<j)][0],mul(gg[s][0],dp[v][j][0]));
                    if((1<<j)<s)cadd(hh[s|(1<<j)][1],mul(gg[s][1],dp[v][j][0]));
                }
            }
            if(dp[v][j][1]){
                f(0,(1<<j)-1,s){
                    cadd(hh[s^(1<<j)][1],mul(gg[s][0],dp[v][j][1]));
                }
            }
        }
        swap(gg,hh);
    }
    f(0,lim+2,i){
        cadd(sum[u],gg[(1<<i)-1][1]);
        ff(0,i-1,j)cadd(dp[u][j][1],gg[(1<<i)-(1<<j)-1][1]);
    }
}

int main(){
    read(n);
    int u,v;
    ff(1,n,i){
        read(u,v);
        e[u].emplace_back(v);
        e[v].emplace_back(u);
    }
    dfs(1,0);
    printf("%d",sum[1]);
    return 0;
}