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 N%j==0: ans*=pow(N//j,mod-2,mod) ans*=pow(N//j,mod-2,mod) ans*=N//j+1 ans*=N//j-1 ans%=mod if N==1: ans=1 print(ans)