#include <bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
int inv(int x){
	int res=1,y=mod-2;
	for(;y;y>>=1,x=1ll*x*x%mod) if(y&1) res=1ll*res*x%mod;
	return res;
}
void solve(){
	int n;
	cin>>n;
	vector<int> a(n),b(n),w(n),f(n),d(n),vis(n);
	for(int&x:a) cin>>x;
	for(int&x:b) cin>>x,x--;
	for(int&x:w) cin>>x;
	for(int i=0;i<n;i++){
		int x=b[i];
		if(a[x]+w[x]<=a[i]) f[i]=1;
		if(a[x]>a[i]) f[i]=2;
	}
	auto calc=[&](auto self,int x)->void{
		if(f[x]||d[x]) return;
		if(vis[x]) {f[x]=1;return;}
		vis[x]=1;
		self(self,b[x]);
		vis[x]=0;
		if(f[b[x]]==1) f[x]=1;
		else d[x]=d[b[x]]+1;
	};
	for(int i=0;i<n;i++){
		calc(calc,i);
		cout<<(a[i]+(f[i]!=1?1ll*w[i]*inv(d[i]+1):0))%mod<<" ";
	}
	cout<<"\n";
}
int main(){
	cin.tie(0)->sync_with_stdio(0);
	int t;
	cin>>t;
	while(t--) solve();
	return 0;
}