#include<bits/stdc++.h>
using namespace std;

using ll = long long;

const int mod = 1e9 + 7;
const int N = 5e5 + 5;

ll fac[N], inv[N];
bool vis[N];
int d[N];
int par[N], a[N], w[N], b[N];

ll binpow(ll a, ll b, ll mod){
    ll res = 1;
    while(b){
        if(b & 1)
            res = res * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return res;
}

ll mul(ll a, ll b){
    a %= mod;
    b %= mod;
    return a * b % mod;
}
ll sub(ll a, ll b){
    a -= b;
    while(a < 0) a += mod;
    return a;
}
ll add(ll a, ll b){
    a += b;
    return a % mod;
}
void init(int n){
    for(int i = 0; i < n;   i++){
        vis[i] = false;
        d[i] = 0;
        par[i] = -1;
    }
}

int solve(int i){
    if(par[i] == -1)return 0;
    if(a[i] < a[b[i]])return 1;
    if(vis[i])return d[i];
    vis[i] = true;
    d[i] = 1 + solve(par[i]);

    return d[i];
}

void solve(){
    int n;
    cin >> n;
    init(n);
    
    vector<int> ans(n, -1);
    for(int i = 0; i < n; i++){
        cin >> a[i];
    }
    for(int i = 0; i < n; i++){
        cin >> b[i]; b[i]--;
        par[i] = b[i];
    }
    for(int i = 0; i < n; i++){
        cin >> w[i];
    }

    for(int i = 0; i < n; i++){
        if(a[i] >= a[b[i]] + w[b[i]] || i == b[i]){
            ans[i] = a[i];
            continue;
        }
        if(a[i] < a[b[i]]){
            ans[i] = (a[i] + w[i]) % mod;
            continue;
        }
    }
    for(int i = 0; i < n; i++){
        solve(i);
    }
    for(int i = 0; i < n; i++){
        if(ans[i] == -1){
            ans[i] = add(mul(inv[d[i] + 1], a[i] + w[i]), mul(sub(1, inv[d[i] + 1]), a[i]));
        }
    }
    for(int i = 0; i < n; i++){
        cout << ans[i] << " \n"[i == n - 1];
    }
}
int main() {
    int t;
    cin >> t;
    fac[0] = 1;
    for(int i = 1; i < N; i++){
        fac[i] = (fac[i - 1] * i) % mod;
        inv[i] = binpow(fac[i], mod - 2, mod);
    }
    while(t--){
        solve();
    }
}