#yuki644

mod=10**9+7
N=10**5+10
n,m=map(int,raw_input().split())
l=[i for i in xrange(N)]
for i in xrange(2,N):
 if l[i]==i:
  j=i
  while j<N:
   l[j]=l[j]/i*(i-1)
   j+=i
res=0
for i in xrange(2,n/m+1):
 res+=2*l[i]
 res%=mod
for i in xrange(1,n-1):
 res*=i
 res%=mod
print res