#include using namespace std; #define fs first #define sc second #define pb push_back #define mp make_pair #define eb emplace_back #define ALL(A) A.begin(),A.end() #define RALL(A) A.rbegin(),A.rend() typedef long long LL; typedef pair P; template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } template T gcd(T a,T b){return b?gcd(b,a%b):a;} const LL mod=1000000007; const LL LINF=1LL<<62; const int INF=1<<30; int dx[]={1,0,-1,0,1,-1,1,-1}; int dy[]={0,1,0,-1,1,-1,-1,1}; int main(){ int n;cin >> n; vector a, b; for (int i = 0; i < n; i++) { int t,x;cin >> t >> x; x--; if(t==1) a.pb(x); else b.pb(x); } sort(ALL(a)); sort(ALL(b)); LL A=a.size(), B=b.size(); vector> dp(n+1,vector (n+1,0)); //dp[i][j] := iまで見てj個保留しているときの通り数 LL pa=0, pb=0; dp[0][0] = 1; for (int i = 0; i < n; i++) { LL cntb = 0; //i == b[k]の数(保留しているものとマッチングしなければならない) while(pa < A && a[pa] <= i) pa++; while(pb < B && b[pb] <= i){ pb++; cntb++; } for (int j = 0; j <= i; j++) { if(dp[i][j] == 0) continue; LL ra = pa - (i - (j - cntb + pb)); //paまで見てマッチングしていないaの数 if(j-cntb>=0 && ra>=0) dp[i+1][j-cntb] = (dp[i+1][j-cntb] + dp[i][j] * ra)%mod; //iをaとマッチング if(j+1-cntb >= 0) dp[i+1][j+1-cntb] = (dp[i+1][j+1-cntb] + (j + 1) * dp[i][j])%mod; //iをbとマッチング } } cout << dp[n][0] << endl; return 0; }