#include #include using namespace std; class ModInt { using ll = long long; public: int value; static int MOD; ModInt(ll value = 0) { this->value = value % MOD; if (this->value < 0) this->value += MOD; } operator int() const noexcept { return this->value; } ModInt& operator=(const ModInt& x) { if (this != &x) { this->value = x.value; } return *this; } bool operator==(const ModInt& x) const { return this->value == x.value; } bool operator!=(const ModInt& x) const { return !(*this == x); } ModInt operator+() const { return value; } ModInt operator-() const { return MOD - value; } ModInt operator~() const { return (*this) ^ ModInt(MOD - 2); } ModInt operator++() { return *this += 1; } ModInt operator--() { return *this -= 1; } ModInt operator++(int) { ModInt before = *this; ++(*this); return before; } ModInt operator--(int) { ModInt before = *this; --(*this); return before; } template ModInt operator+=(const T& x) { int sum = this->value + int(x); return *this = (sum < MOD ? sum : sum - MOD); } template ModInt operator-=(const T& x) { int diff = this->value - int(x); return *this = (diff >= 0 ? diff : diff + MOD); } template ModInt operator*=(const T& x) { return *this = ll(this->value) * ll(int(x)) % MOD; } template ModInt operator/=(const T& x) { return *this = (this->value % int(x) == 0 ? ModInt(this->value / int(x)) : *this * ~x); } template ModInt operator^=(const T& x) { int n = int(x); if (n == 0) return ModInt(1); if (n & 1) { return (*this) = ((*this) ^ ModInt(n - 1)) * (*this); } else { return (*this) = ((*this) * (*this)) ^ ModInt(n / 2); } } template ModInt operator+(const T& b) const { return ModInt(*this) += b; } template ModInt operator-(const T& b) const { return ModInt(*this) -= b; } template ModInt operator*(const T& b) const { return ModInt(*this) *= b; } template ModInt operator/(const T& b) const { return ModInt(*this) /= b; } template ModInt operator^(const T& b) const { return ModInt(*this) ^= b; } }; int ModInt::MOD = 1000000007; ostream& operator<<(ostream& os, const ModInt& x) { return os << x.value; } istream& operator>>(istream& is, ModInt& x) { return is >> x.value; } int main() { int N, M, K; cin >> N >> M >> K; vector L(M), R(M); for (int i = 0; i < M; ++i) cin >> L[i] >> R[i]; vector dp(N + 2, 0); dp[1] = 1; for (int k = 0; k < K; ++k) { for (int i = 1; i <= N; ++i) { dp[i] += dp[i - 1]; } vector ndp(N + 2, 0); for (int i = 0; i < M; ++i) { ModInt add = dp[R[i]] - dp[L[i] - 1]; ndp[L[i]] += add; ndp[R[i] + 1] -= add; } for (int i = 1; i <= N; ++i) { ndp[i] += ndp[i - 1]; } dp = ndp; } cout << dp[N] << endl; return 0; }