#include #include #include #include using namespace std; typedef long long ll; typedef pair P; typedef pair T; long long int INF = 1e18; const ll fact_table = 800008; double Pi = 3.1415926535897932384626; priority_queue pql; priority_queue

pqp; priority_queue

bag; //big priority queue priority_queue ,greater > pqls; priority_queue ,greater

> pqps; //small priority queue //top pop ll dx[8]={1,0,-1,0,1,1,-1,-1}; ll dy[8]={0,1,0,-1,1,-1,-1,1}; char dir[] = "DRUL"; //↓,→,↑,← #define p(x) cout< number2){ swap(number1,number2); } if(number1 == 0 || number1 == number2){ return number2; }else{ return gcd(number2 % number1,number1); } } void YES(bool condition){ if(condition){ p("YES"); }else{ p("NO"); } return; } void Yes(bool condition){ if(condition){ p("Yes"); }else{ p("No"); } return; } /* ll fact[fact_table + 5],rfact[fact_table + 5]; void c3_init(){ fact[0] = rfact[0] = 1; for(ll i=1; i<=fact_table; i++){ fact[i] = (fact[i-1]*i) % mod; } rfact[fact_table] = rui(fact[fact_table],mod - 2); for(ll i=fact_table; i>=1; i--){ rfact[i-1] = rfact[i] * i; rfact[i-1] %= mod; } return;} ll c3(ll n,ll r){ return (((fact[n] * rfact[r]) % mod ) * rfact[n-r]) % mod;} */ bool multicase = false; ll n,m,num,sum,a,b,c,d,e,h,q; ll w, ans; ll k; char s[500005], t[500005]; ll x[500005], y[500005]; bool cond[5005]; ll dp[5005][4005]; void init(){ struct timeval _time; gettimeofday(&_time, NULL); ll usec = _time.tv_usec * 1000000; srand(usec); return; } void solve(){ cin >> n; vector right, left; for(ll i=0;i> a >> x[i]; if(a == 1){ left.pb(x[i]); }else{ right.pb(x[i]); } } sort(left.begin(), left.end()); sort(right.begin(), right.end()); dp[0][0] = 1; ll cand = 0, waste = 0; ll idx1 = 0, idx2 = 0; for(ll i=1;i<=n;i++){ while(idx1 < left.size() && left[idx1] <= i){ idx1 ++; cand ++; } for(ll j=0;j<=i;j++){ // dp[i][j] ? // 捨てる場合 if(j > 0)dp[i][j] += dp[i-1][j-1]; dp[i][j] %= mod; // 捨てない場合 if(i != j)dp[i][j] += dp[i-1][j] * (cand - (i - 1 - j)) % mod; dp[i][j] %= mod; } while(idx2 < right.size() && right[idx2] == i){ for(ll j=0;j<=i;j++){ if(j < waste){ dp[i][j] = 0; }else{ dp[i][j] *= (j - waste) % mod; dp[i][j] %= mod; } } idx2 ++; waste ++; } } ll ans = dp[n][right.size()]; p(ans); return; } int main(){ // init(); ll q; if(multicase){ cin >> q; }else{ q = 1; } while(q--){ solve(); } return 0; }