MOD = 10 ** 9 + 7

n, k = map(int, input().split())
x = k * (k + 1) // 2 % MOD
ans = pow(x + k, n, MOD) - pow(x, n, MOD)
print(ans % MOD)