n, k = map(int, input().split())
x = 1
y = 1
mod = 1000000000 + 7
for _ in range(n):
    x = x * k * (k + 3) // 2 % mod
    y = y * k * (k + 1) // 2 % mod
print((x - y + mod) % mod)