p,k=map(int,input().split())
MOD=10**9+7
dp=[[0,0]]*(k+1)
dp[0][0]=1 
for i in range(1,k+1):
    dp[i][0]=dp[i-1][0]*(p+1)+dp[i-1][1]
    dp[i][0]%=MOD
    dp[i][1]=dp[i-1][0]*(p-1)+dp[i-1][1]*(2*p-1)
    dp[i][1]%=MOD
print(dp[k][0])