n,k= map(int,input().split())
mod = 10**9+7
pn = pow(n,k,mod)
pn1 = pow(n-1,k,mod)
print((n*(pn-pn1))%mod)