N, P = map(int, input().split())

MOD = 10**9 + 7
count = 0
now = 1

total = 1
for i in range(1, N+1):
    total *= i
    total %= MOD

total = pow(total, total, MOD)


while P**now <= N:
    count += N // (P**now)
    count %= MOD
    now += 1

print((count*total) % MOD)