MOD = 10 ** 9 + 7

n, k = map(int, input().split())

m = 100
inv = [1 for _ in range(m + 1)]
for x in range(2, m + 1):
  inv[x] = (- (MOD // x) * inv[MOD % x]) % MOD

def rep_comb(n, r):
  ans = 1
  for x in range(r):
    ans = (ans * (n + r - 1 - x)) % MOD
    ans = (ans * inv[x + 1]) % MOD
  return ans

cur = n
divs = {}
d = 2
while d * d <= n:
  while cur % d == 0:
    divs[d] = divs.get(d, 0) + 1
    cur //= d
  d += 1

if cur > 1:
  divs[cur] = divs.get(cur, 0) + 1

ans = 1
for _, v in divs.items():
  cnt = 0
  for r in range(v + 1):
    cnt = (cnt + rep_comb(k, r)) % MOD
  ans = (ans * cnt) % MOD

print(ans)