N,M= map(int, input().split()) mod=int(1e9+7) for i in range(M): j=i+1 ans=pow(N//j+1,N%j,mod)*pow(N//j,j-N%j,mod) ans%=mod if j==N: ans=0 print(ans)