#include<bits/stdc++.h>
#define int long long
using namespace std;
int TT,n,m,fa[300003],a[300003],sxt[300003][3],lft[300003],rgt[300003],dep[300003],idx;
int mxdep[300003],mxson[300003];
vector<int>lst[300003];
vector<int>son[300003];
struct SegT{
	int mk;
	int mxv;
}T[2000003];
void pushup(int now){T[now].mxv=max(T[now*2].mxv,T[now*2+1].mxv);return;}
void addmk(int now,int v){
	T[now].mxv+=v;
	T[now].mk+=v;
	return;
}
void pushdown(int now){
	addmk(now*2,T[now].mk);
	addmk(now*2+1,T[now].mk);
	T[now].mk=0;
	return;
}
int Query(int now,int l,int r,int ql,int qr){
	if(ql>qr)return 0;
	if(l>=ql&&r<=qr)return T[now].mxv;
	pushdown(now);
	if(ql>((l+r)>>1))return Query(now*2+1,((l+r)>>1)+1,r,ql,qr);
	if(qr<=((l+r)>>1))return Query(now*2,l,((l+r)>>1),ql,qr);
	return max(Query(now*2,l,((l+r)>>1),ql,qr),Query(now*2+1,((l+r)>>1)+1,r,ql,qr));
}
void modify(int now,int l,int r,int wz,int v){
	if(l==r){T[now].mxv=v;return;}
	pushdown(now);
	if(wz<=((l+r)>>1))modify(now*2,l,((l+r)>>1),wz,v);
	else modify(now*2+1,((l+r)>>1)+1,r,wz,v);
	pushup(now);
	return;
}
void m0dify(int now,int l,int r,int ql,int qr,int v){
	if(l>qr||r<ql)return;
	if(l>=ql&&r<=qr){addmk(now,v);return;}
	pushdown(now);
	m0dify(now*2,l,((l+r)>>1),ql,qr,v);
	m0dify(now*2+1,((l+r)>>1)+1,r,ql,qr,v);
	pushup(now);
	return;
}
void dfs(int now){
	lft[now]=rgt[now]=0;
	dep[now]=mxdep[now]=dep[fa[now]]+1;
	mxson[now]=0;
	for(auto i:son[now]){
		dfs(i);
		if(mxdep[i]>mxdep[mxson[now]])mxson[now]=i;
		mxdep[now]=max(mxdep[now],mxdep[i]+1);
	}
	return;
}
void dfs2(int now){
	lft[now]=rgt[now]=++idx;
	modify(1,1,n,idx,a[now]);
	if(mxson[now]){
		dfs2(mxson[now]);
		rgt[now]=rgt[mxson[now]];
		m0dify(1,1,n,lft[now],lft[now],max(Query(1,1,n,lft[now]+1,rgt[now]),0ll));
	}
	for(auto i:son[now]){
		if(i==mxson[now])continue;
		dfs2(i);
		m0dify(1,1,n,lft[now],lft[now],max(Query(1,1,n,lft[i],rgt[i]),0ll));
		for(int j=1;j<=rgt[i]-lft[i]+1;j++){
			modify(1,1,n,lft[now]+j,max(0ll,Query(1,1,n,lft[i]+j-1,rgt[i]))+max(0ll,Query(1,1,n,lft[now]+j,rgt[now])));
		}
	}
	for(auto i:lst[now])m0dify(1,1,n,lft[now],min(rgt[now],lft[now]+sxt[i][1]),-sxt[i][2]);
	return;
}
void sol(){
	cin>>n>>m;
	for(int i=1;i<=n;i++){
		son[i].clear();
		son[i].shrink_to_fit();
		lst[i].clear();
		lst[i].shrink_to_fit();
	}
	for(int i=1;i<=4*n;i++)T[i].mk=T[i].mxv=0;
	for(int i=2;i<=n;i++)cin>>fa[i];
	for(int i=2;i<=n;i++)son[fa[i]].emplace_back(i);
	for(int i=1;i<=n;i++)cin>>a[i];
	for(int i=1;i<=m;i++)cin>>sxt[i][0]>>sxt[i][1]>>sxt[i][2];
	for(int i=1;i<=m;i++)lst[sxt[i][0]].emplace_back(i);
	dfs(1);
	idx=0;
	dfs2(1);
	cout<<max(Query(1,1,n,lft[1],rgt[1]),0ll)<<'\n';
	return;
}
signed main(){
	ios::sync_with_stdio(false);
	cin>>TT;
	while(TT--)sol();
	return 0;
}