N,M = map(int,input().split())

P = 10 ** 9 + 7

import sys
if N < 2 * M:
    print(0)
    exit()

C = N + 5
fact = [1] * C
for i in range(2,C):
    fact[i] = fact[i-1] * i % P
    
Sum = 0
ans = 0
for i in range(N // M,0,-1):
    m = i * M
    if N < 2 * m:continue
    j = N // m
    tmp = j * (j - 1) * fact[N-2] % P
    ans = tmp - Sum
    Sum += ans
print(ans)