N,P=map(int,input().split())
mod=10**9+7
ans=0
NN=N
while NN:
    ans+=NN//P
    NN//=P
f=1
for i in range(1,N+1):
    f*=i
    f%=mod*(mod-1)
ans*=pow(f,f,mod*(mod-1))
ans%=mod
print(ans)