Q = 10 ** 9 + 7
N, P = map(int,input().split())
A = 0
B = C = 1
for i in range(1, N + 1):
	B, C = B * i % Q, C * i % (Q - 1)
while N:
	N //= P
	A += N
print(A * pow(B, C, Q) % Q)