#include <bits/stdc++.h>

 #define mid ((l+r)>>1)

 #define mod 201314

 using namespace std;

 int n,m,N,M,an,p,q,o;

 int size[],fa[],top[],pos[],son[],bro[],l[],r[],ans[];

 int sum[],flag[];

 struct quer

 {

     int r,x,bel;

 } que[];

 bool operator<(quer a,quer b){    return a.r<b.r;}

 void ins(int a,int b,int c)

 {

     que[++N].r=a;

     que[N].x=b;

     que[N].bel=c;

 }

 int build(int now)

 {

     size[now]=;

     for(int i=son[now];i;i=bro[i])

         size[now]+=build(i);

     return size[now];

 }

 void pou(int now,int to)

 {

     int ma=,id=;top[now]=to;pos[now]=++M;

     for(int i=son[now];i;i=bro[i])

         if(size[i]>ma) ma=size[i],id=i;

     if(id) pou(id,to);

     for(int i=son[now];i;i=bro[i])

         if(i!=id) pou(i,i);

 }

 void push(int now,int l,int r)

 {

     if(flag[now] && l!=r)

     {

         sum[now*]+=(mid-l+)*flag[now];

         sum[now*]%=mod;

         sum[now*+]+=(r-mid)*flag[now];

         sum[now*]%=mod;

         flag[now*]+=flag[now];

         flag[now*]%=mod;

         flag[now*+]+=flag[now];

         flag[now*+]%=mod;

         flag[now]=;

     }

 }

 void add(int now,int l,int r,int x,int y)

 {

     if(l==x && r==y)

     {

         sum[now]+=r-l+;flag[now]+=;

         sum[now]%=mod;

         return;

     }

     push(now,l,r);

     if(x<=mid) add(now*,l,mid,x,min(mid,y));

     if(y>mid) add(now*+,mid+,r,max(x,mid+),y);

     sum[now]=sum[now*]+sum[now*+];

 }

 int query(int now,int l,int r,int x,int y)

 {

     if(l==x && r==y)

         return sum[now];

     push(now,l,r);

     int ans=;

     if(x<=mid) ans+=query(now*,l,mid,x,min(mid,y));

     if(y>mid) ans+=query(now*+,mid+,r,max(x,mid+),y);

     return ans;

 }

 void link(int now)

 {

     for(;now;now=fa[top[now]])

         add(,,M,pos[top[now]],pos[now]);

 }

 int qu(int now)

 {

     for(an=;now;now=fa[top[now]])

         an=(an+query(,,M,pos[top[now]],pos[now]))%mod;

     return an;

 }

 int main()

 {

     scanf("%d%d",&n,&m);

     for(int i=;i<=n;i++)

         scanf("%d",&fa[i]),++fa[i],

         bro[i]=son[fa[i]],son[fa[i]]=i;

     build();pou(,);

     for(int i=;i<=m;i++)

     {

         scanf("%d%d%d",&p,&q,&o);

         l[i]=p;r[i]=++q;++o;

         if(p)ins(p,o,i);ins(q,o,i);

     }

     sort(que+,que+N+);

     int now=;

     for(int i=;i<=N;i++)

     {

         while(que[i].r>now) link(++now);

         int ret=qu(que[i].x);

         if(now==l[que[i].bel])

             ans[que[i].bel]-=ret;

         else

             ans[que[i].bel]+=ret;

     }

     for(int i=;i<=m;i++)

         printf("%d\n",(ans[i]+mod)%mod);

     return ;

 }