#include using namespace std; using ll = long long; constexpr char newl = '\n'; // https://noshi91.hatenablog.com/entry/2019/03/31/174006 template struct ModInt { using u64 = std::uint_fast64_t; static constexpr u64 MOD = Modulus; u64 val; constexpr ModInt(const u64 x = 0) noexcept : val(x % MOD) {} constexpr ModInt operator+() const noexcept { return ModInt(*this); } constexpr ModInt operator-() const noexcept { ModInt res(*this); if (res.val != 0) res.val = MOD - res.val; return res; } // prefix increment/decrement constexpr ModInt& operator++() noexcept { return *this += ModInt(1); } constexpr ModInt& operator--() noexcept { return *this -= ModInt(1); } // postfix increment/decrement constexpr ModInt& operator++(int) noexcept { ModInt tmp(*this); ++*this; return tmp; } constexpr ModInt& operator--(int) noexcept { ModInt tmp(*this); --*this; return tmp; } constexpr ModInt operator+(const ModInt& rhs) const noexcept { return ModInt(*this) += rhs; } constexpr ModInt operator-(const ModInt& rhs) const noexcept { return ModInt(*this) -= rhs; } constexpr ModInt operator*(const ModInt& rhs) const noexcept { return ModInt(*this) *= rhs; } constexpr ModInt operator/(const ModInt& rhs) const noexcept { return ModInt(*this) /= rhs; } constexpr ModInt& operator+=(const ModInt& rhs) noexcept { val += rhs.val; if (val >= MOD) val -= MOD; return *this; } constexpr ModInt& operator-=(const ModInt& rhs) noexcept { if (val < rhs.val) val += MOD; val -= rhs.val; return *this; } constexpr ModInt& operator*=(const ModInt& rhs) noexcept { val = val * rhs.val % MOD; return *this; } // prime Modulus only constexpr ModInt& operator/=(const ModInt& rhs) noexcept { return *this *= rhs.inv(); } // prime Modulus only constexpr ModInt inv() const noexcept { return pow(*this, MOD - 2); } }; template constexpr ModInt pow(ModInt x, std::uint_fast64_t n) { ModInt res(1); while (n) { if (n & 1) res *= x; x *= x; n >>= 1; } return res; } template istream& operator>>(istream& is, ModInt& x) { std::uint_fast64_t val; is >> val; x = ModInt(val); return is; } template ostream& operator<<(ostream& os, const ModInt& x) { return os << x.val; } using mint = ModInt<1000000007>; int main() { cin.tie(nullptr); ios::sync_with_stdio(false); mint p; int K; cin >> p >> K; vector dp(2, 0); dp[0] = 1; for (int i = 0; i < K; i++) { vector ndp(2, 0); ndp[0] = dp[0] * (p + 1) + dp[1] * 2; ndp[1] = dp[0] * (p - 1) + dp[1] * (p * 2 - 2); dp = move(ndp); } cout << dp[0] << newl; return 0; }