import io
import sys

_INPUT = """\
6
2 2
6 3
"""

def solve(test):
  N,P=map(int,input().split())
  mod=10**9+7
  tmp=0
  for i in range(N):
    now=i+1
    while now%P==0:
      tmp+=1
      now//=P
  m=1
  for i in range(N):
    m*=i+1
    m%=mod
  n=1
  for i in range(N):
    n*=i+1
    n%=mod-2
  ans=pow(m,n,mod)*tmp%mod
  if test==0:
    print(ans)
  else:
    return None

def random_input():
  from random import randint,shuffle
  N=randint(1,10)
  M=randint(1,N)
  A=list(range(1,M+1))+[randint(1,M) for _ in range(N-M)]
  shuffle(A)
  return (" ".join(map(str, [N,M]))+"\n"+" ".join(map(str, A))+"\n")*3

def simple_solve():
  return []

def main(test):
  if test==0:
    solve(0)
  elif test==1:
    sys.stdin = io.StringIO(_INPUT)
    case_no=int(input())
    for _ in range(case_no):
      solve(0)
  else:
    for i in range(1000):
      sys.stdin = io.StringIO(random_input())
      x=solve(1)
      y=simple_solve()
      if x!=y:
        print(i,x,y)
        print(*[line for line in sys.stdin],sep='')
        break

#0:提出用、1:与えられたテスト用、2:ストレステスト用
main(0)