#include<bits/stdc++.h>
typedef long long ll;
#define pb push_back
#define fr first
#define sc second
#define endl '\n'
using namespace std;
#define mid ((left+right)>>1)
const ll mod=1e9+7;
ll sum(ll x,ll y){
return (x+y)%mod;
}
ll sub(ll x,ll y){
return sum(x,mod-(y%mod));
}
ll mul(ll x,ll y){
x%=mod;y%=mod;
return x*y%mod;
}
ll fpow(ll x,ll y){
ll res=1;
while(y>0){
if(y&1){
res=mul(res,x);
}
x=mul(x,x);
y>>=1;
}
return res;
}
int n;
ll d;
vector<int>komsu[100023];
int dp[100023],cnt[100023];
int dp2[100023],cnt2[100023];
int dp3[100023],cnt3[100023];
int par[100023];
ll ans=0;
void dfs2(int pos){
dp[pos]=dp2[pos];
cnt[pos]=cnt2[pos];
if(par[pos]){
dp3[pos]=dp[par[pos]];
cnt3[pos]=cnt[par[pos]];
if(dp2[pos]==0){
dp3[pos]--;
if(dp3[pos]==0){
cnt3[pos]=0;
for(int x:komsu[par[pos]]){
if(x==pos)continue;
if(x==par[par[pos]])continue;
cnt3[pos]+=cnt2[x];
}
if(par[par[pos]]){
cnt3[pos]+=cnt3[par[pos]];
}
}
else if(dp3[pos]==1){
for(int x:komsu[par[pos]]){
if(x==pos)continue;
if(x==par[par[pos]])continue;
if(dp2[x]==0){
cnt3[pos]+=cnt2[x];
}
}
if(par[par[pos]]&&dp3[par[pos]]==0){
cnt3[pos]+=cnt3[par[pos]];
}
}
}
else if(dp3[pos]==0){
cnt3[pos]-=cnt2[pos];
}
dp[pos]=!dp3[pos];
cnt[pos]=0;
for(int x:komsu[pos]){
if(x==par[pos])continue;
dp[pos]+=!dp2[x];
}
if(dp[pos]==0){
cnt[pos]++;
}
for(int x:komsu[pos]){
if(x==par[pos]){
if(dp[pos]==1){
if(dp3[pos]==0){
cnt[pos]+=cnt3[pos];
}
}
else if(dp[pos]==0){
cnt[pos]+=cnt3[pos];
}
continue;
}
if(dp[pos]==1){
if(dp2[x]==0){
cnt[pos]+=cnt2[x];
}
}
else if(dp[pos]==0){
cnt[pos]+=cnt2[x];
}
}
}
for(int x:komsu[pos]){
if(x==par[pos])continue;
dfs2(x);
}
}
void dfs(int pos){
for(int x:komsu[pos]){
if(x==par[pos])continue;
par[x]=pos;
dfs(x);
dp2[pos]+=!dp2[x];
}
if(dp2[pos]>1)return;
if(dp2[pos]==0)cnt2[pos]++;
for(int x:komsu[pos]){
if(x==par[pos])continue;
if(dp2[pos]==1){
if(dp2[x]==0){
cnt2[pos]+=cnt2[x];
}
}
else{
cnt2[pos]+=cnt2[x];
}
}
}
const bool deb=false;
int main(){
ios_base::sync_with_stdio(23^23);cin.tie(NULL);
cin>>n>>d;
for(int i=1;i<n;i++){
int x,y;cin>>x>>y;
komsu[x].pb(y);
komsu[y].pb(x);
}
dfs(1);
dfs2(1);
int ez=0;
ans=0;
ll kazan[2]={0,0};
ll devam[2]={0,0};
ll iht[2]={0,0};
if(dp[1]){
ans=mul(n-cnt[1],fpow(n,2*d-1));
iht[1]=cnt[1];
}
else{
iht[0]=cnt[1];
}
for(int i=1;i<=n;i++){
if(deb)cout<<dp[i]<<"-"<<cnt[i]<<endl;
if(dp[i]){
ez++;
kazan[1]=sum(kazan[1],n-cnt[i]);
devam[0]=sum(devam[0],cnt[i]);
}
else{
kazan[0]=sum(kazan[0],n-cnt[i]);
devam[1]=sum(devam[1],cnt[i]);
}
}
if(deb)cout<<ans<<endl;
if(deb)cout<<"devam: "<<devam[0]<<" "<<devam[1]<<endl;
if(deb)cout<<"kazan: "<<kazan[0]<<" "<<kazan[1]<<endl;
for(int i=1;i<d;i++){
if(deb)cout<<"iht: "<<iht[0]<<" "<<iht[1]<<endl;
ans=sum(ans,mul(fpow(n,(d-i)*2-1),
sum(mul(kazan[1],iht[1]),mul(kazan[0],iht[0]))));
ll iht2[2]={0,0};
iht2[0]=sum(mul(iht[0],devam[0]),mul(iht[1],devam[1]));
iht2[1]=sum(mul(iht[0],devam[1]),mul(iht[1],devam[0]));
iht[0]=iht2[0];
iht[1]=iht2[1];
if(deb)cout<<ans<<endl;
}
if(deb)cout<<"iht: "<<iht[0]<<" "<<iht[1]<<endl;
ans=sum(ans,sum(mul(iht[0],n-ez),mul(iht[1],ez)));
cout<<ans<<endl;
}
