#include <bits/stdc++.h>
using namespace std;

const int64_t MOD = 1e9+7;
void add(int64_t& a, int64_t b){
    a = (a+b) % MOD;
}
void mul(int64_t& a, int64_t b){
    a = a*b % MOD;
}

int main(){
    int64_t P, K;
    cin >> P >> K;

    static int64_t dp[1000001][2];
    dp[0][0] = 1;
    for(int i=0; i<K; i++){
        add(dp[i+1][0], dp[i][0]+dp[i][1]);
        add(dp[i+1][1], (P-1)*(dp[i][0]+dp[i][1]));
        add(dp[i+1][0], P*dp[i][0] + dp[i][1]);
        add(dp[i+1][1], (P-1)*dp[i][1]);
    }
    cout << dp[K][0] << endl;
    return 0;
}