#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 5e5 + 5;
int a[N], b[N], w[N];
int p[N];
const int mod = 1e9 + 7;
int f[N];
int n; 
LL qmi(LL x, LL k)
{
    LL res = 1;
    while(k)
    {
        if(k & 1) res = res * x % mod;
        k >>= 1;
        x = x * x % mod;
    }
    return res % mod;
}
bool st[N];
int dfs(int u)
{
    st[u] = true;
    int j = b[u];
    if(a[u] < a[j]) p[u] = 1;
    else if(a[u] >= a[j] + w[j]) p[u] = 0;
    else
    {
        if(st[j]) return p[u] = 0;
        int t = dfs(j);
        if(!t)p[u] = 0;
        else p[u] = t + 1;
    }
    return p[u];
}
void solve()
{
    cin >> n;
    for(int i = 1; i <= n; i++) cin >> a[i], st[i] = false;
    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(!st[i]) dfs(i);
        // cout << p[i] << " ";
        int t = 0;
        if(p[i]) t = qmi(f[p[i]], mod - 2);
        cout << (a[i] + t * w[i] % mod) % mod << " ";
    }

    cout << "\n";
}

int main()
{
    f[0] = 1;
    for(int i = 1; i < N; i++) f[i] = f[i - 1] * i % mod;
    int t; cin >> t;
    while(t--)
        solve();
}