import sys input = sys.stdin.readline from collections import * N, K = map(int, input().split()) MOD = 10**9+7 print((pow(K*(K+3)//2, N, MOD)-pow(K*(K+1)//2, N, MOD))%MOD)