import sys
input = lambda : sys.stdin.readline().rstrip()

sys.setrecursionlimit(2*10**5+10)
write = lambda x: sys.stdout.write(x+"\n")
debug = lambda x: sys.stderr.write(x+"\n")
writef = lambda x: print("{:.12f}".format(x))


n,m = list(map(int, input().split()))
M = 10**9+7
vals = [0]*(n+1)
v = 0
v2 = 0
if m>n:
    ans = 0
else:
    for i in range(1,n+1)[::-1]:
        v = (n//i)*(n//i-1)%M
        for j in range(2*i,n+1,i):
            v -= vals[j]
            v %= M
        vals[i] = v
    ans = vals[m] % M
    for i in range(2, n-1):
        ans *= i
        ans %= M
print(ans%M)