#include<iostream>
#include<vector>
using namespace std;
const long mod=1e9+7;
int N,K;
vector<int>G[1010];
long dp[1010][1010];
void dfs1(int u,int p)
{
	for(int i=0;i<K;i++)dp[u][i]=1;
	for(int v:G[u])if(v!=p)
	{
		dfs1(v,u);
		int t=0;
		for(int i=0;i<K;i++)
		{
			t+=dp[v][i];
			if(t>=mod)t-=mod;
			(dp[u][i]*=t)%=mod;
		}
	}
}
long ans;
void dfs2(int u,int p)
{
	vector<int>cs;
	for(int v:G[u])if(v!=p)cs.push_back(v);
	vector<vector<long> >L(cs.size()+1,vector<long>(K,1));
	vector<vector<long> >R(cs.size()+1,vector<long>(K,1));
	long t=p==-1;
	for(int i=0;i<K;i++)
	{
		(ans+=dp[u][i]*t)%=mod;
		if(p!=-1)
		{
			t+=dp[p][i];
			if(t>=mod)t-=mod;
		}
		L[0][i]=t;
	}
	for(int i=0;i<cs.size();i++)
	{
		long t=0;
		for(int j=0;j<K;j++)
		{
			t+=dp[cs[i]][j];
			if(t>=mod)t-=mod;
			L[i+1][j]=L[i][j]*t%mod;
		}
	}
	for(int i=cs.size();i--;)
	{
		long t=0;
		for(int j=0;j<K;j++)
		{
			t+=dp[cs[i]][j];
			if(t>=mod)t-=mod;
			R[i][j]=R[i+1][j]*t%mod;
		}
	}
	for(int i=0;i<cs.size();i++)
	{
		for(int j=0;j<K;j++)
		{
			dp[u][j]=L[i][j]*R[i+1][j]%mod;
		}
		dfs2(cs[i],u);
	}
}
main()
{
	cin>>N>>K;
	for(int i=1;i<N;i++)
	{
		int a,b;cin>>a>>b;
		a--,b--;
		G[a].push_back(b);
		G[b].push_back(a);
	}
	dfs1(0,-1);
	dfs2(0,-1);
	cout<<ans<<endl;
}