def solve():
  N, K = map(int, input().split())
  mod = 10**9+7
  ans = pow(K*(K+3)//2,N,mod)-pow(K*(K+1)//2,N,mod)
  return ans % mod
print(solve())