#include using namespace std; using ll = long long; const ll MOD = 1e9 + 7; template struct ModInt { ll value; ModInt() : value(0) {} ModInt(const ll& v) : value(v % MOD) { if ( value < 0 ) value += MOD; } ModInt& operator=(const ModInt& o) { value = o.value; return *this; } inline bool operator==(const ModInt& o) { return value == o.value; } inline bool operator!=(const ModInt& o) { return value != o.value; } ModInt& operator+=(const ModInt& o) { value += o.value; if ( value >= MOD ) value -= MOD; return *this; } ModInt& operator-=(const ModInt& o) { value += MOD - o.value; if (value >= MOD ) value -= MOD; return *this; } ModInt& operator*=(const ModInt& o) { value *= o.value; value %= MOD; return *this; } ModInt& operator/=(const ModInt& o) { value *= inverse(o.value); value %= MOD; return *this; } ModInt operator-() { return ModInt(-value); } ModInt operator+(const ModInt& o) { ModInt r(*this); r += o; return r; } ModInt operator-(const ModInt& o) { ModInt r(*this); r -= o; return r; } ModInt operator*(const ModInt& o) { ModInt r(*this); r *= o; return r; } ModInt operator/(const ModInt& o) { ModInt r(*this); r /= o; return r; } ModInt pow(ll n) { ModInt t(*this), r = 1; while ( n > 0 ) { if ( n & 1 ) r *= t; t *= t; n >>= 1; } return r; } private: ll inverse(ll a) { ModInt t(a); return t.pow(MOD - 2).value; } }; template std::ostream& operator<< (std::ostream& os, const ModInt& o) { os << o.value; return os; } template std::istream& operator>> (std::istream& is, ModInt& o) { is >> o.value; return is; } //--------------------------- // 二項係数 template struct Combination { ll n_max; vector fac, finv, inv; Combination(ll n_max) { fac.resize(n_max+1); finv.resize(n_max+1); inv.resize(n_max+1); fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[0] = 0, inv[1] = 1; for (int i = 2; i <= n_max; i++){ fac[i] = fac[i - 1] * i; inv[i] = - inv[MOD%i] * (MOD / i); finv[i] = finv[i - 1] * inv[i]; } } T C(ll n, ll k) { if (n < k) return 0; if (n < 0 || k < 0) return 0; return fac[n] * finv[k] * finv[n - k]; } T P(ll n, ll k) { if (n < k) return 0; if (n < 0 || k < 0) return 0; return fac[n] * finv[n - k]; } T H(ll n, ll k) { return C(n+k-1, k); } }; using Int = ModInt; ll solve() { ll N, M; cin >> N >> M; vector T(M), X(M), Y(M); for ( int i = 0; i < M; i++ ) { cin >> T[i] >> X[i] >> Y[i]; } Combination cb(N*2+10); Int ans = Int(2*N) * cb.C(2*N,N); for ( int i = 0; i < M; i++ ) { ll nx = X[i]; if ( T[i] == 1 ) nx++; ans -= cb.C(X[i]+Y[i],X[i]) * cb.C(2*N-(X[i]+Y[i]+1),N-nx); } return ans.value; } int main() { auto ans = solve(); cout << ans << "\n"; return 0; }