#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=5e5+10;
const int mod=1e9+7;
int num,q[N],ne[N],hd[N],d[N];
void addnum(int x,int y){
    q[++num]=y;
    ne[num]=hd[x];
    hd[x]=num;
    d[y]++;
}
int inf[N],fac[N];
int fp(int x,int y){
    if(y==0)return 1;
    if(y==1)return x;
    int sum=fp(x,y/2);
    if(y%2==0)return sum*sum%mod;
    else return sum*sum%mod*x%mod;
}
int n,vis[N],a[N],b[N],w[N],ans[N];
void dfs(int t,int len){
    vis[t]=1;
    if(ans[t]==0){
        ans[t]=(a[t]+w[t]*inf[len]%mod)%mod;
        // cout<<t<<' '<<len<<' '<<ans[t]<<endl;
    }
    for(int i=hd[t]; i; i=ne[i]){
        int v=q[i];
        dfs(v,len+1);
    }
}
signed main(){
    int t;
    cin>>t;
    fac[0]=1;
    for(int i=1; i<=N-5; i++)fac[i]=fac[i-1]*i%mod;
    inf[N-5]=fp(fac[N-5],mod-2);
    for(int i=N-6; i>=0; i--)inf[i]=inf[i+1]*(i+1)%mod;
    while(t--){
        cin>>n;
        num=0;
        for(int i=1; i<=n; i++)hd[i]=0,ans[i]=0,vis[i]=0,d[i]=0;
        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++){
            if(a[i]<a[b[i]])ans[i]=a[i]+w[i];
            else if(a[i]>=a[b[i]]+w[b[i]])ans[i]=a[i];
            else addnum(b[i],i);
        }
        for(int i=1; i<=n; i++){
            if(!vis[i]&&d[i]==0){
                dfs(i,1);
            }
        }
        for(int i=1; i<=n; i++)
            if(ans[i]==0)cout<<a[i]<<' ';
            else cout<<ans[i]<<' ';
        cout<<endl;
    }
}