#define _USE_MATH_DEFINES #include #include using namespace std; /* #include using namespace atcoder; */ /* #include #include using bll = boost::multiprecision::cpp_int; using bdouble = boost::multiprecision::number>; using namespace boost::multiprecision; */ #if defined(LOCAL_TEST) || defined(LOCAL_DEV) #define BOOST_STACKTRACE_USE_ADDR2LINE #define BOOST_STACKTRACE_ADDR2LINE_LOCATION /usr/local/opt/binutils/bin/addr2line #define _GNU_SOURCE 1 #include #endif #ifdef LOCAL_TEST namespace std { template class dvector : public std::vector { public: dvector() : std::vector() {} explicit dvector(size_t n, const T& value = T()) : std::vector(n, value) {} dvector(const std::vector& v) : std::vector(v) {} dvector(const std::initializer_list il) : std::vector(il) {} template dvector(const Itr first, const Itr last) : std::vector(first, last) {} template , std::nullptr_t> = nullptr> std::vector::reference operator[](size_t n) { if (this->size() <= n) { std::cerr << boost::stacktrace::stacktrace() << '\n' << "vector::_M_range_check: __n (which is " << n << ") >= this->size() (which is " << this->size() << ")" << '\n'; } return this->at(n); } template , std::nullptr_t> = nullptr> const T_ operator[](size_t n) const { if (this->size() <= n) { std::cerr << boost::stacktrace::stacktrace() << '\n' << "vector::_M_range_check: __n (which is " << n << ") >= this->size() (which is " << this->size() << ")" << '\n'; } return this->at(n); } template , std::nullptr_t> = nullptr> T_& operator[](size_t n) { if (this->size() <= n) { std::cerr << boost::stacktrace::stacktrace() << '\n' << "vector::_M_range_check: __n (which is " << n << ") >= this->size() (which is " << this->size() << ")" << '\n'; } return this->at(n); } template , std::nullptr_t> = nullptr> const T_& operator[](size_t n) const { if (this->size() <= n) { std::cerr << boost::stacktrace::stacktrace() << '\n' << "vector::_M_range_check: __n (which is " << n << ") >= this->size() (which is " << this->size() << ")" << '\n'; } return this->at(n); } }; template, typename Allocator = std::allocator> class dmultiset : public std::multiset { public: dmultiset() : std::multiset() {} dmultiset(const std::multiset& m) : std::multiset(m) {} dmultiset(const std::initializer_list il) : std::multiset(il) {} dmultiset(const Compare& comp) : std::multiset(comp) {} const typename std::multiset::iterator erase(const typename std::multiset::iterator it) { return std::multiset::erase(it); } size_t erase([[maybe_unused]] const T& x) { std::cerr << boost::stacktrace::stacktrace() << '\n'; assert(false); } size_t erase_all_elements(const T& x) { return std::multiset::erase(x); } }; } #define vector dvector #define multiset dmultiset class SIGFPE_exception : std::exception {}; class SIGSEGV_exception : std::exception {}; void catch_SIGFPE([[maybe_unused]] int e) { std::cerr << boost::stacktrace::stacktrace() << '\n'; throw SIGFPE_exception(); } void catch_SIGSEGV([[maybe_unused]] int e) { std::cerr << boost::stacktrace::stacktrace() << '\n'; throw SIGSEGV_exception(); } signed convertedmain(); signed main() { signal(SIGFPE, catch_SIGFPE); signal(SIGSEGV, catch_SIGSEGV); return convertedmain(); } #define main() convertedmain() #else #define erase_all_elements erase #endif #ifdef LOCAL_DEV template std::ostream& operator<<(std::ostream& s, const std::pair& p) { return s << "(" << p.first << ", " << p.second << ")"; } template std::ostream& operator<<(std::ostream& s, const std::array& a) { s << "{ "; for (size_t i = 0; i < N; ++i){ s << a[i] << "\t"; } s << "}"; return s; } template std::ostream& operator<<(std::ostream& s, const std::set& se) { s << "{ "; for (auto itr = se.begin(); itr != se.end(); ++itr){ s << (*itr) << "\t"; } s << "}"; return s; } template std::ostream& operator<<(std::ostream& s, const std::multiset& se) { s << "{ "; for (auto itr = se.begin(); itr != se.end(); ++itr){ s << (*itr) << "\t"; } s << "}"; return s; } template std::ostream& operator<<(std::ostream& s, const std::map& m) { s << "{\n"; for (auto itr = m.begin(); itr != m.end(); ++itr){ s << "\t" << (*itr).first << " : " << (*itr).second << "\n"; } s << "}"; return s; } template std::ostream& operator<<(std::ostream& s, const std::deque& v) { for (size_t i = 0; i < v.size(); ++i){ s << v[i]; if (i < v.size() - 1) s << "\t"; } return s; } template std::ostream& operator<<(std::ostream& s, const std::vector& v) { for (size_t i = 0; i < v.size(); ++i){ s << v[i]; if (i < v.size() - 1) s << "\t"; } return s; } template std::ostream& operator<<(std::ostream& s, const std::vector>& vv) { s << "\\\n"; for (size_t i = 0; i < vv.size(); ++i){ s << vv[i] << "\n"; } return s; } void debug_impl() { std::cerr << '\n'; } template void debug_impl(Head head, Tail... tail) { std::cerr << " " << head << (sizeof...(tail) ? "," : ""); debug_impl(tail...); } #define debug(...) do { std::cerr << ":" << __LINE__ << " (" << #__VA_ARGS__ << ") ="; debug_impl(__VA_ARGS__); } while (false) constexpr inline long long prodlocal([[maybe_unused]] long long prod, [[maybe_unused]] long long local) { return local; } #else #define debug(...) do {} while (false) constexpr inline long long prodlocal([[maybe_unused]] long long prod, [[maybe_unused]] long long local) { return prod; } #endif //#define int long long using ll = long long; //INT_MAX = (1<<31)-1 = 2147483647, INT64_MAX = (1LL<<63)-1 = 9223372036854775807 constexpr ll INF = numeric_limits::max() == INT_MAX ? (ll)1e9 + 7 : (ll)1e18; constexpr ll MOD = (ll)1e9 + 7; //primitive root = 5 //constexpr ll MOD = 998244353; //primitive root = 3 constexpr double EPS = 1e-9; constexpr ll dx[4] = {1, 0, -1, 0}; constexpr ll dy[4] = {0, 1, 0, -1}; constexpr ll dx8[8] = {1, 0, -1, 0, 1, 1, -1, -1}; constexpr ll dy8[8] = {0, 1, 0, -1, 1, -1, 1, -1}; #define rep(i, n) for(ll i=0, i##_length=(n); i< i##_length; ++i) #define repeq(i, n) for(ll i=1, i##_length=(n); i<=i##_length; ++i) #define rrep(i, n) for(ll i=(n)-1; i>=0; --i) #define rrepeq(i, n) for(ll i=(n) ; i>=1; --i) #define all(v) (v).begin(), (v).end() #define rall(v) (v).rbegin(), (v).rend() void p() { std::cout << '\n'; } template void p(Head head, Tail... tail) { std::cout << head << (sizeof...(tail) ? " " : ""); p(tail...); } template inline void pv(std::vector& v) { for(ll i=0, N=v.size(); i inline bool chmax(T& a, T b) { return a < b && (a = b, true); } template inline bool chmin(T& a, T b) { return a > b && (a = b, true); } template inline void uniq(std::vector& v) { v.erase(std::unique(v.begin(), v.end()), v.end()); } template inline ll sz(T& v) { return v.size(); } template std::vector make_vector_impl(std::vector& sizes, typename std::enable_if<(N==1), T const&>::type x) { return std::vector(sizes.front(),x); } template auto make_vector_impl(std::vector& sizes, typename std::enable_if<(N>1), T const&>::type x) { ll size=sizes.back(); sizes.pop_back(); return std::vector(sizes,x))>(size,make_vector_impl(sizes,x)); } template auto make_vector(ll const (&sizes)[N], T const& x=T()) { std::vector s(N); for(size_t i=0; i(s,x); } /*-----8<-----template-----8<-----*/ //[lib]modint.cpp #include constexpr int ModintConstexprSize = prodlocal(35000, 10); template struct ModintConstexpr { unsigned int inv[N + 1]; constexpr inline long long modpow(long long a, long long n) { long long res = 1; while (n) { if (n & 1) res = res * a % MOD; a = a * a % MOD; n >>= 1; } return res; } constexpr ModintConstexpr() : inv() { inv[1] = 1; for (int i = 2; i <= N; ++i) inv[i] = modpow(i, MOD - 2); } }; #ifdef __INTELLISENSE__ #pragma diag_suppress 28 #endif constexpr auto modintconstexpr = ModintConstexpr<>(); namespace atcoder { template * = nullptr> struct static_modint_custom : internal::static_modint_base { using mint = static_modint_custom; #if __has_include() using inv_cache_type = __gnu_pbds::gp_hash_table; #else using inv_cache_type = map; #endif static inv_cache_type& inv_cache() { static inv_cache_type inv_cache; return inv_cache; } public: static constexpr int mod() { return m; } static mint raw(int v) { mint x; x._v = v; return x; } static_modint_custom() : _v(0) {} template * = nullptr> static_modint_custom(T v) { long long x = (long long)(v % (long long)(umod())); if (x < 0) x += umod(); _v = (unsigned int)(x); } template * = nullptr> static_modint_custom(T v) { _v = (unsigned int)(v % umod()); } unsigned int val() const { return _v; } mint& operator++() { _v++; if (_v == umod()) _v = 0; return *this; } mint& operator--() { if (_v == 0) _v = umod(); _v--; return *this; } mint operator++(int) { mint result = *this; ++*this; return result; } mint operator--(int) { mint result = *this; --*this; return result; } mint& operator+=(const mint& rhs) { _v += rhs._v; if (_v >= umod()) _v -= umod(); return *this; } mint& operator-=(const mint& rhs) { _v -= rhs._v; if (_v >= umod()) _v += umod(); return *this; } mint& operator*=(const mint& rhs) { unsigned long long z = _v; z *= rhs._v; _v = (unsigned int)(z % umod()); return *this; } mint& operator/=(const mint& rhs) { return *this = *this * rhs.inv(); } mint operator+() const { return *this; } mint operator-() const { return mint() - *this; } mint pow(long long n) const { mint x = *this, r = 1; if (n < 0) { x = inv(), n = -n; } while (n) { if (n & 1) r *= x; x *= x; n >>= 1; } return r; } mint inv() const { if (umod() == MOD && _v <= ModintConstexprSize) return raw(modintconstexpr.inv[_v]); auto& cachemap = inv_cache(); unsigned int cached = cachemap[_v]; if (cached) return raw(cached); mint x = inv_calc(); cachemap[_v] = x._v; return x; } mint inv_calc() const { if (prime) { assert(_v); return pow(umod() - 2); } else { auto eg = internal::inv_gcd(_v, m); assert(eg.first == 1); return eg.second; } } friend mint operator+(const mint& lhs, const mint& rhs) { return mint(lhs) += rhs; } friend mint operator-(const mint& lhs, const mint& rhs) { return mint(lhs) -= rhs; } friend mint operator*(const mint& lhs, const mint& rhs) { return mint(lhs) *= rhs; } friend mint operator/(const mint& lhs, const mint& rhs) { return mint(lhs) /= rhs; } friend bool operator==(const mint& lhs, const mint& rhs) { return lhs._v == rhs._v; } friend bool operator!=(const mint& lhs, const mint& rhs) { return lhs._v != rhs._v; } friend bool operator<(const mint& lhs, const mint& rhs) { return lhs._v < rhs._v; } private: unsigned int _v; static constexpr unsigned int umod() { return m; } static constexpr bool prime = internal::is_prime; }; } // namespace atcoder using Modint = atcoder::static_modint_custom; std::istream& operator>>(std::istream& is, Modint& x){ ll a; is>>a; x = a; return is; } std::ostream& operator<<(std::ostream& os, const Modint& x){ return os << x.val(); } inline Modint operator""_M(unsigned long long x) noexcept { return Modint(x); } std::vector fac_vec(1, 1), inv_vec(1, 1); inline void reserve(size_t a){ if(fac_vec.size() >= a) return; if(a < fac_vec.size() * 2) a = fac_vec.size() * 2; if(a >= (size_t)Modint::mod()) a = Modint::mod(); fac_vec.reserve(a); while(fac_vec.size() < a) fac_vec.push_back(fac_vec.back() * Modint(fac_vec.size())); inv_vec.resize(fac_vec.size()); inv_vec.back() = fac_vec.back().inv(); for(ll i = inv_vec.size() - 1; !inv_vec[i - 1].val(); i--) inv_vec[i - 1] = inv_vec[i] * i; } inline Modint factorial(ll n){ if(n < 0) return 0; reserve(n + 1); return fac_vec[n]; } inline Modint nPk(ll n, ll r){ if(r < 0 || n < r) return 0; if(n >> 24){ Modint ans = 1; for(ll i = 0; i < r; i++) ans *= n--; return ans; } reserve(n + 1); return fac_vec[n] * inv_vec[n - r]; } inline Modint nCk(ll n, ll r){ if(r < 0 || n < r) return 0; r = min(r, n - r); reserve(r + 1); return inv_vec[r] * nPk(n, r); } inline Modint nHk(ll n, ll r){ return nCk(n + r - 1, n - 1); } //n種類のものから重複を許してr個選ぶ=玉r個と仕切りn-1個 inline Modint catalan(ll n){ reserve(n * 2 + 1); return fac_vec[n * 2] * inv_vec[n] * inv_vec[n + 1]; } /*-----8<-----library-----8<-----*/ void solve() { ll N,M; cin>>N>>M; Modint ans = nCk(N + N, N)*(N*2); rep(i, M) { ll type,x,y; cin >> type >> x >> y; if(type==1){ ans -= nCk(x+y,x) * nCk(N-x+N-y-1, N-y); } else { ans -= nCk(x+y,x) * nCk(N-x+N-y-1, N-x); } } p(ans); } signed main() { std::cin.tie(nullptr); std::ios::sync_with_stdio(false); //ll Q; cin >> Q; while(Q--)solve(); solve(); return 0; }