#include #include #include #include #include using namespace std; #include template struct combination{ vectorfac,ifac; combination(size_t N=0):fac(1,1),ifac(1,1) { make_table(N); } void make_table(size_t N) { if(fac.size()>N)return; size_t now=fac.size(); N=max(N,now*2); fac.resize(N+1); ifac.resize(N+1); for(size_t i=now;i<=N;i++)fac[i]=fac[i-1]*i; ifac[N]=1/fac[N]; for(size_t i=N;i-->now;)ifac[i]=ifac[i+1]*(i+1); } T factorial(size_t n) { make_table(n); return fac[n]; } T invfac(size_t n) { make_table(n); return ifac[n]; } T P(size_t n,size_t k) { if(nC; int N; int A[2<<17]; vectorG[2<<17]; pairpr[2<<17]; int ord[2<<17]; pairdfs(int u) { int sz=0; mint ret=1; for(int v:G[u]) { auto t=dfs(v); ret*=t.second; ret*=C.C(sz+t.first,t.first); sz+=t.first; } return make_pair(sz+1,ret); } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cin>>N; vectorQ; setE; for(int i=0;i>A[i]; if(A[i]==0)Q.push_back(i); E.insert(i); pr[i]=make_pair(-1,-1); } for(int i=0;i