p,k=map(int,input().split())
M=10**9+7
q1,q2=1,0
for i in range(k):
  nq1=q1+q1*p+q2+q2
  nq2=q1*(p-1)+q2*(p-1)+q2*(p-1)
  q1=nq1%M
  q2=nq2%M
print(q1)