n, k = map(int, input().split())
mod = 1000000000 + 7
ans = pow(k * (k + 3) // 2, n, mod) - pow(k * (k + 1) // 2, n, mod)
print((ans + mod) % mod)