#include <bits/stdc++.h>
using namespace std;
#define int long long
const int N=5e5+10;
const int mod=1e9+7;
#define ll long long
struct node
{
    int to;
};
vector<node> h[N];
int dis[N];
int n;
int a[N],b[N],w[N];
int d[N];
ll fact[N],infact[N];
int ans=0;
int kk[N];
ll qmi(ll a,ll b,ll p)
{
    ll cnt=1;
    while(b)
    {
        if(b&1) cnt=(cnt*a)%mod;
        b>>=1;
        a=(a*a)%mod;
    }
    return cnt;
}
void add(int a,int b)
{
    h[a].push_back({b});
    d[b]++;
}
void bfs(int pos)
{
    queue<int> q;
    q.push(pos);
    dis[pos]=1;
    while(q.size())
    {
        auto t=q.front();
        q.pop();
        for(int i=0;i<h[t].size();i++)
        {
            int j=h[t][i].to;
            if(dis[j]!=0) continue;
            dis[j]=dis[t]+1;
            d[j]--;
            kk[j]=(infact[dis[j]]*(a[j]+w[j])%mod+(1-infact[dis[j]]+mod)%mod*a[j]%mod)%mod;
          //  ans=(ans+infact[dis[j]]*w[j]%mod)%mod;
            q.push(j);
        }
    }

}
void solve()
{
    for(int i=1;i<=n;i++)
    {
        h[i].clear();
        d[i]=0;
        dis[i]=0;
    }
    ans=0;
    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++)
    {
        if(a[b[i]]<=a[i]&&a[b[i]]+w[b[i]]>a[i])
        {
           // cout<<b[i]<<"aa"<<i<<endl;
            add(b[i],i);
        }
        else if(a[b[i]]>a[i])
        {
            kk[i]=(a[i]+w[i])%mod;
        }
        else kk[i]=a[i];
    }
    for(int i=1;i<=n;i++)
    {
        if(dis[i]==0&&d[i]==0) 
        {
            bfs(i);
            //cout<<i<<"aa"<<endl;
        }
    }
    for(int i=1;i<=n;i++)
    {
        if(d[i]!=0)
        {
            kk[i]=a[i];
        }
    }
    for(int i=1;i<=n;i++)
    {
        cout<<kk[i]<<' ';
    }
    cout<<endl;

}
signed main()
{
    //cout<<qmi(2,mod-2,mod)<<endl;
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int T;
    cin>>T;
    fact[0]=infact[0]=1;
    for(ll i=1;i<N;i++)  //此处不能取到N,否则数组越界
    {
        fact[i]=(fact[i-1]*i)%mod;
        infact[i]=(infact[i-1]*qmi(i,mod-2,mod))%mod;
    }
    while(T--)
    {
        solve();
    }
}