# include <bits/stdc++.h>
using namespace std;
using ll=long long;
constexpr int mod=1e9+7;
int fac[100010],finv[100010];
inline void add(int &x,const auto &y){x=(x+y)%mod;}
inline int C(int n,int m){return (ll)fac[n]*finv[m]%mod*finv[n-m]%mod;}
int power(int a,int b)
{
	int ans=1;
	for(;b;b>>=1,a=(ll)a*a%mod)
		if(b&1) ans=(ll)ans*a%mod;
	return ans;
}
int main()
{
	int n,m,ans=0;
	cin>>n>>m;m=n-m;
	fac[0]=fac[1]=1;
	for(int i=2;i<=n+1;i++) fac[i]=(ll)fac[i-1]*i%mod;
	finv[n+1]=power(fac[n+1],mod-2);
	for(int i=n;i>=0;i--) finv[i]=(ll)finv[i+1]*(i+1)%mod;
	for(int i=0;i<=n-m;i++) add(ans,(ll)(n+m+i&1?mod-C(n+1,i+m+1):C(n+1,i+m+1))*power(i,n));
	cout<<ans<<endl;
	return 0;
}