N,K=map(int,input().split())
mod=10**9+7


x=N*(N-1)//2%mod

#import math
#A=math.factorial(N)//math.factorial(N-K)

A=1
for i in range(K):
    A=A*(N-i)%mod

B=x*A*((K-1)*K//2)*pow(N,mod-2,mod)*pow((N-1),mod-2,mod)%mod
C=N*(N+1)//2*A*(N-K)*pow(N,mod-2,mod)%mod

print((B+C)%mod)