P = 10 ** 9 + 7
N, K = map(int, input().split())
a = K * (K + 1) // 2
print((pow(a + K, N, P) - pow(a, N, P)) % P)