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