#include <bits/stdc++.h>
using namespace std;
#define For(Ti,Ta,Tb) for(auto Ti=(Ta);Ti<=(Tb);++Ti)
#define Dec(Ti,Ta,Tb) for(auto Ti=(Ta);Ti>=(Tb);--Ti)
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define range(Tx) begin(Tx),end(Tx)
using ll=long long;
const int N=4.5e5+5;
int n,q,fa[N],n_,poi[N];
ll a[N];
vector<int> e[N];
struct Rect{
	int xl,xr,yl,yr;
	bool include(const Rect &x) const{
		return xl<=x.xl&&x.xr<=xr&&yl<=x.yl&&x.yr<=yr;
	}
	bool out_of(const Rect &x) const{
		return xl>x.xr||xr<x.xl||yl>x.yr||yr<x.yl;
	}
};
struct Point{
	int x,y,i;
	ll v;
}ins[N];
struct KDT{
	struct Node{
		Rect p;
		int cnt;
		ll add,sum;
	}t[N*4];
	void apply(int p,ll x){
		t[p].add+=x;
		t[p].sum+=t[p].cnt*x;
	}
	void push(int p){
		if(t[p].add){
			apply(p*2,t[p].add);
			apply(p*2+1,t[p].add);
			t[p].add=0;
		}
	}
	void pull(int p){
		t[p].sum=t[p*2].sum+t[p*2+1].sum;
	}
	void build(int p,int l,int r,bool axis){
		t[p].add=0;
		if(l==r){
			t[p].p.xl=t[p].p.xr=ins[l].x;
			t[p].p.yl=t[p].p.yr=ins[l].y;
			t[p].cnt=1;
			t[p].sum=ins[l].v;
			poi[ins[l].i]=p;
			return;
		}
		int mid=(l+r)/2;
		nth_element(ins+l,ins+mid,ins+r+1,[axis](const Point& p1,const Point& p2){
			return axis?p1.y<p2.y:p1.x<p2.x;
		});
		build(p*2,l,mid,!axis);
		build(p*2+1,mid+1,r,!axis);
		t[p].p.xl=min(t[p*2].p.xl,t[p*2+1].p.xl);
		t[p].p.xr=max(t[p*2].p.xr,t[p*2+1].p.xr);
		t[p].p.yl=min(t[p*2].p.yl,t[p*2+1].p.yl);
		t[p].p.yr=max(t[p*2].p.yr,t[p*2+1].p.yr);
		t[p].cnt=t[p*2].cnt+t[p*2+1].cnt;
		pull(p);
	}
	void range_add(int p,const Rect& x,int v){
		if(x.out_of(t[p].p)){
			return;
		}
		if(x.include(t[p].p)){
			apply(p,v);
			return;
		}
		push(p);
		range_add(p*2,x,v);
		range_add(p*2+1,x,v);
		pull(p);
	}
	ll query(int p,const Rect& x){
		if(x.out_of(t[p].p)){
			return 0;
		}
		if(x.include(t[p].p)){
			return t[p].sum;
		}
		push(p);
		return query(p*2,x)+query(p*2+1,x);
	}
	void iterate(int p,int l,int r){
		if(l>=r) return;
		push(p);
		int mid=(l+r)/2;
		iterate(p*2,l,mid);
		iterate(p*2+1,mid+1,r);
	}
}kdt;
int bel[N],dep[N],dfn[N],dfx,siz[N];
void dfs(int u){
	dep[u]=dep[fa[u]]+1;
	dfn[u]=++dfx;
	siz[u]=1;
	for(int v:e[u]){
		dfs(v);
		siz[u]+=siz[v];
	}
}
void rebuild(){
	kdt.iterate(1,1,n);
	For(i,1,n) a[i]=kdt.t[poi[i]].sum;
	dfx=0;
	dfs(1);
	n=n_;
	For(i,1,n) ins[i]={dfn[i],dep[i],i,a[i]};
	kdt.build(1,1,n,0);
}
bool subtree(int u,int v){
	return dfn[u]<=dfn[v]&&dfn[v]<dfn[u]+siz[u];
}
void dfs_add(int u,int d,int x){
	if(dep[u]>d) return;
	a[u]+=x;
	for(int v:e[u]){
		dfs_add(v,d,x);
	}
}
ll dfs_query(int u,int d){
	if(dep[u]>d) return 0;
	ll res=a[u];
	for(int v:e[u]){
		res+=dfs_query(v,d);
	}
	return res;
}
int main(){
	cin.tie(nullptr)->sync_with_stdio(false);
	#ifdef zyz
		assert(freopen("D.in","r",stdin));
		assert(freopen("D.out","w",stdout));
	#endif
	cin>>n>>q;
	For(i,1,n) cin>>a[i];
	For(i,2,n){
		cin>>fa[i];
		e[fa[i]].push_back(i);
	}
	n_=n;
	n=0;
	rebuild();
	For(_,1,q){
		int tp;
		cin>>tp;
		if(tp==1){
			int u,d,x;
			cin>>u>>d>>x;
			if(u<=n){
				kdt.range_add(1,{dfn[u],dfn[u]+siz[u]-1,dep[u],dep[u]+d},x);
				For(i,n+1,n_){
					if(dep[i]<=dep[u]+d&&subtree(u,bel[i])){
						a[i]+=x;
					}
				}
			}else{
				dfs_add(u,dep[u]+d,x);
			}
		}else if(tp==2){
			int u,d;
			ll ans;
			cin>>u>>d;
			if(u<=n){
				ans=kdt.query(1,{dfn[u],dfn[u]+siz[u]-1,dep[u],dep[u]+d});
				For(i,n+1,n_){
					if(dep[i]<=dep[u]+d&&subtree(u,bel[i])){
						ans+=a[i];
					}
				}
			}else{
				ans=dfs_query(u,dep[u]+d);
			}
			cout<<ans<<'\n';
		}else{
			++n_;
			cin>>fa[n_]>>a[n_];
			e[fa[n_]].push_back(n_);
			dep[n_]=dep[fa[n_]]+1;
			bel[n_]=(fa[n_]<=n?fa[n_]:bel[fa[n_]]);
			if((n_-n)%5000==0){
				rebuild();
			}
		}
	}
	return 0;
}