#include <bits/stdc++.h>
using namespace std;

#ifdef LOCAL
#include "dump.hpp"
#else
#define dump(...)
#define dumpv(...)
#endif

#define rep(i, n) for (int i = 0; i < (n); i++)
#define mins(x, y) (x = min(x, y))
#define maxs(x, y) (x = max(x, y))
using ll = long long;
using vi = vector<int>;
using vl = vector<ll>;
using vvi = vector<vi>;
using vvl = vector<vl>;
using P = pair<int, int>;
const int MOD = 1e9 + 7;
const int INF = 1001001001;
const ll LINF = 1001002003004005006ll;

struct mint {
  ll v;
  mint() : v(0) {}
  mint(ll v) : v((v % MOD + MOD) % MOD) {}
  // operator ll() const { return v; }
  mint pow(ll t) const;
  mint inv() const;
};
// mint& operator+=(mint& a, mint b) { return a = a.v + b.v; }
// mint& operator-=(mint& a, mint b) { return a = a.v - b.v; }
// mint& operator*=(mint& a, mint b) { return a = a.v * b.v; }
mint& operator+=(mint& a, mint b) {
  if ((a.v += b.v) >= MOD) a.v -= MOD;
  return a;
}
mint& operator-=(mint& a, mint b) {
  if ((a.v += MOD - b.v) >= MOD) a.v -= MOD;
  return a;
}
mint& operator*=(mint& a, mint b) {
  (a.v *= b.v) %= MOD;
  return a;
}
mint operator+(mint a, mint b) { return a += b; }
mint operator-(mint a, mint b) { return a -= b; }
mint operator*(mint a, mint b) { return a *= b; }
mint operator-(mint a) { return (mint)0 - a; }
mint& operator/=(mint& a, mint b) { return a *= b.inv(); }
mint operator/(mint a, mint b) { return a /= b; }
mint mint::pow(ll t) const {
  if (!t) return 1;
  mint a = pow(t >> 1);
  a *= a;
  if (t & 1) a *= *this;
  return a;
}
mint mint::inv() const { return pow(MOD - 2); }
bool operator==(mint a, mint b) { return a.v == b.v; }
bool operator!=(mint a, mint b) { return a.v != b.v; }
istream& operator>>(istream& os, mint& a) { return (os >> a.v); }
ostream& operator<<(ostream& os, const mint& a) { return (os << a.v); }

// struct mint {
//   ll x;
//   mint(ll x = 0) : x((x % MOD + MOD) % MOD) {}
//   mint operator-() const { return mint(-x); }
//   mint& operator+=(const mint a) {
//     if ((x += a.x) >= MOD) x -= MOD;
//     return *this;
//   }
//   mint& operator-=(const mint a) {
//     if ((x += MOD - a.x) >= MOD) x -= MOD;
//     return *this;
//   }
//   mint& operator*=(const mint a) {
//     (x *= a.x) %= MOD;
//     return *this;
//   }
//   mint operator+(const mint a) const {
//     mint res(*this);
//     return res += a;
//   }
//   mint operator-(const mint a) const {
//     mint res(*this);
//     return res -= a;
//   }
//   mint operator*(const mint a) const {
//     mint res(*this);
//     return res *= a;
//   }
//   mint pow(ll t) const {
//     if (!t) return 1;
//     mint a = pow(t >> 1);
//     a *= a;
//     if (t & 1) a *= *this;
//     return a;
//   }
//   mint inv() const { return pow(MOD - 2); }
//   mint& operator/=(const mint a) { return (*this) *= a.inv(); }
//   mint operator/(const mint a) const {
//     mint res(*this);
//     return res /= a;
//   }
// };
// ostream& operator<<(ostream& os, const mint& a) { return (os << a.x); }

mint mfact(int n) { // n! (mod MOD)
  static vector<mint> fact(1, 1);
  int pn = fact.size();
  if (pn < n + 1) {
    fact.resize(n + 1);
    for (int i = pn; i <= n; i++) fact[i] = fact[i - 1] * i;
  }
  return fact[n];
}
mint mifact(int n) { // n! (mod MOD)
  static vector<mint> ifact(1, 1);
  int pn = ifact.size();
  if (pn < n + 1) {
    ifact.resize(n + 1);
    ifact[n] = mfact(n).inv();
    for (int i = n; i >= pn; --i) ifact[i - 1] = ifact[i] * i;
  }
  return ifact[n];
}
mint nPk(int n, int k) {
  if (n < 0 || k < 0 || n < k) return 0;
  if (k == 0) return 1;
  return mfact(n) * mfact(n - k).inv();
}
mint nCk(int n, int k) {
  if (n < 0 || k < 0 || n < k) return 0;
  if (k == 0) return 1;
  // return mfact(n) * mfact(k).inv() * mfact(n - k).inv();
  return mfact(n) * mifact(k) * mifact(n - k);
}

void solve() {
  int n, k;
  cin >> n >> k;
  k = min(k, n - 1);
  mint ans = 0;
  mifact(n);
  dump(mfact(1));
  for (int i = 0; i <= k; i++) {
    ans += nCk(n, i) * nCk(n - i - 1 + i, i);
  }
  cout << ans << endl;
}

int main() {
  cin.tie(0);
  ios::sync_with_stdio(false);
  cout << fixed << setprecision(15);
  // freopen("temp.1", "r", stdin);
  solve();
  return 0;
}

Measured memory_bytes: 3328