#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define cmin(a,b) a=min(a,b)
#define cmax(a,b) a=max(a,b)
#define modd(a,b) a=a%b
#define pb push_back
#define db double
#define lowbit(x) x&(-x)
#define cerr(x) cout<<#x<<"="<<x<<endl
#define fi first
#define se second
const ll mod=1e9+7;
#define N 500010
ll n,a[N],b[N],w[N],ans[N];
ll fact[N],inv[N];
ll p[N],dist[N],flag[N],col[N];
bool vis[N];
ll ksm(ll a,ll b) {ll ans=1,bs=a;while(b) {if(b&1) ans=ans*bs%mod;bs=bs*bs%mod;b>>=1;}return ans;}
void init(int n=500000) {
    fact[0]=1;
    for(int i=1;i<=n;i++) fact[i]=fact[i-1]*i%mod;
    inv[n]=ksm(fact[n],mod-2);
    for(int i=n-1;~i;i--) inv[i]=inv[i+1]*(i+1)%mod;
}
void dfs(int u) {
    if(vis[u]) return;
    vis[u]=1;
    dfs(b[u]);
    col[u]=col[b[u]];
    if(flag[u]) {
        if(p[u]==0) dist[u]=0,col[u]=0;
        else dist[u]=1,col[u]=1;
    }
    else {
        dist[u]=dist[b[u]]+1;
        if(col[b[u]]) ans[u]=(ans[u]+inv[dist[u]]*w[u]%mod)%mod;
    }
}
void solve() {
    cin>>n;
    for(int i=1;i<=n;i++) cin>>a[i];
    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++) flag[i]=0,p[i]=0,ans[i]=a[i],vis[i]=0;
    for(int i=1;i<=n;i++) {
        if(i==b[i] || a[i]>=a[b[i]]+w[b[i]]) p[i]=0,ans[i]=a[i],flag[i]=1,col[i]=0;
        else if(a[i]<a[b[i]]) p[i]=1,ans[i]=a[i]+w[i],flag[i]=1,col[i]=1;
    }
    for(int i=1;i<=n;i++) {
        if(vis[i]) continue;
        dfs(i);
    }
    for(int i=1;i<=n;i++) cout<<ans[i]<<" ";
    cout<<endl;
}
int main() {
    ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    init();
    int T;
    cin>>T;
    while(T--) 
    solve();
    return 0;
}
/*
虽然说对于全排列而言a在b前面和b在a前面的概率是一样的,但是如果固定了a在b前面,
来考虑b和c的关系的话,那b在c前面和c在b前面的概率是不同的,因为a和b的相对位置
固定了的情况下.所以对于某种情况而言应该是一种1/dis(i)!阶乘的形式
*/