#include <iostream>
#include <algorithm>
#include <vector>

#define int long long

using namespace std;

const int mod=1e9+7;

int a[500001],b[500001],w[500001],c[500001],root[500001],depth[500001],vis[500001],pre_inv[500001],root_num;
vector<int> son[500001];

int quick_pow(int base,int index)
{
	int ans=1,tmp=base;
	
	while(index){
		if(index&1){
			ans*=tmp;	ans%=mod;
		}
		tmp*=tmp;	tmp%=mod;
		index>>=1;
	}
	
	return ans;
}

int get_inv(int x)
{
	return quick_pow(x,mod-2);
}

void dfs(int u,int d)
{
	depth[u]=d;
	
	for(auto v:son[u])
	{
		dfs(v,d+1);
	}
	
	return;
}

void solve()
{
	int n;
	
	cin>>n;
	root_num=0;
	for(int i=1;i<=n;i++){
		cin>>a[i];
		
		vector<int>().swap(son[i]);
		depth[i]=c[i]=root[i]=vis[i]=0;
	}
	for(int i=1;i<=n;i++){
		cin>>b[i];
	}
	for(int i=1;i<=n;i++){
		cin>>w[i];
	}
	
	for(int i=1;i<=n;i++)
	{
		for(int now=i;!vis[now];now=b[now])
		{
			vis[now]=1;
			
			if(a[now]<a[b[now]])
			{
				c[now]=a[now]+w[now]%mod;
				root[++root_num]=now;
				break;
			}
			else if(a[now]>=a[b[now]]+w[b[now]])
			{
				c[now]=a[now];
				root[++root_num]=now;
				break;
			}
			
			c[now]=a[now]+w[now]%mod;
			son[b[now]].push_back(now);
			
			if(b[now]==i)
			{
				now=i;
				while(1)
				{
					c[now]=a[now];
					son[now].pop_back();
					root[++root_num]=now;
					now=b[now];
					if(now==i) break;
				}
				break;
			}
		}
	}
	
	for(int i=1;i<=root_num;i++)
	{
		dfs(root[i],1);
	}
	
	for(int i=1;i<=n;i++)
	{
		cout<<(pre_inv[depth[i]]*c[i]+(1-pre_inv[depth[i]]+mod)*a[i])%mod<<' ';
	}
	cout<<'\n';
	
	return;
}

signed main()
{
	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
	
	pre_inv[1]=1;
	for(int i=2;i<=500000;i++){
		pre_inv[i]=pre_inv[i-1]*get_inv(i);	pre_inv[i]%=mod;
	}
	
	int t;
	cin>>t;
	while(t--){
		solve();
	}
	
	return 0;
}