#pragma GCC optimize("Ofast","inline")
#include<bits/stdc++.h>
using namespace std;
const int N=2.5e5+500,mod=998244353;

void ad(int&x,int y)
{
    x+=y;
x-=x>=mod?mod:0;
}

int n,q,S=555,a[N],fa[N],b[N],ro[N],f[N][20],ex[N],dep[N],son[N],mx[N],hb[N],las,la[N],T,val[N],vis[N];
vector<int>up[N],down[N];

struct block
{
    int*w,*fir;
}B[N/500+10];

int find(int x,int k)
{
    if(k<=0)return x;
    int p=f[x][hb[k]],t=ro[p],d=dep[p]-dep[t];
    k-=1<<hb[k];
    if(d>=k)return down[t][d-k];
    return up[t][k-d];
}

void insert(int x,int&ans,int d,int pre)
{
int &w=val[x];
    if(vis[x]!=T)vis[x]=T,w=pre;
    ad(ans,d*(d+2ll*w)%mod);
    ad(w,d);
}

main()
{
ios::sync_with_stdio(0);
    cin.tie(0) , cout.tie(0);
    cin>>n>>q;
    for(int i=1;i<=n;i++)cin>>a[i];dep[1]=1;
    for(int i=2;i<=n;i++)cin>>fa[i],dep[i]=dep[fa[i]]+1,f[i][0]=fa[i];
    for(int i=1,pos=1;i<=n;i++)
    {
        while(pos<=n&&fa[pos]<i)++pos;
        la[i]=pos-1;
    }
    for(int i=2;i<=n;i++)for(int j=1;(1<<j)<dep[i];j++)f[i][j]=f[f[i][j-1]][j-1];
    for(int i=0;(1<<i)<dep[n];i++)for(int j=(1<<i);j<(1<<i+1)&&j<dep[n];j++)hb[j]=i;
    for(int i=n;i>=1;i--)
    {
        mx[i]=max(mx[i],dep[i]);mx[fa[i]]=max(mx[fa[i]],mx[i]);
        if(mx[son[fa[i]]]<mx[i])son[fa[i]]=i;
    }
    son[0]=0;
    for(int i=1;i<=n;i++)ro[i]=son[fa[i]]==i?ro[fa[i]]:i;
    for(int i=1;i<=n;i++)down[ro[i]].push_back(i);
    for(int i=1;i<=n;i++)if(i==ro[i])
    {
        int x=i,re=mx[i]-dep[i]+1;
        while(re--)up[i].push_back(x),x=fa[x];
    }
    for(int i=1;i<=n;i++)b[i]=(i-1)/S+1;
    for(int i=1;i<=b[n];i++)
    {
        int l=(i-1)*S+1,r=min(n,i*S);
        B[i].w=new int[r+10];B[i].fir=new int[r+10];
        for(int j=0;j<r+10;j++)B[i].fir[j]=B[i].w[j]=0;
        for(int j=1;j<l;j++)B[i].w[j]=B[i-1].w[j],ex[j]=0;
        for(int j=l;j<=r;j++)ex[j]=a[j];
        for(int j=r;j>=1;j--)ad(ex[fa[j]],ex[j]),ad(B[i].w[j],ex[j]);
int now=0;
        for(int j=1;j<=r;j++)ad(now,1ll*B[i].w[j]*B[i].w[j]%mod),B[i].fir[j]=now;
    }
    while(q--)
    {
        int l,r;cin>>l>>r;l^=las;r^=las;T++;
        int rr=min(r,la[l]),k=b[r]-1,w=0;
        if(l<=k*S)w=(B[k].fir[min(rr,k*S)]-B[k].fir[l-1]+mod)%mod;
        for(int i=max(l,k*S+1);i<=r;i++)
        {
            int d=find(i,dep[i]-dep[rr]);
            if(d>rr)d=fa[d];
            insert(d,w,a[i],d<=k*S?B[k].w[d]:0);
        }
        cout<<(las=(w+mod)%mod)<<'\n';
    }
}