#include<bits/stdc++.h>
#define REP(i,n) for(int i=0;i<(n);i++)
#define ALL(v) (v).begin(),(v).end()
#define int long long
using namespace std;

//-----------------------------------------------------------------------

const int mod=1e9+7;
const int N=10010;
int fact[N],finv[N],inv[N];

int init()
{
    fact[0]=fact[1]=finv[0]=finv[1]=inv[1]=1;
    for(int i=2;i<N;i++){
	inv[i]=inv[mod%i]*(mod-mod/i)%mod;
	fact[i]=(fact[i-1]*i)%mod;
	finv[i]=(finv[i-1]*inv[i])%mod;
    }
}
int comb(int n,int k)
{
    if(k<0 or n<k) return 0;
    return fact[n]*finv[k]%mod*finv[n-k]%mod;
}

int mod_pow(int x,int n)
{
    int res=1;
    while(n>0){
	if(n&1) res=res*x%mod;
	x=x*x%mod; n>>=1;
    }
    return res;
}

signed main()
{
    init();
    int N,K; cin>>N>>K;
    vector<int> b(K+1,1);
    for(int i=1;i<=K;i++){
	int sum=0;
	for(int j=0;j<i;j++){
	    int tmp=b[j]*comb(i+1,j);
	    if(j&1) tmp=mod-tmp;
	    sum=(sum+tmp)%mod;
	}
	if(!(i&1)) sum=mod-sum;
	b[i]=sum*inv[i+1]%mod;
    }

    int ans=0;
    for(int i=0;i<=K;i++){
	int tmp=comb(K+1,i)*b[i]%mod*mod_pow(N%mod,K+1-i);
	ans=(ans+tmp)%mod;
    }
    cout<<ans*inv[K+1]%mod<<endl;

}