#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
#include<random>
#include<functional>
using namespace std;
#define fi first
#define se second
#define mp make_pair
#define pb emplace_back
#define dbg(x) cerr<<"In the line "<<__LINE__<<" the "<<#x<<" = "<<x<<'\n'
#define dpi(x,y) cerr<<"In the line "<<__LINE__<<" the "<<#x<<" = "<<x<<" the "<<#y<<" = "<<y<<'\n';
typedef long long ll;
typedef vector<int> vi;
template<typename T>T &read(T &r){
	r=0;bool w=0;char ch=getchar();
	while(ch<'0'||ch>'9')w=ch=='-'?1:0,ch=getchar();
	while(ch>='0'&&ch<='9')r=r*10+ch-'0',ch=getchar();
	return r=w?-r:r;
}
inline int lowbit(int x){return x&(-x);}
const int N=300010;
const int B=100;//
int n,m;
int fa[N];
vi eg[N];
int dep[N],L[N],R[N],dft,mx;
int rk[N];
vi vec[N],tree[N];
vi pos[201][201],tr[201][201];
void modify(int o,int x,int v){
	int k=tree[o].size();--k;
	for(;x<=k;x+=lowbit(x))tree[o][x]+=v;
}
int query(int o,int x){
	int s=0;for(;x;x-=lowbit(x))s+=tree[o][x];
	return s;
}
void modify(int u,int v,int x,int w){
	int k=tr[u][v].size();--k;
	for(;x<=k;x+=lowbit(x))tr[u][v][x]+=w;
}
int query(int u,int v,int x){
	int s=0;for(;x;x-=lowbit(x))s+=tr[u][v][x];
	return s;
}
void dfs(int x){
	L[x]=++dft;mx=max(mx,dep[x]);rk[dft]=x;
	vec[dep[x]].pb(dft);
	for(int i=1;i<=B;i++)pos[i][dep[x]%i].pb(dft);
	for(auto v:eg[x]){
		dep[v]=dep[x]+1;
		dfs(v);
	}
	R[x]=dft;
}
void solve(){
	read(n);read(m);
	for(int i=1;i<=n;i++)vi().swap(eg[i]),vi().swap(vec[i]);
	for(int i=2;i<=n;i++)read(fa[i]),eg[fa[i]].pb(i);
	for(int i=1;i<=B;i++)
		for(int j=0;j<i;j++)
			vi().swap(pos[i][j]);
	mx=0;
	dep[1]=1;dfs(1);
	for(int i=1;i<=mx;i++){
		vi().swap(tree[i]);
		tree[i].resize(vec[i].size()+1);
	}
	for(int i=1;i<=B;i++)
		for(int j=0;j<i;j++){
			vi().swap(tr[i][j]);
			tr[i][j].resize(pos[i][j].size()+1);
		}
	auto LB=[&](const vi &ve,int &x){
		return lower_bound(ve.begin(),ve.end(),x)-ve.begin();
	};
	auto UB=[&](const vi &ve,int &x){
		return upper_bound(ve.begin(),ve.end(),x)-ve.begin();
	};
	for(int o=1;o<=m;o++){
		int op;read(op);
		if(o%1000==0)cerr<<o<<' ';
		if(op==1){
			int a,x,y,z;read(a);read(x);read(y);read(z);
			if(x>B){
				for(int i=dep[a]+y;i<=mx;i+=x){
						int l=LB(vec[i],L[a]);
						int r=UB(vec[i],R[a])-1;
						++l;++r;
						if(l<=r){
							modify(i,l,z);
							modify(i,r+1,-z);
						}
						else break;
					}
			}
			else{
				y=(y+dep[a])%x;
				int l=LB(pos[x][y],L[a]);
				int r=UB(pos[x][y],R[a])-1;
				++l;++r;
				if(l<=r){
//					dpi(l,r);
//					dpi(rk[pos[x][y][l-1]],rk[pos[x][y][r-1]]);
					modify(x,y,l,z);
					modify(x,y,r+1,-z);
				}
			}
		}
		else{
			int a;read(a);
			int p=LB(vec[dep[a]],L[a]);
			int s=query(dep[a],p+1);
			for(int i=1;i<=B;i++){
				s+=query(i,dep[a]%i,LB(pos[i][dep[a]%i],L[a])+1);
			}
			cout << s << '\n';
		}
	}
}
signed main(){
//	freopen("data.in","r",stdin); 
//	freopen("data.out","w",stdout); 
	int T;read(T);
	while(T--)solve();
	return 0;
}