local mod = 1000000007
local mfl = math.floor

local function bmul(x, y)
  local x1, y1 = mfl(x / 31623), mfl(y / 31623)
  local x0, y0 = x - x1 * 31623, y - y1 * 31623
  return (x1 * y1 * 14122 + (x1 * y0 + x0 * y1) * 31623 + x0 * y0) % mod
end

local function badd(x, y)
  return (x + y) % mod
end

local p, k = io.read("*n", "*n")
local dp1, dp2 = 1, 0
local dst1, dst2 = 0, 0
for i = 1, k do
  dst1 = badd(bmul(dp1, p + 1), bmul(dp2, 2))
  dst2 = badd(bmul(dp1, p - 1), bmul(2, bmul(dp2, p - 1)))
  dp1, dp2 = dst1, dst2
end
print(dp1)