#include<bits/stdc++.h>
using namespace std;
const int N=100005;
const long long mod=998244353;
long long f[N],g[N],dis[N],ans;
vector<int>e[N],r[N];
int fa[N][21],dep[N],n,q,K;
void dfs(const int k)
{
	r[k].resize(e[k].size()+2),f[k]=1;
	for(int i:e[k])
		dfs(i),f[k]=f[k]*(f[i]+1)%mod;
	r[k][e[k].size()]=1;
	for(int i=e[k].size()-1;i>=0;--i)
		r[k][i]=r[k][i+1]*(f[e[k][i]]+1)%mod;
}
void dfs2(const int k)
{
	dis[k]=(dis[fa[k][0]]+f[k])%mod,dep[k]=dep[fa[k][0]]+1;
	long long s=1;
	for(int i=0;i<e[k].size();++i)
		g[e[k][i]]=(g[k]+1)*s%mod*r[k][i+1]%mod,s=s*(f[e[k][i]]+1)%mod,dfs2(e[k][i]);
}
int LCA(int u,int v)
{
	if(dep[u]<dep[v])
		swap(u,v);
	for(int i=K;i>=0;--i)
		if(dep[fa[u][i]]>=dep[v])
			u=fa[u][i];
	if(u==v)
		return u;
	for(int i=K;i>=0;--i)
		if(fa[u][i]!=fa[v][i])
			u=fa[u][i],v=fa[v][i];
	return fa[u][0];
}
int main()
{
	ios::sync_with_stdio(0),cin.tie(0);
	cin>>n>>q;
	int u,v,w;
	for(int i=2;i<=n;++i)
		cin>>u,e[u].push_back(i),fa[i][0]=u;
	fa[1][0]=dep[1]=1;
	dfs(1);
	dfs2(1);
	K=__lg(n);
	for(int i=1;i<=K;++i)
		for(int j=1;j<=n;++j)
			fa[j][i]=fa[fa[j][i-1]][i-1];
	while(q--)
	{
		cin>>u>>v,w=LCA(u,v);
		ans=dis[u]+dis[v]-2*dis[w];
		ans=(ans+(g[w]+1)*f[w])%mod;
		cout<<(ans+mod)%mod<<'\n';
	}
	return 0;
}