#include<bits/stdc++.h>
#define ll long long
const int p=1e9+7;
using namespace std;

const int N=5e5+5;
ll inv[N],ans[N];
int a[N],b[N],w[N],d[N];
bool vis[N];

int ksm(ll a,int b) {
	ll ret=1;
	while(b) {
		if(b&1) ret=ret*a%p;
		a=a*a%p,b>>=1;
	}
	return ret;
}
int dfs(int u) {
	if(vis[u]) {
		return d[u];
	}
	vis[u]=1;
	int dep=dfs(b[u])+1; 
	if(a[u]<a[b[u]]) ans[u]=1,dep=0;
		else if(a[u]>=a[b[u]]+w[b[u]]) ans[u]=0,dep=0;
			else {
				if(ans[b[u]]) {
					ans[u]=inv[dep+1];
				} else ans[u]=0;
			}
	return d[u]=(ans[u]==0?0:dep);
}
int main() {
	ll jc=1;
	for(int i=2;i<=5e5;++i) jc=jc*i%p;
	inv[(int)5e5]=ksm(jc,p-2);
	for(int i=(int)5e5-1;i;--i) {
		inv[i]=inv[i+1]*(i+1)%p;
	}
	int T; scanf("%d",&T);
	while(T--) {
		int n; scanf("%d",&n);
		for(int i=1;i<=n;++i) {
			scanf("%d",&a[i]);
		}
		for(int i=1;i<=n;++i) {
			scanf("%d",&b[i]);
		}
		for(int i=1;i<=n;++i) {
			scanf("%d",&w[i]);
		}
		for(int i=1;i<=n;++i) {
			vis[i]=0; ans[i]=0,d[i]=0;
		}
		for(int i=1;i<=n;++i) {
			if(!vis[i]) {
				dfs(i);
			}
			printf("%lld ",(a[i]+w[i]*ans[i])%p);
		}
		puts("");
	}
	return 0;
}