#include<bits/stdc++.h>
#define int long long
#define mod 1000000007
using namespace std;
int n,v,d[1000006],p[1000006],a[1000006],b[1000006],dp[1000006][4][2],f[2][5][2];
vector<int>e[1000006];
void dfs(int x){
	memset(dp[x],0,sizeof(dp[x]));
	for(int i:e[x]) dfs(i);
	memset(f,0,sizeof(f));
	int t=0;
	f[0][0][0]=1; 
	for(int i:e[x]){
		memset(f[t^1],0,sizeof(f[t]));
		for(int j=0;j<2;j++){
			for(int k=0;k+j<2;k++){
				f[t^1][0][j+k]+=f[t][0][j]*dp[i][0][k];
				for(int l=1;l<4;l++) f[t^1][l][j+k]+=f[t][l][j]*dp[i][0][k]+f[t][0][j]*dp[i][l][k];
				f[t^1][4][j+k]+=f[t][2][j]*dp[i][3][k]+f[t][3][j]*dp[i][2][k]+f[t][4][j]*dp[i][0][k];
			}
		}
		t^=1;
		for(int j=0;j<2;j++) for(int k=0;k<4;k++) f[t][k][j]%=mod;
	}
	for(int i=0;i<2;i++){
		dp[x][0][i]-=f[t][0][i]*d[x]+f[t][1][i];
		dp[x][1][i]+=f[t][0][i]*a[x]%mod*b[x];
		dp[x][2][i]+=(f[t][0][i]+f[t][2][i])*a[x];
		dp[x][3][i]+=(f[t][0][i]+f[t][3][i])*b[x];
	}
	dp[x][0][1]-=f[t][0][0]+f[t][2][0]+f[t][3][0]+f[t][4][0];
	for(int i=0;i<2;i++) for(int j=0;j<4;j++) dp[x][j][i]%=mod;
}
void doing(){
	scanf("%lld %lld",&n,&v);
	for(int i=1;i<=n;i++) scanf("%lld",&d[i]),d[i]-=v;
	for(int i=2;i<=n;i++) scanf("%lld %lld %lld",&p[i],&a[i],&b[i]),a[i]-=v,b[i]-=v;
	for(int i=1;i<=n;i++) e[i].clear();
	for(int i=2;i<=n;i++) e[p[i]].push_back(i);
	dfs(1);
	int ans=(dp[1][0][1]*v%mod+dp[1][0][0]+mod*2)%mod;
	if(n&1) ans=(mod-ans)%mod;
	printf("%lld\n",ans);
}
signed main(){
	int T; 
	cin>>T;
	while(T--) doing();
}