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