def comb(n,r):
  if r>n:
    return 0
  tmp1=1
  tmp2=1
  for i in range(r):
    tmp1*=n-i
    tmp2*=i+1
    tmp1%=mod
    tmp2%=mod
  return (tmp1*pow(tmp2,mod-2,mod)%mod)

def fact(n):
  #fact(18)=[[2, 1],[3, 2]] (18 = 2**1 * 3**2)
  arr = []
  temp = n
  for i in range(2, int(-(-n**0.5//1))+1):
    if temp%i==0:
      cnt=0
      while temp%i==0:
        cnt+=1
        temp //= i
      arr.append([i, cnt])
  if temp!=1:
    arr.append([temp, 1])
  if arr==[]:
    arr.append([n, 1])
  return arr

m,k=map(int,input().split())
mod = 10**9+7

if m==1:
  print(1)
  exit()
ans=1

for i,j in fact(m):
  ans*=comb(k+j,j)
  ans%=mod
print(ans)