#include<bits/stdc++.h>
using namespace std;
#include "circuit.h"
//#include "stub.cpp"
const int MAX=200005,mod=(1e9)+2022;
int cont[MAX],val[MAX],state[MAX],n,m;
vector<vector<int>>adj;
void dfs(int v){
    if(v>=n){
        val[v]=1;
        return;
    }
    int sz=adj[v].size();
    val[v]=sz;
    for(auto&u:adj[v]){
        dfs(u);
        val[v]=(val[v]*1LL*val[u])%mod;
    }
}
void dfs(int v,int cur){
    if(v>=n){
        cont[v-n]=cur;
        return;
    }
    int sz=adj[v].size();
    int pre[sz],suf[sz];
    for(int i=0;i<sz;i++)suf[i]=pre[i]=val[adj[v][i]];
    for(int i=1;i<sz;i++)pre[i]=(pre[i]*1LL*pre[i-1])%mod;
    for(int i=sz-2;i>=0;i--)suf[i]=(suf[i]*1LL*suf[i+1])%mod;
    for(int i=0;i<sz;i++){
        int nw=cur;
        if(i)nw=(nw*1LL*pre[i-1])%mod;
        if(i+1<sz)nw=(nw*1LL*suf[i+1])%mod;
        dfs(adj[v][i],nw);
    }
}
void init(int N,int M,vector<int>P,vector<int>A){
    n=N;m=M;
    adj.resize(N+M);
    for(int i=1;i<N+M;i++)adj[P[i]].push_back(i);
    for(int i=0;i<M;i++)state[i]=A[i];
    dfs(0);
    dfs(0,1);
}
int count_ways(int L, int R) {
    for(int i=L;i<=R;i++)state[i-n]^=1;
    int res=0;
    for(int i=0;i<m;i++)res=(res+state[i]*cont[i])%mod;
    return res;
}