#include<bits/stdc++.h>
using namespace std ;
typedef long long ll ;
typedef unsigned long long ull ;
typedef pair < int , int > pii ;
typedef vector < int > vi ;
#define fi first
#define se second
mt19937 rng(chrono::high_resolution_clock::now().time_since_epoch().count());

#define rep(i, a, b) for(int i = a; i < (b); ++i)
#define all(x) begin(x), end(x)
#define sz(x) (int)(x).size()


const int MAXN = 3e5 + 7 ;
const int MOD = 1e9 + 7 ;
const int LOG = 22 ;

int n ;
vector < int > v[ MAXN ] ;
int subtree[ MAXN ] ;


ll dp[ MAXN ][ LOG ] , wtf[ MAXN ][ LOG ] ; //
ll root[ MAXN ][ LOG ] ; //
ll complete[ MAXN ] ;

int mxpw[ MAXN ] ;
ll f[ ( 1 << LOG ) ][ LOG ] ;
int bit_rv[ 2 * ( 1 << LOG ) ] ;

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

void dfs ( int x , int prv ) {
    vector < int > children ;
    for ( auto y : v[ x ] ) {
        if ( y == prv ) { continue ; }
        dfs ( y , x ) ;
        children.push_back ( y ) ;
    }
    sort ( children.begin ( ) , children.end ( ) , cmp ) ;
    mxpw[ x ] = 0 ;
    subtree[ x ] = 1 ;
    for ( int i = 0 ; i < LOG ; ++ i ) {
        f[ 0 ][ i ] = 0 ;
    }
    f[ 0 ][ LOG - 1 ] = 1 ;
    for ( auto y : children ) {
        subtree[ x ] += subtree[ y ] ;
        int rem = mxpw[ x ] ;
        // mxpw[ x ] = max ( mxpw[ x ] , min ( mxpw[ x ] + 1 , mxpw[ y ] + 1 ) ) ;
        while ( mxpw[ x ] < LOG - 1 && ( 1 << mxpw[ x ] ) < 2 * subtree[ x ] ) { ++ mxpw[ x ] ; }
        if ( rem < mxpw[ x ] ) {
            int lo = ( 1 << rem ) ;
            int hi = ( 1 << mxpw[ x ] ) - 1 ;
            for ( int i = lo ; i <= hi ; ++ i ) {
                for ( int j = 0 ; j < LOG ; ++ j ) {
                    f[ i ][ j ] = 0 ;
                }
            }
        }
        for ( int i = ( 1 << mxpw[ x ] ) - 1 ; i >= 0 ; -- i ) {
            for ( int dig = 0 ; dig <= mxpw[ x ] ; ++ dig ) {
                wtf[ i ][ dig ] = 0 ;
            }
            wtf[ i ][ LOG - 1 ] = 0 ;
        }
        for ( int i = ( 1 << mxpw[ x ] ) - 1 ; i >= 0 ; -- i ) {
            for ( int j = 0 ; j <= mxpw[ y ] ; ++ j ) { 
                if ( ( i & ( 1 << j ) ) == 0 ) { 
                    for ( int dig = 0 ; dig <= mxpw[ x ] ; ++ dig ) {
                        if ( j >= dig ) { continue ; }
                        wtf[ ( i + ( 1 << j ) ) ][ dig ] = ( wtf[ ( i + ( 1 << j ) ) ][ dig ] + f[ i ][ dig ] * root[ y ][ j ] ) % MOD ;
                    }
                    wtf[ ( i + ( 1 << j ) ) ][ LOG - 1 ] = ( wtf[ ( i + ( 1 << j ) ) ][ LOG - 1 ] + f[ i ][ LOG - 1 ] * root[ y ][ j ] ) % MOD ;
                    wtf[ i ][ j ] = ( wtf[ i ][ j ] + f[ i ][ LOG - 1 ] * dp[ y ][ j ] ) % MOD ; 
                }
                wtf[ i ][ j ] = ( wtf[ i ][ j ] + f[ i ][ j ] * complete[ y ] ) % MOD ;
            }
            wtf[ i ][ LOG - 1 ] = ( wtf[ i ][ LOG - 1 ] + f[ i ][ LOG - 1 ] * complete[ y ] ) % MOD ;
        }
        for ( int i = ( 1 << mxpw[ x ] ) - 1 ; i >= 0 ; -- i ) {
            for ( int dig = 0 ; dig <= mxpw[ x ] ; ++ dig ) {
                f[ i ][ dig ] = wtf[ i ][ dig ] ;
            }
            f[ i ][ LOG - 1 ] = wtf[ i ][ LOG - 1 ] ;
        }
    }
    /**
    for ( int i = 0 ; i < ( 1 << mxpw[ x ] ) ; ++ i ) {
        for ( int j = 0 ; j <= mxpw[ x ] ; ++ j )  {
            printf ( "vertex %d --> f[ %d ][ %d ] = %lld\n" , x , i , j , f[ i ][ j ] ) ;
        }
        printf ( "vertex %d --> f[ %d ][ -1 ] = %lld\n" , x , i , f[ i ][ LOG - 1 ] ) ;
    }
    **/
    for ( int mask = 0 ; mask < ( 1 << mxpw[ x ] ) ; ++ mask ) {
        int lw = bit_rv[ ( ( ~mask ) & -( ~mask ) ) ] ; // smallest 0 bit id
        for ( int p1 = 0 ; p1 < LOG ; ++ p1 ) {
            if ( ( mask & ( 1 << p1 ) ) > 0 ) { continue ; }
            if ( p1 == LOG - 1 ) { //
                if ( ( ( mask + ( 1 << lw ) ) & ( mask + ( 1 << lw ) + 1 ) ) == 0 ) { 
                    dp[ x ][ lw ] = ( dp[ x ][ lw ] + f[ mask ][ p1 ] ) % MOD ;
                }
                if ( ( mask & ( mask + 1 ) ) == 0 ) {
                    root[ x ][ lw ] = ( root[ x ][ lw ] + f[ mask ][ p1 ] ) % MOD ;
                    complete[ x ] = ( complete[ x ] + f[ mask ][ p1 ] ) % MOD ;
                }
            }
            else {
                if ( ( mask + ( 1 << lw ) ) == ( 1 << p1 ) - 1 ) {
                    dp[ x ][ lw ] = ( dp[ x ][ lw ] + f[ mask ][ p1 ] ) % MOD ;
                }
                if ( mask == ( 1 << p1 ) - 1 ) {
                    complete[ x ] = ( complete[ x ] + f[ mask ][ p1 ] ) % MOD ;
                }
            }
        }
    }
    /**
    printf ( "----- vertex %d\n" , x ) ;
    printf ( "complete[ %d ] = %lld\n" , x , complete[ x ] ) ;
    for ( int i = 0 ; i <= mxpw[ x ] ; ++ i ) {
        printf ( "dp[ %d ][ %d ] = %lld\n" , x , i , dp[ x ][ i ] ) ;
        printf ( "root[ %d ][ %d ] = %lld\n" , x , i , root[ x ][ i ] ) ;
    }
    **/
}

void solve ( ) {
    cin >> n ;
    for ( int i = 1 , x , y ; i < n ; ++ i ) {
        cin >> x >> y ;
        v[ x ].push_back ( y ) ;
        v[ y ].push_back ( x ) ;
    }
    for ( int i = 0 ; i < 2 * ( 1 << LOG ) ; ++ i ) {
        bit_rv[ i ] = -1 ;
    }
    for ( int i = 0 ; i <= LOG ; ++ i ) {
        bit_rv[ ( 1 << i ) ] = i ;
    }
    dfs ( 1 , -1 ) ;
    cout << complete[ 1 ] << "\n" ;
}

int main ( ) {
    ios_base :: sync_with_stdio ( false ) ;
    cin.tie ( NULL ) ;
    int t = 1 ; // cin >> t ;
    while ( t -- ) { solve ( ) ; }
    return 0 ;
}