n, k = map(int, input().split())
mod = 10 ** 9 + 7
ans = 0
coma = 1
comb = 1
for i in range(min(k+1, n)):
    ans += coma * comb
    ans %= mod
    coma *= (n-i) * pow(i+1, mod-2, mod)
    coma %= mod
    comb *= (n-i-1) * pow(i+1, mod-2, mod)
    comb %= mod

print(ans) n, k = map(int, input().split()) n map(int, input().split()) map int input().split() ().split () input split k mod = 10 ** 9 + 7 mod 10 ** 9 + 7 10 ** 9 10 9 7 ans = 0 ans 0 coma = 1 coma 1 comb = 1 comb 1 for i in range(min(k+1, n)):
    ans += coma * comb
    ans %= mod
    coma *= (n-i) * pow(i+1, mod-2, mod)
    coma %= mod
    comb *= (n-i-1) * pow(i+1, mod-2, mod)
    comb %= mod i range(min(k+1, n)) range min(k+1, n) min k+1 k 1 n ans += coma * comb ans coma * comb coma comb ans %= mod ans mod coma *= (n-i) * pow(i+1, mod-2, mod) coma (n-i) * pow(i+1, mod-2, mod) n-i n i pow(i+1, mod-2, mod) pow i+1 i 1 mod-2 mod 2 mod coma %= mod coma mod comb *= (n-i-1) * pow(i+1, mod-2, mod) comb (n-i-1) * pow(i+1, mod-2, mod) n-i-1 n-i n i 1 pow(i+1, mod-2, mod) pow i+1 i 1 mod-2 mod 2 mod comb %= mod comb mod print(ans) print ans k = map(int, input().split()) map(int, input().split()) k comb = 1 1 comb coma %= mod mod coma coma = 1 1 coma ans %= mod mod ans comb %= mod mod comb mod = 10 ** 9 + 7 10 ** 9 + 7 mod n, k = map(int, input().split()) map(int, input().split()) n comb *= (n-i-1) * pow(i+1, mod-2, mod) (n-i-1) * pow(i+1, mod-2, mod) comb coma *= (n-i) * pow(i+1, mod-2, mod) (n-i) * pow(i+1, mod-2, mod) coma ans = 0 0 ans ans += coma * comb coma * comb ans
n, k = map(int, input().split())
mod = 10 ** 9 + 7
ans = 0
coma = 1
comb = 1
for i in range(min(k+1, n)):
    ans += coma * comb
    ans %= mod
    coma *= (n-i) * pow(i+1, mod-2, mod)
    coma %= mod
    comb *= (n-i-1) * pow(i+1, mod-2, mod)
    comb %= mod

print(ans)