#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=2e5+5;
vector<int> G[maxn];
int n,d,a[maxn];
vector<int> D[maxn];//每个数的因子
vector< pair<int,int> > D2[maxn];
int g[maxn],prime[maxn],mu[maxn],tot; 
int siz[maxn];
inline void dfs1(int x,int fa){
	if(fa) G[x].erase(find(G[x].begin(),G[x].end(),fa));
	siz[x]=1;
	for(auto &y:G[x]){
		dfs1(y,x);
		siz[x]+=siz[y];
		if(siz[y]>siz[G[x][0]]) swap(y,G[x][0]);
	}
}
int f[maxn],czc[maxn];
int ans;
#define ll long long
const ll INF=4e10;
ll pw[maxn];
inline void add(int x,int s,int op){
	int v=a[x];
	if(op==1){
		for(auto l:D[v]) czc[l]+=g[v]*op;
		for(auto t:D2[v]){
			int fl=g[v]*mu[t.first]*op;
			ans+=fl*czc[t.second];
		}
	}
	else{
		for(auto t:D2[v]){
			int fl=g[v]*mu[t.first]*op;
			ans+=fl*czc[t.second];
		}
		for(auto l:D[v]) czc[l]+=g[v]*op;
	}
	for(auto y:G[x]){
		if(y==s) continue;
		add(y,s,op);
	}
}
inline void dfs2(int x,int op){ 
	for(int i=1;i<(int)G[x].size();i++)
		dfs2(G[x][i],0);
	if(G[x].size()) dfs2(G[x][0],1);
	add(x,G[x].size()?G[x][0]:0,1);
	f[x]=ans;
	if(!op){
		add(x,0,-1);
	}
}
signed main(){
	scanf("%lld%lld",&n,&d);
	for(int i=1;i<=n;i++){
		ll v=1;
		for(int j=1;j<=d+1;j++) v*=(ll)i;
		if(v>INF) break;
		pw[i]=v;
	}
	for(int i=1,x,y;i<n;i++){
		scanf("%lld%lld",&x,&y);
		G[x].push_back(y);
		G[y].push_back(x);
	}
	for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
	for(int i=1;i<=n;i++) for(int j=1;j*i<=n;j++) D[i*j].push_back(i);
	g[1]=mu[1]=1;
	for(int i=2;i<=n;i++){
		if(!prime[i]){
			prime[++tot]=i;
			mu[i]=g[i]=-1;
		}
		for(int j=1;j<=tot && i*prime[j]<=n;j++){
			if(i%prime[j]==0){
				mu[i*prime[j]]=0;
				g[i*prime[j]]=g[i]*g[prime[j]];//完全积性函数的优越性 
				prime[i*prime[j]]=1;
				break;
			}
			else{
				mu[i*prime[j]]=mu[i]*mu[prime[j]];
				g[i*prime[j]]=g[i]*g[prime[j]];
				prime[i*prime[j]]=1;
			}
		}
		for(int k=2;pw[k] && pw[k]<=i;k++){
			if(i%pw[k]==0){
				g[i]=0;
				break;
			}
		}
	}
	for(int i=1;i<=n;i++){
		for(auto y:D[i])
			if(pw[y] && pw[y]/__gcd(pw[y],(ll)i)<=n)
				D2[i].push_back({y,pw[y]/__gcd(pw[y],(ll)i)});
	}
	dfs1(1,0);
	dfs2(1,1);
	for(int i=1;i<=n;i++) printf("%lld\n",f[i]);
	return 0;
}