#include <iostream>
#include <vector>
using namespace std;
typedef long long ll;
typedef vector<ll> vl;

const int MAX_N=100005;
const ll mod=1e9+7;

vl Euler_Phi(){
	vl e(MAX_N);
	for(int i=0;i<MAX_N;i++) e[i]=i;
	for(int i=2;i<MAX_N;i++) if(e[i]==i){
		for(int j=i;j<MAX_N;j+=i) e[j]=e[j]/i*(i-1);
	}
	return e;
}

int N,M;

int main(){
	cin>>N>>M;
	vl e=Euler_Phi();
	ll res=0;
	for(int i=2;i<=N/M;i++) (res+=2*e[i])%=mod;
	for(int i=1;i<N-1;i++) (res*=i)%=mod;
	cout<<res<<endl;
}