#include using namespace std; //a^b mod p O(logP) int64_t ModPow(int64_t a, int64_t b, int64_t p) { int64_t ans = 1; while(b > 0) { if(b & 1) ans = ans * a % p; //bの最下位bitが1ならa^(2^i)をかける a = a * a % p; b >>= 1; //bを1bit右にシフト } return ans; } vector fact; // 要素数(N + 1) vector fact_inv; // 要素数(N + 1) //a! mod pと(a!)^-1 mod pを配列に記録 O(NlogP) (やや遅い まだ高速化可能) //ModComb実行前にこれを実行すること void CombInit(int64_t n, int64_t p) { fact.resize(n + 1); fact_inv.resize(n + 1); fact.at(0) = 1; fact.at(1) = 1; for(int64_t i = 2; i <= n; i++) { fact.at(i) = fact.at(i - 1) * i; fact.at(i) %= p; } fact_inv.at(0) = 1; fact_inv.at(1) = 1; for(int64_t i = 2; i <= n; i++) fact_inv.at(i) = ModPow(fact.at(i), p - 2, p); } //nCk mod p O(1) int64_t ModComb(int64_t n, int64_t k, int64_t p) { return ((fact.at(n) * fact_inv.at(k)) % p * fact_inv.at(n - k)) % p; } int main() { int64_t n, m; cin >> n >> m; const int64_t mod = 1000000007; CombInit(m, mod); int64_t ans = 0; for(int64_t i = n; i <= m; i++) { ans += ModComb(i, n, mod); ans %= mod; } cout << ans << endl; return 0; }