#include using namespace std; typedef long long ll; //a^n (mod m) const long long mod = 1e9+7; long long modpow(long long a, long long n) { long long res = 1; while (n > 0) { if (n & 1) res = res * a % mod; a = a * a % mod; n >>= 1; } return res; } //mod mでの a の逆元を求める long long modinv(long long a) { long long b = mod, u = 1, v = 0; while (b) { long long t = a / b; a -= t * b; swap(a, b); u -= t * v; swap(u, v); } u %= mod; if (u < 0) u += mod; return u; } signed main(){ //cout << setprecision(12) ; ios::sync_with_stdio(false); cin.tie(0); ll n,p; cin>>n>>p; ll ans=0; ans += 2*(n-1)*(n-p) + p*(p-1); ans %= mod; ans *= modinv(2*n*(n-1)); ans %= mod; ans *= (n*(n-1)/2)%mod; ans%=mod; ans += n-p; ans%=mod; for(int i=n;i>n-p;i--){ ans *= i; ans %=mod; } cout<