#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int MAXN = 2e5 + 5;
const ll linf = 0x3f3f3f3f3f3f3f3f;

inline void chk_max(ll &a,ll b){ if(a<b) a=b;}

struct Edge
{
	int next,to;
}e[MAXN<<1];
int head[MAXN],etot=0;
inline void add(int u,int v)
{
	e[++etot] = (Edge){ head[u],v};
	head[u] = etot;
}

int a[MAXN];
ll f[MAXN];

int mx, dep[MAXN], anc[MAXN];
vector<int> nodes[MAXN];
void dfs_dep(int u,int fa)
{
	mx = max(mx, dep[u]);
	anc[u] = fa;
	nodes[dep[u]].push_back(u);
	for(int i=head[u]; i; i=e[i].next)
	{
		int v = e[i].to;
		if(v == fa) continue;
		dep[v] = dep[u] + 1;
		dfs_dep(v,u);
	}
}

void solve(void)
{
	int n;
	scanf("%d",&n);
	for(int i=1; i<=n; ++i) head[i] = 0;
	etot = 0;
	for(int i=2; i<=n; ++i)
	{
		int v;
		scanf("%d",&v);
		add(i,v); add(v,i);
	}
	for(int i=2; i<=n; ++i) scanf("%d",&a[i]);
	
	mx = 0;
	for(int i=0; i<=n; ++i) nodes[i].clear();
	dfs_dep(1,0);
	
	ll ans = 0;
	for(int k=1; k<=mx; ++k)
	{
		vector<int> &vec = nodes[k];
		
		static int dsc[MAXN];
		int dcnt = 0;
		for(int i=0; i<(int)vec.size(); ++i)
			dsc[++dcnt] = a[vec[i]];
		sort(dsc+1,dsc+dcnt+1);
		dcnt = unique(dsc+1,dsc+dcnt+1) - dsc - 1;
		for(int i=0; i<(int)vec.size(); ++i)
			a[vec[i]] = lower_bound(dsc+1,dsc+dcnt+1,a[vec[i]]) - dsc;
		
		static ll pre[MAXN], suf[MAXN];
		for(int i=0; i<=dcnt+1; ++i) pre[i] = suf[i] = -linf;
		
		for(int i=0; i<(int)vec.size(); ++i)
		{
			int u = vec[i];
			chk_max(pre[a[u]], f[anc[u]] - dsc[a[u]]);
			chk_max(suf[a[u]], f[anc[u]] + dsc[a[u]]);
		}
		for(int i=1; i<=dcnt; ++i) chk_max(pre[i], pre[i-1]);
		for(int i=dcnt; i>=1; --i) chk_max(suf[i], suf[i+1]);
		
		for(int i=0; i<(int)vec.size(); ++i)
		{
			int u = vec[i];
			f[u] = f[anc[u]] + max(dsc[dcnt] - dsc[a[u]], dsc[a[u]] - dsc[1]);
			chk_max(f[u], pre[a[u]] + dsc[a[u]]);
			chk_max(f[u], suf[a[u]] - dsc[a[u]]);
			chk_max(ans, f[u]);
		}
	}
	printf("%lld\n",ans);
}

int main(void)
{
	int T;
	scanf("%d",&T);
	while(T--) solve();
	return 0;
}
