N,M = map(int,input().split()) mod = 10**9+7 for i in range(1,M+1): Q = (N+1)//i r = (N+1)%i res = pow(Q,i-r-1,mod) * pow(Q+1,r,mod) * (Q-1) print(res%mod)