#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn = 1e5+5;
const int mod = 1e9+7;

struct node{
    int c=0;
    array<int,2> a[2];
    node(){a[0]=a[1]={0,0};}
    friend node operator+(node &x,node y){
        node res;
        res.c=x.c|y.c;
        for(int i=0;i<=1;i++){
            int ni=i|y.c;
            (res.a[ni][0]+=x.a[i][0])%=mod;
            (res.a[ni][1]+=x.a[i][1])%=mod;

            ni=x.c|i;
            (res.a[ni][0]+=y.a[i][0])%=mod;
            (res.a[ni][1]+=y.a[i][1])%=mod;
        }      
        return res;
    }
    void add(){
        if(c) (a[1][0]+=1)%=mod,(a[1][1]+=1)%=mod;
        else (a[1][0]+=1)%=mod,(a[0][1]+=1)%=mod;
    }
}dp[maxn],sum[maxn];

int n;
ll D;
vector<int> edge[maxn];
node rev(node cur){
    cur.c^=1;
    swap(cur.a[0],cur.a[1]);
    return cur;
}

void dfs(int u,int p){
    for(int v:edge[u]){
        if(v==p) continue;
        dfs(v,u);
        sum[u]=sum[u]+rev(dp[v]);
    }
    dp[u]=sum[u];dp[u].add();
}

node val[maxn];
void redfs(int u,int p){
    sum[u]=sum[u]+rev(val[u]);sum[u].add();

    node cur;cur=cur+rev(val[u]);
    for(int v:edge[u]){
        if(v==p) continue;
        val[v]=cur;cur=cur+rev(dp[v]);
    }
    cur=node();
    reverse(edge[u].begin(),edge[u].end());
    for(int v:edge[u]){
        if(v==p) continue;
        val[v]=val[v]+cur;
        cur=cur+rev(dp[v]);
        val[v].add();
        redfs(v,u);
    }
}

array<int,2> mul0(array<int,2> val,array<array<int,2>,2> ss){
    array<int,2> res={0,0};
    for(int i=0;i<=1;i++) for(int j=0;j<=1;j++) (res[j]+=1LL*val[i]*ss[i][j]%mod)%=mod;
    return res;
}
array<array<int,2>,2> mul1(array<array<int,2>,2> val,array<array<int,2>,2> ss){
    array<array<int,2>,2> res={0,0};
    for(int i=0;i<=1;i++) for(int j=0;j<=1;j++) for(int k=0;k<=1;k++) (res[i][j]+=1LL*val[i][k]*ss[k][j]%mod)%=mod;
    return res;
}

signed main(){
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);cout.tie(NULL);
    cin >> n >> D;
    for(int i=1;i<n;i++){
        int u,v;cin >> u >> v;
        edge[u].push_back(v);
        edge[v].push_back(u);
    }
    dfs(1,0);
    val[1].c=1;
    redfs(1,0);
    array<array<int,2>,2> ss;
    ss[0]=ss[1]={0,0};
    for(int u=1;u<=n;u++){
        for(int i=0;i<=1;i++) for(int j=0;j<=1;j++) (ss[j][i]+=sum[u].a[i][j])%=mod;
    }
    array<int,2> val={0,0};
    for(int u=1;u<=n;u++) val[sum[u].c]++;
    D--;
    while(D){
        if(D&1) val=mul0(val,ss);
        ss=mul1(ss,ss);D>>=1;
    }
    cout << (1LL*sum[1].a[1][0]*val[0]+1LL*sum[1].a[1][1]*val[1])%mod << '\n';
}