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