#include<bits/stdc++.h>
using namespace std;
#define FOR(i,s,t) for(int i=(s),_t=(t); i<=_t; ++i)
#define DOR(i,s,t) for(int i=(s),_t=(t); i>=_t; --i)
typedef long long ll;
const int N=4e5+50;
const int Mod=998244353;

int X[N];

struct node {
    int L,R,ls,rs;
} T[N];
int tid,ntid;
void build(int L,int R,int &p) {
    if(L+1==R) p=++ntid;
    else p=++tid;
    T[p].L=L;
    T[p].R=R;
    if(L+1==R) return ;
    // cerr << L << ' ' << R << ' ' << X[p] << endl;
    build(L,X[p],T[p].ls);
    build(X[p],R,T[p].rs);
}
vector<int> lp[N],rp[N];
vector<int> g[N];
int n,m;
int sum[N];
void clear() {
    FOR(i,0,n) sum[i]=0;
}
void add(int x) {
    for(; x<=n; x+=x&-x) {
        sum[x]++;
    }
}
int qry(int x) {
    int s=0;
    for(; x; x^=x&-x) {
        s+=sum[x];
    }
    return s;
}
int mark[N];
ll f[N][3];
void dfs(int x) {
    if(!T[x].ls) {
        if(mark[x]) {
            f[x][0]=1;
            f[x][1]=1;
            f[x][2]=0;
        } else {
            f[x][0]=1;
            f[x][1]=0;
            f[x][2]=1;
        }
        return ;
    }
    int l=T[x].ls,r=T[x].rs;
    dfs(l);
    dfs(r);
    FOR(u,0,1) FOR(i,0,2) FOR(j,0,2) {
        int a=u || (i==0 && j==0);
        int b=(i==0) || (u && j==0);
        int c=(j==0) || (u && i==0);
        // if(u==0 && i==2 && j==2) cerr << mark[x] << endl;
        if(a) {
            if((!b && i==1) || (!c && j==1)) ;
            else {
                f[x][0]=(f[x][0]+f[l][i]*f[r][j])%Mod;
            }
        } else {
            if(mark[x] || (!b && i==1) || (!c && j==1)) {
                int nb=(i==0) || (1 && j==0);
                int nc=(j==0) || (1 && i==0);
                if((!nb && i==1) || (!nc && j==1)) ;
                else {
                    // if(u==0 && i==2 && j==2) cerr << f[l][i]*f[r][j] << endl;
                    f[x][1]=(f[x][1]+f[l][i]*f[r][j])%Mod;
                }
            } else {
                f[x][2]=(f[x][2]+f[l][i]*f[r][j])%Mod;
            }
        }
    }
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cin >> n >> m;
    FOR(i,1,n-1) {
        cin >> X[i];
    }
    FOR(i,1,m) {
        int l,r;
        cin >> l >> r;
        rp[l].push_back(r);
        lp[r].push_back(l);
        if(l==0 && r==n) mark[1]=1;
    }
    
    int rt=0;
    ntid=n-1;
    build(0,n,rt);
    tid=ntid;

    // cerr << " -- " << endl;
    FOR(i,1,tid) {
        if(T[i].ls) {
            g[T[i].L].push_back(i);
        }
    }
    FOR(i,0,n-1) {
        for(auto r: rp[i]) {
            add(r);
        }
        for(auto j: g[i]) {
            if(qry(T[j].R-1)-qry(T[T[j].ls].R-1)) mark[T[j].ls]=1;
        }
    }
    // cerr << " -- " << endl;
    clear();
    FOR(i,0,n) g[i].clear();
    FOR(i,1,tid) {
        if(T[i].rs) {
            g[T[i].R].push_back(i);
        }
    }
    DOR(i,n,1) {
        for(auto l: lp[i]) {
            add(l+1);
        }
        for(auto j: g[i]) {
            if(qry(T[T[j].rs].L+1)-qry(T[j].L+1)) mark[T[j].rs]=1;
        }
    }
    // FOR(i,1,tid) if(mark[i]) cerr << i << endl;
    dfs(1);
    // cerr << f[2][0] << ' ' << f[2][1] << ' ' << f[2][2] << endl;
    ll res=(f[1][0]+f[1][2])%Mod;
    cout << res << '\n';
    return 0;
}