p,k = map(int,input().split())
mod = 10**9+7

zero = 1
notzero = 0

for i in range(k):
    zero_old = zero
    notzero_old = notzero
    zero = zero_old * (p+1) + notzero_old * 2
    zero %= mod
    notzero = zero_old * (p-1) + notzero_old * (2*p-2)
    notzero %= mod
    #print(zero,notzero)

print(zero%mod)