#include <bits/stdc++.h>
#define rep(i, n) for(int i=0; i<(n); ++i)
#define rep2(i, s, n) for(int i=s; i<(n); ++i)
#define ALL(v) (v).begin(), (v).end()

using namespace std;

typedef long long ll;

static const int INTINF = (2147483647 >> 1); // 10^9 + 5000
static const ll LLINF = (9223372036854775807 >> 1);

namespace Kunitaka{
    template<
        typename TYPE,
        std::size_t SIZE
    >
    std::size_t array_length(const TYPE (&array)[SIZE])
    {
        return SIZE;
    }
}
using namespace Kunitaka;
namespace Printer{
    void br(){
        cout << endl;
    }

    void loop_cnt(int i){
        cout << "i = " << i << endl;
    }

    void loop_cnt(int i, int j){
        cout << "i = " << i << ", " << "j = " << j << endl;
    }

    template<class T>
    void print(T x){
        cout << x << endl;
    }

    template<class T, class S>
    void print(T x, S y){
        cout << x << "," << y << endl;
    }

    template<class T, class S, class U>
    void print(T x, S y, U z){
        cout << x << "," << y << "," << z << endl;
    }

    template<class T, class S>
    void print(pair<T, S> p){
        cout << p.first << ", " << p.second << endl;
    }

    template<
        typename TYPE,
        std::size_t SIZE
    >
    void print(const TYPE (&array)[SIZE]){
        int lim = array_length(array);
        for(int i=0; i<lim; i++){
            if(i) cout << " ";
            cout << array[i];
        }
        cout << endl;
    }

    template<class S>
    void print(vector<S> v){
        for(int i=0; i<v.size(); i++){
            if(i) cout << " ";
            cout<< v[i];
        }
        cout << endl;
    }

    template<class S>
    void print(vector<vector<S> > vv){
        for(int i=0; i<vv.size(); i++){
            for(int j=0; j<vv[i].size(); j++){
                if(j) cout << " ";
                cout << vv[i][j];
            }
            cout << endl;
        }
    }

    void yesno(bool x){
        if(x) cout << "Yes" << endl;
        else cout << "No" << endl;
    }
};
using namespace Printer;

typedef long long ll;
static const int MAX = 5000001;
static const ll MOD = 1000000007;

ll FACT[MAX];
ll IFACT[MAX];

// ax + by = gcd(a, b) となるような (x, y) を求める
// 多くの場合 a と b は互いに素として ax + by = 1 となる (x, y) を求める
long long extGCD(long long a, long long b, long long &x, long long &y) {
    if (b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    long long d = extGCD(b, a%b, y, x); // 再帰的に解く
    y -= a / b * x;
    return d;
}

// 負の数にも対応した mod (a = -11 とかでも OK)
inline long long mod(long long a, long long m) {
    return (a % m + m) % m;
}

// 逆元計算 (ここでは a と m が互いに素であることが必要)
long long modinv(long long a, long long m) {
    long long x, y;
    extGCD(a, m, x, y);
    return mod(x, m); // 気持ち的には x % m だが、x が負かもしれないので
}

ll nCr(ll n, ll r)
{
    if(n < 0 || r < 0 || r > n)return 0;
    if(r > n / 2)r = n - r;
    return FACT[n]*IFACT[n-r]%MOD*IFACT[r]%MOD;
}

void ex(){
    int N = MAX;
    FACT[0] = 1;
    IFACT[0] = 1;
    for(int i = 1;i <= N;i++){
        FACT[i] = FACT[i-1] * i % MOD;
        IFACT[i] = modinv(FACT[i], MOD);
    }
    
    ll n= 100;
    ll r = 10;
    ll res = nCr(n, r);
}


int main(int argc, const char * argv[]) {
    //提出時、消す----//
    //--------------//


    // input
    cout << fixed << setprecision(10);
    ll n, k; cin >> n >> k;


    int N = MAX;
    FACT[0] = 1;
    IFACT[0] = 1;
    for(int i = 1;i <= N;i++){
        FACT[i] = FACT[i-1] * i % MOD;
        IFACT[i] = modinv(FACT[i], MOD);
    }
    

    ll ans = 0;
    rep(i, min(k+1, n)){
        ll m = i; // 0の数
        // print(nCr(n-m+m-1, m));
        ans += nCr(n, m) % MOD * nCr(n-m+m-1, m) % MOD;
        ans %= MOD;
    }

    print(ans);

    
    
    return 0;
}


Measured memory_bytes: 78336