/** * @FileName a.cpp * @Author kanpurin * @Created 2020.06.08 21:44:55 **/ #include "bits/stdc++.h" using namespace std; typedef long long ll; constexpr int MOD = 1e9 + 7; struct mint { private: long long x; public: mint(long long x = 0) : x((MOD + x) % MOD) {} mint(std::string &s) { long long z = 0; for (int i = 0; i < s.size(); i++) { z *= 10; z += s[i] - '0'; z %= MOD; } this->x = z; } mint &operator+=(const mint &a) { if ((x += a.x) >= MOD) x -= MOD; return *this; } mint &operator-=(const mint &a) { if ((x += MOD - a.x) >= MOD) x -= MOD; return *this; } mint &operator*=(const mint &a) { (x *= a.x) %= MOD; return *this; } mint &operator/=(const mint &a) { long long n = MOD - 2; mint u = 1, b = a; while (n > 0) { if (n & 1) { u *= b; } b *= b; n >>= 1; } return *this *= u; } mint operator+(const mint &a) const { mint res(*this); return res += a; } mint operator-(const mint &a) const { mint res(*this); return res -= a; } mint operator*(const mint &a) const { mint res(*this); return res *= a; } mint operator/(const mint &a) const { mint res(*this); return res /= a; } friend std::ostream &operator<<(std::ostream &os, const mint &n) { return os << n.x; } bool operator==(const mint &a) const { return this->x == a.x; } }; int main() { int n, m; cin >> n >> m; mint ans = 0; int k = n / m; for (int i = 2; i <= k; i++) { for (int j = i + i; j <= k; j += i) { ans -= 1; } } ans += mint(k - 1) * (k - 2) / 2 + k - 1; ans *= 2; for (int i = 0; i < n - 2; i++) { ans *= i + 1; } cout << ans << endl; return 0; }