#!/usr/bin/ python3.8
import sys
read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines
from functools import lru_cache

N, M = map(int, read().split())
MOD = 10**9 + 7


@lru_cache(None)
def f(N):
    """ count 1 <= a < b <= N s.t. gcd(a,b) == 1"""
    if N <= 1:
        return 0
    ret = N * (N - 1) // 2
    sqN = int(N ** .5)
    for d in range(2, N + 1):
        n = N // d
        if n <= sqN:
            break
        ret -= f(n)
    for n in range(1, sqN + 1):
        cnt = (N // n) - (N // (n + 1))
        ret -= f(n) * cnt
    return ret


x = 2 * f(N // M) % MOD
for i in range(1, N - 1):
    x *= i
    x %= MOD
print(x)