import sys
input = sys.stdin.readline
from collections import *

P, K = map(int, input().split())
MOD = 10**9+7
dp = [0]*(K+1)
dp[0] = 1

for i in range(1, K+1):
    dp[i] = ((P+1)*dp[i-1])%MOD
    
    if i>=2:
        dp[i] += (pow(2*P, i-1, MOD)-dp[i-1])*2
        dp[i] %= MOD

print(dp[K])