#include <iostream>
using namespace std;
typedef long long ll;

const ll mod=1e9+7;
ll N,K;

ll f(ll x){
	return x*(x+1)/2%mod;
}

int main(){
	cin>>N>>K;
	ll A=f(N),B=f(N-1)*f(K-1)%mod;
	for(int i=0;i<K;i++) (A*=N-i-1)%=mod;
	for(int i=0;i<K-2;i++) (B*=N-i-2)%=mod;
	cout<<(A+B)%mod<<endl;
}