MOD = 10 **9 + 7
def main():
    p,K = map(int,input().split())
    dp = [0] * (K + 1)
    dp[0] = 1
    for i in range(K):
        dp[i + 1] = (p - 1)*dp[i] + 2*pow(2*p,i,MOD)
        dp[i + 1] %= MOD
    print(dp[-1])
if __name__ == '__main__':
    main()