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
    
q = N // M
dat = [0] * (q + 1)
for i in range(1,q+1):
    m = i * M
    j = N // m
    dat[i] = j * (j - 1) * fact[N-2] % P
for i in range(q,0,-1):
    for j in range(2 * i,q + 1,i):
        dat[i] -= dat[j]
    dat[i] %= P
        
print(dat[1])