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