#include #include #include #include #include #include using namespace std; using default_counter_t = int64_t; // macro #define let auto const & #define overload4(a, b, c, d, name, ...) name #define rep1(n) \ for (default_counter_t i = 0, end_i = default_counter_t(n); i < end_i; ++i) #define rep2(i, n) \ for (default_counter_t i = 0, end_##i = default_counter_t(n); i < end_##i; \ ++i) #define rep3(i, a, b) \ for (default_counter_t i = default_counter_t(a), \ end_##i = default_counter_t(b); \ i < end_##i; ++i) #define rep4(i, a, b, c) \ for (default_counter_t i = default_counter_t(a), \ end_##i = default_counter_t(b); \ i < end_##i; i += default_counter_t(c)) #define rep(...) overload4(__VA_ARGS__, rep4, rep3, rep2, rep1)(__VA_ARGS__) #define rrep1(n) \ for (default_counter_t i = default_counter_t(n) - 1; i >= 0; --i) #define rrep2(i, n) \ for (default_counter_t i = default_counter_t(n) - 1; i >= 0; --i) #define rrep3(i, a, b) \ for (default_counter_t i = default_counter_t(b) - 1, \ begin_##i = default_counter_t(a); \ i >= begin_##i; --i) #define rrep4(i, a, b, c) \ for (default_counter_t \ i = (default_counter_t(b) - default_counter_t(a) - 1) / \ default_counter_t(c) * default_counter_t(c) + \ default_counter_t(a), \ begin_##i = default_counter_t(a); \ i >= begin_##i; i -= default_counter_t(c)) #define rrep(...) \ overload4(__VA_ARGS__, rrep4, rrep3, rrep2, rrep1)(__VA_ARGS__) #define ALL(f, c, ...) \ (([&](decltype((c)) cccc) { \ return (f)(std::begin(cccc), std::end(cccc), ##__VA_ARGS__); \ })(c)) // function template constexpr C &Sort(C &a) { std::sort(std::begin(a), std::end(a)); return a; } template constexpr auto &Min(C const &a) { return *std::min_element(std::begin(a), std::end(a)); } template constexpr auto &Max(C const &a) { return *std::max_element(std::begin(a), std::end(a)); } template constexpr auto Total(C const &c) { return std::accumulate(std::begin(c), std::end(c), C(0)); } template auto CumSum(std::vector const &v) { std::vector a(v.size() + 1, T(0)); for (std::size_t i = 0; i < a.size() - size_t(1); ++i) a[i + 1] = a[i] + v[i]; return a; } template constexpr bool ChMax(T &a, T const &b) { if (a < b) { a = b; return true; } return false; } template constexpr bool ChMin(T &a, T const &b) { if (b < a) { a = b; return true; } return false; } void In(void) { return; } template void In(First &first, Rest &... rest) { cin >> first; In(rest...); return; } template void VectorIn(vector &v, const I n) { v.resize(size_t(n)); rep(i, v.size()) cin >> v[i]; } void Out(void) { cout << "\n"; return; } template void Out(First first, Rest... rest) { cout << first << " "; Out(rest...); return; } constexpr auto yes(const bool c) { return c ? "yes" : "no"; } constexpr auto Yes(const bool c) { return c ? "Yes" : "No"; } constexpr auto YES(const bool c) { return c ? "YES" : "NO"; } // http://koturn.hatenablog.com/entry/2018/08/01/010000 template inline std::vector MakeNdVector(T n, U val) noexcept { static_assert(std::is_integral::value, "[MakeNdVector] The 1st argument must be an integer"); return std::vector(std::forward(n), std::forward(val)); } template inline decltype(auto) MakeNdVector(T n, Args &&... args) noexcept { static_assert(std::is_integral::value, "[MakeNdVector] The 1st argument must be an integer"); return std::vector(args)...))>( std::forward(n), MakeNdVector(std::forward(args)...)); } template 0), std::nullptr_t>::type = nullptr> struct NdvectorImpl { using type = std::vector::type>; }; // struct ndvector_impl template struct NdvectorImpl { using type = std::vector; }; // struct ndvector_impl template using NdVector = typename NdvectorImpl::type; #ifdef USE_STACK_TRACE_LOGGER #ifdef __clang__ #pragma clang diagnostic push #pragma clang diagnostic ignored "-Weverything" #include #pragma clang diagnostic pop #endif //__clang__ #endif // USE_STACK_TRACE_LOGGER template struct ModInt { int64_t x; ModInt() : x(0) {} ModInt(int64_t y) : x(y >= 0 ? y % Mod : (Mod - (-y) % Mod) % Mod) {} ModInt &operator+=(const ModInt &p) { if ((x += p.x) >= Mod) x -= Mod; return *this; } ModInt &operator-=(const ModInt &p) { if ((x += Mod - p.x) >= Mod) x -= Mod; return *this; } ModInt &operator*=(const ModInt &p) { x = int64_t(1LL * x * p.x % Mod); return *this; } ModInt &operator/=(const ModInt &p) { *this *= p.Inverse(); return *this; } ModInt operator-() const { return ModInt(-x); } ModInt operator+(const ModInt &p) const { return ModInt(*this) += p; } ModInt operator-(const ModInt &p) const { return ModInt(*this) -= p; } ModInt operator*(const ModInt &p) const { return ModInt(*this) *= p; } ModInt operator/(const ModInt &p) const { return ModInt(*this) /= p; } bool operator==(const ModInt &p) const { return x == p.x; } bool operator!=(const ModInt &p) const { return x != p.x; } ModInt Inverse() const { int64_t a = x, b = Mod, u = 1, v = 0, t; while (b > 0) { t = a / b; swap(a -= t * b, b); swap(u -= t * v, v); } return ModInt(u); } ModInt Pow(int64_t n) const { ModInt ret(1), mul(x); while (n > 0) { if (n & 1) ret *= mul; mul *= mul; n >>= 1; } return ret; } friend ostream &operator<<(ostream &os, const ModInt &p) { return os << p.x; } friend istream &operator>>(istream &is, ModInt &a) { int64_t t; is >> t; a = ModInt(t); return (is); } static int64_t GetMod() { return Mod; } }; using Event = pair; signed main(int argc, char *argv[]) { (void)argc; #ifdef USE_STACK_TRACE_LOGGER google::InitGoogleLogging(argv[0]); google::InstallFailureSignalHandler(); #else (void)argv; #endif // USE_STACK_TRACE_LOGGER int64_t N, M, K; In(N, M, K); vector L(M), R(M); rep(i, M) { In(L[i], R[i]); } using Mi = ModInt<1000000007>; auto dp = MakeNdVector(K + 1, N + 1, Mi(0)); dp[0][1] = 1; rep(k, K) { auto cum = MakeNdVector(N + 2, Mi(0)); rep(i, 1, N + 1) { cum[i + 1] += cum[i] + dp[k][i]; } rep(i, M) { dp[k + 1][L[i]] += cum[R[i] + 1] - cum[L[i]]; if (R[i] + 1 <= N) { dp[k + 1][R[i] + 1] -= cum[R[i] + 1] - cum[L[i]]; } } rep(i, 1, N) { dp[k + 1][i + 1] += dp[k + 1][i]; } } cout << dp[K][N] << endl; return EXIT_SUCCESS; }