#include<bits/stdc++.h>

using namespace std;
using ll = long long;
const ll mod = 1e9 + 7;



int main()
{
    int p, k; cin >> p >> k;
    vector<vector<ll>> dp(k + 1, vector<ll>(2));   
    dp[0][0] = 1;

    for(int i = 0; i < k; i++) {
        dp[i + 1][0] = dp[i][0] * (p + 1) + dp[i][1] * 2 * (p - 1); // dp[i + 1][0] は dp[i][0]から0から掛けて遷移する通りがP個 + 足し算が1個 非零から足し算と掛け算
        dp[i + 1][0] %= mod;
        dp[i + 1][1] = dp[i][0] + dp[i][1] * 2 * (p - 1); //dp[i + 1][1]はdp[i][0]から足し算で一個, dp[i][1]から掛け算と割り算で(p-1)こずつ
    }
    cout << (dp[k][0]) % mod << endl;
}

//本質は非零に共通点があること