#yuki1049
mod=10**9+7
p,k=map(int,input().split())
a=1
b=0
for i in range(k):
 a,b=a*(p+1)+2*b,(a+2*b)*(p-1)
 a%=mod
 b%=mod
print(a)