#include<bits/stdc++.h>
#define For(i, a, b) for(int i = (a); i <= (b); i++)
#define Rof(i, a, b) for(int i = (a); i >= (b); i--)
using namespace std;
using pii = pair<int, int>;
const int N = 1e6 + 5;
struct node{
	int x, y, z;
	int op, w, id;
}p[N << 3];
int n, q, tot, fa[N], dep[N], siz[N], son[N], ans[N];
int dft, L[N], R[N], dfu[N], tax[N];
vector<int> e[N];
void dfs(int now){
	siz[now] = 1; dep[now] = dep[fa[now]] + 1;
	L[now] = ++dft; dfu[dft] = now;
	for(int to : e[now]){
		dfs(to); siz[now] += siz[to];
		if(siz[to] > siz[son[now]]) son[now] = to;
	}
	R[now] = dft;
}
int lst[N]; set<int> S[N];
void dot(int now, bool keep){
	for(int to : e[now]){
		if(to == son[now]) continue;
		dot(to, 0);
	}
	if(son[now]) dot(son[now], 1);
	S[dep[now]].insert(now);
	for(int to : e[now]){
		if(to == son[now]) continue;
		For(i, L[to], R[to]){
			int x = dfu[i], d = dep[x];
			auto it = S[d].upper_bound(x);
			if(it != S[d].end()){
				int y = *it, tim = dep[x] - dep[now];
				if(x > lst[y]){
					p[++tot] = {y, lst[y], tim, 0, -1, 0};
					p[++tot] = {y, lst[y] = x, tim, 0, 1, 0};
				}
			}
			if(it != S[d].begin()){
				int y = *(--it), tim = dep[x] - dep[now];
				if(y > lst[x]){
					p[++tot] = {x, lst[x], tim, 0, -1, 0};
					p[++tot] = {x, lst[x] = y, tim, 0, 1, 0};
				}
			}
		}
		For(i, L[to], R[to]) S[dep[dfu[i]]].insert(dfu[i]);
	}
	if(!keep){
		For(i, L[now], R[now]) S[dep[dfu[i]]].erase(dfu[i]);
	}
}
struct BIT{
	int t[N];
	int lowbit(int x){return x & (-x);}
	void add(int x, int y){x++; for(; x <= n + 1; x += lowbit(x)) t[x] += y;}
	int query(int x){x++; int res = 0; for(; x; x -= lowbit(x)) res += t[x]; return res;}
}bit;
bool cmp1(const node &a, const node &b){return a.x != b.x ? a.x < b.x : (a.y != b.y ? a.y < b.y : (a.z != b.z ? a.z < b.z : a.op < b.op));}
bool cmp2(const node &a, const node &b){return a.y != b.y ? a.y < b.y : (a.z != b.z ? a.z < b.z : a.op < b.op);}
void cdq(int l, int r){
	if(l == r) return;
	int mid = (l + r) >> 1;
	cdq(l, mid); cdq(mid + 1, r);
	int j = l;
	For(i, mid + 1, r){
		while(j <= mid && p[j].y <= p[i].y){
			if(!p[j].op) bit.add(p[j].z, p[j].w);
			j++;
		}
		if(!p[i].op) continue;
		ans[p[i].id] += p[i].w * bit.query(p[i].z);
	}
	while(j > l) if(!p[--j].op) bit.add(p[j].z, -p[j].w);
	inplace_merge(p + l, p + mid + 1, p + r + 1, cmp2);
}
void Solve(){
	cin >> n >> q; int rt = -1;
	For(i, 1, n){
		cin >> fa[i]; if(!fa[i]) rt = i;
		if(fa[i]) e[fa[i]].push_back(i);
	}
	dfs(rt); dot(rt, 1);
	For(i, 1, q){
		int l, r, x; cin >> l >> r >> x;
		p[++tot] = {r, l - 1, x, 1, 1, i};
		p[++tot] = {l - 1, l - 1, x, 1, -1, i};
	}
	For(i, 1, n) p[++tot] = {i, 0, 0, 0, 1, 0};
	For(i, 1, n){
		p[++tot] = {i, lst[i], dep[i], 0, -1, 0};
		p[++tot] = {i, n, dep[i], 0, 1, 0};
	}
	sort(p + 1, p + tot + 1, cmp1); cdq(1, tot);
	For(i, 1, q) cout << ans[i] << '\n';
}
signed main(){
	cin.tie(0)->sync_with_stdio(0);
	int T = 1; //cin >> T;
	while(T--) Solve();
	return 0;
}