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

const int N=1e6+5,T=20;

int n;
int f[N][T+1];
int v[N];
vector<int> ed[N];
int dp[N];
int cf[N];
int ans;

void dfsi(int now){
    dp[now]=dp[f[now][0]]+1;
    for(auto nxt:ed[now]){
        dfsi(nxt);
    }
}

void binit(){
    for(int t=1;t<=T;t++){
        for(int i=1;i<=n;i++){
            f[i][t]=f[f[i][t-1]][t-1];
        }
    }
}

int lca(int a,int b){
    if(dp[a]>dp[b])swap(a,b);
    for(int t=T;t>=0;t--){
        if(dp[f[b][t]]>=dp[a])b=f[b][t];
    }
    if(a==b)return a;
    for(int t=T;t>=0;t--){
        if(f[a][t]!=f[b][t]){
            a=f[a][t];
            b=f[b][t];
        }
    }
    return a;
}

void dfss(int now){
    for(auto nxt:ed[now]){
        dfss(nxt);
        cf[now]+=cf[nxt];
        if(cf[nxt]>1&&v[nxt]!=0)ans++;
    }
}

void solve(){
    cin>>n;
    f[1][0]=1;
    for(int i=2;i<=n;i++)cin>>f[i][0];
    for(int i=1;i<n;i++)ed[i].clear();
    for(int i=2;i<=n;i++){
        cin>>v[i];
        ed[f[i][0]].push_back(i);
    }
    dfsi(1);
    binit();
    for(int i=2,j;i<=n;i++){
        cin>>j;
        int k=lca(i,j);
        cf[k]-=2;
        cf[i]++;cf[j]++;
    }
    ans=0;
    dfss(1);
    cout<<ans<<"\n";
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);cout.tie(0);
    int t;cin>>t;
    while(t--)solve();
    return 0;
}