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