#include<iostream>
using namespace std;
long mod=1e9+7;
long power(long a,long b){return b?power(a*a%mod,b/2)*(b%2?a:1)%mod:1;}
long F[10101],I[10101];
long C(int N,int K){return F[N]*I[K]%mod*I[N-K]%mod;}
long N,K;
long S[10101];
main()
{
	F[0]=1;
	for(int i=1;i<10101;i++)F[i]=F[i-1]*i%mod;
	I[10100]=power(F[10100],mod-2);
	for(int i=10100;i--;)I[i]=I[i+1]*(i+1)%mod;
	cin>>N>>K;
	S[0]=N%mod;
	for(int i=1;i<=K;i++)
	{
		long now=power((N+1)%mod,i+1)+mod-1;
		for(int j=0;j<i;j++)
		{
			now=(now+mod-C(i+1,j)*S[j]%mod)%mod;
		}
		now=now*power(i+1,mod-2)%mod;
		S[i]=now;
	}
	cout<<S[K]<<endl;
}