N, M = map(int, input().split()) mod = int(1e9+7) for X in range(1, M+1): q = (N+1) // X r = (N+1) % X ans = pow(q, X-r-1, mod) * pow(q+1, r, mod) * (q-1) ans %= mod print(ans)